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