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

Last change on this file since 3582 was 3582, checked in by Mathieu Morlighem, 15 years ago

Added discontinuous Galerkin for Balanced thickness solution
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
8class Object;
9class Node;
10class Hook;
11class ElementProperties;
12class DataSet;
13struct IoModel;
14
15
16#include "./Object.h"
17#include "./Element.h"
18#include "./Matpar.h"
19#include "./Numpar.h"
20#include "./Matice.h"
21#include "./Tria.h"
22#include "./Hook.h"
23#include "./ElementProperties.h"
24#include "../ModelProcessorx/IoModel.h"
25#include "./ParameterInputs.h"
26#include "./Node.h"
27
28class Penta: public Element{
29
30 private:
31
32 int id;
33 Hook hnodes; //hook to 6 nodes
34 Hook hmatice; //hook to 1 matice
35 Hook hmatpar; //hook to 1 matpar
36 Hook hnumpar; //hook to 1 numpar
37
38 ElementProperties properties;
39
40 public:
41
42 /*FUNCTION constructors, destructors {{{1*/
43 Penta();
44 Penta(int penta_id,int* penta_node_ids, int penta_matice_id, int penta_matpar_id, int penta_numpar_id, ElementProperties* properties);
45 Penta(int penta_id,Hook* penta_hnodes, Hook* penta_hmatice, Hook* penta_hmatpar, Hook* penta_hnumpar, ElementProperties* penta_properties);
46 Penta(int i, IoModel* iomodel);
47 ~Penta();
48 /*}}}*/
49 /*FUNCTION object management {{{1*/
50 void Configure(void* loads,void* nodes,void* materials,void* parameters);
51 Object* copy();
52 void DeepEcho();
53 void Demarshall(char** pmarshalled_dataset);
54 void Echo();
55 int Enum();
56 int GetId();
57 char* GetName();
58 void Marshall(char** pmarshalled_dataset);
59 int MarshallSize();
60 int MyRank();
61 void* SpawnTria(int g0, int g1, int g2);
62 void UpdateFromDakota(void* inputs);
63 void UpdateFromInputs(void* inputs);
64 void SetClone(int* minranks);
65
66 /*}}}*/
67 /*FUNCTION element numerical routines {{{1*/
68 void CreateKMatrix(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
69 void CreateKMatrixDiagnosticHoriz( Mat Kgg, void* inputs, int analysis_type,int sub_analysis_type);
70 void CreateKMatrixDiagnosticVert( Mat Kgg, void* inputs, int analysis_type,int sub_analysis_type);
71 void CreatePVector(Vec pg, void* inputs, int analysis_type,int sub_analysis_type);
72 void GetDofList(int* doflist,int* pnumberofdofs);
73 void GetDofList1(int* doflist);
74 void* GetMatPar();
75 int GetShelf();
76 void GetNodes(void** nodes);
77 int GetOnBed();
78 void Du(Vec du_g,void* inputs,int analysis_type,int sub_analysis_type);
79 void Gradj(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type,char* control_type);
80 void GradjDrag(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type);
81 void GradjB(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type);
82 double Misfit(void* inputs,int analysis_type,int sub_analysis_type);
83 double SurfaceArea(void* inputs,int analysis_type,int sub_analysis_type);
84 double CostFunction(void* inputs,int analysis_type,int sub_analysis_type);
85
86 void GetThicknessList(double* thickness_list);
87 void GetBedList(double* bed_list);
88
89 void GetStrainRate(double* epsilon, double* velocity, double* xyz_list, double* gauss_coord);
90 void GetB(double* pB, double* xyz_list, double* gauss_coord);
91 void GetBPrime(double* B, double* xyz_list, double* gauss_coord);
92 void GetB_vert(double* B, double* xyz_list, double* gauss_coord);
93 void GetBPrime_vert(double* B, double* xyz_list, double* gauss_coord);
94 void GetJacobianDeterminant(double* Jdet, double* xyz_list,double* gauss_coord);
95 void GetNodalFunctionsDerivatives(double* dh1dh6,double* xyz_list, double* gauss_coord);
96 void GetJacobian(double* J, double* xyz_list,double* gauss_coord);
97 void GetNodalFunctionsDerivativesReference(double* dl1dl6,double* gauss_coord);
98 void GetJacobianInvert(double* Jinv, double* xyz_list,double* gauss_coord);
99 void CreatePVectorDiagnosticHoriz( Vec pg, void* inputs,int analysis_type,int sub_analysis_type);
100 void CreatePVectorDiagnosticVert( Vec pg, void* inputs,int analysis_type,int sub_analysis_type);
101 void GetParameterValue(double* pvalue, double* v_list,double* gauss_coord);
102 void GetParameterDerivativeValue(double* p, double* p_list,double* xyz_list, double* gauss_coord);
103 void GetNodalFunctions(double* l1l6, double* gauss_coord);
104 void FieldExtrude(Vec field,double* field_serial,char* field_name, int iscollapsed);
105 void ComputeBasalStress(Vec sigma_b,void* inputs,int analysis_type,int sub_analysis_type);
106 void ComputePressure(Vec p_g,void* inputs,int analysis_type,int sub_analysis_type);
107 void ComputeStrainRate(Vec eps,void* inputs,int analysis_type,int sub_analysis_type);
108 void CreateKMatrixSlopeCompute(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
109 void CreatePVectorSlopeCompute( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
110 void CreateKMatrixPrognostic(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
111 void CreatePVectorPrognostic( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
112
113 void CreateKMatrixDiagnosticStokes( Mat Kgg, void* vinputs, int analysis_type,int sub_analysis_type);
114 void CreatePVectorDiagnosticStokes( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
115 void ReduceMatrixStokes(double* Ke_reduced, double* Ke_temp);
116 void GetMatrixInvert(double* Ke_invert, double* Ke);
117 void SurfaceNormal(double* surface_normal, double xyz_list[3][3]);
118 void GetStrainRateStokes(double* epsilon, double* velocity, double* xyz_list, double* gauss_coord);
119 void GetBStokes(double* B, double* xyz_list, double* gauss_coord);
120 void GetBprimeStokes(double* B_prime, double* xyz_list, double* gauss_coord);
121 void GetLStokes(double* LStokes, double* gauss_coord_tria);
122 void GetLprimeStokes(double* LprimeStokes, double* xyz_list, double* gauss_coord_tria, double* gauss_coord);
123 void GetNodalFunctionsDerivativesStokes(double* dh1dh7,double* xyz_list, double* gauss_coord);
124 void GetNodalFunctionsDerivativesReferenceStokes(double* dl1dl7,double* gauss_coord);
125 void ReduceVectorStokes(double* Pe_reduced, double* Ke_temp, double* Pe_temp);
126 void GetNodalFunctionsStokes(double* l1l7, double* gauss_coord);
127 void CreateKMatrixThermal(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
128 void GetB_conduct(double* B_conduct, double* xyz_list, double* gauss_coord);
129 void GetB_advec(double* B_advec, double* xyz_list, double* gauss_coord);
130 void GetBprime_advec(double* Bprime_advec, double* xyz_list, double* gauss_coord);
131 void GetB_artdiff(double* B_artdiff, double* xyz_list, double* gauss_coord);
132 void CreateKMatrixMelting(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
133 void CreatePVectorThermal( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
134 void CreatePVectorMelting( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
135 void GetPhi(double* phi, double* epsilon, double viscosity);
136 double MassFlux(double* segment,double* ug);
137 /*}}}*/
138
139
140};
141#endif /* _PENTA_H */
Note: See TracBrowser for help on using the repository browser.