source: issm/trunk-jpl/src/c/classes/Loads/Numericalflux.h@ 15372

Last change on this file since 15372 was 15372, checked in by Mathieu Morlighem, 12 years ago

NEW:merged PentaHook and TriaHook into ElementHook
File size: 4.1 KB
Line 
1/*!\file Numericalflux.h
2 * \brief: header file for icefront object
3 */
4
5#ifndef _NUMERICALFLUX_H_
6#define _NUMERICALFLUX_H_
7
8/*Headers:*/
9/*{{{*/
10#include "./Load.h"
11class Hook;
12class Parameters;
13class Inputs;
14class IoModel;
15class ElementMatrix;
16class ElementVector;
17/*}}}*/
18
19class Numericalflux: public Load {
20
21 public:
22 int id;
23 int analysis_type;
24
25 /*Hooks*/
26 Hook *helement;
27 Hook *hnodes;
28 Hook *hvertices;
29
30 /*Corresponding fields*/
31 Element *element;
32 Vertex **vertices;
33 Node **nodes;
34 Parameters *parameters;
35 Inputs *inputs;
36
37 /*Numericalflux constructors,destructors {{{*/
38 Numericalflux();
39 Numericalflux(int numericalflux_id,int i,int i1,int i2,int e1,int e2,IoModel* iomodel,int analysis_type);
40 ~Numericalflux();
41 /*}}}*/
42 /*Object virtual functions definitions:{{{ */
43 void Echo();
44 void DeepEcho();
45 int Id();
46 int ObjectEnum();
47 Object *copy();
48 /*}}}*/
49 /*Update virtual functions resolution: {{{*/
50 void InputUpdateFromVector(IssmDouble* vector, int name, int type){/*Do nothing*/}
51 void InputUpdateFromVector(int* vector, int name, int type){_error_("Not implemented yet!");}
52 void InputUpdateFromVector(bool* vector, int name, int type){_error_("Not implemented yet!");}
53 void InputUpdateFromMatrixDakota(IssmDouble* matrix, int nrows, int ncols, int name, int type){/*Do nothing*/}
54 void InputUpdateFromVectorDakota(IssmDouble* vector, int name, int type){/*Do nothing*/}
55 void InputUpdateFromVectorDakota(int* vector, int name, int type){_error_("Not implemented yet!");}
56 void InputUpdateFromVectorDakota(bool* vector, int name, int type){_error_("Not implemented yet!");}
57 void InputUpdateFromConstant(IssmDouble constant, int name){/*Do nothing*/};
58 void InputUpdateFromConstant(int constant, int name){/*Do nothing*/};
59 void InputUpdateFromConstant(bool constant, int name){_error_("Not implemented yet!");}
60 void InputUpdateFromSolution(IssmDouble* solution){_error_("Not implemented yet!");}
61 void InputUpdateFromIoModel(int index, IoModel* iomodel){_error_("not implemented yet");};
62 /*}}}*/
63 /*Load virtual functions definitions: {{{*/
64 void Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters);
65 void SetCurrentConfiguration(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters);
66 void CreateKMatrix(Matrix<IssmDouble>* Kff, Matrix<IssmDouble>* Kfs);
67 void CreatePVector(Vector<IssmDouble>* pf);
68 void GetNodesSidList(int* sidlist);
69 int GetNumberOfNodes(void);
70 void CreateJacobianMatrix(Matrix<IssmDouble>* Jff){_error_("Not implemented yet");};
71 bool IsPenalty(void);
72 void PenaltyCreateJacobianMatrix(Matrix<IssmDouble>* Jff,IssmDouble kmax){_error_("Not implemented yet");};
73 void PenaltyCreateKMatrix(Matrix<IssmDouble>* Kff, Matrix<IssmDouble>* kfs, IssmDouble kmax);
74 void PenaltyCreatePVector(Vector<IssmDouble>* pf, IssmDouble kmax);
75 void SetwiseNodeConnectivity(int* d_nz,int* o_nz,Node* node,bool* flags,int set1_enum,int set2_enum);
76 bool InAnalysis(int analysis_type);
77 /*}}}*/
78 /*Numericalflux management:{{{*/
79 void GetNormal(IssmDouble* normal,IssmDouble xyz_list[4][3]);
80 ElementMatrix* CreateKMatrixPrognostic(void);
81 ElementMatrix* CreateKMatrixPrognosticInternal(void);
82 ElementMatrix* CreateKMatrixPrognosticBoundary(void);
83 ElementMatrix* CreateKMatrixBalancethickness(void);
84 ElementMatrix* CreateKMatrixBalancethicknessInternal(void);
85 ElementMatrix* CreateKMatrixBalancethicknessBoundary(void);
86 ElementMatrix* CreateKMatrixAdjointBalancethickness(void);
87 ElementMatrix* CreateKMatrixAdjointBalancethicknessInternal(void);
88 ElementMatrix* CreateKMatrixAdjointBalancethicknessBoundary(void);
89 ElementVector* CreatePVectorPrognostic(void);
90 ElementVector* CreatePVectorPrognosticInternal(void);
91 ElementVector* CreatePVectorPrognosticBoundary(void);
92 ElementVector* CreatePVectorBalancethickness(void);
93 ElementVector* CreatePVectorBalancethicknessInternal(void);
94 ElementVector* CreatePVectorBalancethicknessBoundary(void);
95 ElementVector* CreatePVectorAdjointBalancethickness(void);
96 /*}}}*/
97
98};
99
100#endif /* _NUMERICALFLUX_H_ */
Note: See TracBrowser for help on using the repository browser.