| 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:*/
|
|---|
| 11 | /*{{{*/
|
|---|
| 12 | #include "../../datastructures/datastructures.h"
|
|---|
| 13 | #include "../../toolkits/toolkits.h"
|
|---|
| 14 | #include "../Update.h"
|
|---|
| 15 | class DataSet;
|
|---|
| 16 | class Parameters;
|
|---|
| 17 | class Parameter;
|
|---|
| 18 | class Elements;
|
|---|
| 19 | class Loads;
|
|---|
| 20 | class Nodes;
|
|---|
| 21 | class Node;
|
|---|
| 22 | class Vertices;
|
|---|
| 23 | class Vertex;
|
|---|
| 24 | class Materials;
|
|---|
| 25 | class Material;
|
|---|
| 26 | class Matpar;
|
|---|
| 27 | class Inputs;
|
|---|
| 28 | class Input;
|
|---|
| 29 | class Gauss;
|
|---|
| 30 | class ElementVector;
|
|---|
| 31 | template <class doublematrix> class Matrix;
|
|---|
| 32 | template <class doubletype> class Vector;
|
|---|
| 33 | class ElementMatrix;
|
|---|
| 34 | class ElementVector;
|
|---|
| 35 | /*}}}*/
|
|---|
| 36 |
|
|---|
| 37 | class Element: public Object,public Update{
|
|---|
| 38 |
|
|---|
| 39 | public:
|
|---|
| 40 | Inputs *inputs;
|
|---|
| 41 | Node **nodes;
|
|---|
| 42 | Vertex **vertices;
|
|---|
| 43 | Material *material;
|
|---|
| 44 | Matpar *matpar;
|
|---|
| 45 | Parameters *parameters;
|
|---|
| 46 |
|
|---|
| 47 | public:
|
|---|
| 48 | /*Constructors/Destructores*/
|
|---|
| 49 | Element();
|
|---|
| 50 | ~Element();
|
|---|
| 51 |
|
|---|
| 52 | /*Functions*/
|
|---|
| 53 | void AddInput(Input* input_in);
|
|---|
| 54 | bool AllActive(void);
|
|---|
| 55 | bool AnyActive(void);
|
|---|
| 56 | void CoordinateSystemTransform(IssmDouble** ptransform,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 57 | void DeleteMaterials(void);
|
|---|
| 58 | void FindParam(bool* pvalue,int paramenum);
|
|---|
| 59 | void FindParam(int* pvalue,int paramenum);
|
|---|
| 60 | void FindParam(IssmDouble* pvalue,int paramenum);
|
|---|
| 61 | void FindParam(int** pvalues,int* psize,int paramenum);
|
|---|
| 62 | void GetDofList(int** pdoflist,int approximation_enum,int setenum);
|
|---|
| 63 | void GetDofListVelocity(int** pdoflist,int setenum);
|
|---|
| 64 | void GetDofListPressure(int** pdoflist,int setenum);
|
|---|
| 65 | Input* GetInput(int inputenum);
|
|---|
| 66 | void GetInputListOnNodes(IssmDouble* pvalue,int enumtype);
|
|---|
| 67 | void GetInputListOnNodes(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue);
|
|---|
| 68 | void GetInputListOnNodesVelocity(IssmDouble* pvalue,int enumtype);
|
|---|
| 69 | void GetInputListOnVertices(IssmDouble* pvalue,int enumtype);
|
|---|
| 70 | void GetInputListOnVertices(IssmDouble* pvalue,int enumtype,IssmDouble defaultvalue);
|
|---|
| 71 | void GetInputValue(bool* pvalue,int enum_type);
|
|---|
| 72 | void GetInputValue(int* pvalue,int enum_type);
|
|---|
| 73 | void GetInputValue(IssmDouble* pvalue,int enum_type);
|
|---|
| 74 | void GetInputValue(IssmDouble* pvalue,Gauss* gauss,int enum_type);
|
|---|
| 75 | IssmDouble GetMaterialParameter(int enum_in);
|
|---|
| 76 | void GetPhi(IssmDouble* phi, IssmDouble* epsilon, IssmDouble viscosity);
|
|---|
| 77 | void GetVerticesSidList(int* sidlist);
|
|---|
| 78 | void GetVerticesConnectivityList(int* connectivitylist);
|
|---|
| 79 | bool IsInput(int name);
|
|---|
| 80 | bool IsFloating();
|
|---|
| 81 | void ResultInterpolation(int* pinterpolation,int*nodesperelement,int output_enum);
|
|---|
| 82 | void ResultToVector(Vector<IssmDouble>* vector,int output_enum);
|
|---|
| 83 | void ResultToPatch(IssmDouble* values,int nodesperelement,int output_enum);
|
|---|
| 84 | void StrainRateSSA(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 85 | void StrainRateSSA1d(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input);
|
|---|
| 86 | void StrainRateHO(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 87 | void StrainRateHO2dvertical(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 88 | void StrainRateFS(IssmDouble* epsilon,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input);
|
|---|
| 89 | IssmDouble TMeltingPoint(IssmDouble pressure);
|
|---|
| 90 | void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,int cs_enum);
|
|---|
| 91 | void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 92 | void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 93 | void TransformLoadVectorCoord(ElementVector* pe,int cs_enum);
|
|---|
| 94 | void TransformLoadVectorCoord(ElementVector* pe,int* cs_array);
|
|---|
| 95 | void TransformLoadVectorCoord(ElementVector* pe,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 96 | void TransformLoadVectorCoord(ElementVector* pe,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 97 | void TransformSolutionCoord(IssmDouble* solution,int cs_enum);
|
|---|
| 98 | void TransformSolutionCoord(IssmDouble* solution,int* cs_array);
|
|---|
| 99 | void TransformSolutionCoord(IssmDouble* solution,int numnodes,int cs_enum);
|
|---|
| 100 | void TransformSolutionCoord(IssmDouble* solution,int numnodes,int* cs_array);
|
|---|
| 101 | void TransformSolutionCoord(IssmDouble* solution,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 102 | void TransformSolutionCoord(IssmDouble* solution,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 103 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int cs_enum);
|
|---|
| 104 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int* cs_array);
|
|---|
| 105 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int cs_enum);
|
|---|
| 106 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int* cs_array);
|
|---|
| 107 | void ViscousHeatingCreateInput(void);
|
|---|
| 108 | void ViscosityFS(IssmDouble* pviscosity,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input);
|
|---|
| 109 | void ViscosityHO(IssmDouble* pviscosity,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 110 | void ViscosityL1L2(IssmDouble* pviscosity,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* surf);
|
|---|
| 111 | void ViscositySSA(IssmDouble* pviscosity,int dim,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input);
|
|---|
| 112 | void ViscositySSADerivativeEpsSquare(IssmDouble* pmu_prime,IssmDouble* epsilon);
|
|---|
| 113 | void ViscosityHODerivativeEpsSquare(IssmDouble* pmu_prime,IssmDouble* epsilon);
|
|---|
| 114 | void ViscosityFSDerivativeEpsSquare(IssmDouble* pmu_prime,IssmDouble* epsilon);
|
|---|
| 115 | void TransformLoadVectorCoord(ElementVector* pe,int numnodes,int transformenum){_error_("not implemented yet");};/*Tiling only*/
|
|---|
| 116 | void TransformLoadVectorCoord(ElementVector* pe,int numnodes,int* transformenum_list){_error_("not implemented yet");};/*Tiling only*/
|
|---|
| 117 | void TransformStiffnessMatrixCoord(ElementMatrix* Ke,int numnodes,int* transformenum_list){_error_("not implemented yet");};/*Tiling only*/
|
|---|
| 118 |
|
|---|
| 119 | /*Virtual functions*/
|
|---|
| 120 | virtual void AddBasalInput(int input_enum, IssmDouble* values, int interpolation_enum)=0;
|
|---|
| 121 | virtual void AddInput(int input_enum, IssmDouble* values, int interpolation_enum)=0;
|
|---|
| 122 | virtual IssmDouble CharacteristicLength(void)=0;
|
|---|
| 123 | virtual void Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters)=0;
|
|---|
| 124 | virtual void SetCurrentConfiguration(Elements* elements,Loads* loads,Nodes* nodes,Materials* materials,Parameters* parameters)=0;
|
|---|
| 125 | virtual void SetwiseNodeConnectivity(int* d_nz,int* o_nz,Node* node,bool* flags,int* flagsindices,int set1_enum,int set2_enum)=0;
|
|---|
| 126 | virtual void ElementSizes(IssmDouble* phx,IssmDouble* phy,IssmDouble* phz)=0;
|
|---|
| 127 | virtual void ThermalToEnthalpy(IssmDouble* penthalpy,IssmDouble temperature,IssmDouble waterfraction,IssmDouble pressure)=0;
|
|---|
| 128 | virtual void EnthalpyToThermal(IssmDouble* ptemperature,IssmDouble* pwaterfraction,IssmDouble enthalpy,IssmDouble pressure)=0;
|
|---|
| 129 | virtual IssmDouble EnthalpyDiffusionParameter(IssmDouble enthalpy,IssmDouble pressure)=0;
|
|---|
| 130 | virtual IssmDouble EnthalpyDiffusionParameterVolume(int numvertices,IssmDouble* enthalpy,IssmDouble* pressure)=0;
|
|---|
| 131 | virtual int FiniteElement(void)=0;
|
|---|
| 132 | virtual IssmDouble MinEdgeLength(IssmDouble* xyz_list)=0;
|
|---|
| 133 | virtual void NodalFunctions(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| 134 | virtual void NodalFunctionsP1(IssmDouble* basis,Gauss* gauss)=0;
|
|---|
| 135 | virtual void NodalFunctionsVelocity(IssmDouble* basis, Gauss* gauss)=0;
|
|---|
| 136 | virtual void NodalFunctionsPressure(IssmDouble* basis, Gauss* gauss)=0;
|
|---|
| 137 | virtual void NodalFunctionsDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 138 | virtual void NodalFunctionsP1Derivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 139 | virtual void NodalFunctionsMINIDerivatives(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 140 | virtual void NodalFunctionsDerivativesVelocity(IssmDouble* dbasis,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 141 | virtual void NormalSection(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| 142 | virtual void NormalTop(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| 143 | virtual void NormalBase(IssmDouble* normal,IssmDouble* xyz_list)=0;
|
|---|
| 144 | virtual IssmDouble PureIceEnthalpy(IssmDouble pressure)=0;
|
|---|
| 145 |
|
|---|
| 146 | virtual Element* GetUpperElement(void)=0;
|
|---|
| 147 | virtual Element* GetLowerElement(void)=0;
|
|---|
| 148 | virtual Element* GetSurfaceElement(void)=0;
|
|---|
| 149 | virtual Element* GetBasalElement(void)=0;
|
|---|
| 150 |
|
|---|
| 151 | virtual void JacobianDeterminant(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 152 | virtual void JacobianDeterminantLine(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 153 | virtual void JacobianDeterminantSurface(IssmDouble* Jdet, IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 154 | virtual void JacobianDeterminantBase(IssmDouble* Jdet,IssmDouble* xyz_list_base,Gauss* gauss)=0;
|
|---|
| 155 | virtual void JacobianDeterminantTop(IssmDouble* Jdet,IssmDouble* xyz_list_base,Gauss* gauss)=0;
|
|---|
| 156 | virtual void GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution,int solutionenum)=0;
|
|---|
| 157 | virtual int GetNodeIndex(Node* node)=0;
|
|---|
| 158 | virtual int GetNumberOfNodes(void)=0;
|
|---|
| 159 | virtual int GetNumberOfNodesVelocity(void)=0;
|
|---|
| 160 | virtual int GetNumberOfNodesPressure(void)=0;
|
|---|
| 161 | virtual int GetNumberOfVertices(void)=0;
|
|---|
| 162 | virtual void GetNodesSidList(int* sidlist)=0;
|
|---|
| 163 | virtual void GetNodesLidList(int* sidlist)=0;
|
|---|
| 164 |
|
|---|
| 165 | virtual int Id()=0;
|
|---|
| 166 | virtual int Sid()=0;
|
|---|
| 167 | virtual bool IsNodeOnShelfFromFlags(IssmDouble* flags)=0;
|
|---|
| 168 | virtual bool IsOnBed()=0;
|
|---|
| 169 | virtual bool IsOnSurface()=0;
|
|---|
| 170 | virtual bool IsIceInElement()=0;
|
|---|
| 171 | virtual void GetGroundedPart(int* point1,IssmDouble* fraction1,IssmDouble* fraction2, bool* mainlyfloating)=0;
|
|---|
| 172 | virtual IssmDouble GetGroundedPortion(IssmDouble* xyz_list)=0;
|
|---|
| 173 | virtual void GetInputValue(IssmDouble* pvalue,Node* node,int enumtype)=0;
|
|---|
| 174 | virtual Node* GetNode(int node_number)=0;
|
|---|
| 175 | virtual void GetVerticesCoordinates(IssmDouble** xyz_list)=0;
|
|---|
| 176 | virtual void GetVerticesCoordinatesBase(IssmDouble** xyz_list)=0;
|
|---|
| 177 | virtual void GetVerticesCoordinatesTop(IssmDouble** xyz_list)=0;
|
|---|
| 178 | virtual IssmDouble GetXcoord(Gauss* gauss)=0;
|
|---|
| 179 | virtual IssmDouble GetYcoord(Gauss* gauss)=0;
|
|---|
| 180 | virtual IssmDouble GetZcoord(Gauss* gauss)=0;
|
|---|
| 181 | virtual int GetElementType(void)=0;
|
|---|
| 182 |
|
|---|
| 183 | virtual IssmDouble SurfaceArea(void)=0;
|
|---|
| 184 | virtual void InputDepthAverageAtBase(int enum_type,int average_enum_type)=0;
|
|---|
| 185 | virtual void InputChangeName(int enum_type,int enum_type_old)=0;
|
|---|
| 186 | virtual void ComputeBasalStress(Vector<IssmDouble>* sigma_b)=0;
|
|---|
| 187 | virtual void ComputeStrainRate(Vector<IssmDouble>* eps)=0;
|
|---|
| 188 | virtual void ComputeStressTensor(void)=0;
|
|---|
| 189 |
|
|---|
| 190 | virtual void Update(int index, IoModel* iomodel,int analysis_counter,int analysis_type,int finite_element)=0;
|
|---|
| 191 | virtual void InputDuplicate(int original_enum,int new_enum)=0;
|
|---|
| 192 | virtual void InputCreate(IssmDouble* vector,IoModel* iomodel,int M,int N,int vector_type,int vector_enum,int code)=0;
|
|---|
| 193 | virtual void InputUpdateFromSolutionOneDofCollapsed(IssmDouble* solution,int inputenum)=0;
|
|---|
| 194 | virtual void InputUpdateFromSolutionOneDof(IssmDouble* solution,int inputenum)=0;
|
|---|
| 195 |
|
|---|
| 196 | virtual int NodalValue(IssmDouble* pvalue, int index, int natureofdataenum)=0;
|
|---|
| 197 | virtual int NumberofNodesVelocity(void)=0;
|
|---|
| 198 | virtual int NumberofNodesPressure(void)=0;
|
|---|
| 199 | virtual Gauss* NewGauss(void)=0;
|
|---|
| 200 | virtual Gauss* NewGauss(int order)=0;
|
|---|
| 201 | virtual Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order)=0;
|
|---|
| 202 | virtual Gauss* NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order_horiz,int order_vert)=0;
|
|---|
| 203 | virtual Gauss* NewGauss(int point1,IssmDouble fraction1,IssmDouble fraction2,bool mainlyfloating,int order)=0;
|
|---|
| 204 | virtual Gauss* NewGaussBase(int order)=0;
|
|---|
| 205 | virtual Gauss* NewGaussLine(int vertex1,int vertex2,int order)=0;
|
|---|
| 206 | virtual Gauss* NewGaussTop(int order)=0;
|
|---|
| 207 | virtual ElementVector* NewElementVector(int approximation_enum=NoneApproximationEnum)=0;
|
|---|
| 208 | virtual ElementMatrix* NewElementMatrix(int approximation_enum=NoneApproximationEnum)=0;
|
|---|
| 209 | virtual ElementMatrix* NewElementMatrixCoupling(int number_nodes,int approximation_enum=NoneApproximationEnum)=0;
|
|---|
| 210 | virtual void InputScale(int enum_type,IssmDouble scale_factor)=0;
|
|---|
| 211 | virtual void GetVectorFromInputs(Vector<IssmDouble>* vector, int name_enum)=0;
|
|---|
| 212 | virtual IssmDouble TimeAdapt()=0;
|
|---|
| 213 | virtual void PositiveDegreeDay(IssmDouble* pdds,IssmDouble* pds,IssmDouble signorm)=0;
|
|---|
| 214 | virtual void Delta18oParameterization(void)=0;
|
|---|
| 215 | virtual void SmbGradients(void)=0;
|
|---|
| 216 | virtual Element* SpawnBasalElement(void)=0;
|
|---|
| 217 | virtual Element* SpawnTopElement(void)=0;
|
|---|
| 218 | virtual void ReduceMatrices(ElementMatrix* Ke,ElementVector* pe)=0;
|
|---|
| 219 | virtual void ResetFSBasalBoundaryCondition()=0;
|
|---|
| 220 | virtual void SetTemporaryElementType(int element_type_in)=0;
|
|---|
| 221 | virtual IssmDouble StabilizationParameter(IssmDouble u, IssmDouble v, IssmDouble w, IssmDouble diameter, IssmDouble kappa)=0;
|
|---|
| 222 | virtual void ValueP1OnGauss(IssmDouble* pvalue,IssmDouble* values,Gauss* gauss)=0;
|
|---|
| 223 | virtual void ValueP1DerivativesOnGauss(IssmDouble* dvalue,IssmDouble* values,IssmDouble* xyz_list,Gauss* gauss)=0;
|
|---|
| 224 | virtual void ViscousHeating(IssmDouble* pphi,IssmDouble* xyz_list,Gauss* gauss,Input* vx_input,Input* vy_input,Input* vz_input)=0;
|
|---|
| 225 |
|
|---|
| 226 | virtual int VelocityInterpolation()=0;
|
|---|
| 227 | virtual int PressureInterpolation()=0;
|
|---|
| 228 | virtual bool IsZeroLevelset(int levelset_enum)=0;
|
|---|
| 229 | virtual bool IsIcefront(void)=0;
|
|---|
| 230 | virtual void ZeroLevelsetCoordinates(IssmDouble** pxyz_zero,IssmDouble* xyz_list,int levelsetenum)=0;
|
|---|
| 231 | virtual void GetIcefrontCoordinates(IssmDouble** pxyz_front,IssmDouble* xyz_list,int levelsetenum)=0;
|
|---|
| 232 | virtual void GetNormalFromLSF(IssmDouble *pnormal)=0;
|
|---|
| 233 |
|
|---|
| 234 | virtual void AverageOntoPartition(Vector<IssmDouble>* partition_contributions,Vector<IssmDouble>* partition_areas,IssmDouble* vertex_response,IssmDouble* qmu_part)=0;
|
|---|
| 235 | virtual void MinVel(IssmDouble* pminvel)=0;
|
|---|
| 236 | virtual void MaxVel(IssmDouble* pmaxvel)=0;
|
|---|
| 237 | virtual void MinVx(IssmDouble* pminvx)=0;
|
|---|
| 238 | virtual void MaxVx(IssmDouble* pmaxvx)=0;
|
|---|
| 239 | virtual void MaxAbsVx(IssmDouble* pmaxabsvx)=0;
|
|---|
| 240 | virtual void MinVy(IssmDouble* pminvy)=0;
|
|---|
| 241 | virtual void MaxVy(IssmDouble* pmaxvy)=0;
|
|---|
| 242 | virtual void MaxAbsVy(IssmDouble* pmaxabsvy)=0;
|
|---|
| 243 | virtual void MinVz(IssmDouble* pminvz)=0;
|
|---|
| 244 | virtual void MaxVz(IssmDouble* pmaxvz)=0;
|
|---|
| 245 | virtual void MaxAbsVz(IssmDouble* pmaxabsvz)=0;
|
|---|
| 246 | virtual IssmDouble MassFlux(IssmDouble* segment)=0;
|
|---|
| 247 | virtual IssmDouble MassFlux(IssmDouble x1,IssmDouble y1, IssmDouble x2, IssmDouble y2,int segment_id)=0;
|
|---|
| 248 | virtual void ElementResponse(IssmDouble* presponse,int response_enum)=0;
|
|---|
| 249 | virtual IssmDouble IceVolume(void)=0;
|
|---|
| 250 | virtual IssmDouble IceVolumeAboveFloatation(void)=0;
|
|---|
| 251 | virtual IssmDouble TotalSmb(void)=0;
|
|---|
| 252 | virtual IssmDouble Misfit(int modelenum,int observationenum,int weightsenum)=0;
|
|---|
| 253 | virtual IssmDouble MisfitArea(int weightsenum)=0;
|
|---|
| 254 | virtual int VertexConnectivity(int vertexindex)=0;
|
|---|
| 255 | virtual void VerticalSegmentIndices(int** pindices,int* pnumseg)=0;
|
|---|
| 256 |
|
|---|
| 257 | #ifdef _HAVE_GIA_
|
|---|
| 258 | virtual void GiaDeflection(Vector<IssmDouble>* wg,Vector<IssmDouble>* dwgdt,IssmDouble* x,IssmDouble* y)=0;
|
|---|
| 259 | #endif
|
|---|
| 260 |
|
|---|
| 261 | virtual void Gradj(Vector<IssmDouble>* gradient,int control_type,int control_index)=0;
|
|---|
| 262 | virtual IssmDouble ThicknessAbsMisfit(void)=0;
|
|---|
| 263 | virtual IssmDouble SurfaceAbsVelMisfit(void)=0;
|
|---|
| 264 | virtual IssmDouble SurfaceRelVelMisfit(void)=0;
|
|---|
| 265 | virtual IssmDouble SurfaceLogVelMisfit(void)=0;
|
|---|
| 266 | virtual IssmDouble SurfaceLogVxVyMisfit(void)=0;
|
|---|
| 267 | virtual IssmDouble SurfaceAverageVelMisfit(void)=0;
|
|---|
| 268 | virtual IssmDouble ThicknessAbsGradient(void)=0;
|
|---|
| 269 | virtual IssmDouble ThicknessAlongGradient(void)=0;
|
|---|
| 270 | virtual IssmDouble ThicknessAcrossGradient(void)=0;
|
|---|
| 271 | virtual IssmDouble BalancethicknessMisfit(void)=0;
|
|---|
| 272 | virtual IssmDouble RheologyBbarAbsGradient(void)=0;
|
|---|
| 273 | virtual IssmDouble DragCoefficientAbsGradient(void)=0;
|
|---|
| 274 | virtual void ControlInputGetGradient(Vector<IssmDouble>* gradient,int enum_type,int control_index)=0;
|
|---|
| 275 | virtual void ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index)=0;
|
|---|
| 276 | virtual void ControlInputScaleGradient(int enum_type, IssmDouble scale)=0;
|
|---|
| 277 | virtual void ControlToVectors(Vector<IssmPDouble>* vector_control, Vector<IssmPDouble>* vector_gradient,int control_enum)=0;
|
|---|
| 278 | virtual void GetVectorFromControlInputs(Vector<IssmDouble>* gradient,int control_enum,int control_index,const char* data)=0;
|
|---|
| 279 | virtual void SetControlInputsFromVector(IssmDouble* vector,int control_enum,int control_index)=0;
|
|---|
| 280 | virtual void InputControlUpdate(IssmDouble scalar,bool save_parameter)=0;
|
|---|
| 281 |
|
|---|
| 282 | virtual void UpdateConstraintsExtrudeFromBase(void)=0;
|
|---|
| 283 | virtual void UpdateConstraintsExtrudeFromTop(void)=0;
|
|---|
| 284 |
|
|---|
| 285 | // virtual void GetHydrologyDCInefficientHmax(IssmDouble* ph_max,int index)=0;
|
|---|
| 286 | virtual void GetHydrologyDCInefficientHmax(IssmDouble* ph_max, Node* innode)=0;
|
|---|
| 287 | virtual void GetHydrologyTransfer(Vector<IssmDouble>* transfer)=0;
|
|---|
| 288 | virtual void HydrologyEPLGetMask(Vector<IssmDouble>* mask)=0;
|
|---|
| 289 | virtual void HydrologyEPLGetActive(Vector<IssmDouble>* active)=0;
|
|---|
| 290 | virtual void ComputeEPLThickness(void)=0;
|
|---|
| 291 |
|
|---|
| 292 | virtual void MigrateGroundingLine(IssmDouble* sheet_ungrounding)=0;
|
|---|
| 293 | virtual void PotentialUngrounding(Vector<IssmDouble>* potential_sheet_ungrounding)=0;
|
|---|
| 294 | virtual int UpdatePotentialUngrounding(IssmDouble* potential_sheet_ungrounding,Vector<IssmDouble>* vec_nodes_on_iceshelf,IssmDouble* nodes_on_iceshelf)=0;
|
|---|
| 295 |
|
|---|
| 296 | };
|
|---|
| 297 | #endif
|
|---|
![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
