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

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

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