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

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

Vector parallel fetch not supported. Had to go around it
File size: 3.8 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 "../Model.h"
14
15void CreateParametersControl(DataSet* parameters,Model* model,ConstDataHandle model_handle,int* pcount){
16
17 int i;
18
19 Param* param = NULL;
20 int count;
21 int analysis_type;
22 int numberofdofspernode;
23
24 double* fit=NULL;
25 double* optscal=NULL;
26 double* maxiter=NULL;
27 double* control_parameter=NULL;
28 double* vx_obs=NULL;
29 double* vy_obs=NULL;
30
31 /*Get counter : */
32 count=*pcount;
33
34 /*control_type: */
35 count++;
36 param= new Param(count,"control_type",STRING);
37 param->SetString(model->control_type);
38 parameters->AddObject(param);
39
40
41 /*nsteps: */
42 count++;
43 param= new Param(count,"nsteps",INTEGER);
44 param->SetInteger(model->nsteps);
45 parameters->AddObject(param);
46
47 /*tolx: */
48 count++;
49 param= new Param(count,"tolx",DOUBLE);
50 param->SetDouble(model->tolx);
51 parameters->AddObject(param);
52
53 /*mincontrolconstraint: */
54 count++;
55 param= new Param(count,"mincontrolconstraint",DOUBLE);
56 param->SetDouble(model->mincontrolconstraint);
57 parameters->AddObject(param);
58
59 /*maxcontrolconstraint: */
60 count++;
61 param= new Param(count,"maxcontrolconstraint",DOUBLE);
62 param->SetDouble(model->maxcontrolconstraint);
63 parameters->AddObject(param);
64
65 /*epsvel: */
66 count++;
67 param= new Param(count,"epsvel",DOUBLE);
68 param->SetDouble(model->epsvel);
69 parameters->AddObject(param);
70
71 /*meanvel: */
72 count++;
73 param= new Param(count,"meanvel",DOUBLE);
74 param->SetDouble(model->meanvel);
75 parameters->AddObject(param);
76
77 /*Now, recover fit, optscal and maxiter as vectors: */
78 ModelFetchData((void**)&model->fit,NULL,NULL,model_handle,"fit","Matrix","Mat");
79 ModelFetchData((void**)&model->optscal,NULL,NULL,model_handle,"optscal","Matrix","Mat");
80 ModelFetchData((void**)&model->maxiter,NULL,NULL,model_handle,"maxiter","Matrix","Mat");
81
82 count++;
83 param= new Param(count,"fit",DOUBLEVEC);
84 param->SetDoubleVec(model->fit,model->nsteps);
85 parameters->AddObject(param);
86
87 count++;
88 param= new Param(count,"optscal",DOUBLEVEC);
89 param->SetDoubleVec(model->optscal,model->nsteps);
90 parameters->AddObject(param);
91
92 count++;
93 param= new Param(count,"maxiter",DOUBLEVEC);
94 param->SetDoubleVec(model->maxiter,model->nsteps);
95 parameters->AddObject(param);
96
97 xfree((void**)&model->fit);
98 xfree((void**)&model->optscal);
99 xfree((void**)&model->maxiter);
100
101 /*Get vx_obs, vy_obs, and the parameter value: */
102 ModelFetchData((void**)&vx_obs,NULL,NULL,model_handle,"vx_obs","Matrix","Mat");
103 ModelFetchData((void**)&vy_obs,NULL,NULL,model_handle,"vy_obs","Matrix","Mat");
104 ModelFetchData((void**)&control_parameter,NULL,NULL,model_handle,model->control_type,"Matrix","Mat");
105
106 if(vx_obs)for(i=0;i<model->numberofnodes;i++)vx_obs[i]=vx_obs[i]/model->yts;
107 if(vy_obs)for(i=0;i<model->numberofnodes;i++)vy_obs[i]=vy_obs[i]/model->yts;
108
109 count++;
110 param= new Param(count,"vx_obs",DOUBLEVEC);
111 if(vx_obs) param->SetDoubleVec(vx_obs,model->numberofnodes);
112 else param->SetDoubleVec(vx_obs,0);
113 parameters->AddObject(param);
114
115 count++;
116 param= new Param(count,"vy_obs",DOUBLEVEC);
117 if(vy_obs) param->SetDoubleVec(vy_obs,model->numberofnodes);
118 else param->SetDoubleVec(vy_obs,0);
119 parameters->AddObject(param);
120
121 count++;
122 param= new Param(count,"control_parameter",DOUBLEVEC);
123 if(control_parameter) param->SetDoubleVec(control_parameter,model->numberofnodes);
124 else param->SetDoubleVec(control_parameter,0);
125 parameters->AddObject(param);
126
127
128 xfree((void**)&vx_obs);
129 xfree((void**)&vy_obs);
130 xfree((void**)&control_parameter);
131
132
133 /*Assign output pointer: */
134 *pcount=count;
135}
136
137
Note: See TracBrowser for help on using the repository browser.