1 | /*!\file ControlInput.c
|
---|
2 | * \brief: implementation of the ControlInput object
|
---|
3 | */
|
---|
4 |
|
---|
5 | #ifdef HAVE_CONFIG_H
|
---|
6 | #include <config.h>
|
---|
7 | #else
|
---|
8 | #error "Cannot compile with HAVE_CONFIG_H symbol! run configure first!"
|
---|
9 | #endif
|
---|
10 |
|
---|
11 | #include <stdio.h>
|
---|
12 | #include <string.h>
|
---|
13 | #include "../../classes.h"
|
---|
14 | #include "../../../EnumDefinitions/EnumDefinitions.h"
|
---|
15 | #include "../../../shared/shared.h"
|
---|
16 | #include "../../../Container/Container.h"
|
---|
17 | #include "../../../include/include.h"
|
---|
18 |
|
---|
19 | /*ControlInput constructors and destructor*/
|
---|
20 | /*FUNCTION ControlInput::ControlInput(){{{*/
|
---|
21 | ControlInput::ControlInput(){
|
---|
22 | control_id = 0;
|
---|
23 | values = NULL;
|
---|
24 | savedvalues = NULL;
|
---|
25 | minvalues = NULL;
|
---|
26 | maxvalues = NULL;
|
---|
27 | gradient = NULL;
|
---|
28 | }
|
---|
29 | /*}}}*/
|
---|
30 | /*FUNCTION ControlInput::ControlInput(int enum_type,int enum_input,IssmDouble* pvalues,IssmDouble* pmin,IssmDouble* pmax,int id){{{*/
|
---|
31 | ControlInput::ControlInput(int in_enum_type,int enum_input,IssmDouble* pvalues,IssmDouble* pmin,IssmDouble* pmax,int id){
|
---|
32 |
|
---|
33 | control_id=id;
|
---|
34 | enum_type=in_enum_type;
|
---|
35 |
|
---|
36 | switch(enum_input){
|
---|
37 | case TriaP1InputEnum:
|
---|
38 | values =new TriaP1Input(enum_type,pvalues);
|
---|
39 | savedvalues=new TriaP1Input(enum_type,pvalues);
|
---|
40 | minvalues =new TriaP1Input(enum_type,pmin);
|
---|
41 | maxvalues =new TriaP1Input(enum_type,pmax);
|
---|
42 | break;
|
---|
43 | case PentaP1InputEnum:
|
---|
44 | values =new PentaP1Input(enum_type,pvalues);
|
---|
45 | savedvalues=new PentaP1Input(enum_type,pvalues);
|
---|
46 | minvalues =new PentaP1Input(enum_type,pmin);
|
---|
47 | maxvalues =new PentaP1Input(enum_type,pmax);
|
---|
48 | break;
|
---|
49 | default:
|
---|
50 | _error_("Input of Enum " << EnumToStringx(enum_input) << " not supported yet by ControlInput");
|
---|
51 | }
|
---|
52 | gradient =NULL;
|
---|
53 | }
|
---|
54 | /*}}}*/
|
---|
55 | /*FUNCTION ControlInput::~ControlInput(){{{*/
|
---|
56 | ControlInput::~ControlInput(){
|
---|
57 | delete values;
|
---|
58 | delete savedvalues;
|
---|
59 | delete minvalues;
|
---|
60 | delete maxvalues;
|
---|
61 | delete gradient;
|
---|
62 | }
|
---|
63 | /*}}}*/
|
---|
64 |
|
---|
65 | /*Object virtual functions definitions:*/
|
---|
66 | /*FUNCTION ControlInput::Echo {{{*/
|
---|
67 | void ControlInput::Echo(void){
|
---|
68 | this->DeepEcho();
|
---|
69 | }
|
---|
70 | /*}}}*/
|
---|
71 | /*FUNCTION ControlInput::DeepEcho{{{*/
|
---|
72 | void ControlInput::DeepEcho(void){
|
---|
73 |
|
---|
74 | _printLine_("ControlInput:");
|
---|
75 | _printLine_(" enum: " << this->enum_type << " (" << EnumToStringx(this->enum_type) << ")");
|
---|
76 | _printLine_("---values: "); if (values) values->Echo();
|
---|
77 | _printLine_("---savedvalues: ");if (savedvalues) savedvalues->Echo();
|
---|
78 | _printLine_("---minvalues: "); if (minvalues) minvalues->Echo();
|
---|
79 | _printLine_("---maxvalues: "); if (maxvalues) maxvalues->Echo();
|
---|
80 | _printLine_("---gradient: "); if (gradient) gradient->Echo();
|
---|
81 | }
|
---|
82 | /*}}}*/
|
---|
83 | /*FUNCTION ControlInput::Id{{{*/
|
---|
84 | int ControlInput::Id(void){ return -1; }
|
---|
85 | /*}}}*/
|
---|
86 | /*FUNCTION ControlInput::MyRank{{{*/
|
---|
87 | int ControlInput::MyRank(void){
|
---|
88 | extern int my_rank;
|
---|
89 | return my_rank;
|
---|
90 | }
|
---|
91 | /*}}}*/
|
---|
92 | /*FUNCTION ControlInput::ObjectEnum{{{*/
|
---|
93 | int ControlInput::ObjectEnum(void){
|
---|
94 |
|
---|
95 | return ControlInputEnum;
|
---|
96 |
|
---|
97 | }
|
---|
98 | /*}}}*/
|
---|
99 | /*FUNCTION ControlInput::copy{{{*/
|
---|
100 | Object* ControlInput::copy() {
|
---|
101 |
|
---|
102 | ControlInput* output=NULL;
|
---|
103 |
|
---|
104 | output = new ControlInput();
|
---|
105 | output->enum_type=this->enum_type;
|
---|
106 | output->control_id=this->control_id;
|
---|
107 |
|
---|
108 | if(values) output->values=(Input*)this->values->copy();
|
---|
109 | if(savedvalues) output->savedvalues=(Input*)this->savedvalues->copy();
|
---|
110 | if(minvalues) output->minvalues=(Input*)this->minvalues->copy();
|
---|
111 | if(maxvalues) output->maxvalues=(Input*)this->maxvalues->copy();
|
---|
112 | if(gradient) output->gradient=(Input*)this->gradient->copy();
|
---|
113 |
|
---|
114 | return output;
|
---|
115 | }
|
---|
116 | /*}}}*/
|
---|
117 |
|
---|
118 | /*ControlInput management*/
|
---|
119 | /*FUNCTION ControlInput::InstanceEnum{{{*/
|
---|
120 | int ControlInput::InstanceEnum(void){
|
---|
121 |
|
---|
122 | return this->enum_type;
|
---|
123 |
|
---|
124 | }
|
---|
125 | /*}}}*/
|
---|
126 |
|
---|
127 | /*Object functions*/
|
---|
128 | /*FUNCTION ControlInput::Constrain(){{{*/
|
---|
129 | void ControlInput::Constrain(void){
|
---|
130 |
|
---|
131 | Input* newvalues=NULL;
|
---|
132 |
|
---|
133 | newvalues=this->values->PointwiseMin(maxvalues);
|
---|
134 | delete values; this->values=newvalues;
|
---|
135 | newvalues=this->values->PointwiseMax(minvalues);
|
---|
136 | delete values; this->values=newvalues;
|
---|
137 | }/*}}}*/
|
---|
138 | /*FUNCTION ControlInput::Constrain(IssmDouble min, IssmDouble max){{{*/
|
---|
139 | void ControlInput::Constrain(IssmDouble min, IssmDouble max){
|
---|
140 | values->Constrain(min,max);
|
---|
141 | }/*}}}*/
|
---|
142 | /*FUNCTION ControlInput::Extrude{{{*/
|
---|
143 | void ControlInput::Extrude(void){
|
---|
144 | values->Extrude();
|
---|
145 | savedvalues->Extrude();
|
---|
146 | //gradient->Extrude();
|
---|
147 | }/*}}}*/
|
---|
148 | /*FUNCTION ControlInput::GetGradient{{{*/
|
---|
149 | void ControlInput::GetGradient(Vector* gradient_vec,int* doflist){
|
---|
150 | if(gradient) gradient->GetVectorFromInputs(gradient_vec,doflist);
|
---|
151 | }/*}}}*/
|
---|
152 | /*FUNCTION ControlInput::ScaleGradient{{{*/
|
---|
153 | void ControlInput::ScaleGradient(IssmDouble scaling_factor){
|
---|
154 | if(!gradient) _error_("Gradient of ControlInput " << EnumToStringx(enum_type) << " not found");
|
---|
155 | gradient->Scale(scaling_factor);
|
---|
156 | }/*}}}*/
|
---|
157 | /*FUNCTION ControlInput::SetGradient{{{*/
|
---|
158 | void ControlInput::SetGradient(Input* gradient_in){
|
---|
159 |
|
---|
160 | /*Get enum for current gradient*/
|
---|
161 | switch(this->control_id){
|
---|
162 | case 1:
|
---|
163 | gradient_in->ChangeEnum(Gradient1Enum);
|
---|
164 | break;
|
---|
165 | case 2:
|
---|
166 | gradient_in->ChangeEnum(Gradient2Enum);
|
---|
167 | break;
|
---|
168 | case 3:
|
---|
169 | gradient_in->ChangeEnum(Gradient3Enum);
|
---|
170 | break;
|
---|
171 | default:
|
---|
172 | _error_("more than 3 controls not implemented yet (Gradient " << this->control_id << " was requested). EnumDefinitions.h needs to be updated.");
|
---|
173 | }
|
---|
174 |
|
---|
175 | /*Delete old gradient and assign new gradient*/
|
---|
176 | if(gradient) delete gradient;
|
---|
177 | gradient=gradient_in;
|
---|
178 |
|
---|
179 | }/*}}}*/
|
---|
180 | /*FUNCTION ControlInput::SetInput{{{*/
|
---|
181 | void ControlInput::SetInput(Input* in_input){
|
---|
182 |
|
---|
183 | delete values; this->values=in_input;
|
---|
184 | this->SaveValue(); //because this is what SpawnResult saves FIXME
|
---|
185 |
|
---|
186 | }/*}}}*/
|
---|
187 | /*FUNCTION ControlInput::SpawnResult{{{*/
|
---|
188 | ElementResult* ControlInput::SpawnResult(int step, IssmDouble time){
|
---|
189 | return savedvalues->SpawnResult(step,time);
|
---|
190 | }/*}}}*/
|
---|
191 | /*FUNCTION ControlInput::SpawnTriaInput{{{*/
|
---|
192 | Input* ControlInput::SpawnTriaInput(int* indices){
|
---|
193 | return values->SpawnTriaInput(indices);
|
---|
194 | }/*}}}*/
|
---|
195 | /*FUNCTION ControlInput::SpawnGradient{{{*/
|
---|
196 | ElementResult* ControlInput::SpawnGradient(int step, IssmDouble time){
|
---|
197 | _assert_(gradient);
|
---|
198 | return gradient->SpawnResult(step,time);
|
---|
199 | }/*}}}*/
|
---|
200 | /*FUNCTION ControlInput::GetVectorFromInputs(Vector* vector,int* doflist){{{*/
|
---|
201 | void ControlInput::GetVectorFromInputs(Vector* vector,int* doflist){
|
---|
202 | values->GetVectorFromInputs(vector,doflist);
|
---|
203 | }/*}}}*/
|
---|
204 | /*FUNCTION ControlInput::GetVectorFromInputs(Vector* vector,int* doflist,const char* data){{{*/
|
---|
205 | void ControlInput::GetVectorFromInputs(Vector* vector,int* doflist,const char* data){
|
---|
206 | if(strcmp(data,"value")==0){
|
---|
207 | _assert_(values);
|
---|
208 | values->GetVectorFromInputs(vector,doflist);
|
---|
209 | }
|
---|
210 | else if (strcmp(data,"lowerbound")==0){
|
---|
211 | _assert_(minvalues);
|
---|
212 | minvalues->GetVectorFromInputs(vector,doflist);
|
---|
213 | }
|
---|
214 | else if (strcmp(data,"upperbound")==0){
|
---|
215 | _assert_(maxvalues);
|
---|
216 | maxvalues->GetVectorFromInputs(vector,doflist);
|
---|
217 | }
|
---|
218 | else if (strcmp(data,"gradient")==0){
|
---|
219 | _assert_(gradient);
|
---|
220 | gradient->GetVectorFromInputs(vector,doflist);
|
---|
221 | }
|
---|
222 | else{
|
---|
223 | _error_("Data " << data << " not supported yet");
|
---|
224 | }
|
---|
225 | }/*}}}*/
|
---|
226 | /*FUNCTION ControlInput::GetInputAverage(IssmDouble* pvalue){{{*/
|
---|
227 | void ControlInput::GetInputAverage(IssmDouble* pvalue){
|
---|
228 | values->GetInputAverage(pvalue);
|
---|
229 | }/*}}}*/
|
---|
230 | /*FUNCTION ControlInput::GetInputValue(bool* pvalue){{{*/
|
---|
231 | void ControlInput::GetInputValue(bool* pvalue){
|
---|
232 | values->GetInputValue(pvalue);
|
---|
233 | }/*}}}*/
|
---|
234 | /*FUNCTION ControlInput::GetInputValue(int* pvalue){{{*/
|
---|
235 | void ControlInput::GetInputValue(int* pvalue){
|
---|
236 | values->GetInputValue(pvalue);
|
---|
237 | }/*}}}*/
|
---|
238 | /*FUNCTION ControlInput::GetInputValue(IssmDouble* pvalue){{{*/
|
---|
239 | void ControlInput::GetInputValue(IssmDouble* pvalue){
|
---|
240 | values->GetInputValue(pvalue);
|
---|
241 | }/*}}}*/
|
---|
242 | /*FUNCTION ControlInput::GetInputValue(IssmDouble* pvalue){{{*/
|
---|
243 | void ControlInput::GetInputValue(IssmDouble* pvalue,GaussTria* gauss){
|
---|
244 | values->GetInputValue(pvalue,gauss);
|
---|
245 | }/*}}}*/
|
---|
246 | /*FUNCTION ControlInput::GetInputValue(IssmDouble* pvalue){{{*/
|
---|
247 | void ControlInput::GetInputValue(IssmDouble* pvalue,GaussPenta* gauss){
|
---|
248 | values->GetInputValue(pvalue,gauss);
|
---|
249 | }/*}}}*/
|
---|
250 | /*FUNCTION ControlInput::GetInputDerivativeValue(IssmDouble* derivativevalues, IssmDouble* xyz_list, GaussTria* gauss){{{*/
|
---|
251 | void ControlInput::GetInputDerivativeValue(IssmDouble* derivativevalues, IssmDouble* xyz_list, GaussTria* gauss){
|
---|
252 | values->GetInputDerivativeValue(derivativevalues,xyz_list,gauss);
|
---|
253 | }/*}}}*/
|
---|
254 | /*FUNCTION ControlInput::GetInputDerivativeValue(IssmDouble* derivativevalues, IssmDouble* xyz_list, GaussPenta* gauss){{{*/
|
---|
255 | void ControlInput::GetInputDerivativeValue(IssmDouble* derivativevalues, IssmDouble* xyz_list, GaussPenta* gauss){
|
---|
256 | values->GetInputDerivativeValue(derivativevalues,xyz_list,gauss);
|
---|
257 | }/*}}}*/
|
---|
258 | /*FUNCTION ControlInput::SaveValue{{{*/
|
---|
259 | void ControlInput::SaveValue(void){
|
---|
260 | if(!values) _error_("Values of " << EnumToStringx(this->enum_type) << " not found");
|
---|
261 |
|
---|
262 | if(savedvalues) delete this->savedvalues;
|
---|
263 | this->savedvalues=(Input*)this->values->copy();
|
---|
264 | }/*}}}*/
|
---|
265 | /*FUNCTION ControlInput::UpdateValue{{{*/
|
---|
266 | void ControlInput::UpdateValue(IssmDouble scalar){
|
---|
267 | if(!gradient) _error_("Gradient of " << EnumToStringx(this->enum_type) << " not found");
|
---|
268 | if(!savedvalues) _error_("Values of " << EnumToStringx(this->enum_type) << " not found");
|
---|
269 |
|
---|
270 | if(values) delete this->values;
|
---|
271 | this->values=(Input*)this->savedvalues->copy();
|
---|
272 | this->values->AXPY(gradient,scalar);
|
---|
273 | }/*}}}*/
|
---|
274 | /*FUNCTION ControlInput::VerticallyIntegrate{{{*/
|
---|
275 | void ControlInput::VerticallyIntegrate(Input* thickness_input){
|
---|
276 | values->VerticallyIntegrate(thickness_input);
|
---|
277 | }/*}}}*/
|
---|
278 | /*FUNCTION ControlInput::Configure{{{*/
|
---|
279 | void ControlInput::Configure(Parameters* parameters){
|
---|
280 | /*do nothing: */
|
---|
281 | }
|
---|
282 | /*}}}*/
|
---|