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

Last change on this file since 2892 was 2892, checked in by Eric.Larour, 15 years ago

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