source: issm/trunk/src/c/analyses/HydrologyGlaDSAnalysis.h@ 24313

/*! \file HydrologyGlaDSAnalysis.h 
 *  \brief: header file for generic external result object
 */

#ifndef _HydrologyGlaDSAnalysis_
#define _HydrologyGlaDSAnalysis_

/*Headers*/
#include "./Analysis.h"

class HydrologyGlaDSAnalysis: public Analysis{

        public:
                /*Model processing*/
                void CreateConstraints(Constraints* constraints,IoModel* iomodel);
                void CreateLoads(Loads* loads, IoModel* iomodel);
                void CreateNodes(Nodes* nodes,IoModel* iomodel,bool isamr=false);
                int  DofsPerNode(int** doflist,int domaintype,int approximation);
                void UpdateElements(Elements* elements,IoModel* iomodel,int analysis_counter,int analysis_type);
                void UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum);

                /*Finite element Analysis*/
                void           Core(FemModel* femmodel);
                ElementVector* CreateDVector(Element* element);
                ElementMatrix* CreateJacobianMatrix(Element* element);
                ElementMatrix* CreateKMatrix(Element* element);
                ElementVector* CreatePVector(Element* element);
                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);

                /*Specific to GlaDS*/
                void SetChannelCrossSectionOld(FemModel* femmodel);
                void UpdateSheetThickness(FemModel* femmodel);
                void UpdateSheetThickness(Element*  element);
                void UpdateChannelCrossSection(FemModel* femmodel);
                void UpdateEffectivePressure(FemModel* femmodel);
                void UpdateEffectivePressure(Element* element);
};
#endif