1 | /*! \file AdjointHorizAnalysis.h
|
---|
2 | * \brief: header file for generic external result object
|
---|
3 | */
|
---|
4 |
|
---|
5 | #ifndef _AdjointHorizAnalysis_
|
---|
6 | #define _AdjointHorizAnalysis_
|
---|
7 |
|
---|
8 | /*Headers*/
|
---|
9 | #include "./Analysis.h"
|
---|
10 |
|
---|
11 | class AdjointHorizAnalysis: public Analysis{
|
---|
12 |
|
---|
13 | public:
|
---|
14 | /*Model processing*/
|
---|
15 | int DofsPerNode(int** doflist,int domaintype,int approximation);
|
---|
16 | void UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum);
|
---|
17 | void UpdateElements(Elements* elements,IoModel* iomodel,int analysis_counter,int analysis_type);
|
---|
18 | void CreateNodes(Nodes* nodes,IoModel* iomodel);
|
---|
19 | void CreateConstraints(Constraints* constraints,IoModel* iomodel);
|
---|
20 | void CreateLoads(Loads* loads, IoModel* iomodel);
|
---|
21 |
|
---|
22 | /*Finite element Analysis*/
|
---|
23 | void Core(FemModel* femmodel);
|
---|
24 | ElementVector* CreateDVector(Element* element);
|
---|
25 | ElementMatrix* CreateJacobianMatrix(Element* element);
|
---|
26 | ElementMatrix* CreateKMatrix(Element* element);
|
---|
27 | ElementMatrix* CreateKMatrixSSA(Element* element);
|
---|
28 | ElementMatrix* CreateKMatrixHO(Element* element);
|
---|
29 | ElementMatrix* CreateKMatrixFS(Element* element);
|
---|
30 | ElementVector* CreatePVector(Element* element);
|
---|
31 | ElementVector* CreatePVectorSSA(Element* element);
|
---|
32 | ElementVector* CreatePVectorHO(Element* element);
|
---|
33 | ElementVector* CreatePVectorFS(Element* element);
|
---|
34 | void GetSolutionFromInputs(Vector<IssmDouble>* solution,Element* element);
|
---|
35 | void GradientJ(Vector<IssmDouble>* gradient,Element* element,int control_type,int control_index);
|
---|
36 | void GradientJDragGradient(Element* element,Vector<IssmDouble>* gradient,int control_index);
|
---|
37 | void GradientJBGradient(Element* element,Vector<IssmDouble>* gradient,int control_index);
|
---|
38 | void GradientJDragSSA(Element* element,Vector<IssmDouble>* gradient,int control_index);
|
---|
39 | void GradientJDragHO(Element* element,Vector<IssmDouble>* gradient,int control_index);
|
---|
40 | void GradientJDragFS(Element* element,Vector<IssmDouble>* gradient,int control_index);
|
---|
41 | void GradientJBbarSSA(Element* element,Vector<IssmDouble>* gradient,int control_index);
|
---|
42 | void GradientJBbarHO(Element* element,Vector<IssmDouble>* gradient,int control_index);
|
---|
43 | void GradientJBbarFS(Element* element,Vector<IssmDouble>* gradient,int control_index);
|
---|
44 | void GradientJDSSA(Element* element,Vector<IssmDouble>* gradient,int control_index);
|
---|
45 | void InputUpdateFromSolution(IssmDouble* solution,Element* element);
|
---|
46 | void UpdateConstraints(FemModel* femmodel);
|
---|
47 | void InputUpdateFromSolutionHoriz(IssmDouble* solution,Element* element);
|
---|
48 | void InputUpdateFromSolutionFS(IssmDouble* solution,Element* element);
|
---|
49 | };
|
---|
50 | #endif
|
---|