| [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"
|
|---|
| [15012] | 13 | #include "../Update.h"
|
|---|
| [3683] | 14 |
|
|---|
| 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;
|
|---|
| [13216] | 26 | template <class doublematrix> class Matrix;
|
|---|
| 27 | template <class doubletype> class Vector;
|
|---|
| [3683] | 28 |
|
|---|
| [15012] | 29 | #include "../../toolkits/toolkits.h"
|
|---|
| [3683] | 30 | /*}}}*/
|
|---|
| 31 |
|
|---|
| [4244] | 32 | class Element: public Object,public Update{
|
|---|
| [3683] | 33 |
|
|---|
| 34 | public:
|
|---|
| [13623] | 35 |
|
|---|
| [3683] | 36 | virtual ~Element(){};
|
|---|
| [13623] | 37 |
|
|---|
| [14761] | 38 | virtual void Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters)=0;
|
|---|
| [14953] | 39 | virtual void SetCurrentConfiguration(Elements* elements,Loads* loads,Nodes* nodes,Materials* materials,Parameters* parameters)=0;
|
|---|
| [16042] | 40 | virtual void SetwiseNodeConnectivity(int* d_nz,int* o_nz,Node* node,bool* flags,int* flagsindices,int set1_enum,int set2_enum)=0;
|
|---|
| [16693] | 41 | virtual void AddInput(int input_enum, IssmDouble* values, int interpolation_enum)=0;
|
|---|
| [14435] | 42 | virtual void CreateKMatrix(Matrix<IssmDouble>* Kff, Matrix<IssmDouble>* Kfs)=0;
|
|---|
| 43 | virtual void CreateDVector(Vector<IssmDouble>* df)=0;
|
|---|
| [13216] | 44 | virtual void CreatePVector(Vector<IssmDouble>* pf)=0;
|
|---|
| 45 | virtual void CreateJacobianMatrix(Matrix<IssmDouble>* Jff)=0;
|
|---|
| [16675] | 46 | virtual void FindParam(int* pvalue,int paramenum)=0;
|
|---|
| 47 | virtual void FindParam(IssmDouble* pvalue,int paramenum)=0;
|
|---|
| 48 | virtual int FiniteElement(void)=0;
|
|---|
| [16693] | 49 | virtual Matpar* GetMatparPointer(void)=0;
|
|---|
| 50 | virtual void TransformSolutionCoord(IssmDouble* values,int transformenum)=0;
|
|---|
| [16675] | 51 | virtual void GetDofList(int** pdoflist,int approximation_enum,int setenum)=0;
|
|---|
| 52 | virtual void GetDofListVelocity(int** pdoflist,int setenum)=0;
|
|---|
| 53 | virtual void GetDofListPressure(int** pdoflist,int setenum)=0;
|
|---|
| 54 | virtual void GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution,int solutionenum)=0;
|
|---|
| [5719] | 55 | virtual int GetNodeIndex(Node* node)=0;
|
|---|
| [13886] | 56 | virtual int GetNumberOfNodes(void)=0;
|
|---|
| [13892] | 57 | virtual void GetNodesSidList(int* sidlist)=0;
|
|---|
| [16042] | 58 | virtual void GetNodesLidList(int* sidlist)=0;
|
|---|
| [13892] | 59 |
|
|---|
| [7089] | 60 | virtual int Sid()=0;
|
|---|
| [10143] | 61 | virtual bool IsFloating()=0;
|
|---|
| [12471] | 62 | virtual bool IsNodeOnShelfFromFlags(IssmDouble* flags)=0;
|
|---|
| [5745] | 63 | virtual bool IsOnBed()=0;
|
|---|
| [16693] | 64 | virtual void GetInputListOnNodes(IssmDouble* pvalue,int enumtype)=0;
|
|---|
| 65 | virtual void GetInputListOnNodes(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue)=0;
|
|---|
| [12471] | 66 | virtual void GetInputListOnVertices(IssmDouble* pvalue,int enumtype)=0;
|
|---|
| 67 | virtual void GetInputListOnVertices(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue)=0;
|
|---|
| [16675] | 68 | virtual Input* GetInput(int inputenum)=0;
|
|---|
| [12471] | 69 | virtual void GetInputValue(IssmDouble* pvalue,Node* node,int enumtype)=0;
|
|---|
| [16675] | 70 | virtual void GetInputValue(bool* pvalue,int enum_type)=0;
|
|---|
| 71 | virtual void GetInputValue(int* pvalue,int enum_type)=0;
|
|---|
| 72 | virtual void GetInputValue(IssmDouble* pvalue,int enum_type)=0;
|
|---|
| [16188] | 73 | virtual void GetMaterialInputValue(IssmDouble* pvalue,Node* node,int enumtype)=0;
|
|---|
| [13623] | 74 |
|
|---|
| [12471] | 75 | virtual IssmDouble SurfaceArea(void)=0;
|
|---|
| [4931] | 76 | virtual void InputDepthAverageAtBase(int enum_type,int average_enum_type,int object_enum)=0;
|
|---|
| [16693] | 77 | virtual void InputChangeName(int enum_type,int enum_type_old)=0;
|
|---|
| [13216] | 78 | virtual void ComputeBasalStress(Vector<IssmDouble>* sigma_b)=0;
|
|---|
| 79 | virtual void ComputeStrainRate(Vector<IssmDouble>* eps)=0;
|
|---|
| [16461] | 80 | virtual void ResultInterpolation(int* pinterpolation,int output_enum)=0;
|
|---|
| [16589] | 81 | virtual void ResultToVector(Vector<IssmDouble>* vector,int output_enum)=0;
|
|---|
| [15613] | 82 | virtual void Update(int index, IoModel* iomodel,int analysis_counter,int analysis_type,int finite_element)=0;
|
|---|
| [8967] | 83 | virtual void InputDuplicate(int original_enum,int new_enum)=0;
|
|---|
| [16208] | 84 | virtual void InputCreate(IssmDouble* vector,IoModel* iomodel,int M,int N,int vector_type,int vector_enum,int code)=0;
|
|---|
| [13623] | 85 |
|
|---|
| [15130] | 86 | virtual int NodalValue(IssmDouble* pvalue, int index, int natureofdataenum)=0;
|
|---|
| [16675] | 87 | virtual int NumberofNodesVelocity(void)=0;
|
|---|
| 88 | virtual int NumberofNodesPressure(void)=0;
|
|---|
| 89 | virtual Gauss* NewGauss(void)=0;
|
|---|
| [12471] | 90 | virtual void InputScale(int enum_type,IssmDouble scale_factor)=0;
|
|---|
| [13216] | 91 | virtual void GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum)=0;
|
|---|
| [12471] | 92 | virtual IssmDouble TimeAdapt()=0;
|
|---|
| 93 | virtual void PositiveDegreeDay(IssmDouble* pdds,IssmDouble* pds,IssmDouble signorm)=0;
|
|---|
| [12748] | 94 | virtual void Delta18oParameterization(void)=0;
|
|---|
| [12677] | 95 | virtual void SmbGradients()=0;
|
|---|
| [10440] | 96 | virtual void ResetCoordinateSystem()=0;
|
|---|
| [16675] | 97 | virtual int VelocityInterpolation()=0;
|
|---|
| 98 | virtual int PressureInterpolation()=0;
|
|---|
| [9817] | 99 |
|
|---|
| 100 | #ifdef _HAVE_RESPONSES_
|
|---|
| [15387] | 101 | virtual void AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part)=0;
|
|---|
| [15130] | 102 | virtual void MinVel(IssmDouble* pminvel)=0;
|
|---|
| 103 | virtual void MaxVel(IssmDouble* pmaxvel)=0;
|
|---|
| 104 | virtual void MinVx(IssmDouble* pminvx)=0;
|
|---|
| 105 | virtual void MaxVx(IssmDouble* pmaxvx)=0;
|
|---|
| 106 | virtual void MaxAbsVx(IssmDouble* pmaxabsvx)=0;
|
|---|
| 107 | virtual void MinVy(IssmDouble* pminvy)=0;
|
|---|
| 108 | virtual void MaxVy(IssmDouble* pmaxvy)=0;
|
|---|
| 109 | virtual void MaxAbsVy(IssmDouble* pmaxabsvy)=0;
|
|---|
| 110 | virtual void MinVz(IssmDouble* pminvz)=0;
|
|---|
| 111 | virtual void MaxVz(IssmDouble* pmaxvz)=0;
|
|---|
| 112 | virtual void MaxAbsVz(IssmDouble* pmaxabsvz)=0;
|
|---|
| 113 | virtual IssmDouble MassFlux(IssmDouble* segment)=0;
|
|---|
| [16388] | 114 | virtual IssmDouble MassFlux(IssmDouble x1,IssmDouble y1, IssmDouble x2, IssmDouble y2,int segment_id)=0;
|
|---|
| [15129] | 115 | virtual void ElementResponse(IssmDouble* presponse,int response_enum)=0;
|
|---|
| [12471] | 116 | virtual IssmDouble IceVolume(void)=0;
|
|---|
| [16275] | 117 | virtual IssmDouble IceVolumeAboveFloatation(void)=0;
|
|---|
| [12773] | 118 | virtual IssmDouble TotalSmb(void)=0;
|
|---|
| [9817] | 119 | #endif
|
|---|
| [4055] | 120 |
|
|---|
| [14589] | 121 | #ifdef _HAVE_GIA_
|
|---|
| [14807] | 122 | virtual void GiaDeflection(Vector<IssmDouble>* wg,Vector<IssmDouble>* dwgdt,IssmDouble* x,IssmDouble* y)=0;
|
|---|
| [14589] | 123 | #endif
|
|---|
| 124 |
|
|---|
| [9817] | 125 | #ifdef _HAVE_CONTROL_
|
|---|
| [13216] | 126 | virtual void Gradj(Vector<IssmDouble>* gradient,int control_type,int control_index)=0;
|
|---|
| [16314] | 127 | virtual IssmDouble ThicknessAbsMisfit(void)=0;
|
|---|
| 128 | virtual IssmDouble SurfaceAbsVelMisfit(void)=0;
|
|---|
| 129 | virtual IssmDouble SurfaceRelVelMisfit(void)=0;
|
|---|
| 130 | virtual IssmDouble SurfaceLogVelMisfit(void)=0;
|
|---|
| 131 | virtual IssmDouble SurfaceLogVxVyMisfit(void)=0;
|
|---|
| 132 | virtual IssmDouble SurfaceAverageVelMisfit(void)=0;
|
|---|
| 133 | virtual IssmDouble ThicknessAbsGradient(void)=0;
|
|---|
| 134 | virtual IssmDouble ThicknessAlongGradient(void)=0;
|
|---|
| 135 | virtual IssmDouble ThicknessAcrossGradient(void)=0;
|
|---|
| 136 | virtual IssmDouble BalancethicknessMisfit(void)=0;
|
|---|
| 137 | virtual IssmDouble RheologyBbarAbsGradient(void)=0;
|
|---|
| 138 | virtual IssmDouble DragCoefficientAbsGradient(void)=0;
|
|---|
| [13216] | 139 | virtual void ControlInputGetGradient(Vector<IssmDouble>* gradient,int enum_type,int control_index)=0;
|
|---|
| [12471] | 140 | virtual void ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index)=0;
|
|---|
| 141 | virtual void ControlInputScaleGradient(int enum_type, IssmDouble scale)=0;
|
|---|
| [16494] | 142 | virtual void ControlToVectors(Vector<IssmPDouble>* vector_control, Vector<IssmPDouble>* vector_gradient,int control_enum)=0;
|
|---|
| [13216] | 143 | virtual void GetVectorFromControlInputs(Vector<IssmDouble>* gradient,int control_enum,int control_index,const char* data)=0;
|
|---|
| [12471] | 144 | virtual void SetControlInputsFromVector(IssmDouble* vector,int control_enum,int control_index)=0;
|
|---|
| 145 | virtual void InputControlUpdate(IssmDouble scalar,bool save_parameter)=0;
|
|---|
| [9817] | 146 | #endif
|
|---|
| [14796] | 147 |
|
|---|
| [16442] | 148 | virtual void UpdateConstraintsExtrudeFromBase(void)=0;
|
|---|
| 149 | virtual void UpdateConstraintsExtrudeFromTop(void)=0;
|
|---|
| 150 |
|
|---|
| [15853] | 151 | #ifdef _HAVE_THERMAL_
|
|---|
| [16272] | 152 | virtual void UpdateBasalConstraintsEnthalpy(void)=0;
|
|---|
| [16026] | 153 | virtual void ComputeBasalMeltingrate(void)=0;
|
|---|
| [16442] | 154 | virtual void DrainWaterfraction(void)=0;
|
|---|
| [15853] | 155 | #endif
|
|---|
| 156 |
|
|---|
| [14796] | 157 | #ifdef _HAVE_HYDROLOGY_
|
|---|
| [14967] | 158 | virtual void GetHydrologyDCInefficientHmax(IssmDouble* ph_max, Node* innode)=0;
|
|---|
| [15015] | 159 | virtual void GetHydrologyTransfer(Vector<IssmDouble>* transfer)=0;
|
|---|
| [15249] | 160 | virtual void HydrologyEPLGetMask(Vector<IssmDouble>* mask)=0;
|
|---|
| 161 | virtual void HydrologyEPLGetActive(Vector<IssmDouble>* active)=0;
|
|---|
| [16600] | 162 | virtual void ComputeEPLThickness(void)=0;
|
|---|
| [16685] | 163 | virtual void UpdateConstraintsL2ProjectionEPL(void)=0;
|
|---|
| [14796] | 164 | #endif
|
|---|
| [15643] | 165 |
|
|---|
| [15386] | 166 | #ifdef _HAVE_GROUNDINGLINE_
|
|---|
| [15962] | 167 | virtual void MigrateGroundingLine(IssmDouble* sheet_ungrounding)=0;
|
|---|
| [15386] | 168 | virtual void PotentialUngrounding(Vector<IssmDouble>* potential_sheet_ungrounding)=0;
|
|---|
| 169 | virtual int UpdatePotentialUngrounding(IssmDouble* potential_sheet_ungrounding,Vector<IssmDouble>* vec_nodes_on_iceshelf,IssmDouble* nodes_on_iceshelf)=0;
|
|---|
| 170 | #endif
|
|---|
| 171 |
|
|---|
| [3683] | 172 | };
|
|---|
| 173 | #endif
|
|---|
![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
