source: issm/trunk/src/c/ModelProcessorx/Control/CreateParametersControl.cpp@ 3332

Last change on this file since 3332 was 3332, checked in by Mathieu Morlighem, 15 years ago

Do not use throw ErrorException -> ISSMERROR macro
Removed all FUNCT definitions (now useless)
File size: 6.3 KB
Line 
1/*!\file: CreateParametersControl.cpp
2 * \brief driver for creating parameters dataset, for control analysis.
3 */
4
5#include "../../DataSet/DataSet.h"
6#include "../../toolkits/toolkits.h"
7#include "../../EnumDefinitions/EnumDefinitions.h"
8#include "../../objects/objects.h"
9#include "../../shared/shared.h"
10#include "../../include/macros.h"
11#include "../IoModel.h"
12
13void CreateParametersControl(DataSet** pparameters,IoModel* iomodel,ConstDataHandle iomodel_handle){
14
15 int i;
16
17 DataSet* parameters=NULL;
18 Param* param = NULL;
19 int count;
20
21 double* fit=NULL;
22 double* weights=NULL;
23 double* cm_jump=NULL;
24 double* optscal=NULL;
25 double* maxiter=NULL;
26 double* control_parameter=NULL;
27 double* param_g=NULL;
28
29 double* vx_obs=NULL;
30 double* vy_obs=NULL;
31 double* u_g_obs=NULL;
32
33 double* vx=NULL;
34 double* vy=NULL;
35 double* vz=NULL;
36
37 /*Get parameters: */
38 parameters=*pparameters;
39 count=parameters->Size();
40
41 //control analysis?
42 count++;
43 param= new Param(count,"control_analysis",DOUBLE);
44 param->SetDouble(iomodel->control_analysis);
45 parameters->AddObject(param);
46
47 if(iomodel->control_analysis){
48
49 /*control_type: */
50 count++;
51 param= new Param(count,"control_type",STRING);
52 param->SetString(iomodel->control_type);
53 parameters->AddObject(param);
54
55 /*extrude_param: */
56 count++;
57 param= new Param(count,"extrude_param",DOUBLE);
58 if (strcmp(iomodel->control_type,"drag")==0) param->SetDouble(0);
59 else if (strcmp(iomodel->control_type,"B")==0) param->SetDouble(1);
60 else ISSMERROR(exprintf("control_type %s not supported yet!",iomodel->control_type));
61 parameters->AddObject(param);
62
63 /*control_steady: */
64 count++;
65 param= new Param(count,"control_steady",DOUBLE);
66 param->SetDouble(0);
67 parameters->AddObject(param);
68
69 /*nsteps: */
70 count++;
71 param= new Param(count,"nsteps",DOUBLE);
72 param->SetDouble(iomodel->nsteps);
73 parameters->AddObject(param);
74
75 /*tolx: */
76 count++;
77 param= new Param(count,"tolx",DOUBLE);
78 param->SetDouble(iomodel->tolx);
79 parameters->AddObject(param);
80
81 /*eps_cm: */
82 count++;
83 param= new Param(count,"eps_cm",DOUBLE);
84 param->SetDouble(iomodel->eps_cm);
85 parameters->AddObject(param);
86
87 /*epsvel: */
88 count++;
89 param= new Param(count,"epsvel",DOUBLE);
90 param->SetDouble(iomodel->epsvel);
91 parameters->AddObject(param);
92
93 /*meanvel: */
94 count++;
95 param= new Param(count,"meanvel",DOUBLE);
96 param->SetDouble(iomodel->meanvel);
97 parameters->AddObject(param);
98
99 /*cm_noisedmp: */
100 count++;
101 param= new Param(count,"cm_noisedmp",DOUBLE);
102 param->SetDouble(iomodel->cm_noisedmp);
103 parameters->AddObject(param);
104
105 /*cm_mindmp_value: */
106 count++;
107 param= new Param(count,"cm_mindmp_value",DOUBLE);
108 param->SetDouble(iomodel->cm_mindmp_value);
109 parameters->AddObject(param);
110
111 /*cm_mindmp_slope: */
112 count++;
113 param= new Param(count,"cm_mindmp_slope",DOUBLE);
114 param->SetDouble(iomodel->cm_mindmp_slope);
115 parameters->AddObject(param);
116
117 /*cm_maxdmp_value: */
118 count++;
119 param= new Param(count,"cm_maxdmp_value",DOUBLE);
120 param->SetDouble(iomodel->cm_maxdmp_value);
121 parameters->AddObject(param);
122
123 /*cm_maxdmp_slope: */
124 count++;
125 param= new Param(count,"cm_maxdmp_slope",DOUBLE);
126 param->SetDouble(iomodel->cm_maxdmp_slope);
127 parameters->AddObject(param);
128
129 /*cm_min: */
130 count++;
131 param= new Param(count,"cm_min",DOUBLE);
132 param->SetDouble(iomodel->cm_min);
133 parameters->AddObject(param);
134
135 /*cm_max: */
136 count++;
137 param= new Param(count,"cm_max",DOUBLE);
138 param->SetDouble(iomodel->cm_max);
139 parameters->AddObject(param);
140
141 /*cm_gradient: */
142 count++;
143 param= new Param(count,"cm_gradient",DOUBLE);
144 param->SetDouble(iomodel->cm_gradient);
145 parameters->AddObject(param);
146
147 /*Now, recover fit, optscal and maxiter as vectors: */
148 IoModelFetchData(&iomodel->fit,NULL,NULL,iomodel_handle,"fit");
149 IoModelFetchData(&iomodel->cm_jump,NULL,NULL,iomodel_handle,"cm_jump");
150 IoModelFetchData(&iomodel->optscal,NULL,NULL,iomodel_handle,"optscal");
151 IoModelFetchData(&iomodel->maxiter,NULL,NULL,iomodel_handle,"maxiter");
152
153 count++;
154 param= new Param(count,"fit",DOUBLEVEC);
155 param->SetDoubleVec(iomodel->fit,iomodel->nsteps);
156 parameters->AddObject(param);
157
158 count++;
159 param= new Param(count,"cm_jump",DOUBLEVEC);
160 param->SetDoubleVec(iomodel->cm_jump,iomodel->nsteps);
161 parameters->AddObject(param);
162
163 count++;
164 param= new Param(count,"optscal",DOUBLEVEC);
165 param->SetDoubleVec(iomodel->optscal,iomodel->nsteps);
166 parameters->AddObject(param);
167
168 count++;
169 param= new Param(count,"maxiter",DOUBLEVEC);
170 param->SetDoubleVec(iomodel->maxiter,iomodel->nsteps);
171 parameters->AddObject(param);
172
173 xfree((void**)&iomodel->fit);
174 xfree((void**)&iomodel->cm_jump);
175 xfree((void**)&iomodel->optscal);
176 xfree((void**)&iomodel->maxiter);
177
178 /*Now, recover weights as vectors: */
179 IoModelFetchData(&iomodel->weights,NULL,NULL,iomodel_handle,"weights");
180
181 count++;
182 param= new Param(count,"weights",DOUBLEVEC);
183 param->SetDoubleVec(iomodel->weights,iomodel->numberofnodes);
184 parameters->AddObject(param);
185
186 xfree((void**)&iomodel->weights);
187
188 /*Get vx, vx_obs, vy, vy_obs, and the parameter value: */
189 IoModelFetchData(&vx,NULL,NULL,iomodel_handle,"vx");
190 IoModelFetchData(&vy,NULL,NULL,iomodel_handle,"vy");
191 IoModelFetchData(&vz,NULL,NULL,iomodel_handle,"vz");
192 IoModelFetchData(&vx_obs,NULL,NULL,iomodel_handle,"vx_obs");
193 IoModelFetchData(&vy_obs,NULL,NULL,iomodel_handle,"vy_obs");
194 IoModelFetchData(&control_parameter,NULL,NULL,iomodel_handle,iomodel->control_type);
195
196 u_g_obs=(double*)xcalloc(iomodel->numberofnodes*2,sizeof(double));
197 if(vx_obs)for(i=0;i<iomodel->numberofnodes;i++)u_g_obs[2*i+0]=vx_obs[i]/iomodel->yts;
198 if(vy_obs)for(i=0;i<iomodel->numberofnodes;i++)u_g_obs[2*i+1]=vy_obs[i]/iomodel->yts;
199
200 count++;
201 param= new Param(count,"u_g_obs",DOUBLEVEC);
202 param->SetDoubleVec(u_g_obs,2*iomodel->numberofnodes,2);
203 parameters->AddObject(param);
204
205 param_g=(double*)xcalloc(iomodel->numberofnodes,sizeof(double));
206 for(i=0;i<iomodel->numberofnodes;i++)param_g[i]=control_parameter[i];
207
208 count++;
209 param= new Param(count,"param_g",DOUBLEVEC);
210 param->SetDoubleVec(param_g,iomodel->numberofnodes,1);
211 parameters->AddObject(param);
212
213 xfree((void**)&vx);
214 xfree((void**)&vy);
215 xfree((void**)&vz);
216 xfree((void**)&vx_obs);
217 xfree((void**)&vy_obs);
218 xfree((void**)&u_g_obs);
219 xfree((void**)&param_g);
220 xfree((void**)&control_parameter);
221 }
222
223 /*Assign output pointer: */
224 *pparameters=parameters;
225}
Note: See TracBrowser for help on using the repository browser.