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;
|
---|
17 | class Patch;
|
---|
18 |
|
---|
19 | #include "../../toolkits/toolkits.h"
|
---|
20 | /*}}}*/
|
---|
21 |
|
---|
22 | class Element: public Object,public Update{
|
---|
23 |
|
---|
24 | public:
|
---|
25 |
|
---|
26 | virtual ~Element(){};
|
---|
27 |
|
---|
28 | virtual void Configure(Elements* elements,Loads* loads,DataSet* nodes,Materials* materials,Parameters* parameters)=0;
|
---|
29 | virtual void SetCurrentConfiguration(Elements* elements,Loads* loads,DataSet* nodes,Materials* materials,Parameters* parameters)=0;
|
---|
30 | virtual void CreateKMatrix(Mat Kff, Mat Kfs,Vec df)=0;
|
---|
31 | virtual void CreatePVector(Vec pf)=0;
|
---|
32 | virtual void CreateJacobianMatrix(Mat Jff)=0;
|
---|
33 | virtual void GetSolutionFromInputs(Vec solution)=0;
|
---|
34 | virtual int GetNodeIndex(Node* node)=0;
|
---|
35 | virtual int Sid()=0;
|
---|
36 | virtual bool IsFloating()=0;
|
---|
37 | virtual bool IsNodeOnShelf()=0;
|
---|
38 | virtual bool IsNodeOnShelfFromFlags(double* flags)=0;
|
---|
39 | virtual bool IsOnBed()=0;
|
---|
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;
|
---|
43 |
|
---|
44 | virtual double SurfaceArea(void)=0;
|
---|
45 | virtual void InputDepthAverageAtBase(int enum_type,int average_enum_type,int object_enum)=0;
|
---|
46 | virtual void ComputeBasalStress(Vec sigma_b)=0;
|
---|
47 | virtual void ComputeStrainRate(Vec eps)=0;
|
---|
48 | virtual void PatchSize(int* pnumrows, int* pnumvertices,int* pnumnodes)=0;
|
---|
49 | virtual void PatchFill(int* pcount, Patch* patch)=0;
|
---|
50 | virtual void ListResultsInfo(int** results_enums,int** results_size,double** results_times,int** results_steps,int* num_results)=0;
|
---|
51 | virtual void DeleteResults(void)=0;
|
---|
52 | virtual void Update(int index, IoModel* iomodel,int analysis_counter,int analysis_type)=0;
|
---|
53 | virtual void InputToResult(int enum_type,int step,double time)=0;
|
---|
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;
|
---|
57 | virtual void ProcessResultsUnits(void)=0;
|
---|
58 | virtual void RequestedOutput(int output_enum,int step,double time)=0;
|
---|
59 |
|
---|
60 | virtual void InputScale(int enum_type,double scale_factor)=0;
|
---|
61 | virtual void GetVectorFromInputs(Vec vector, int name_enum)=0;
|
---|
62 | virtual void GetVectorFromResults(Vec vector,int id,int interp)=0;
|
---|
63 | virtual void InputArtificialNoise(int enum_type,double min,double max)=0;
|
---|
64 | virtual bool InputConvergence(double* eps, int* enums,int num_enums,int* criterionenums,double* criterionvalues,int num_criterionenums)=0;
|
---|
65 | virtual void AverageOntoPartition(Vec partition_contributions,Vec partition_areas,double* vertex_response,double* qmu_part)=0;
|
---|
66 | virtual int* GetHorizontalNeighboorSids(void)=0;
|
---|
67 | virtual double TimeAdapt()=0;
|
---|
68 | virtual void MigrateGroundingLine(double* old_floating_ice,double* sheet_ungrounding)=0;
|
---|
69 | virtual void PotentialSheetUngrounding(Vec potential_sheet_ungrounding)=0;
|
---|
70 | virtual void PositiveDegreeDay(void)=0;
|
---|
71 | virtual int UpdatePotentialSheetUngrounding(double* potential_sheet_ungrounding,Vec vec_nodes_on_iceshelf,double* nodes_on_iceshelf)=0;
|
---|
72 | virtual void ResetCoordinateSystem()=0;
|
---|
73 | virtual void SmearFunction(Vec smearedvector,double (*WeightFunction)(double distance,double radius),double radius)=0;
|
---|
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;
|
---|
88 | virtual void ElementResponse(double* presponse,int response_enum,bool process_units)=0;
|
---|
89 | virtual double IceVolume(void)=0;
|
---|
90 | virtual int NodalValue(double* pvalue, int index, int natureofdataenum,bool process_units)=0;
|
---|
91 | #endif
|
---|
92 |
|
---|
93 | #ifdef _HAVE_CONTROL_
|
---|
94 | virtual void Gradj(Vec gradient,int control_type,int control_index)=0;
|
---|
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;
|
---|
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;
|
---|
106 | virtual void ControlInputScaleGradient(int enum_type, double scale)=0;
|
---|
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;
|
---|
109 | virtual void InputControlUpdate(double scalar,bool save_parameter)=0;
|
---|
110 | #endif
|
---|
111 | };
|
---|
112 | #endif
|
---|