source: issm/trunk/src/c/objects/Penta.h@ 2112

Last change on this file since 2112 was 2112, checked in by Eric.Larour, 15 years ago

Heavy debugging of Mass flux capability
File size: 6.7 KB
Line 
1/*! \file Penta.h
2 * \brief: header file for penta object
3 */
4
5#ifndef _PENTA_H
6#define _PENTA_H
7
8#include "./Object.h"
9#include "./Element.h"
10#include "./Matpar.h"
11#include "./Matice.h"
12#include "./Node.h"
13#include "./Tria.h"
14#include "./ParameterInputs.h"
15
16class Penta: public Element{
17
18 private:
19 int id;
20
21 /*nodes: */
22 int node_ids[6]; //node ids
23 Node* nodes[6]; //node pointers
24 int node_offsets[6]; //node offsets in nodes dataset
25
26 /*materials: */
27 int mid;
28 Matice* matice;
29 int matice_offset;
30
31 int mparid;
32 Matpar* matpar;
33 int matpar_offset;
34
35
36 double h[6];
37 double s[6];
38 double b[6];
39 double k[6];
40 double melting[6];
41 double accumulation[6];
42 int friction_type;
43 double p;
44 double q;
45 int shelf;
46 int onbed;
47 int onsurface;
48 bool onwater;
49 double meanvel;/*!scaling ratio for velocities*/
50 double epsvel; /*!minimum velocity to avoid infinite velocity ratios*/
51 int collapse;
52 double geothermalflux[6];
53 int artdiff;
54 int thermal_steadystate;
55 double viscosity_overshoot;
56 double stokesreconditioning;
57
58 public:
59
60 Penta();
61 Penta( int id, int mid, int mparid, int node_ids[6], double h[6], double s[6], double b[6], double k[6], int friction_type,
62 double p, double q, int shelf, int onbed, int onsurface, double meanvel,double epsvel,
63 int collapse, double melting[6], double accumulation[6], double geothermalflux[6],
64 int artdiff, int thermal_steadystate,double viscosity_overshoot,double stokesreconditioning,bool onwater);
65 ~Penta();
66
67 void Echo();
68 void DeepEcho();
69 void Marshall(char** pmarshalled_dataset);
70 int MarshallSize();
71 char* GetName();
72 void Demarshall(char** pmarshalled_dataset);
73 int Enum();
74 int GetId();
75 int MyRank();
76 void Configure(void* loads,void* nodes,void* materials);
77 void CreateKMatrix(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
78 void CreateKMatrixDiagnosticHoriz( Mat Kgg, void* inputs, int analysis_type,int sub_analysis_type);
79 void CreateKMatrixDiagnosticVert( Mat Kgg, void* inputs, int analysis_type,int sub_analysis_type);
80 void CreatePVector(Vec pg, void* inputs, int analysis_type,int sub_analysis_type);
81 void UpdateFromInputs(void* inputs);
82 void GetDofList(int* doflist,int* pnumberofdofs);
83 void GetDofList1(int* doflist);
84 void* GetMatPar();
85 int GetShelf();
86 void GetNodes(void** nodes);
87 int GetOnBed();
88 void Du(Vec du_g,void* inputs,int analysis_type,int sub_analysis_type);
89 void Gradj(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type,char* control_type);
90 void GradjDrag(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type);
91 void GradjB(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type);
92 double Misfit(void* inputs,int analysis_type,int sub_analysis_type);
93
94 void GetThicknessList(double* thickness_list);
95 void GetBedList(double* bed_list);
96 Object* copy();
97 void* SpawnTria(int g0, int g1, int g2);
98
99 void GetStrainRate(double* epsilon, double* velocity, double* xyz_list, double* gauss_l1l2l3l4);
100 void GetB(double* pB, double* xyz_list, double* gauss_l1l2l3l4);
101 void GetBPrime(double* B, double* xyz_list, double* gauss_l1l2l3l4);
102 void GetB_vert(double* B, double* xyz_list, double* gauss_l1l2l3l4);
103 void GetBPrime_vert(double* B, double* xyz_list, double* gauss_l1l2l3l4);
104 void GetJacobianDeterminant(double* Jdet, double* xyz_list,double* gauss_l1l2l3l4);
105 void GetNodalFunctionsDerivativesBasic(double* dh1dh2dh3dh4dh5dh6_basic,double* xyz_list, double* gauss_l1l2l3l4);
106 void GetJacobian(double* J, double* xyz_list,double* gauss_l1l2l3l4);
107 void GetNodalFunctionsDerivativesParams(double* dl1dl2dl3dl4dl5dl6,double* gauss_l1l2l3l4);
108 void GetJacobianInvert(double* Jinv, double* xyz_list,double* gauss_l1l2l3l4);
109 void CreatePVectorDiagnosticHoriz( Vec pg, void* inputs,int analysis_type,int sub_analysis_type);
110 void CreatePVectorDiagnosticVert( Vec pg, void* inputs,int analysis_type,int sub_analysis_type);
111 void GetParameterValue(double* pvalue, double* v_list,double* gauss_l1l2l3l4);
112 void GetParameterDerivativeValue(double* p, double* p_list,double* xyz_list, double* gauss_l1l2l3l4);
113 void GetNodalFunctions(double* l1l2l3l4l5l6, double* gauss_l1l2l3l4);
114 void FieldExtrude(Vec field,double* field_serial,char* field_name, int iscollapsed);
115 void ComputePressure(Vec p_g);
116 void CreateKMatrixSlopeCompute(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
117 void CreatePVectorSlopeCompute( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
118 void CreateKMatrixPrognostic(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
119 void CreatePVectorPrognostic( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
120
121 void CreateKMatrixDiagnosticStokes( Mat Kgg, void* vinputs, int analysis_type,int sub_analysis_type);
122 void CreatePVectorDiagnosticStokes( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
123 void ReduceMatrixStokes(double* Ke_reduced, double* Ke_temp);
124 void GetMatrixInvert(double* Ke_invert, double* Ke);
125 void SurfaceNormal(double* surface_normal, double xyz_list[3][3]);
126 void GetStrainRateStokes(double* epsilon, double* velocity, double* xyz_list, double* gauss_coord);
127 void GetBStokes(double* B, double* xyz_list, double* gauss_coord);
128 void GetBprimeStokes(double* B_prime, double* xyz_list, double* gauss_coord);
129 void GetLStokes(double* LStokes, double* gauss_coord_tria);
130 void GetLprimeStokes(double* LprimeStokes, double* xyz_list, double* gauss_coord_tria, double* gauss_coord);
131 void GetNodalFunctionsDerivativesBasicStokes(double* dh1dh7_basic,double* xyz_list, double* gauss_coord);
132 void GetNodalFunctionsDerivativesParamsStokes(double* dl1dl7,double* gauss_coord);
133 void ReduceVectorStokes(double* Pe_reduced, double* Ke_temp, double* Pe_temp);
134 void GetNodalFunctionsStokes(double* l1l7, double* gauss_coord);
135 void CreateKMatrixThermal(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
136 void GetB_conduct(double* B_conduct, double* xyz_list, double* gauss_coord);
137 void GetB_advec(double* B_advec, double* xyz_list, double* gauss_coord);
138 void GetBprime_advec(double* Bprime_advec, double* xyz_list, double* gauss_coord);
139 void GetB_artdiff(double* B_artdiff, double* xyz_list, double* gauss_coord);
140 void CreateKMatrixMelting(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
141 void CreatePVectorThermal( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
142 void CreatePVectorMelting( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
143 void GetPhi(double* phi, double* epsilon, double viscosity);
144 double MassFlux(double* segment,double* ug);
145
146
147};
148#endif /* _PENTA_H */
149
150
Note: See TracBrowser for help on using the repository browser.