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

Last change on this file since 20696 was 20696, checked in by Mathieu Morlighem, 9 years ago

BUG: some errors for AD and dakota using the new marshall system
File size: 5.7 KB
Line 
1/*!\file: CreateParametersDakota.cpp
2 * \brief general driver for creating parameters dataset
3 */
4
5#include "../../../toolkits/toolkits.h"
6#include "../../../classes/classes.h"
7#include "../../../shared/shared.h"
8#include "../../MeshPartitionx/MeshPartitionx.h"
9#include "../ModelProcessorx.h"
10
11void CreateParametersDakota(Parameters* parameters,IoModel* iomodel,char* rootpath){
12
13 /*variable declarations*/
14 int i;
15 int *part = NULL;
16 double *dpart = NULL;
17 char **responsedescriptors = NULL;
18 int numresponsedescriptors;
19 char **variabledescriptors = NULL;
20 int numvariabledescriptors;
21 char *descriptor = NULL;
22 double *dakota_parameter = NULL;
23
24 //qmu files
25 char *qmuinname = NULL;
26 char *qmuerrname = NULL;
27 char *qmuoutname = NULL;
28
29 //descriptors:
30 char tag[50];
31
32 bool dakota_analysis = false;
33 char *name = NULL;
34 int numberofresponses;
35 int nrows,ncols;
36
37 /*recover parameters: */
38 iomodel->FindConstant(&dakota_analysis,"md.qmu.isdakota");
39
40 if(dakota_analysis){
41
42 iomodel->FindConstant(&name,"md.miscellaneous.name");
43 iomodel->FindConstant(&numberofresponses,"md.qmu.numberofresponses");
44
45 /*name of qmu input, error and output files*/
46 qmuinname=xNew<char>((strlen(rootpath)+strlen(name)+strlen(".qmu.in")+1));
47 sprintf(qmuinname,"%s%s%s",rootpath,name,".qmu.in");
48 parameters->AddObject(new StringParam(QmuInNameEnum,qmuinname));
49
50 qmuoutname=xNew<char>((strlen(rootpath)+strlen(name)+strlen(".qmu.out")+1));
51 sprintf(qmuoutname,"%s%s%s",rootpath,name,".qmu.out");
52 parameters->AddObject(new StringParam(QmuOutNameEnum,qmuoutname));
53
54 qmuerrname=xNew<char>((strlen(rootpath)+strlen(name)+strlen(".qmu.err")+1));
55 sprintf(qmuerrname,"%s%s%s",rootpath,name,".qmu.err");
56 parameters->AddObject(new StringParam(QmuErrNameEnum,qmuerrname));
57
58 /*Fetch variable descriptors*/
59 iomodel->FindConstant(&variabledescriptors,&numvariabledescriptors,"md.qmu.variabledescriptors");
60
61 /*Ok, we have all the variable descriptors. Build a parameter with it: */
62 parameters->AddObject(new StringArrayParam(QmuVariabledescriptorsEnum,variabledescriptors,numvariabledescriptors));
63
64 /*Fetch response descriptors*/
65 iomodel->FindConstant(&responsedescriptors,&numresponsedescriptors,"md.qmu.responsedescriptors");
66
67 /*Ok, we have all the response descriptors. Build a parameter with it: */
68 parameters->AddObject(new StringArrayParam(QmuResponsedescriptorsEnum,responsedescriptors,numresponsedescriptors));
69 parameters->AddObject(new IntParam(QmuNumberofresponsesEnum,numberofresponses));
70
71 /*Deal with partitioning*/
72 /*partition vertices in iomodel->qmu_npart parts, unless a partition is already present: */
73 parameters->AddObject(iomodel->CopyConstantObject("md.qmu.numberofpartitions",QmuNumberofpartitionsEnum));
74 iomodel->FetchData(&dpart,NULL,NULL,"md.qmu.partition");
75 if(!dpart){
76 /*Partition elements and vertices and nodes: */
77 ElementsAndVerticesPartitioning(&iomodel->my_elements,&iomodel->my_vertices,iomodel);
78
79 dpart=xNew<double>(iomodel->numberofvertices);
80 for(i=0;i<iomodel->numberofvertices;i++)dpart[i]=iomodel->my_vertices[i];
81 }
82 parameters->AddObject(new DoubleVecParam(QmuPartitionEnum,dpart,iomodel->numberofvertices));
83
84 /*Deal with data needed because of qmu variables*/
85 for(i=0;i<numvariabledescriptors;i++){
86 printf("-------------- file: CreateParametersDakota.cpp line: %i\n",__LINE__);
87 if (strncmp(variabledescriptors[i],"scaled_",7)==0){
88 /*Ok, we are dealing with a variable that is distributed over nodes. Recover the name of the variable (ex: scaled_Thickness): */
89 sscanf(variabledescriptors[i],"scaled_%s",tag);
90 printf("STRING: %s\n",tag);
91
92 /*Convert tag to field name*/
93 char* fieldname = NULL;
94 int param_enum = -1;
95 if(strcmp(tag,"Thickness")==0){
96 const char* field = "md.geometry.thickness";
97 param_enum = ThicknessEnum;
98 fieldname=xNew<char>((strlen(field)+1)); xMemCpy<char>(fieldname,field,(strlen(field)+1));
99 }
100 else if(strcmp(tag,"MaterialsRheologyB")==0){
101 const char* field = "md.materials.rheology_B";
102 param_enum = MaterialsRheologyBEnum;
103 fieldname=xNew<char>((strlen(field)+1)); xMemCpy<char>(fieldname,field,(strlen(field)+1));
104 }
105 else if(strcmp(tag,"SmbMassBalance")==0){
106 const char* field = "md.smb.mass_balance";
107 param_enum = SmbMassBalanceEnum;
108 fieldname=xNew<char>((strlen(field)+1)); xMemCpy<char>(fieldname,field,(strlen(field)+1));
109 }
110 else if(strcmp(tag,"FrictionCoefficient")==0){
111 const char* field = "md.friction.coefficient";
112 param_enum = FrictionCoefficientEnum;
113 fieldname=xNew<char>((strlen(field)+1)); xMemCpy<char>(fieldname,field,(strlen(field)+1));
114 }
115 else{
116 _error_("Tag \""<<tag<<"\" not supported yet");
117 }
118
119 /*Recover data: */
120 iomodel->FetchData(&dakota_parameter,&nrows,&ncols,fieldname);
121 if(nrows==iomodel->numberofvertices){
122 parameters->AddObject(new DoubleMatParam(param_enum,dakota_parameter,nrows,ncols));
123 }
124 else{
125 parameters->AddObject(new DoubleTransientMatParam(param_enum,dakota_parameter,nrows,ncols));
126 }
127 xDelete<double>(dakota_parameter);
128 xDelete<char>(fieldname);
129 }
130 }
131
132 /*clean-up*/
133 for(i=0;i<numresponsedescriptors;i++){
134 descriptor=responsedescriptors[i];
135 xDelete<char>(descriptor);
136 }
137 xDelete<char*>(responsedescriptors);
138 for(i=0;i<numvariabledescriptors;i++){
139 descriptor=variabledescriptors[i];
140 xDelete<char>(descriptor);
141 }
142 xDelete<char*>(variabledescriptors);
143 xDelete<int>(part);
144 xDelete<double>(dpart);
145 xDelete<char>(qmuinname);
146 xDelete<char>(qmuerrname);
147 xDelete<char>(qmuoutname);
148 }
149
150 /*Free data*/
151 xDelete<char>(name);
152}
Note: See TracBrowser for help on using the repository browser.