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
RevLine 
[1]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;
[202]21
22 /*diagnostic: */
23 double* vx=NULL;
24 double* vy=NULL;
25 double* vz=NULL;
[586]26
[202]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
[1]34 double* vx_obs=NULL;
35 double* vy_obs=NULL;
36 double* control_parameter=NULL;
[758]37 double* param_g=NULL;
[202]38
39 double* u_g_obs=NULL;
40
[596]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
[586]50 /*prognostic: */
51 double* a_g=NULL;
52 double* h_g=NULL;
53 double* accumulation=NULL;
54 double* thickness=NULL;
55
[202]56 /*diverse: */
[1]57 int numberofnodes;
58 int analysis_type;
[613]59 int sub_analysis_type;
[1]60 int count;
61
62 parameters->FindParam((void*)&analysis_type,"analysis_type");
[613]63 parameters->FindParam((void*)&sub_analysis_type,"sub_analysis_type");
[1]64 count=parameters->Size();
[202]65
66 parameters->FindParam((void*)&numberofnodes,"numberofnodes");
[1]67
[202]68 /*First serialize partition vector: */
[304]69 if(part)VecToMPISerial(&partition,part);
[202]70
[586]71 if ( (analysis_type==ControlAnalysisEnum()) ||
72 (analysis_type==DiagnosticAnalysisEnum()) ||
73 (analysis_type==ThermalAnalysisEnum()) ||
74 (analysis_type==PrognosticAnalysisEnum())
75 ){
[202]76
[209]77 parameters->FindParam((void*)&vx,"vx");
78 parameters->FindParam((void*)&vy,"vy");
79 parameters->FindParam((void*)&vz,"vz");
[202]80
[304]81 if(numberofnodes){
82 u_g=(double*)xcalloc(numberofnodes*3,sizeof(double));
[202]83
[304]84 if(vx){
85 for(i=0;i<numberofnodes;i++){
86 u_g[(int)(3*partition[i]+0)]=vx[i];
87 }
[202]88 }
[304]89 if(vy){
90 for(i=0;i<numberofnodes;i++){
91 u_g[(int)(3*partition[i]+1)]=vy[i];
92 }
[202]93 }
[304]94 if(vz){
95 for(i=0;i<numberofnodes;i++){
96 u_g[(int)(3*partition[i]+2)]=vz[i];
97 }
[202]98 }
99
[304]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);
[202]105
[304]106 /*erase old parameters: */
107 param=(Param*)parameters->FindParamObject("vx");
108 parameters->DeleteObject((Object*)param);
[202]109
[304]110 param=(Param*)parameters->FindParamObject("vy");
111 parameters->DeleteObject((Object*)param);
[202]112
[304]113 param=(Param*)parameters->FindParamObject("vz");
114 parameters->DeleteObject((Object*)param);
115 }
[202]116 }
117
[1]118 if(analysis_type==ControlAnalysisEnum()){
119
[209]120 parameters->FindParam((void*)&vx_obs,"vx_obs");
121 parameters->FindParam((void*)&vy_obs,"vy_obs");
122 parameters->FindParam((void*)&control_parameter,"control_parameter");
[1]123
124 /*Now, from vx_obs and vy_obs, build u_g_obs, correctly partitioned: */
[202]125
[1]126 u_g_obs=(double*)xcalloc(numberofnodes*2,sizeof(double));
[758]127 param_g=(double*)xcalloc(numberofnodes*2,sizeof(double));
[1]128
129 for(i=0;i<numberofnodes;i++){
[202]130 u_g_obs[(int)(2*partition[i]+0)]=vx_obs[i];
131 u_g_obs[(int)(2*partition[i]+1)]=vy_obs[i];
[758]132 param_g[(int)(2*partition[i]+0)]=control_parameter[i];
[1]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++;
[758]142 param= new Param(count,"param_g",DOUBLEVEC);
143 param->SetDoubleVec(param_g,2*numberofnodes);
[1]144 parameters->AddObject(param);
[202]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);
[1]155 }
156
[503]157 if(analysis_type==ThermalAnalysisEnum()){
158
159 parameters->FindParam((void*)&pressure,"pressure");
[596]160
161 /*Now, from temperature and melting, build t_g and m_g, correctly partitioned: */
[503]162 p_g=(double*)xcalloc(numberofnodes,sizeof(double));
163
[596]164 for(i=0;i<numberofnodes;i++){
[503]165 p_g[(int)(partition[i])]= pressure[i];
[596]166 }
[503]167
[596]168 /*Now, create new parameter: */
[503]169 count++;
170 param= new Param(count,"p_g",DOUBLEVEC);
171 param->SetDoubleVec(p_g,numberofnodes);
172 parameters->AddObject(param);
173
[596]174 /*erase old parameter: */
[613]175 param=(Param*)parameters->FindParamObject("pressure");
176 parameters->DeleteObject((Object*)param);
[596]177
[613]178 if (sub_analysis_type==TransientAnalysisEnum()){
[596]179
[613]180 parameters->FindParam((void*)&temperature,"temperature");
181
[713]182 /*Now, from temperature , build t_g , correctly partitioned: */
[613]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 }
[586]199 }
[548]200
[713]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
[586]223 if(analysis_type==PrognosticAnalysisEnum()){
[548]224
[586]225 parameters->FindParam((void*)&melting,"melting");
[596]226 parameters->FindParam((void*)&thickness,"thickness");
227 parameters->FindParam((void*)&accumulation,"accumulation");
[586]228
[596]229 /*Now, from melting accumulation and thickness, build m_g, h_g, correctly partitioned: */
[586]230 m_g=(double*)xcalloc(numberofnodes,sizeof(double));
[596]231 h_g=(double*)xcalloc(numberofnodes,sizeof(double));
232 a_g=(double*)xcalloc(numberofnodes,sizeof(double));
[548]233
[586]234 for(i=0;i<numberofnodes;i++){
235 m_g[(int)(partition[i])]= melting[i];
[596]236 h_g[(int)(partition[i])]= thickness[i];
237 a_g[(int)(partition[i])]= accumulation[i];
[586]238 }
[548]239
[586]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);
[548]245
[586]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: */
[596]257 param=(Param*)parameters->FindParamObject("melting");
[586]258 parameters->DeleteObject((Object*)param);
259
[596]260 param=(Param*)parameters->FindParamObject("thickness");
261 parameters->DeleteObject((Object*)param);
[586]262
[596]263 param=(Param*)parameters->FindParamObject("accumulation");
264 parameters->DeleteObject((Object*)param);
[503]265 }
[1]266 xfree((void**)&partition);
[202]267
268 xfree((void**)&vx);
269 xfree((void**)&vy);
270 xfree((void**)&vz);
271 xfree((void**)&u_g);
272
[1]273 xfree((void**)&vx_obs);
274 xfree((void**)&vy_obs);
[758]275 xfree((void**)&param_g);
[503]276 xfree((void**)&pressure);
[1]277 xfree((void**)&control_parameter);
278 xfree((void**)&u_g_obs);
279 xfree((void**)&p_g);
[586]280 xfree((void**)&a_g);
281 xfree((void**)&h_g);
[552]282 xfree((void**)&m_g);
[1]283}
Note: See TracBrowser for help on using the repository browser.