source: issm/trunk-jpl/src/c/modules/ModelProcessorx/Dakota/CreateParametersDakota.cpp@ 14996

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

CHG: integrated Container/ directory into src/c/classes/objects directory. No reason to have the containers
and the objects that they contain defined in different places.
File size: 4.8 KB
Line 
1/*!\file: CreateParametersDakota.cpp
2 * \brief general driver for creating parameters dataset
3 */
4
5#include "../../../toolkits/toolkits.h"
6#include "../../../shared/Enum/Enum.h"
7#include "../../../classes/objects/objects.h"
8#include "../../../shared/shared.h"
9#include "../../../shared/io/io.h"
10#include "../../MeshPartitionx/MeshPartitionx.h"
11#include "../ModelProcessorx.h"
12
13void CreateParametersDakota(Parameters** pparameters,IoModel* iomodel,char* rootpath,int solution_type,int analysis_type){
14
15 /*variable declarations: {{{*/
16 int i;
17 Parameters* parameters = NULL;
18
19 int* part=NULL;
20 double* dpart=NULL;
21
22 char** responsedescriptors=NULL;
23 int numresponsedescriptors;
24 char** variabledescriptors=NULL;
25 int numvariabledescriptors;
26 char* descriptor=NULL;
27 double* dakota_parameter=NULL;
28
29 //qmu files
30 char* qmuinname=NULL;
31 char* qmuerrname=NULL;
32 char* qmuoutname=NULL;
33
34 //descriptors:
35 char tag[50];
36
37 bool dakota_analysis=false;
38 char* name=NULL;
39 int numberofresponses;
40 int numberofvertices;
41 int nrows;
42 int ncols;
43
44 /*}}}*/
45
46 /*recover parameters : */
47 parameters=*pparameters;
48
49 /*recover parameters: */
50 iomodel->Constant(&dakota_analysis,QmuIsdakotaEnum);
51
52 if(dakota_analysis){
53
54 iomodel->Constant(&name,MiscellaneousNameEnum);
55 iomodel->Constant(&numberofresponses,QmuNumberofresponsesEnum);
56 iomodel->Constant(&numberofvertices,MeshNumberofverticesEnum);
57
58 /*name of qmu input, error and output files:{{{*/
59 qmuinname=xNew<char>((strlen(rootpath)+strlen(name)+strlen(".qmu.in")+1));
60 sprintf(qmuinname,"%s%s%s",rootpath,name,".qmu.in");
61 parameters->AddObject(new StringParam(QmuInNameEnum,qmuinname));
62
63 qmuoutname=xNew<char>((strlen(rootpath)+strlen(name)+strlen(".qmu.out")+1));
64 sprintf(qmuoutname,"%s%s%s",rootpath,name,".qmu.out");
65 parameters->AddObject(new StringParam(QmuOutNameEnum,qmuoutname));
66
67 qmuerrname=xNew<char>((strlen(rootpath)+strlen(name)+strlen(".qmu.err")+1));
68 sprintf(qmuerrname,"%s%s%s",rootpath,name,".qmu.err");
69 parameters->AddObject(new StringParam(QmuErrNameEnum,qmuerrname));
70 /*}}}*/
71 /*Fetch variable descriptors: {{{*/
72 iomodel->FetchData(&variabledescriptors,&numvariabledescriptors,QmuVariabledescriptorsEnum);
73
74 /*Ok, we have all the variable descriptors. Build a parameter with it: */
75 parameters->AddObject(new StringArrayParam(QmuVariabledescriptorsEnum,variabledescriptors,numvariabledescriptors));
76
77 /*}}}*/
78 /*Fetch response descriptors: {{{*/
79 iomodel->FetchData(&responsedescriptors,&numresponsedescriptors,QmuResponsedescriptorsEnum);
80
81 /*Ok, we have all the response descriptors. Build a parameter with it: */
82 parameters->AddObject(new StringArrayParam(QmuResponsedescriptorsEnum,responsedescriptors,numresponsedescriptors));
83 parameters->AddObject(new IntParam(QmuNumberofresponsesEnum,numberofresponses));
84 /*}}}*/
85 /*Deal with partitioning: {{{*/
86 /*partition vertices in iomodel->qmu_npart parts, unless a partition is already present: */
87
88 parameters->AddObject(iomodel->CopyConstantObject(QmuNumberofpartitionsEnum));
89 iomodel->FetchData(&dpart,NULL,NULL,QmuPartitionEnum);
90
91 if(!dpart){
92
93 /*Partition elements and vertices and nodes: */
94 ElementsAndVerticesPartitioning(&iomodel->my_elements,&iomodel->my_vertices,iomodel);
95
96 dpart=xNew<double>(numberofvertices);
97 for(i=0;i<numberofvertices;i++)dpart[i]=iomodel->my_vertices[i];
98 }
99 parameters->AddObject(new DoubleVecParam(QmuPartitionEnum,dpart,numberofvertices));
100 /*}}}*/
101 /*Deal with data needed because of qmu variables: {{{*/
102
103 for(i=0;i<numvariabledescriptors;i++){
104
105 if (strncmp(variabledescriptors[i],"scaled_",7)==0){
106 /*Ok, we are dealing with a variable that is distributed over nodes. Recover the name of the variable (ex: scaled_Thickness): */
107 sscanf(variabledescriptors[i],"scaled_%s",tag);
108
109 /*Recover data: */
110 iomodel->FetchData(&dakota_parameter,&nrows,&ncols,StringToEnumx(tag));
111
112 /*Add to parameters: */
113 if(nrows==numberofvertices){
114 parameters->AddObject(new DoubleMatParam(StringToEnumx(tag),dakota_parameter,nrows,ncols));
115 }
116 else{
117 parameters->AddObject(new DoubleTransientMatParam(StringToEnumx(tag),dakota_parameter,nrows,ncols));
118 }
119
120 /*Free ressources:*/
121 xDelete<double>(dakota_parameter);
122 }
123 }
124 /*}}}*/
125 /*Free data: {{{*/
126 for(i=0;i<numresponsedescriptors;i++){
127 descriptor=responsedescriptors[i];
128 xDelete<char>(descriptor);
129 }
130 xDelete<char*>(responsedescriptors);
131
132 for(i=0;i<numvariabledescriptors;i++){
133 descriptor=variabledescriptors[i];
134 xDelete<char>(descriptor);
135 }
136 xDelete<char*>(variabledescriptors);
137 xDelete<int>(part);
138 xDelete<double>(dpart);
139 xDelete<char>(qmuinname);
140 xDelete<char>(qmuerrname);
141 xDelete<char>(qmuoutname);
142 /*}}}*/
143 } //if(dakota_analysis)
144
145 /*Free data*/
146 xDelete<char>(name);
147
148 /*Assign output pointer: */
149 *pparameters=parameters;
150}
Note: See TracBrowser for help on using the repository browser.