HydrologyDCEfficientAnalysis.h

Go to the documentation of this file.

 
#ifndef _HydrologyDCEfficientAnalysis_
#define _HydrologyDCEfficientAnalysis_
 
/*Headers*/
#include "./Analysis.h"
class Node;
class Input2;
class HydrologyDCEfficientAnalysis: public Analysis{
 
   public:
      /*Model processing*/
      int  DofsPerNode(int** doflist,int domaintype,int approximation);
      void UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum);
      void UpdateElements(Elements* elements,Inputs2* inputs2,IoModel* iomodel,int analysis_counter,int analysis_type);
      void CreateNodes(Nodes* nodes,IoModel* iomodel,bool isamr=false);
      void CreateConstraints(Constraints* constraints,IoModel* iomodel);
      void CreateLoads(Loads* loads, IoModel* iomodel);
      void InitZigZagCounter(FemModel* femmodel);
      void ResetCounter(FemModel* femmodel);
 
      /*Finite element Analysis*/
      void           Core(FemModel* femmodel);
      ElementVector* CreateDVector(Element* element);
      ElementMatrix* CreateJacobianMatrix(Element* element);
      ElementMatrix* CreateKMatrix(Element* element);
      ElementVector* CreatePVector(Element* element);
      void GetB(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss);
      void GetSolutionFromInputs(Vector<IssmDouble>* solution,Element* element);
      void GradientJ(Vector<IssmDouble>* gradient,Element* element,int control_type,int control_index);
      void InputUpdateFromSolution(IssmDouble* solution,Element* element);
      void UpdateConstraints(FemModel* femmodel);
 
      /*Intermediaries*/
      IssmDouble EplStoring(Element* element,Gauss* gauss, Input2* epl_thick_input, Input2* epl_head_input, Input2* base_input);
      IssmDouble EplTransmitivity(Element* element,Gauss* gauss, Input2* epl_thick_input, Input2* epl_head_input, Input2* base_input);
      void GetHydrologyDCInefficientHmax(IssmDouble* ph_max,Element* element, Node* innode);
      IssmDouble GetHydrologyKMatrixTransfer(Element* element);
      IssmDouble GetHydrologyPVectorTransfer(Element* element, Gauss* gauss, Input2* sed_head_input);
      void ComputeEPLThickness(FemModel* femmodel);
      void HydrologyEPLGetMask(Vector<IssmDouble>* vec_mask, Vector<IssmDouble>* recurence, Element* element);
      void HydrologyEPLGetActive(Vector<IssmDouble>* active_vec, Element* element);
};
#endif