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

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

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