[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 |
|
---|
| 8 | #ifndef _ELEMENT_H_
|
---|
| 9 | #define _ELEMENT_H_
|
---|
| 10 |
|
---|
| 11 | /*Headers:*/
|
---|
| 12 | /*{{{1*/
|
---|
| 13 | #include "../Object.h"
|
---|
| 14 |
|
---|
| 15 | class DataSet;
|
---|
| 16 | class Parameters;
|
---|
[4039] | 17 | class Patch;
|
---|
[3683] | 18 |
|
---|
| 19 | #include "../../toolkits/toolkits.h"
|
---|
| 20 | /*}}}*/
|
---|
| 21 |
|
---|
[4244] | 22 | class Element: public Object,public Update{
|
---|
[3683] | 23 |
|
---|
| 24 | public:
|
---|
| 25 |
|
---|
| 26 | virtual ~Element(){};
|
---|
[4248] | 27 |
|
---|
[4218] | 28 | virtual void Configure(Elements* elements,Loads* loads,DataSet* nodes,Materials* materials,Parameters* parameters)=0;
|
---|
[4575] | 29 | virtual void SetCurrentConfiguration(Elements* elements,Loads* loads,DataSet* nodes,Materials* materials,Parameters* parameters)=0;
|
---|
[8800] | 30 | virtual void CreateKMatrix(Mat Kff, Mat Kfs,Vec df)=0;
|
---|
| 31 | virtual void CreatePVector(Vec pf)=0;
|
---|
[11527] | 32 | virtual void CreateJacobianMatrix(Mat Jff)=0;
|
---|
[4043] | 33 | virtual void GetSolutionFromInputs(Vec solution)=0;
|
---|
[5719] | 34 | virtual int GetNodeIndex(Node* node)=0;
|
---|
[7089] | 35 | virtual int Sid()=0;
|
---|
[10143] | 36 | virtual bool IsFloating()=0;
|
---|
[7089] | 37 | virtual bool IsNodeOnShelf()=0;
|
---|
[7323] | 38 | virtual bool IsNodeOnShelfFromFlags(double* flags)=0;
|
---|
[5745] | 39 | virtual bool IsOnBed()=0;
|
---|
[10135] | 40 | virtual void GetInputListOnVertices(double* pvalue,int enumtype)=0;
|
---|
| 41 | virtual void GetInputListOnVertices(double* pvalue,int enumtype,double defaultvalue)=0;
|
---|
| 42 | virtual void GetInputValue(double* pvalue,Node* node,int enumtype)=0;
|
---|
[9817] | 43 |
|
---|
[4043] | 44 | virtual double SurfaceArea(void)=0;
|
---|
[4931] | 45 | virtual void InputDepthAverageAtBase(int enum_type,int average_enum_type,int object_enum)=0;
|
---|
[4043] | 46 | virtual void ComputeBasalStress(Vec sigma_b)=0;
|
---|
| 47 | virtual void ComputeStrainRate(Vec eps)=0;
|
---|
[4039] | 48 | virtual void PatchSize(int* pnumrows, int* pnumvertices,int* pnumnodes)=0;
|
---|
| 49 | virtual void PatchFill(int* pcount, Patch* patch)=0;
|
---|
[11027] | 50 | virtual void ListResultsInfo(int** results_enums,int** results_size,double** results_times,int** results_steps,int* num_results)=0;
|
---|
[4873] | 51 | virtual void DeleteResults(void)=0;
|
---|
[4025] | 52 | virtual void Update(int index, IoModel* iomodel,int analysis_counter,int analysis_type)=0;
|
---|
[4039] | 53 | virtual void InputToResult(int enum_type,int step,double time)=0;
|
---|
[8967] | 54 | virtual void InputDuplicate(int original_enum,int new_enum)=0;
|
---|
| 55 | virtual void InputCreate(double scalar,int name,int code)=0;
|
---|
| 56 | virtual void InputCreate(double* vector, int index,IoModel* iomodel,int M,int N,int vector_type,int vector_enum,int code)=0;
|
---|
[9775] | 57 | virtual void ProcessResultsUnits(void)=0;
|
---|
[9110] | 58 | virtual void RequestedOutput(int output_enum,int step,double time)=0;
|
---|
[9817] | 59 |
|
---|
[4076] | 60 | virtual void InputScale(int enum_type,double scale_factor)=0;
|
---|
[9874] | 61 | virtual void GetVectorFromInputs(Vec vector, int name_enum)=0;
|
---|
[11027] | 62 | virtual void GetVectorFromResults(Vec vector,int id,int interp)=0;
|
---|
[5578] | 63 | virtual void InputArtificialNoise(int enum_type,double min,double max)=0;
|
---|
[4880] | 64 | virtual bool InputConvergence(double* eps, int* enums,int num_enums,int* criterionenums,double* criterionvalues,int num_criterionenums)=0;
|
---|
[5518] | 65 | virtual void AverageOntoPartition(Vec partition_contributions,Vec partition_areas,double* vertex_response,double* qmu_part)=0;
|
---|
[7089] | 66 | virtual int* GetHorizontalNeighboorSids(void)=0;
|
---|
[6130] | 67 | virtual double TimeAdapt()=0;
|
---|
[10377] | 68 | virtual void MigrateGroundingLine(double* old_floating_ice,double* sheet_ungrounding)=0;
|
---|
[7323] | 69 | virtual void PotentialSheetUngrounding(Vec potential_sheet_ungrounding)=0;
|
---|
[11527] | 70 | virtual void PositiveDegreeDay(void)=0;
|
---|
[7323] | 71 | virtual int UpdatePotentialSheetUngrounding(double* potential_sheet_ungrounding,Vec vec_nodes_on_iceshelf,double* nodes_on_iceshelf)=0;
|
---|
[10440] | 72 | virtual void ResetCoordinateSystem()=0;
|
---|
[11527] | 73 | virtual void SmearFunction(Vec smearedvector,double (*WeightFunction)(double distance,double radius),double radius)=0;
|
---|
[9817] | 74 |
|
---|
| 75 | #ifdef _HAVE_RESPONSES_
|
---|
| 76 | virtual void MinVel(double* pminvel, bool process_units)=0;
|
---|
| 77 | virtual void MaxVel(double* pmaxvel, bool process_units)=0;
|
---|
| 78 | virtual void MinVx(double* pminvx, bool process_units)=0;
|
---|
| 79 | virtual void MaxVx(double* pmaxvx, bool process_units)=0;
|
---|
| 80 | virtual void MaxAbsVx(double* pmaxabsvx, bool process_units)=0;
|
---|
| 81 | virtual void MinVy(double* pminvy, bool process_units)=0;
|
---|
| 82 | virtual void MaxVy(double* pmaxvy, bool process_units)=0;
|
---|
| 83 | virtual void MaxAbsVy(double* pmaxabsvy, bool process_units)=0;
|
---|
| 84 | virtual void MinVz(double* pminvz, bool process_units)=0;
|
---|
| 85 | virtual void MaxVz(double* pmaxvz, bool process_units)=0;
|
---|
| 86 | virtual void MaxAbsVz(double* pmaxabsvz, bool process_units)=0;
|
---|
| 87 | virtual double MassFlux(double* segment,bool process_units)=0;
|
---|
[10703] | 88 | virtual void ElementResponse(double* presponse,int response_enum,bool process_units)=0;
|
---|
[9880] | 89 | virtual double IceVolume(void)=0;
|
---|
[9206] | 90 | virtual int NodalValue(double* pvalue, int index, int natureofdataenum,bool process_units)=0;
|
---|
[9817] | 91 | #endif
|
---|
[4055] | 92 |
|
---|
[9817] | 93 | #ifdef _HAVE_CONTROL_
|
---|
[11527] | 94 | virtual void Gradj(Vec gradient,int control_type,int control_index)=0;
|
---|
[9817] | 95 | virtual double ThicknessAbsMisfit(bool process_units ,int weight_index)=0;
|
---|
| 96 | virtual double SurfaceAbsVelMisfit(bool process_units ,int weight_index)=0;
|
---|
| 97 | virtual double SurfaceRelVelMisfit(bool process_units ,int weight_index)=0;
|
---|
| 98 | virtual double SurfaceLogVelMisfit(bool process_units ,int weight_index)=0;
|
---|
| 99 | virtual double SurfaceLogVxVyMisfit(bool process_units,int weight_index)=0;
|
---|
| 100 | virtual double SurfaceAverageVelMisfit(bool process_units,int weight_index)=0;
|
---|
| 101 | virtual double ThicknessAbsGradient(bool process_units,int weight_index)=0;
|
---|
| 102 | virtual double RheologyBbarAbsGradient(bool process_units,int weight_index)=0;
|
---|
| 103 | virtual double DragCoefficientAbsGradient(bool process_units,int weight_index)=0;
|
---|
[11527] | 104 | virtual void ControlInputGetGradient(Vec gradient,int enum_type,int control_index)=0;
|
---|
| 105 | virtual void ControlInputSetGradient(double* gradient,int enum_type,int control_index)=0;
|
---|
[9817] | 106 | virtual void ControlInputScaleGradient(int enum_type, double scale)=0;
|
---|
[11527] | 107 | virtual void GetVectorFromControlInputs(Vec gradient,int control_enum,int control_index,const char* data)=0;
|
---|
| 108 | virtual void SetControlInputsFromVector(double* vector,int control_enum,int control_index)=0;
|
---|
[9817] | 109 | virtual void InputControlUpdate(double scalar,bool save_parameter)=0;
|
---|
| 110 | #endif
|
---|
[3683] | 111 | };
|
---|
| 112 | #endif
|
---|