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

Last change on this file since 2956 was 2956, checked in by seroussi, 15 years ago

changed Params to Reference and Basic to nothing
File size: 6.8 KB
RevLine 
[1]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"
[96]13#include "./Tria.h"
[246]14#include "./ParameterInputs.h"
[1]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;
[2333]34
35 int numparid;
36 Numpar* numpar;
37 int numpar_offset;
[1]38
39
40 double h[6];
41 double s[6];
42 double b[6];
43 double k[6];
[128]44 double melting[6];
45 double accumulation[6];
[1]46 int friction_type;
47 double p;
48 double q;
49 int shelf;
50 int onbed;
51 int onsurface;
[2333]52 bool onwater;
[1]53 int collapse;
54 double geothermalflux[6];
55 int thermal_steadystate;
56
57 public:
58
59 Penta();
[2333]60 Penta( int id, int mid, int mparid, int numparid, int node_ids[6], double h[6], double s[6], double b[6], double k[6], int friction_type,
61 double p, double q, int shelf, int onbed, int onsurface, int collapse, double melting[6], double accumulation[6], double geothermalflux[6],
62 int thermal_steadystate,bool onwater);
[1]63 ~Penta();
64
65 void Echo();
[803]66 void DeepEcho();
[1]67 void Marshall(char** pmarshalled_dataset);
68 int MarshallSize();
69 char* GetName();
70 void Demarshall(char** pmarshalled_dataset);
71 int Enum();
72 int GetId();
73 int MyRank();
[2333]74 void Configure(void* loads,void* nodes,void* materials,void* parameters);
[465]75 void CreateKMatrix(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
76 void CreateKMatrixDiagnosticHoriz( Mat Kgg, void* inputs, int analysis_type,int sub_analysis_type);
77 void CreateKMatrixDiagnosticVert( Mat Kgg, void* inputs, int analysis_type,int sub_analysis_type);
78 void CreatePVector(Vec pg, void* inputs, int analysis_type,int sub_analysis_type);
[246]79 void UpdateFromInputs(void* inputs);
[96]80 void GetDofList(int* doflist,int* pnumberofdofs);
[301]81 void GetDofList1(int* doflist);
[246]82 void* GetMatPar();
[1]83 int GetShelf();
[246]84 void GetNodes(void** nodes);
[1]85 int GetOnBed();
[1185]86 void Du(Vec du_g,void* inputs,int analysis_type,int sub_analysis_type);
[1188]87 void Gradj(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type,char* control_type);
88 void GradjDrag(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type);
89 void GradjB(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type);
[1185]90 double Misfit(void* inputs,int analysis_type,int sub_analysis_type);
[1]91
92 void GetThicknessList(double* thickness_list);
93 void GetBedList(double* bed_list);
94 Object* copy();
[246]95 void* SpawnTria(int g0, int g1, int g2);
[1]96
[2956]97 void GetStrainRate(double* epsilon, double* velocity, double* xyz_list, double* gauss_coord);
98 void GetB(double* pB, double* xyz_list, double* gauss_coord);
99 void GetBPrime(double* B, double* xyz_list, double* gauss_coord);
100 void GetB_vert(double* B, double* xyz_list, double* gauss_coord);
101 void GetBPrime_vert(double* B, double* xyz_list, double* gauss_coord);
102 void GetJacobianDeterminant(double* Jdet, double* xyz_list,double* gauss_coord);
103 void GetNodalFunctionsDerivatives(double* dh1dh6,double* xyz_list, double* gauss_coord);
104 void GetJacobian(double* J, double* xyz_list,double* gauss_coord);
105 void GetNodalFunctionsDerivativesReference(double* dl1dl6,double* gauss_coord);
106 void GetJacobianInvert(double* Jinv, double* xyz_list,double* gauss_coord);
[465]107 void CreatePVectorDiagnosticHoriz( Vec pg, void* inputs,int analysis_type,int sub_analysis_type);
108 void CreatePVectorDiagnosticVert( Vec pg, void* inputs,int analysis_type,int sub_analysis_type);
[2956]109 void GetParameterValue(double* pvalue, double* v_list,double* gauss_coord);
110 void GetParameterDerivativeValue(double* p, double* p_list,double* xyz_list, double* gauss_coord);
111 void GetNodalFunctions(double* l1l6, double* gauss_coord);
[848]112 void FieldExtrude(Vec field,double* field_serial,char* field_name, int iscollapsed);
[301]113 void ComputePressure(Vec p_g);
[465]114 void CreateKMatrixSlopeCompute(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
115 void CreatePVectorSlopeCompute( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
[607]116 void CreateKMatrixPrognostic(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
117 void CreatePVectorPrognostic( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
[2714]118 void CreateKMatrixBalancedthickness(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
119 void CreatePVectorBalancedthickness( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
[2722]120 void CreateKMatrixBalancedvelocities(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
121 void CreatePVectorBalancedvelocities( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
[1]122
[465]123 void CreateKMatrixDiagnosticStokes( Mat Kgg, void* vinputs, int analysis_type,int sub_analysis_type);
124 void CreatePVectorDiagnosticStokes( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
[394]125 void ReduceMatrixStokes(double* Ke_reduced, double* Ke_temp);
126 void GetMatrixInvert(double* Ke_invert, double* Ke);
127 void SurfaceNormal(double* surface_normal, double xyz_list[3][3]);
128 void GetStrainRateStokes(double* epsilon, double* velocity, double* xyz_list, double* gauss_coord);
129 void GetBStokes(double* B, double* xyz_list, double* gauss_coord);
130 void GetBprimeStokes(double* B_prime, double* xyz_list, double* gauss_coord);
131 void GetLStokes(double* LStokes, double* gauss_coord_tria);
132 void GetLprimeStokes(double* LprimeStokes, double* xyz_list, double* gauss_coord_tria, double* gauss_coord);
[2956]133 void GetNodalFunctionsDerivativesStokes(double* dh1dh7,double* xyz_list, double* gauss_coord);
134 void GetNodalFunctionsDerivativesReferenceStokes(double* dl1dl7,double* gauss_coord);
[394]135 void ReduceVectorStokes(double* Pe_reduced, double* Ke_temp, double* Pe_temp);
136 void GetNodalFunctionsStokes(double* l1l7, double* gauss_coord);
[483]137 void CreateKMatrixThermal(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
138 void GetB_conduct(double* B_conduct, double* xyz_list, double* gauss_coord);
139 void GetB_advec(double* B_advec, double* xyz_list, double* gauss_coord);
140 void GetBprime_advec(double* Bprime_advec, double* xyz_list, double* gauss_coord);
[1676]141 void GetB_artdiff(double* B_artdiff, double* xyz_list, double* gauss_coord);
[483]142 void CreateKMatrixMelting(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
143 void CreatePVectorThermal( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
144 void CreatePVectorMelting( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
145 void GetPhi(double* phi, double* epsilon, double viscosity);
[2112]146 double MassFlux(double* segment,double* ug);
[358]147
[394]148
[1]149};
150#endif /* _PENTA_H */
[96]151
152
Note: See TracBrowser for help on using the repository browser.