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

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

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