source: issm/trunk/src/c/ProcessParamsx/ProcessParamsx.cpp@ 758

Last change on this file since 758 was 758, checked in by Mathieu Morlighem, 16 years ago

fixed control method output
File size: 7.4 KB
Line 
1/*!\file ProcessParamsx
2 * \brief: process parameters using partitioning vector
3 */
4
5#include "./ProcessParamsx.h"
6
7#undef __FUNCT__
8#define __FUNCT__ "ProcessParamsx"
9
10#include "../shared/shared.h"
11#include "../include/macros.h"
12#include "../toolkits/toolkits.h"
13#include "../EnumDefinitions/EnumDefinitions.h"
14
15void ProcessParamsx( DataSet* parameters, Vec part){
16
17 int i;
18
19 double* partition=NULL;
20 Param* param=NULL;
21
22 /*diagnostic: */
23 double* vx=NULL;
24 double* vy=NULL;
25 double* vz=NULL;
26
27 double* u_g=NULL;
28
29 /*control: */
30 Vec vec_vx_obs=NULL;
31 Vec vec_vy_obs=NULL;
32 Vec vec_control_parameter=NULL;
33
34 double* vx_obs=NULL;
35 double* vy_obs=NULL;
36 double* control_parameter=NULL;
37 double* param_g=NULL;
38
39 double* u_g_obs=NULL;
40
41 /*thermal*/
42 double* pressure=NULL;
43 double* temperature=NULL;
44 double* melting=NULL;
45
46 double* p_g=NULL;
47 double* t_g=NULL;
48 double* m_g=NULL;
49
50 /*prognostic: */
51 double* a_g=NULL;
52 double* h_g=NULL;
53 double* accumulation=NULL;
54 double* thickness=NULL;
55
56 /*diverse: */
57 int numberofnodes;
58 int analysis_type;
59 int sub_analysis_type;
60 int count;
61
62 parameters->FindParam((void*)&analysis_type,"analysis_type");
63 parameters->FindParam((void*)&sub_analysis_type,"sub_analysis_type");
64 count=parameters->Size();
65
66 parameters->FindParam((void*)&numberofnodes,"numberofnodes");
67
68 /*First serialize partition vector: */
69 if(part)VecToMPISerial(&partition,part);
70
71 if ( (analysis_type==ControlAnalysisEnum()) ||
72 (analysis_type==DiagnosticAnalysisEnum()) ||
73 (analysis_type==ThermalAnalysisEnum()) ||
74 (analysis_type==PrognosticAnalysisEnum())
75 ){
76
77 parameters->FindParam((void*)&vx,"vx");
78 parameters->FindParam((void*)&vy,"vy");
79 parameters->FindParam((void*)&vz,"vz");
80
81 if(numberofnodes){
82 u_g=(double*)xcalloc(numberofnodes*3,sizeof(double));
83
84 if(vx){
85 for(i=0;i<numberofnodes;i++){
86 u_g[(int)(3*partition[i]+0)]=vx[i];
87 }
88 }
89 if(vy){
90 for(i=0;i<numberofnodes;i++){
91 u_g[(int)(3*partition[i]+1)]=vy[i];
92 }
93 }
94 if(vz){
95 for(i=0;i<numberofnodes;i++){
96 u_g[(int)(3*partition[i]+2)]=vz[i];
97 }
98 }
99
100 /*Now, create new parameters: */
101 count++;
102 param= new Param(count,"u_g",DOUBLEVEC);
103 param->SetDoubleVec(u_g,3*numberofnodes);
104 parameters->AddObject(param);
105
106 /*erase old parameters: */
107 param=(Param*)parameters->FindParamObject("vx");
108 parameters->DeleteObject((Object*)param);
109
110 param=(Param*)parameters->FindParamObject("vy");
111 parameters->DeleteObject((Object*)param);
112
113 param=(Param*)parameters->FindParamObject("vz");
114 parameters->DeleteObject((Object*)param);
115 }
116 }
117
118 if(analysis_type==ControlAnalysisEnum()){
119
120 parameters->FindParam((void*)&vx_obs,"vx_obs");
121 parameters->FindParam((void*)&vy_obs,"vy_obs");
122 parameters->FindParam((void*)&control_parameter,"control_parameter");
123
124 /*Now, from vx_obs and vy_obs, build u_g_obs, correctly partitioned: */
125
126 u_g_obs=(double*)xcalloc(numberofnodes*2,sizeof(double));
127 param_g=(double*)xcalloc(numberofnodes*2,sizeof(double));
128
129 for(i=0;i<numberofnodes;i++){
130 u_g_obs[(int)(2*partition[i]+0)]=vx_obs[i];
131 u_g_obs[(int)(2*partition[i]+1)]=vy_obs[i];
132 param_g[(int)(2*partition[i]+0)]=control_parameter[i];
133 }
134
135 /*Now, create new parameters: */
136 count++;
137 param= new Param(count,"u_g_obs",DOUBLEVEC);
138 param->SetDoubleVec(u_g_obs,2*numberofnodes);
139 parameters->AddObject(param);
140
141 count++;
142 param= new Param(count,"param_g",DOUBLEVEC);
143 param->SetDoubleVec(param_g,2*numberofnodes);
144 parameters->AddObject(param);
145
146 /*erase old parameter: */
147 param=(Param*)parameters->FindParamObject("vx_obs");
148 parameters->DeleteObject((Object*)param);
149
150 param=(Param*)parameters->FindParamObject("vy_obs");
151 parameters->DeleteObject((Object*)param);
152
153 param=(Param*)parameters->FindParamObject("control_parameter");
154 parameters->DeleteObject((Object*)param);
155 }
156
157 if(analysis_type==ThermalAnalysisEnum()){
158
159 parameters->FindParam((void*)&pressure,"pressure");
160
161 /*Now, from temperature and melting, build t_g and m_g, correctly partitioned: */
162 p_g=(double*)xcalloc(numberofnodes,sizeof(double));
163
164 for(i=0;i<numberofnodes;i++){
165 p_g[(int)(partition[i])]= pressure[i];
166 }
167
168 /*Now, create new parameter: */
169 count++;
170 param= new Param(count,"p_g",DOUBLEVEC);
171 param->SetDoubleVec(p_g,numberofnodes);
172 parameters->AddObject(param);
173
174 /*erase old parameter: */
175 param=(Param*)parameters->FindParamObject("pressure");
176 parameters->DeleteObject((Object*)param);
177
178 if (sub_analysis_type==TransientAnalysisEnum()){
179
180 parameters->FindParam((void*)&temperature,"temperature");
181
182 /*Now, from temperature , build t_g , correctly partitioned: */
183 t_g=(double*)xcalloc(numberofnodes,sizeof(double));
184
185 for(i=0;i<numberofnodes;i++){
186 t_g[(int)(partition[i])]=temperature[i];
187 }
188
189 /*Now, create new parameter: */
190 count++;
191 param= new Param(count,"t_g",DOUBLEVEC);
192 param->SetDoubleVec(t_g,numberofnodes);
193 parameters->AddObject(param);
194
195 /*erase old parameter: */
196 param=(Param*)parameters->FindParamObject("temperature");
197 parameters->DeleteObject((Object*)param);
198 }
199 }
200
201 if(analysis_type==MeltingAnalysisEnum() && sub_analysis_type==TransientAnalysisEnum()){
202
203 parameters->FindParam((void*)&melting,"melting");
204
205 /*Now, from melting , build m_g , correctly partitioned: */
206 m_g=(double*)xcalloc(numberofnodes,sizeof(double));
207
208 for(i=0;i<numberofnodes;i++){
209 m_g[(int)(partition[i])]=melting[i];
210 }
211
212 /*Now, create new parameter: */
213 count++;
214 param= new Param(count,"m_g",DOUBLEVEC);
215 param->SetDoubleVec(m_g,numberofnodes);
216 parameters->AddObject(param);
217
218 /*erase old parameter: */
219 param=(Param*)parameters->FindParamObject("melting");
220 parameters->DeleteObject((Object*)param);
221 }
222
223 if(analysis_type==PrognosticAnalysisEnum()){
224
225 parameters->FindParam((void*)&melting,"melting");
226 parameters->FindParam((void*)&thickness,"thickness");
227 parameters->FindParam((void*)&accumulation,"accumulation");
228
229 /*Now, from melting accumulation and thickness, build m_g, h_g, correctly partitioned: */
230 m_g=(double*)xcalloc(numberofnodes,sizeof(double));
231 h_g=(double*)xcalloc(numberofnodes,sizeof(double));
232 a_g=(double*)xcalloc(numberofnodes,sizeof(double));
233
234 for(i=0;i<numberofnodes;i++){
235 m_g[(int)(partition[i])]= melting[i];
236 h_g[(int)(partition[i])]= thickness[i];
237 a_g[(int)(partition[i])]= accumulation[i];
238 }
239
240 /*Now, create new parameter: */
241 count++;
242 param= new Param(count,"m_g",DOUBLEVEC);
243 param->SetDoubleVec(m_g,numberofnodes);
244 parameters->AddObject(param);
245
246 count++;
247 param= new Param(count,"h_g",DOUBLEVEC);
248 param->SetDoubleVec(h_g,numberofnodes);
249 parameters->AddObject(param);
250
251 count++;
252 param= new Param(count,"a_g",DOUBLEVEC);
253 param->SetDoubleVec(a_g,numberofnodes);
254 parameters->AddObject(param);
255
256 /*erase old parameter: */
257 param=(Param*)parameters->FindParamObject("melting");
258 parameters->DeleteObject((Object*)param);
259
260 param=(Param*)parameters->FindParamObject("thickness");
261 parameters->DeleteObject((Object*)param);
262
263 param=(Param*)parameters->FindParamObject("accumulation");
264 parameters->DeleteObject((Object*)param);
265 }
266 xfree((void**)&partition);
267
268 xfree((void**)&vx);
269 xfree((void**)&vy);
270 xfree((void**)&vz);
271 xfree((void**)&u_g);
272
273 xfree((void**)&vx_obs);
274 xfree((void**)&vy_obs);
275 xfree((void**)&param_g);
276 xfree((void**)&pressure);
277 xfree((void**)&control_parameter);
278 xfree((void**)&u_g_obs);
279 xfree((void**)&p_g);
280 xfree((void**)&a_g);
281 xfree((void**)&h_g);
282 xfree((void**)&m_g);
283}
Note: See TracBrowser for help on using the repository browser.