source: issm/trunk/src/c/classes/Elements/Tria.h@ 22758

Last change on this file since 22758 was 22758, checked in by Mathieu Morlighem, 7 years ago

merged trunk-jpl and trunk for revision 22757
File size: 13.0 KB
Line 
1/*! \file Tria.h
2 * \brief: header file for tria object
3 */
4
5#ifndef _TRIA_H_
6#define _TRIA_H_
7
8/*Headers:*/
9/*{{{*/
10#include "./Element.h"
11#include "./ElementHook.h"
12#include "./TriaRef.h"
13class Parameters;
14class Inputs;
15class IoModel;
16class Results;
17class Node;
18class Material;
19class Matpar;
20class Seg;
21class ElementMatrix;
22class ElementVector;
23class Vertex;
24class GaussTria;
25
26#include "../../shared/Exceptions/exceptions.h"
27#include "../../shared/Enum/Enum.h"
28/*}}}*/
29
30class Tria: public Element,public ElementHook,public TriaRef{
31
32 public:
33
34 /*Tria constructors, destructors {{{*/
35 Tria(){};
36 Tria(int tria_id,int tria_sid,IoModel* iomodel,int nummodels);
37 ~Tria();
38 /*}}}*/
39 /*Object virtual functions definitions:{{{ */
40 Object *copy();
41 void Marshall(char** pmarshalled_data,int* pmarshalled_data_size, int marshall_direction);
42 int ObjectEnum();
43 /*}}}*/
44 /*Update virtual functions resolution: {{{*/
45 #ifdef _HAVE_DAKOTA_
46 void InputUpdateFromMatrixDakota(IssmDouble* matrix, int nows, int ncols, int name, int type);
47 void InputUpdateFromVectorDakota(IssmDouble* vector, int name, int type);
48 #endif
49 void InputUpdateFromIoModel(int index, IoModel* iomodel);
50 void InputUpdateFromVector(IssmDouble* vector, int name, int type);
51 /*}}}*/
52 /*Element virtual functions definitions: {{{*/
53 void AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part);
54 void CalvingRateVonmises();
55 void CalvingCrevasseDepth();
56 void CalvingRateLevermann();
57 IssmDouble CharacteristicLength(void);
58 void ComputeBasalStress(Vector<IssmDouble>* sigma_b);
59 void ComputeDeviatoricStressTensor();
60 void ComputeEsaStrainAndVorticity();
61 void ComputeSigmaNN();
62 void ComputeStressTensor();
63 void ComputeSurfaceNormalVelocity();
64 void Configure(Elements* elements,Loads* loads,Nodes* nodesin,Vertices* verticesin,Materials* materials,Parameters* parameters);
65 void ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int N, int M);
66 void ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index);
67 void ControlToVectors(Vector<IssmPDouble>* vector_control, Vector<IssmPDouble>* vector_gradient,int control_enum);
68 void CreateDistanceInputFromSegmentlist(IssmDouble* segments,int numsegments,int distanceenum);
69 int EdgeOnBaseIndex();
70 void EdgeOnBaseIndices(int* pindex1,int* pindex);
71 int EdgeOnSurfaceIndex();
72 void EdgeOnSurfaceIndices(int* pindex1,int* pindex);
73 void ElementResponse(IssmDouble* presponse,int response_enum);
74 void ElementSizes(IssmDouble* hx,IssmDouble* hy,IssmDouble* hz);
75 int FiniteElement(void);
76 IssmDouble FloatingArea(bool scaled);
77 void FSContactMigration(Vector<IssmDouble>* vertexgrounded,Vector<IssmDouble>* vertexfloating);
78 Element* GetBasalElement(void){_error_("not implemented yet");};
79 void GetLevelsetPositivePart(int* point1,IssmDouble* fraction1,IssmDouble* fraction2, bool* mainlynegative,IssmDouble* levelsetvalues);
80 void GetGroundedPart(int* point1,IssmDouble* fraction1, IssmDouble* fraction2,bool* mainlyfloating);
81 IssmDouble GetGroundedPortion(IssmDouble* xyz_list);
82 void GetIcefrontCoordinates(IssmDouble** pxyz_front,IssmDouble* xyz_list,int levelsetenum);
83 void GetLevelCoordinates(IssmDouble** pxyz_front,IssmDouble* xyz_list,int levelsetenum,IssmDouble level);
84 int GetNodeIndex(Node* node);
85 int GetNumberOfNodes(void);
86 int GetNumberOfNodes(int enum_type);
87 int GetNumberOfVertices(void);
88 void GetVectorFromControlInputs(Vector<IssmDouble>* gradient,int control_enum,int control_index,const char* data,int offset,bool onsid);
89 void GetVectorFromControlInputs(Vector<IssmDouble>* gradient,int control_enum,int control_index,const char* data,bool onsid);
90 void GetVerticesCoordinatesBase(IssmDouble** pxyz_list);
91 void GetVerticesCoordinatesTop(IssmDouble** pxyz_list);
92 IssmDouble GroundedArea(bool scaled);
93 bool HasEdgeOnBase();
94 bool HasEdgeOnSurface();
95 IssmDouble IceVolume(bool scaled);
96 IssmDouble IceVolumeAboveFloatation(bool scaled);
97 void InputDepthAverageAtBase(int enum_type,int average_enum_type);
98 void InputExtrude(int enum_type,int start){_error_("not implemented"); /*For penta only*/};
99 bool IsFaceOnBoundary(void);
100 bool IsIcefront(void);
101 bool IsNodeOnShelfFromFlags(IssmDouble* flags);
102 bool IsOnBase();
103 bool IsOnSurface();
104 bool IsZeroLevelset(int levelset_enum);
105 IssmDouble Masscon(IssmDouble* levelset);
106 IssmDouble MassFlux(IssmDouble* segment);
107 IssmDouble MassFlux(IssmDouble x1,IssmDouble y1, IssmDouble x2, IssmDouble y2,int segment_id);
108 void MaterialUpdateFromTemperature(void){_error_("not implemented yet");};
109 IssmDouble Misfit(int modelenum,int observationenum,int weightsenum);
110 IssmDouble MisfitArea(int weightsenum);
111 int NodalValue(IssmDouble* pvalue, int index, int natureofdataenum);
112 int NumberofNodesPressure(void);
113 int NumberofNodesVelocity(void);
114 void PicoUpdateBoxid(int* pmax_boxid_basin);
115 void PicoUpdateBox(int loopboxid);
116 void PotentialUngrounding(Vector<IssmDouble>* potential_sheet_ungrounding);
117 int PressureInterpolation();
118 void ReduceMatrices(ElementMatrix* Ke,ElementVector* pe);
119 void ResetFSBasalBoundaryCondition(void);
120 void ResetHooks();
121 void ResetLevelsetFromSegmentlist(IssmDouble* segments,int numsegments);
122 void SetControlInputsFromVector(IssmDouble* vector,int control_enum,int control_index,int offset,int N,int M);
123 void SetControlInputsFromVector(IssmDouble* vector,int control_enum,int control_index);
124 void SetCurrentConfiguration(Elements* elements,Loads* loads,Nodes* nodes,Materials* materials,Parameters* parameters);
125 Element* SpawnBasalElement(void);
126 Element* SpawnTopElement(void);
127 void StrainRateparallel();
128 void StrainRateperpendicular();
129 void StressIntensityFactor(void){_error_("not implemented yet");};
130 IssmDouble SurfaceArea(void);
131 int TensorInterpolation();
132 IssmDouble TimeAdapt();
133 IssmDouble TotalFloatingBmb(bool scaled);
134 IssmDouble TotalGroundedBmb(bool scaled);
135 IssmDouble TotalSmb(bool scaled);
136 void Update(int index, IoModel* iomodel,int analysis_counter,int analysis_type,int finitelement);
137 int UpdatePotentialUngrounding(IssmDouble* vertices_potentially_ungrounding,Vector<IssmDouble>* vec_nodes_on_iceshelf,IssmDouble* nodes_on_iceshelf);
138 void ValueP1DerivativesOnGauss(IssmDouble* dvalue,IssmDouble* values,IssmDouble* xyz_list,Gauss* gauss);
139 void ValueP1OnGauss(IssmDouble* pvalue,IssmDouble* values,Gauss* gauss);
140 int VelocityInterpolation();
141 int VertexConnectivity(int vertexindex);
142 void VerticalSegmentIndices(int** pindices,int* pnumseg){_error_("not implemented yet");};
143 void VerticalSegmentIndicesBase(int** pindices,int* pnumseg){_error_("not implemented yet");};
144 void WriteFieldIsovalueSegment(DataSet* segments,int fieldenum,IssmDouble fieldvalue);
145 void WriteLevelsetSegment(DataSet* segments);
146 void ZeroLevelsetCoordinates(IssmDouble** pxyz_zero,IssmDouble* xyz_list,int levelsetenum);
147
148 #ifdef _HAVE_GIAIVINS_
149 void GiaDeflection(Vector<IssmDouble>* wg,Vector<IssmDouble>* dwgdt,IssmDouble* x,IssmDouble* y);
150 #endif
151 #ifdef _HAVE_ESA_
152 void EsaGeodetic2D(Vector<IssmDouble>* pUp,Vector<IssmDouble>* pNorth,Vector<IssmDouble>* pEast, Vector<IssmDouble>* pX,Vector<IssmDouble>* pY,IssmDouble* xx,IssmDouble* yy);
153 void EsaGeodetic3D(Vector<IssmDouble>* pUp,Vector<IssmDouble>* pNorth,Vector<IssmDouble>* pEast,IssmDouble* latitude,IssmDouble* longitude,IssmDouble* radius,IssmDouble* xx,IssmDouble* yy,IssmDouble* zz,IssmDouble eartharea);
154 #endif
155 #ifdef _HAVE_SEALEVELRISE_
156 IssmDouble OceanArea(void);
157 IssmDouble OceanAverage(IssmDouble* Sg);
158 void SealevelriseMomentOfInertia(IssmDouble* dI_list,IssmDouble* Sg_old,IssmDouble eartharea);
159 void SealevelriseEustatic(Vector<IssmDouble>* pSgi,IssmDouble* peustatic,IssmDouble* latitude,IssmDouble* longitude,IssmDouble* radius,IssmDouble oceanarea,IssmDouble eartharea);
160 void SealevelriseNonEustatic(Vector<IssmDouble>* pSgo,IssmDouble* Sg_old,IssmDouble* latitude,IssmDouble* longitude,IssmDouble* radius,IssmDouble eartharea);
161 void SealevelriseGeodetic(Vector<IssmDouble>* pUp,Vector<IssmDouble>* pNorth,Vector<IssmDouble>* pEast,IssmDouble* Sg,IssmDouble* latitude,IssmDouble* longitude,IssmDouble* radius,IssmDouble* xx,IssmDouble* yy,IssmDouble* zz,IssmDouble eartharea);
162 #endif
163 /*}}}*/
164 /*Tria specific routines:{{{*/
165 void AddBasalInput(int input_enum, IssmDouble* values, int interpolation_enum);
166 void AddInput(int input_enum, IssmDouble* values, int interpolation_enum);
167 void AddControlInput(int input_enum, IssmDouble* values,IssmDouble* values_min,IssmDouble* values_max, int interpolation_enum,int id);
168 void DatasetInputCreate(IssmDouble* array,int M,int N,int* individual_enums,int num_inputs,IoModel* iomodel,int input_enum);
169 IssmDouble GetArea(void);
170 IssmDouble GetHorizontalSurfaceArea(void);
171 IssmDouble GetArea3D(void);
172 IssmDouble GetAreaIce(void);
173 IssmDouble GetAreaSpherical(void);
174 void GetAreaCoordinates(IssmDouble *area_coordinates,IssmDouble* xyz_zero,IssmDouble* xyz_list,int numpoints);
175 int GetElementType(void);
176 void GetInputValue(IssmDouble* pvalue,Node* node,int enumtype);
177 void GetLevelsetIntersection(int** pindices, int* pnumiceverts, IssmDouble* fraction, int levelset_enum, IssmDouble level);
178 void GetMaterialInputValue(IssmDouble* pvalue,Node* node,int enumtype);
179 Node* GetNode(int node_number);
180 void InputUpdateFromSolutionOneDof(IssmDouble* solution,int enum_type);
181 void InputUpdateFromSolutionOneDofCollapsed(IssmDouble* solution,int enum_type){_error_("not implemented yet");};
182 void JacobianDeterminant(IssmDouble* pJdet, IssmDouble* xyz_list,Gauss* gauss);
183 void JacobianDeterminantBase(IssmDouble* pJdet,IssmDouble* xyz_list_base,Gauss* gauss);
184 void JacobianDeterminantLine(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss){_error_("not implemented yet");};
185 void JacobianDeterminantSurface(IssmDouble* pJdet, IssmDouble* xyz_list,Gauss* gauss);
186 void JacobianDeterminantTop(IssmDouble* pJdet,IssmDouble* xyz_list_base,Gauss* gauss);
187 IssmDouble MinEdgeLength(IssmDouble* xyz_list){_error_("not implemented yet");};
188 Gauss* NewGauss(void);
189 Gauss* NewGauss(int order);
190 Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order);
191 Gauss* NewGauss(int point1,IssmDouble fraction1,IssmDouble fraction2,bool mainlyfloating,int order);
192 Gauss* NewGauss(int point1,IssmDouble fraction1,IssmDouble fraction2,int order);
193 Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order_horiz,int order_vert);
194 Gauss* NewGaussBase(int order);
195 Gauss* NewGaussLine(int vertex1,int vertex2,int order){_error_("not implemented yet");};
196 Gauss* NewGaussTop(int order);
197 void NodalFunctions(IssmDouble* basis,Gauss* gauss);
198 void NodalFunctionsDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss);
199 void NodalFunctionsDerivativesVelocity(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss);
200 void NodalFunctionsMINIDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss){_error_("not implemented yet");};
201 void NodalFunctionsPressure(IssmDouble* basis,Gauss* gauss);
202 void NodalFunctionsP1(IssmDouble* basis,Gauss* gauss);
203 void NodalFunctionsP1Derivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss);
204 void NodalFunctionsP2(IssmDouble* basis,Gauss* gauss);
205 void NodalFunctionsTensor(IssmDouble* basis,Gauss* gauss);
206 void NodalFunctionsVelocity(IssmDouble* basis,Gauss* gauss);
207 void NormalBase(IssmDouble* normal,IssmDouble* xyz_list);
208 void NormalSection(IssmDouble* normal,IssmDouble* xyz_list);
209 void NormalTop(IssmDouble* normal,IssmDouble* xyz_list);
210 void SetClone(int* minranks);
211 void SetTemporaryElementType(int element_type_in){_error_("not implemented yet");};
212 Seg* SpawnSeg(int index1,int index2);
213 IssmDouble StabilizationParameter(IssmDouble u, IssmDouble v, IssmDouble w, IssmDouble diameter, IssmDouble kappa){_error_("not implemented yet");};
214 void UpdateConstraintsExtrudeFromBase(void);
215 void UpdateConstraintsExtrudeFromTop(void);
216 void ViscousHeating(IssmDouble* pphi,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input){_error_("not implemented yet");};
217 /*}}}*/
218
219};
220#endif /* _TRIA_H */
Note: See TracBrowser for help on using the repository browser.