| [3683] | 1 | /*!\file: Element.h
|
|---|
| 2 | * \brief abstract class for Element object
|
|---|
| 3 | * This class is a place holder for the Tria and the Penta elements.
|
|---|
| 4 | * It is derived from Element, so DataSets can contain them.
|
|---|
| 5 | */
|
|---|
| 6 |
|
|---|
| 7 | #ifndef _ELEMENT_H_
|
|---|
| 8 | #define _ELEMENT_H_
|
|---|
| 9 |
|
|---|
| 10 | /*Headers:*/
|
|---|
| [12365] | 11 | /*{{{*/
|
|---|
| [15067] | 12 | #include "../../datastructures/datastructures.h"
|
|---|
| [16789] | 13 | #include "../../toolkits/toolkits.h"
|
|---|
| [15012] | 14 | #include "../Update.h"
|
|---|
| [3683] | 15 | class DataSet;
|
|---|
| 16 | class Parameters;
|
|---|
| [14996] | 17 | class Elements;
|
|---|
| 18 | class Loads;
|
|---|
| 19 | class Node;
|
|---|
| 20 | class Nodes;
|
|---|
| 21 | class Vertices;
|
|---|
| 22 | class Materials;
|
|---|
| [16693] | 23 | class Matpar;
|
|---|
| [16675] | 24 | class Input;
|
|---|
| 25 | class Gauss;
|
|---|
| [16799] | 26 | class ElementVector;
|
|---|
| [13216] | 27 | template <class doublematrix> class Matrix;
|
|---|
| 28 | template <class doubletype> class Vector;
|
|---|
| [16789] | 29 | class ElementMatrix;
|
|---|
| 30 | class ElementVector;
|
|---|
| [3683] | 31 | /*}}}*/
|
|---|
| 32 |
|
|---|
| [4244] | 33 | class Element: public Object,public Update{
|
|---|
| [3683] | 34 |
|
|---|
| 35 | public:
|
|---|
| [13623] | 36 |
|
|---|
| [3683] | 37 | virtual ~Element(){};
|
|---|
| [13623] | 38 |
|
|---|
| [14761] | 39 | virtual void Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters)=0;
|
|---|
| [14953] | 40 | virtual void SetCurrentConfiguration(Elements* elements,Loads* loads,Nodes* nodes,Materials* materials,Parameters* parameters)=0;
|
|---|
| [16042] | 41 | virtual void SetwiseNodeConnectivity(int* d_nz,int* o_nz,Node* node,bool* flags,int* flagsindices,int set1_enum,int set2_enum)=0;
|
|---|
| [16720] | 42 | virtual void AddBasalInput(int input_enum, IssmDouble* values, int interpolation_enum)=0;
|
|---|
| [16693] | 43 | virtual void AddInput(int input_enum, IssmDouble* values, int interpolation_enum)=0;
|
|---|
| [16699] | 44 | virtual void AddMaterialInput(int input_enum, IssmDouble* values, int interpolation_enum)=0;
|
|---|
| [16789] | 45 | virtual ElementMatrix* CreateKMatrix(void)=0;
|
|---|
| [14435] | 46 | virtual void CreateDVector(Vector<IssmDouble>* df)=0;
|
|---|
| [16789] | 47 | virtual ElementVector* CreatePVector(void)=0;
|
|---|
| [13216] | 48 | virtual void CreateJacobianMatrix(Matrix<IssmDouble>* Jff)=0;
|
|---|
| [16778] | 49 | virtual void DeleteMaterials(void)=0;
|
|---|
| [16734] | 50 | virtual void EnthalpyToThermal(IssmDouble* ptemperature,IssmDouble* pwaterfraction,IssmDouble enthalpy,IssmDouble pressure)=0;
|
|---|
| [16812] | 51 | virtual IssmDouble EnthalpyDiffusionParameter(IssmDouble enthalpy,IssmDouble pressure)=0;
|
|---|
| 52 | virtual IssmDouble EnthalpyDiffusionParameterVolume(int numvertices,IssmDouble* enthalpy,IssmDouble* pressure)=0;
|
|---|
| [16675] | 53 | virtual void FindParam(int* pvalue,int paramenum)=0;
|
|---|
| 54 | virtual void FindParam(IssmDouble* pvalue,int paramenum)=0;
|
|---|
| 55 | virtual int FiniteElement(void)=0;
|
|---|
| [16697] | 56 | virtual IssmDouble GetMaterialParameter(int enum_in)=0;
|
|---|
| [16810] | 57 | virtual IssmDouble MinEdgeLength(IssmDouble* xyz_list)=0;
|
|---|
| [16799] | 58 | virtual void NodalFunctions(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| [16810] | 59 | virtual void NodalFunctionsDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| [16805] | 60 | virtual void NodalFunctionsVelocity(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| 61 | virtual void NodalFunctionsPressure(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| [16807] | 62 | virtual void NormalBase(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| [16813] | 63 | virtual IssmDouble PureIceEnthalpy(IssmDouble pressure)=0;
|
|---|
| [16799] | 64 | virtual void TransformLoadVectorCoord(ElementVector* pe,int transformenum)=0;
|
|---|
| 65 | virtual void TransformLoadVectorCoord(ElementVector* pe,int* transformenum_list)=0;
|
|---|
| 66 | virtual void TransformLoadVectorCoord(ElementVector* pe,int numnodes,int transformenum)=0;
|
|---|
| 67 | virtual void TransformLoadVectorCoord(ElementVector* pe,int numnodes,int* transformenum_list)=0;
|
|---|
| [16693] | 68 | virtual void TransformSolutionCoord(IssmDouble* values,int transformenum)=0;
|
|---|
| [16715] | 69 | virtual void TransformSolutionCoord(IssmDouble* values,int* transformenum_list)=0;
|
|---|
| [16775] | 70 | virtual void TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum)=0;
|
|---|
| 71 | virtual void TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list)=0;
|
|---|
| [16762] | 72 | virtual Element* GetBasalElement(void)=0;
|
|---|
| [16675] | 73 | virtual void GetDofList(int** pdoflist,int approximation_enum,int setenum)=0;
|
|---|
| 74 | virtual void GetDofListVelocity(int** pdoflist,int setenum)=0;
|
|---|
| 75 | virtual void GetDofListPressure(int** pdoflist,int setenum)=0;
|
|---|
| [16799] | 76 | virtual void JacobianDeterminant(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| [16803] | 77 | virtual void JacobianDeterminantBase(IssmDouble* Jdet,IssmDouble* xyz_list_base,Gauss* gauss)=0;
|
|---|
| [16675] | 78 | virtual void GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution,int solutionenum)=0;
|
|---|
| [5719] | 79 | virtual int GetNodeIndex(Node* node)=0;
|
|---|
| [13886] | 80 | virtual int GetNumberOfNodes(void)=0;
|
|---|
| [16715] | 81 | virtual int GetNumberOfNodesVelocity(void)=0;
|
|---|
| 82 | virtual int GetNumberOfNodesPressure(void)=0;
|
|---|
| [16727] | 83 | virtual int GetNumberOfVertices(void)=0;
|
|---|
| [13892] | 84 | virtual void GetNodesSidList(int* sidlist)=0;
|
|---|
| [16042] | 85 | virtual void GetNodesLidList(int* sidlist)=0;
|
|---|
| [13892] | 86 |
|
|---|
| [7089] | 87 | virtual int Sid()=0;
|
|---|
| [10143] | 88 | virtual bool IsFloating()=0;
|
|---|
| [12471] | 89 | virtual bool IsNodeOnShelfFromFlags(IssmDouble* flags)=0;
|
|---|
| [5745] | 90 | virtual bool IsOnBed()=0;
|
|---|
| [16693] | 91 | virtual void GetInputListOnNodes(IssmDouble* pvalue,int enumtype)=0;
|
|---|
| 92 | virtual void GetInputListOnNodes(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue)=0;
|
|---|
| [12471] | 93 | virtual void GetInputListOnVertices(IssmDouble* pvalue,int enumtype)=0;
|
|---|
| 94 | virtual void GetInputListOnVertices(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue)=0;
|
|---|
| [16675] | 95 | virtual Input* GetInput(int inputenum)=0;
|
|---|
| [12471] | 96 | virtual void GetInputValue(IssmDouble* pvalue,Node* node,int enumtype)=0;
|
|---|
| [16675] | 97 | virtual void GetInputValue(bool* pvalue,int enum_type)=0;
|
|---|
| 98 | virtual void GetInputValue(int* pvalue,int enum_type)=0;
|
|---|
| 99 | virtual void GetInputValue(IssmDouble* pvalue,int enum_type)=0;
|
|---|
| [16698] | 100 | virtual void GetVerticesCoordinates(IssmDouble** xyz_list)=0;
|
|---|
| [16803] | 101 | virtual void GetVerticesCoordinatesBase(IssmDouble** xyz_list)=0;
|
|---|
| [16188] | 102 | virtual void GetMaterialInputValue(IssmDouble* pvalue,Node* node,int enumtype)=0;
|
|---|
| [13623] | 103 |
|
|---|
| [12471] | 104 | virtual IssmDouble SurfaceArea(void)=0;
|
|---|
| [4931] | 105 | virtual void InputDepthAverageAtBase(int enum_type,int average_enum_type,int object_enum)=0;
|
|---|
| [16693] | 106 | virtual void InputChangeName(int enum_type,int enum_type_old)=0;
|
|---|
| [13216] | 107 | virtual void ComputeBasalStress(Vector<IssmDouble>* sigma_b)=0;
|
|---|
| 108 | virtual void ComputeStrainRate(Vector<IssmDouble>* eps)=0;
|
|---|
| [16461] | 109 | virtual void ResultInterpolation(int* pinterpolation,int output_enum)=0;
|
|---|
| [16589] | 110 | virtual void ResultToVector(Vector<IssmDouble>* vector,int output_enum)=0;
|
|---|
| [15613] | 111 | virtual void Update(int index, IoModel* iomodel,int analysis_counter,int analysis_type,int finite_element)=0;
|
|---|
| [8967] | 112 | virtual void InputDuplicate(int original_enum,int new_enum)=0;
|
|---|
| [16208] | 113 | virtual void InputCreate(IssmDouble* vector,IoModel* iomodel,int M,int N,int vector_type,int vector_enum,int code)=0;
|
|---|
| [16715] | 114 | virtual void InputUpdateFromSolutionOneDofCollapsed(IssmDouble* solution,int inputenum)=0;
|
|---|
| [16717] | 115 | virtual void InputUpdateFromSolutionOneDof(IssmDouble* solution,int inputenum)=0;
|
|---|
| [13623] | 116 |
|
|---|
| [15130] | 117 | virtual int NodalValue(IssmDouble* pvalue, int index, int natureofdataenum)=0;
|
|---|
| [16675] | 118 | virtual int NumberofNodesVelocity(void)=0;
|
|---|
| 119 | virtual int NumberofNodesPressure(void)=0;
|
|---|
| 120 | virtual Gauss* NewGauss(void)=0;
|
|---|
| [16799] | 121 | virtual Gauss* NewGauss(int order)=0;
|
|---|
| [16803] | 122 | virtual Gauss* NewGaussBase(int order)=0;
|
|---|
| 123 | virtual ElementVector* NewElementVector(int approximation_enum=NoneApproximationEnum)=0;
|
|---|
| [12471] | 124 | virtual void InputScale(int enum_type,IssmDouble scale_factor)=0;
|
|---|
| [13216] | 125 | virtual void GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum)=0;
|
|---|
| [12471] | 126 | virtual IssmDouble TimeAdapt()=0;
|
|---|
| 127 | virtual void PositiveDegreeDay(IssmDouble* pdds,IssmDouble* pds,IssmDouble signorm)=0;
|
|---|
| [12748] | 128 | virtual void Delta18oParameterization(void)=0;
|
|---|
| [16721] | 129 | virtual void SmbGradients(void)=0;
|
|---|
| 130 | virtual Element* SpawnBasalElement(void)=0;
|
|---|
| [16699] | 131 | virtual IssmDouble TMeltingPoint(IssmDouble pressure)=0;
|
|---|
| [16789] | 132 | virtual void ReduceMatrices(ElementMatrix* Ke,ElementVector* pe)=0;
|
|---|
| [10440] | 133 | virtual void ResetCoordinateSystem()=0;
|
|---|
| [16810] | 134 | virtual IssmDouble StabilizationParameter(IssmDouble u, IssmDouble v, IssmDouble w, IssmDouble diameter, IssmDouble kappa)=0;
|
|---|
| 135 | virtual void ViscousHeating(IssmDouble* pphi,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input)=0;
|
|---|
| [16675] | 136 | virtual int VelocityInterpolation()=0;
|
|---|
| 137 | virtual int PressureInterpolation()=0;
|
|---|
| [9817] | 138 |
|
|---|
| 139 | #ifdef _HAVE_RESPONSES_
|
|---|
| [15387] | 140 | virtual void AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part)=0;
|
|---|
| [15130] | 141 | virtual void MinVel(IssmDouble* pminvel)=0;
|
|---|
| 142 | virtual void MaxVel(IssmDouble* pmaxvel)=0;
|
|---|
| 143 | virtual void MinVx(IssmDouble* pminvx)=0;
|
|---|
| 144 | virtual void MaxVx(IssmDouble* pmaxvx)=0;
|
|---|
| 145 | virtual void MaxAbsVx(IssmDouble* pmaxabsvx)=0;
|
|---|
| 146 | virtual void MinVy(IssmDouble* pminvy)=0;
|
|---|
| 147 | virtual void MaxVy(IssmDouble* pmaxvy)=0;
|
|---|
| 148 | virtual void MaxAbsVy(IssmDouble* pmaxabsvy)=0;
|
|---|
| 149 | virtual void MinVz(IssmDouble* pminvz)=0;
|
|---|
| 150 | virtual void MaxVz(IssmDouble* pmaxvz)=0;
|
|---|
| 151 | virtual void MaxAbsVz(IssmDouble* pmaxabsvz)=0;
|
|---|
| 152 | virtual IssmDouble MassFlux(IssmDouble* segment)=0;
|
|---|
| [16388] | 153 | virtual IssmDouble MassFlux(IssmDouble x1,IssmDouble y1, IssmDouble x2, IssmDouble y2,int segment_id)=0;
|
|---|
| [15129] | 154 | virtual void ElementResponse(IssmDouble* presponse,int response_enum)=0;
|
|---|
| [12471] | 155 | virtual IssmDouble IceVolume(void)=0;
|
|---|
| [16275] | 156 | virtual IssmDouble IceVolumeAboveFloatation(void)=0;
|
|---|
| [12773] | 157 | virtual IssmDouble TotalSmb(void)=0;
|
|---|
| [16787] | 158 | virtual IssmDouble Misfit(int modelenum,int observationenum,int weightsenum)=0;
|
|---|
| [9817] | 159 | #endif
|
|---|
| [4055] | 160 |
|
|---|
| [14589] | 161 | #ifdef _HAVE_GIA_
|
|---|
| [14807] | 162 | virtual void GiaDeflection(Vector<IssmDouble>* wg,Vector<IssmDouble>* dwgdt,IssmDouble* x,IssmDouble* y)=0;
|
|---|
| [14589] | 163 | #endif
|
|---|
| 164 |
|
|---|
| [9817] | 165 | #ifdef _HAVE_CONTROL_
|
|---|
| [13216] | 166 | virtual void Gradj(Vector<IssmDouble>* gradient,int control_type,int control_index)=0;
|
|---|
| [16314] | 167 | virtual IssmDouble ThicknessAbsMisfit(void)=0;
|
|---|
| 168 | virtual IssmDouble SurfaceAbsVelMisfit(void)=0;
|
|---|
| 169 | virtual IssmDouble SurfaceRelVelMisfit(void)=0;
|
|---|
| 170 | virtual IssmDouble SurfaceLogVelMisfit(void)=0;
|
|---|
| 171 | virtual IssmDouble SurfaceLogVxVyMisfit(void)=0;
|
|---|
| 172 | virtual IssmDouble SurfaceAverageVelMisfit(void)=0;
|
|---|
| 173 | virtual IssmDouble ThicknessAbsGradient(void)=0;
|
|---|
| 174 | virtual IssmDouble ThicknessAlongGradient(void)=0;
|
|---|
| 175 | virtual IssmDouble ThicknessAcrossGradient(void)=0;
|
|---|
| 176 | virtual IssmDouble BalancethicknessMisfit(void)=0;
|
|---|
| 177 | virtual IssmDouble RheologyBbarAbsGradient(void)=0;
|
|---|
| 178 | virtual IssmDouble DragCoefficientAbsGradient(void)=0;
|
|---|
| [13216] | 179 | virtual void ControlInputGetGradient(Vector<IssmDouble>* gradient,int enum_type,int control_index)=0;
|
|---|
| [12471] | 180 | virtual void ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index)=0;
|
|---|
| 181 | virtual void ControlInputScaleGradient(int enum_type, IssmDouble scale)=0;
|
|---|
| [16494] | 182 | virtual void ControlToVectors(Vector<IssmPDouble>* vector_control, Vector<IssmPDouble>* vector_gradient,int control_enum)=0;
|
|---|
| [13216] | 183 | virtual void GetVectorFromControlInputs(Vector<IssmDouble>* gradient,int control_enum,int control_index,const char* data)=0;
|
|---|
| [12471] | 184 | virtual void SetControlInputsFromVector(IssmDouble* vector,int control_enum,int control_index)=0;
|
|---|
| 185 | virtual void InputControlUpdate(IssmDouble scalar,bool save_parameter)=0;
|
|---|
| [9817] | 186 | #endif
|
|---|
| [14796] | 187 |
|
|---|
| [16442] | 188 | virtual void UpdateConstraintsExtrudeFromBase(void)=0;
|
|---|
| 189 | virtual void UpdateConstraintsExtrudeFromTop(void)=0;
|
|---|
| 190 |
|
|---|
| [15853] | 191 | #ifdef _HAVE_THERMAL_
|
|---|
| [16272] | 192 | virtual void UpdateBasalConstraintsEnthalpy(void)=0;
|
|---|
| [16026] | 193 | virtual void ComputeBasalMeltingrate(void)=0;
|
|---|
| [16816] | 194 | virtual void DrainWaterfraction(IssmDouble* drainrate_element)=0;
|
|---|
| [15853] | 195 | #endif
|
|---|
| 196 |
|
|---|
| [14796] | 197 | #ifdef _HAVE_HYDROLOGY_
|
|---|
| [14967] | 198 | virtual void GetHydrologyDCInefficientHmax(IssmDouble* ph_max, Node* innode)=0;
|
|---|
| [16761] | 199 | virtual void GetHydrologyDCInefficientHmax(IssmDouble* ph_max, int index)=0;
|
|---|
| [15015] | 200 | virtual void GetHydrologyTransfer(Vector<IssmDouble>* transfer)=0;
|
|---|
| [15249] | 201 | virtual void HydrologyEPLGetMask(Vector<IssmDouble>* mask)=0;
|
|---|
| 202 | virtual void HydrologyEPLGetActive(Vector<IssmDouble>* active)=0;
|
|---|
| [16600] | 203 | virtual void ComputeEPLThickness(void)=0;
|
|---|
| [14796] | 204 | #endif
|
|---|
| [15643] | 205 |
|
|---|
| [15386] | 206 | #ifdef _HAVE_GROUNDINGLINE_
|
|---|
| [15962] | 207 | virtual void MigrateGroundingLine(IssmDouble* sheet_ungrounding)=0;
|
|---|
| [15386] | 208 | virtual void PotentialUngrounding(Vector<IssmDouble>* potential_sheet_ungrounding)=0;
|
|---|
| 209 | virtual int UpdatePotentialUngrounding(IssmDouble* potential_sheet_ungrounding,Vector<IssmDouble>* vec_nodes_on_iceshelf,IssmDouble* nodes_on_iceshelf)=0;
|
|---|
| 210 | #endif
|
|---|
| 211 |
|
|---|
| [3683] | 212 | };
|
|---|
| 213 | #endif
|
|---|
![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
