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

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

Added pressure computation
File size: 4.1 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;
34
35
36 double h[6];
37 double s[6];
38 double b[6];
39 double k[6];
[128]40 double melting[6];
41 double accumulation[6];
[1]42 int friction_type;
43 double p;
44 double q;
45 int shelf;
46 int onbed;
47 int onsurface;
48 double meanvel;/*!scaling ratio for velocities*/
49 double epsvel; /*!minimum velocity to avoid infinite velocity ratios*/
50 int collapse;
51 double geothermalflux[6];
52 int artdiff;
53 int thermal_steadystate;
[111]54 double viscosity_overshoot;
[1]55
56 public:
57
58 Penta();
[96]59 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,
[46]60 double p, double q, int shelf, int onbed, int onsurface, double meanvel,double epsvel,
61 int collapse, double melting[6], double accumulation[6], double geothermalflux[6],
[111]62 int artdiff, int thermal_steadystate,double viscosity_overshoot);
[1]63 ~Penta();
64
65 void Echo();
66 void Marshall(char** pmarshalled_dataset);
67 int MarshallSize();
68 char* GetName();
69 void Demarshall(char** pmarshalled_dataset);
70 int Enum();
71 int GetId();
72 int MyRank();
73 void Configure(void* loads,void* nodes,void* materials);
[246]74 void CreateKMatrix(Mat Kgg,void* inputs,int analysis_type);
75 void CreateKMatrixDiagnosticHoriz( Mat Kgg, void* inputs, int analysis_type);
76 void CreateKMatrixDiagnosticVert( Mat Kgg, void* inputs, int analysis_type);
77 void CreatePVector(Vec pg, void* inputs, int analysis_type);
78 void UpdateFromInputs(void* inputs);
[96]79 void GetDofList(int* doflist,int* pnumberofdofs);
[301]80 void GetDofList1(int* doflist);
[246]81 void* GetMatPar();
[1]82 int GetShelf();
[246]83 void GetNodes(void** nodes);
[1]84 int GetOnBed();
[246]85 void Du(Vec du_g,double* u_g,double* u_g_obs,void* inputs,int analysis_type);
86 void Gradj(Vec grad_g,double* u_g,double* lambda_g,void* inputs,int analysis_type,char* control_type);
87 void GradjDrag(Vec grad_g,double* u_g,double* lambda_g,void* inputs,int analysis_type);
88 void GradjB(Vec grad_g,double* u_g,double* lambda_g,void* inputs,int analysis_type);
89 double Misfit(double* u_g,double* u_g_obs,void* inputs,int analysis_type);
[1]90
91 void GetThicknessList(double* thickness_list);
92 void GetBedList(double* bed_list);
93 Object* copy();
[246]94 void* SpawnTria(int g0, int g1, int g2);
[1]95
[96]96 void GetStrainRate(double* epsilon, double* velocity, double* xyz_list, double* gauss_l1l2l3l4);
97 void GetB(double* pB, double* xyz_list, double* gauss_l1l2l3l4);
98 void GetBPrime(double* B, double* xyz_list, double* gauss_l1l2l3l4);
[128]99 void GetB_vert(double* B, double* xyz_list, double* gauss_l1l2l3l4);
100 void GetBPrime_vert(double* B, double* xyz_list, double* gauss_l1l2l3l4);
[96]101 void GetJacobianDeterminant(double* Jdet, double* xyz_list,double* gauss_l1l2l3l4);
102 void GetNodalFunctionsDerivativesBasic(double* dh1dh2dh3dh4dh5dh6_basic,double* xyz_list, double* gauss_l1l2l3l4);
103 void GetJacobian(double* J, double* xyz_list,double* gauss_l1l2l3l4);
104 void GetNodalFunctionsDerivativesParams(double* dl1dl2dl3dl4dl5dl6,double* gauss_l1l2l3l4);
105 void GetJacobianInvert(double* Jinv, double* xyz_list,double* gauss_l1l2l3l4);
[246]106 void CreatePVectorDiagnosticHoriz( Vec pg, void* inputs,int analysis_type);
107 void CreatePVectorDiagnosticVert( Vec pg, void* inputs,int analysis_type);
[96]108 void GetParameterValue(double* pvalue, double* v_list,double* gauss_l1l2l3l4);
109 void GetParameterDerivativeValue(double* p, double* p_list,double* xyz_list, double* gauss_l1l2l3l4);
110 void GetNodalFunctions(double* l1l2l3l4l5l6, double* gauss_l1l2l3l4);
111 void VelocityExtrude(Vec ug,double* ug_serial);
[301]112 void ComputePressure(Vec p_g);
[1]113
114};
115#endif /* _PENTA_H */
[96]116
117
Note: See TracBrowser for help on using the repository browser.