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

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

Massive temporary commit: redesign to include vertices (for galerkin discontinous), to simplify ModelProcessor.
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"
[3383]11#include "./Numpar.h"
[1]12#include "./Matice.h"
[96]13#include "./Tria.h"
[3383]14#include "./Hook.h"
15#include "./ElementProperties.h"
[246]16#include "./ParameterInputs.h"
[3417]17#include "./Node.h"
[1]18
[3417]19class Object;
[3383]20class Node;
[3417]21class Hook;
22class ElementProperties;
23class DataSet;
[3383]24
[1]25class Penta: public Element{
26
27 private:
[3383]28
[1]29 int id;
[3383]30 Hook hnodes; //hook to 6 nodes
31 Hook hmatice; //hook to 1 matice
32 Hook hmatpar; //hook to 1 matpar
33 Hook hnumpar; //hook to 1 numpar
[1]34
[3383]35 ElementProperties properties;
[2333]36
[1]37 public:
38
[3383]39 /*FUNCTION constructors, destructors {{{1*/
[1]40 Penta();
[3383]41 Penta(int penta_id,int* penta_node_ids, int penta_matice_id, int penta_matpar_id, int penta_numpar_id, ElementProperties* properties);
[3389]42 Penta(int penta_id,Hook* penta_hnodes, Hook* penta_hmatice, Hook* penta_hmatpar, Hook* penta_hnumpar, ElementProperties* penta_properties);
[3417]43 Penta(int i, IoModel* iomodel);
[1]44 ~Penta();
[3383]45 /*}}}*/
46 /*FUNCTION object management {{{1*/
[3389]47 void Configure(void* loads,void* nodes,void* materials,void* parameters);
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);
59 void UpdateFromDakota(void* inputs);
60 void UpdateFromInputs(void* inputs);
[3383]61 /*}}}*/
62 /*FUNCTION element numerical routines {{{1*/
[465]63 void CreateKMatrix(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
64 void CreateKMatrixDiagnosticHoriz( Mat Kgg, void* inputs, int analysis_type,int sub_analysis_type);
65 void CreateKMatrixDiagnosticVert( Mat Kgg, void* inputs, int analysis_type,int sub_analysis_type);
66 void CreatePVector(Vec pg, void* inputs, 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();
[1]70 int GetShelf();
[246]71 void GetNodes(void** nodes);
[1]72 int GetOnBed();
[1185]73 void Du(Vec du_g,void* inputs,int analysis_type,int sub_analysis_type);
[1188]74 void Gradj(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type,char* control_type);
75 void GradjDrag(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type);
76 void GradjB(Vec grad_g,void* inputs,int analysis_type,int sub_analysis_type);
[3169]77 double Misfit(void* inputs,int analysis_type,int sub_analysis_type);
[3180]78 double SurfaceArea(void* inputs,int analysis_type,int sub_analysis_type);
[3169]79 double CostFunction(void* inputs,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);
[465]94 void CreatePVectorDiagnosticHoriz( Vec pg, void* inputs,int analysis_type,int sub_analysis_type);
95 void CreatePVectorDiagnosticVert( Vec pg, void* inputs,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);
[301]100 void ComputePressure(Vec p_g);
[465]101 void CreateKMatrixSlopeCompute(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
102 void CreatePVectorSlopeCompute( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
[607]103 void CreateKMatrixPrognostic(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
104 void CreatePVectorPrognostic( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
[2714]105 void CreateKMatrixBalancedthickness(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
106 void CreatePVectorBalancedthickness( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
[2722]107 void CreateKMatrixBalancedvelocities(Mat Kgg,void* vinputs,int analysis_type,int sub_analysis_type);
108 void CreatePVectorBalancedvelocities( Vec pg, void* vinputs, int analysis_type,int sub_analysis_type);
[1]109
[465]110 void CreateKMatrixDiagnosticStokes( Mat Kgg, void* vinputs, int analysis_type,int sub_analysis_type);
111 void CreatePVectorDiagnosticStokes( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
[394]112 void ReduceMatrixStokes(double* Ke_reduced, double* Ke_temp);
113 void GetMatrixInvert(double* Ke_invert, double* Ke);
114 void SurfaceNormal(double* surface_normal, double xyz_list[3][3]);
115 void GetStrainRateStokes(double* epsilon, double* velocity, double* xyz_list, double* gauss_coord);
116 void GetBStokes(double* B, double* xyz_list, double* gauss_coord);
117 void GetBprimeStokes(double* B_prime, double* xyz_list, double* gauss_coord);
118 void GetLStokes(double* LStokes, double* gauss_coord_tria);
119 void GetLprimeStokes(double* LprimeStokes, double* xyz_list, double* gauss_coord_tria, double* gauss_coord);
[2956]120 void GetNodalFunctionsDerivativesStokes(double* dh1dh7,double* xyz_list, double* gauss_coord);
121 void GetNodalFunctionsDerivativesReferenceStokes(double* dl1dl7,double* gauss_coord);
[394]122 void ReduceVectorStokes(double* Pe_reduced, double* Ke_temp, double* Pe_temp);
123 void GetNodalFunctionsStokes(double* l1l7, double* gauss_coord);
[483]124 void CreateKMatrixThermal(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
125 void GetB_conduct(double* B_conduct, double* xyz_list, double* gauss_coord);
126 void GetB_advec(double* B_advec, double* xyz_list, double* gauss_coord);
127 void GetBprime_advec(double* Bprime_advec, double* xyz_list, double* gauss_coord);
[1676]128 void GetB_artdiff(double* B_artdiff, double* xyz_list, double* gauss_coord);
[483]129 void CreateKMatrixMelting(Mat Kgg,void* inputs,int analysis_type,int sub_analysis_type);
130 void CreatePVectorThermal( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
131 void CreatePVectorMelting( Vec pg, void* vinputs,int analysis_type,int sub_analysis_type);
132 void GetPhi(double* phi, double* epsilon, double viscosity);
[2112]133 double MassFlux(double* segment,double* ug);
[3383]134 /*}}}*/
[358]135
[394]136
[1]137};
138#endif /* _PENTA_H */
[96]139
140
Note: See TracBrowser for help on using the repository browser.