| [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 |
|
|---|
| [16853] | 39 | virtual bool AllActive(void)=0;
|
|---|
| 40 | virtual bool AnyActive(void)=0;
|
|---|
| [16832] | 41 | virtual IssmDouble CharacteristicLength(void)=0;
|
|---|
| [14761] | 42 | virtual void Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters)=0;
|
|---|
| [14953] | 43 | virtual void SetCurrentConfiguration(Elements* elements,Loads* loads,Nodes* nodes,Materials* materials,Parameters* parameters)=0;
|
|---|
| [16042] | 44 | virtual void SetwiseNodeConnectivity(int* d_nz,int* o_nz,Node* node,bool* flags,int* flagsindices,int set1_enum,int set2_enum)=0;
|
|---|
| [16720] | 45 | virtual void AddBasalInput(int input_enum, IssmDouble* values, int interpolation_enum)=0;
|
|---|
| [16693] | 46 | virtual void AddInput(int input_enum, IssmDouble* values, int interpolation_enum)=0;
|
|---|
| [16864] | 47 | virtual void AddInput(Input* input_in)=0;
|
|---|
| [16699] | 48 | virtual void AddMaterialInput(int input_enum, IssmDouble* values, int interpolation_enum)=0;
|
|---|
| [16789] | 49 | virtual ElementMatrix* CreateKMatrix(void)=0;
|
|---|
| [14435] | 50 | virtual void CreateDVector(Vector<IssmDouble>* df)=0;
|
|---|
| [16789] | 51 | virtual ElementVector* CreatePVector(void)=0;
|
|---|
| [13216] | 52 | virtual void CreateJacobianMatrix(Matrix<IssmDouble>* Jff)=0;
|
|---|
| [16778] | 53 | virtual void DeleteMaterials(void)=0;
|
|---|
| [16840] | 54 | virtual void ElementSizes(IssmDouble* phx,IssmDouble* phy,IssmDouble* phz)=0;
|
|---|
| [16734] | 55 | virtual void EnthalpyToThermal(IssmDouble* ptemperature,IssmDouble* pwaterfraction,IssmDouble enthalpy,IssmDouble pressure)=0;
|
|---|
| [16812] | 56 | virtual IssmDouble EnthalpyDiffusionParameter(IssmDouble enthalpy,IssmDouble pressure)=0;
|
|---|
| 57 | virtual IssmDouble EnthalpyDiffusionParameterVolume(int numvertices,IssmDouble* enthalpy,IssmDouble* pressure)=0;
|
|---|
| [16876] | 58 | virtual void FindParam(bool* pvalue,int paramenum)=0;
|
|---|
| [16675] | 59 | virtual void FindParam(int* pvalue,int paramenum)=0;
|
|---|
| 60 | virtual void FindParam(IssmDouble* pvalue,int paramenum)=0;
|
|---|
| [16847] | 61 | virtual void FindParam(int** pvalues,int* psize,int paramenum)=0;
|
|---|
| [16675] | 62 | virtual int FiniteElement(void)=0;
|
|---|
| [16697] | 63 | virtual IssmDouble GetMaterialParameter(int enum_in)=0;
|
|---|
| [16810] | 64 | virtual IssmDouble MinEdgeLength(IssmDouble* xyz_list)=0;
|
|---|
| [16799] | 65 | virtual void NodalFunctions(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| [16831] | 66 | virtual void NodalFunctionsVelocity(IssmDouble* basis, Gauss* gauss)=0;
|
|---|
| 67 | virtual void NodalFunctionsPressure(IssmDouble* basis, Gauss* gauss)=0;
|
|---|
| [16810] | 68 | virtual void NodalFunctionsDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| [16831] | 69 | virtual void NodalFunctionsDerivativesVelocity(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| [16839] | 70 | virtual void NormalSection(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| [16827] | 71 | virtual void NormalTop(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| [16844] | 72 | virtual void NormalBase(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| [16813] | 73 | virtual IssmDouble PureIceEnthalpy(IssmDouble pressure)=0;
|
|---|
| [16799] | 74 | virtual void TransformLoadVectorCoord(ElementVector* pe,int transformenum)=0;
|
|---|
| 75 | virtual void TransformLoadVectorCoord(ElementVector* pe,int* transformenum_list)=0;
|
|---|
| 76 | virtual void TransformLoadVectorCoord(ElementVector* pe,int numnodes,int transformenum)=0;
|
|---|
| 77 | virtual void TransformLoadVectorCoord(ElementVector* pe,int numnodes,int* transformenum_list)=0;
|
|---|
| [16693] | 78 | virtual void TransformSolutionCoord(IssmDouble* values,int transformenum)=0;
|
|---|
| [16715] | 79 | virtual void TransformSolutionCoord(IssmDouble* values,int* transformenum_list)=0;
|
|---|
| [16775] | 80 | virtual void TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum)=0;
|
|---|
| 81 | virtual void TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list)=0;
|
|---|
| [16818] | 82 | virtual void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int transformenum)=0;
|
|---|
| 83 | virtual void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int* transformenum_list)=0;
|
|---|
| 84 | virtual void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int numnodes,int transformenum)=0;
|
|---|
| 85 | virtual void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int numnodes,int* transformenum_list)=0;
|
|---|
| [16762] | 86 | virtual Element* GetBasalElement(void)=0;
|
|---|
| [16675] | 87 | virtual void GetDofList(int** pdoflist,int approximation_enum,int setenum)=0;
|
|---|
| 88 | virtual void GetDofListVelocity(int** pdoflist,int setenum)=0;
|
|---|
| 89 | virtual void GetDofListPressure(int** pdoflist,int setenum)=0;
|
|---|
| [16799] | 90 | virtual void JacobianDeterminant(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| [16850] | 91 | virtual void JacobianDeterminantLine(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| [16838] | 92 | virtual void JacobianDeterminantSurface(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| [16803] | 93 | virtual void JacobianDeterminantBase(IssmDouble* Jdet,IssmDouble* xyz_list_base,Gauss* gauss)=0;
|
|---|
| [16827] | 94 | virtual void JacobianDeterminantTop(IssmDouble* Jdet,IssmDouble* xyz_list_base,Gauss* gauss)=0;
|
|---|
| [16675] | 95 | virtual void GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution,int solutionenum)=0;
|
|---|
| [5719] | 96 | virtual int GetNodeIndex(Node* node)=0;
|
|---|
| [13886] | 97 | virtual int GetNumberOfNodes(void)=0;
|
|---|
| [16715] | 98 | virtual int GetNumberOfNodesVelocity(void)=0;
|
|---|
| 99 | virtual int GetNumberOfNodesPressure(void)=0;
|
|---|
| [16727] | 100 | virtual int GetNumberOfVertices(void)=0;
|
|---|
| [13892] | 101 | virtual void GetNodesSidList(int* sidlist)=0;
|
|---|
| [16042] | 102 | virtual void GetNodesLidList(int* sidlist)=0;
|
|---|
| [13892] | 103 |
|
|---|
| [7089] | 104 | virtual int Sid()=0;
|
|---|
| [10143] | 105 | virtual bool IsFloating()=0;
|
|---|
| [12471] | 106 | virtual bool IsNodeOnShelfFromFlags(IssmDouble* flags)=0;
|
|---|
| [5745] | 107 | virtual bool IsOnBed()=0;
|
|---|
| [16827] | 108 | virtual bool IsOnSurface()=0;
|
|---|
| [16872] | 109 | virtual void GetGroundedPart(int* point1,IssmDouble* fraction1,IssmDouble* fraction2, bool* mainlyfloating)=0;
|
|---|
| 110 | virtual IssmDouble GetGroundedPortion(IssmDouble* xyz_list)=0;
|
|---|
| [16693] | 111 | virtual void GetInputListOnNodes(IssmDouble* pvalue,int enumtype)=0;
|
|---|
| 112 | virtual void GetInputListOnNodes(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue)=0;
|
|---|
| [12471] | 113 | virtual void GetInputListOnVertices(IssmDouble* pvalue,int enumtype)=0;
|
|---|
| 114 | virtual void GetInputListOnVertices(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue)=0;
|
|---|
| [16675] | 115 | virtual Input* GetInput(int inputenum)=0;
|
|---|
| [16848] | 116 | virtual Input* GetMaterialInput(int inputenum)=0;
|
|---|
| [12471] | 117 | virtual void GetInputValue(IssmDouble* pvalue,Node* node,int enumtype)=0;
|
|---|
| [16675] | 118 | virtual void GetInputValue(bool* pvalue,int enum_type)=0;
|
|---|
| 119 | virtual void GetInputValue(int* pvalue,int enum_type)=0;
|
|---|
| 120 | virtual void GetInputValue(IssmDouble* pvalue,int enum_type)=0;
|
|---|
| [16873] | 121 | virtual void GetInputValue(IssmDouble* pvalue,Gauss* gauss,int enum_type)=0;
|
|---|
| [16698] | 122 | virtual void GetVerticesCoordinates(IssmDouble** xyz_list)=0;
|
|---|
| [16803] | 123 | virtual void GetVerticesCoordinatesBase(IssmDouble** xyz_list)=0;
|
|---|
| [16827] | 124 | virtual void GetVerticesCoordinatesTop(IssmDouble** xyz_list)=0;
|
|---|
| [16839] | 125 | virtual IssmDouble GetZcoord(Gauss* gauss)=0;
|
|---|
| [16866] | 126 | virtual IssmDouble GetYcoord(Gauss* gauss)=0;
|
|---|
| [16188] | 127 | virtual void GetMaterialInputValue(IssmDouble* pvalue,Node* node,int enumtype)=0;
|
|---|
| [13623] | 128 |
|
|---|
| [12471] | 129 | virtual IssmDouble SurfaceArea(void)=0;
|
|---|
| [4931] | 130 | virtual void InputDepthAverageAtBase(int enum_type,int average_enum_type,int object_enum)=0;
|
|---|
| [16693] | 131 | virtual void InputChangeName(int enum_type,int enum_type_old)=0;
|
|---|
| [13216] | 132 | virtual void ComputeBasalStress(Vector<IssmDouble>* sigma_b)=0;
|
|---|
| 133 | virtual void ComputeStrainRate(Vector<IssmDouble>* eps)=0;
|
|---|
| [16848] | 134 | virtual void ComputeStressTensor(void)=0;
|
|---|
| [16461] | 135 | virtual void ResultInterpolation(int* pinterpolation,int output_enum)=0;
|
|---|
| [16589] | 136 | virtual void ResultToVector(Vector<IssmDouble>* vector,int output_enum)=0;
|
|---|
| [15613] | 137 | virtual void Update(int index, IoModel* iomodel,int analysis_counter,int analysis_type,int finite_element)=0;
|
|---|
| [8967] | 138 | virtual void InputDuplicate(int original_enum,int new_enum)=0;
|
|---|
| [16208] | 139 | virtual void InputCreate(IssmDouble* vector,IoModel* iomodel,int M,int N,int vector_type,int vector_enum,int code)=0;
|
|---|
| [16715] | 140 | virtual void InputUpdateFromSolutionOneDofCollapsed(IssmDouble* solution,int inputenum)=0;
|
|---|
| [16717] | 141 | virtual void InputUpdateFromSolutionOneDof(IssmDouble* solution,int inputenum)=0;
|
|---|
| [13623] | 142 |
|
|---|
| [15130] | 143 | virtual int NodalValue(IssmDouble* pvalue, int index, int natureofdataenum)=0;
|
|---|
| [16675] | 144 | virtual int NumberofNodesVelocity(void)=0;
|
|---|
| 145 | virtual int NumberofNodesPressure(void)=0;
|
|---|
| 146 | virtual Gauss* NewGauss(void)=0;
|
|---|
| [16799] | 147 | virtual Gauss* NewGauss(int order)=0;
|
|---|
| [16838] | 148 | virtual Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order)=0;
|
|---|
| [16839] | 149 | virtual Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order_horiz,int order_vert)=0;
|
|---|
| [16872] | 150 | virtual Gauss* NewGauss(int point1,IssmDouble fraction1,IssmDouble fraction2,bool mainlyfloating,int order)=0;
|
|---|
| [16803] | 151 | virtual Gauss* NewGaussBase(int order)=0;
|
|---|
| [16850] | 152 | virtual Gauss* NewGaussLine(int vertex1,int vertex2,int order)=0;
|
|---|
| [16827] | 153 | virtual Gauss* NewGaussTop(int order)=0;
|
|---|
| [16803] | 154 | virtual ElementVector* NewElementVector(int approximation_enum=NoneApproximationEnum)=0;
|
|---|
| [16818] | 155 | virtual ElementMatrix* NewElementMatrix(int approximation_enum=NoneApproximationEnum)=0;
|
|---|
| [12471] | 156 | virtual void InputScale(int enum_type,IssmDouble scale_factor)=0;
|
|---|
| [13216] | 157 | virtual void GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum)=0;
|
|---|
| [12471] | 158 | virtual IssmDouble TimeAdapt()=0;
|
|---|
| 159 | virtual void PositiveDegreeDay(IssmDouble* pdds,IssmDouble* pds,IssmDouble signorm)=0;
|
|---|
| [12748] | 160 | virtual void Delta18oParameterization(void)=0;
|
|---|
| [16721] | 161 | virtual void SmbGradients(void)=0;
|
|---|
| 162 | virtual Element* SpawnBasalElement(void)=0;
|
|---|
| [16699] | 163 | virtual IssmDouble TMeltingPoint(IssmDouble pressure)=0;
|
|---|
| [16789] | 164 | virtual void ReduceMatrices(ElementMatrix* Ke,ElementVector* pe)=0;
|
|---|
| [10440] | 165 | virtual void ResetCoordinateSystem()=0;
|
|---|
| [16810] | 166 | virtual IssmDouble StabilizationParameter(IssmDouble u, IssmDouble v, IssmDouble w, IssmDouble diameter, IssmDouble kappa)=0;
|
|---|
| [16860] | 167 | virtual void ValueP1OnGauss(IssmDouble* pvalue,IssmDouble* values,Gauss* gauss)=0;
|
|---|
| 168 | virtual void ValueP1DerivativesOnGauss(IssmDouble* dvalue,IssmDouble* values,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| [16810] | 169 | virtual void ViscousHeating(IssmDouble* pphi,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input)=0;
|
|---|
| [16831] | 170 | virtual void ViscosityFS(IssmDouble* pviscosity,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input)=0;
|
|---|
| [16829] | 171 | virtual void ViscosityHO(IssmDouble* pviscosity,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input)=0;
|
|---|
| [16818] | 172 | virtual void ViscositySSA(IssmDouble* pviscosity,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input)=0;
|
|---|
| [16876] | 173 | virtual void ViscositySSADerivativeEpsSquare(IssmDouble* pmu_prime,IssmDouble* epsilon)=0;
|
|---|
| [16886] | 174 | virtual void ViscosityHODerivativeEpsSquare(IssmDouble* pmu_prime,IssmDouble* epsilon)=0;
|
|---|
| 175 | virtual void ViscosityFSDerivativeEpsSquare(IssmDouble* pmu_prime,IssmDouble* epsilon)=0;
|
|---|
| [16876] | 176 | virtual void StrainRateSSA(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input)=0;
|
|---|
| [16886] | 177 | virtual void StrainRateHO(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input)=0;
|
|---|
| 178 | virtual void StrainRateFS(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input)=0;
|
|---|
| [16675] | 179 | virtual int VelocityInterpolation()=0;
|
|---|
| 180 | virtual int PressureInterpolation()=0;
|
|---|
| [16833] | 181 | virtual bool IsZeroLevelset(int levelset_enum)=0;
|
|---|
| [16838] | 182 | virtual void ZeroLevelsetCoordinates(IssmDouble** pxyz_zero,IssmDouble* xyz_list,int levelsetenum)=0;
|
|---|
| [9817] | 183 |
|
|---|
| 184 | #ifdef _HAVE_RESPONSES_
|
|---|
| [15387] | 185 | virtual void AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part)=0;
|
|---|
| [15130] | 186 | virtual void MinVel(IssmDouble* pminvel)=0;
|
|---|
| 187 | virtual void MaxVel(IssmDouble* pmaxvel)=0;
|
|---|
| 188 | virtual void MinVx(IssmDouble* pminvx)=0;
|
|---|
| 189 | virtual void MaxVx(IssmDouble* pmaxvx)=0;
|
|---|
| 190 | virtual void MaxAbsVx(IssmDouble* pmaxabsvx)=0;
|
|---|
| 191 | virtual void MinVy(IssmDouble* pminvy)=0;
|
|---|
| 192 | virtual void MaxVy(IssmDouble* pmaxvy)=0;
|
|---|
| 193 | virtual void MaxAbsVy(IssmDouble* pmaxabsvy)=0;
|
|---|
| 194 | virtual void MinVz(IssmDouble* pminvz)=0;
|
|---|
| 195 | virtual void MaxVz(IssmDouble* pmaxvz)=0;
|
|---|
| 196 | virtual void MaxAbsVz(IssmDouble* pmaxabsvz)=0;
|
|---|
| 197 | virtual IssmDouble MassFlux(IssmDouble* segment)=0;
|
|---|
| [16388] | 198 | virtual IssmDouble MassFlux(IssmDouble x1,IssmDouble y1, IssmDouble x2, IssmDouble y2,int segment_id)=0;
|
|---|
| [15129] | 199 | virtual void ElementResponse(IssmDouble* presponse,int response_enum)=0;
|
|---|
| [12471] | 200 | virtual IssmDouble IceVolume(void)=0;
|
|---|
| [16275] | 201 | virtual IssmDouble IceVolumeAboveFloatation(void)=0;
|
|---|
| [12773] | 202 | virtual IssmDouble TotalSmb(void)=0;
|
|---|
| [16787] | 203 | virtual IssmDouble Misfit(int modelenum,int observationenum,int weightsenum)=0;
|
|---|
| [16849] | 204 | virtual int VertexConnectivity(int vertexindex)=0;
|
|---|
| [16850] | 205 | virtual void VerticalSegmentIndices(int** pindices,int* pnumseg)=0;
|
|---|
| [9817] | 206 | #endif
|
|---|
| [4055] | 207 |
|
|---|
| [14589] | 208 | #ifdef _HAVE_GIA_
|
|---|
| [14807] | 209 | virtual void GiaDeflection(Vector<IssmDouble>* wg,Vector<IssmDouble>* dwgdt,IssmDouble* x,IssmDouble* y)=0;
|
|---|
| [14589] | 210 | #endif
|
|---|
| 211 |
|
|---|
| [9817] | 212 | #ifdef _HAVE_CONTROL_
|
|---|
| [13216] | 213 | virtual void Gradj(Vector<IssmDouble>* gradient,int control_type,int control_index)=0;
|
|---|
| [16314] | 214 | virtual IssmDouble ThicknessAbsMisfit(void)=0;
|
|---|
| 215 | virtual IssmDouble SurfaceAbsVelMisfit(void)=0;
|
|---|
| 216 | virtual IssmDouble SurfaceRelVelMisfit(void)=0;
|
|---|
| 217 | virtual IssmDouble SurfaceLogVelMisfit(void)=0;
|
|---|
| 218 | virtual IssmDouble SurfaceLogVxVyMisfit(void)=0;
|
|---|
| 219 | virtual IssmDouble SurfaceAverageVelMisfit(void)=0;
|
|---|
| 220 | virtual IssmDouble ThicknessAbsGradient(void)=0;
|
|---|
| 221 | virtual IssmDouble ThicknessAlongGradient(void)=0;
|
|---|
| 222 | virtual IssmDouble ThicknessAcrossGradient(void)=0;
|
|---|
| 223 | virtual IssmDouble BalancethicknessMisfit(void)=0;
|
|---|
| 224 | virtual IssmDouble RheologyBbarAbsGradient(void)=0;
|
|---|
| 225 | virtual IssmDouble DragCoefficientAbsGradient(void)=0;
|
|---|
| [13216] | 226 | virtual void ControlInputGetGradient(Vector<IssmDouble>* gradient,int enum_type,int control_index)=0;
|
|---|
| [12471] | 227 | virtual void ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index)=0;
|
|---|
| 228 | virtual void ControlInputScaleGradient(int enum_type, IssmDouble scale)=0;
|
|---|
| [16494] | 229 | virtual void ControlToVectors(Vector<IssmPDouble>* vector_control, Vector<IssmPDouble>* vector_gradient,int control_enum)=0;
|
|---|
| [13216] | 230 | virtual void GetVectorFromControlInputs(Vector<IssmDouble>* gradient,int control_enum,int control_index,const char* data)=0;
|
|---|
| [12471] | 231 | virtual void SetControlInputsFromVector(IssmDouble* vector,int control_enum,int control_index)=0;
|
|---|
| 232 | virtual void InputControlUpdate(IssmDouble scalar,bool save_parameter)=0;
|
|---|
| [9817] | 233 | #endif
|
|---|
| [14796] | 234 |
|
|---|
| [16442] | 235 | virtual void UpdateConstraintsExtrudeFromBase(void)=0;
|
|---|
| 236 | virtual void UpdateConstraintsExtrudeFromTop(void)=0;
|
|---|
| 237 |
|
|---|
| [15853] | 238 | #ifdef _HAVE_THERMAL_
|
|---|
| [16272] | 239 | virtual void UpdateBasalConstraintsEnthalpy(void)=0;
|
|---|
| [16026] | 240 | virtual void ComputeBasalMeltingrate(void)=0;
|
|---|
| [16816] | 241 | virtual void DrainWaterfraction(IssmDouble* drainrate_element)=0;
|
|---|
| [15853] | 242 | #endif
|
|---|
| 243 |
|
|---|
| [14796] | 244 | #ifdef _HAVE_HYDROLOGY_
|
|---|
| [14967] | 245 | virtual void GetHydrologyDCInefficientHmax(IssmDouble* ph_max, Node* innode)=0;
|
|---|
| [16761] | 246 | virtual void GetHydrologyDCInefficientHmax(IssmDouble* ph_max, int index)=0;
|
|---|
| [15015] | 247 | virtual void GetHydrologyTransfer(Vector<IssmDouble>* transfer)=0;
|
|---|
| [15249] | 248 | virtual void HydrologyEPLGetMask(Vector<IssmDouble>* mask)=0;
|
|---|
| 249 | virtual void HydrologyEPLGetActive(Vector<IssmDouble>* active)=0;
|
|---|
| [16600] | 250 | virtual void ComputeEPLThickness(void)=0;
|
|---|
| [14796] | 251 | #endif
|
|---|
| [15643] | 252 |
|
|---|
| [15386] | 253 | #ifdef _HAVE_GROUNDINGLINE_
|
|---|
| [15962] | 254 | virtual void MigrateGroundingLine(IssmDouble* sheet_ungrounding)=0;
|
|---|
| [15386] | 255 | virtual void PotentialUngrounding(Vector<IssmDouble>* potential_sheet_ungrounding)=0;
|
|---|
| 256 | virtual int UpdatePotentialUngrounding(IssmDouble* potential_sheet_ungrounding,Vector<IssmDouble>* vec_nodes_on_iceshelf,IssmDouble* nodes_on_iceshelf)=0;
|
|---|
| 257 | #endif
|
|---|
| 258 |
|
|---|
| [3683] | 259 | };
|
|---|
| 260 | #endif
|
|---|
![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
