Changeset 3703 for issm/trunk/src/c/parallel/ProcessResults.cpp
- Timestamp:
- 05/11/10 08:28:39 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/parallel/ProcessResults.cpp
r3588 r3703 112 112 113 113 /*some flags needed: */ 114 model->FindParam(&dim, "dim");115 model->FindParam(&ishutter, "ishutter");116 model->FindParam(&isstokes, "isstokes");117 model->FindParam(&ismacayealpattyn, "ismacayealpattyn");114 model->FindParam(&dim,DimEnum); 115 model->FindParam(&ishutter,IsHutterEnum); 116 model->FindParam(&isstokes,IsStokesEnum); 117 model->FindParam(&ismacayealpattyn,IsMacayealPattynEnum); 118 118 119 119 /*Recover femmodels first: */ … … 158 158 /*ok, 2 dofs, on number of nodes: */ 159 159 if(ismacayealpattyn){ 160 fem_dh->parameters->FindParam(&numberofnodes, "numberofnodes");160 fem_dh->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 161 161 VecToMPISerial(&partition,fem_dh->partition->vector); 162 fem_dh->parameters->FindParam(&yts, "yts");162 fem_dh->parameters->FindParam(&yts,YtsEnum); 163 163 } 164 164 else{ 165 fem_dhu->parameters->FindParam(&numberofnodes, "numberofnodes");165 fem_dhu->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 166 166 VecToMPISerial(&partition,fem_dhu->partition->vector); 167 fem_dhu->parameters->FindParam(&yts, "yts");167 fem_dhu->parameters->FindParam(&yts,YtsEnum); 168 168 } 169 169 vx=(double*)xmalloc(numberofnodes*sizeof(double)); … … 184 184 /*ok, 3 dofs, on number of nodes: */ 185 185 if(ismacayealpattyn){ 186 fem_dh->parameters->FindParam(&numberofnodes, "numberofnodes");186 fem_dh->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 187 187 VecToMPISerial(&partition,fem_dh->partition->vector); 188 fem_dh->parameters->FindParam(&yts, "yts");188 fem_dh->parameters->FindParam(&yts,YtsEnum); 189 189 } 190 190 else{ 191 fem_dhu->parameters->FindParam(&numberofnodes, "numberofnodes");191 fem_dhu->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 192 192 VecToMPISerial(&partition,fem_dhu->partition->vector); 193 fem_dhu->parameters->FindParam(&yts, "yts");193 fem_dhu->parameters->FindParam(&yts,YtsEnum); 194 194 } 195 195 vx=(double*)xmalloc(numberofnodes*sizeof(double)); … … 207 207 else{ 208 208 /* 4 dofs on number of nodes. discard pressure: */ 209 fem_ds->parameters->FindParam(&numberofnodes, "numberofnodes");209 fem_ds->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 210 210 VecToMPISerial(&partition,fem_ds->partition->vector); 211 fem_ds->parameters->FindParam(&yts, "yts");211 fem_ds->parameters->FindParam(&yts,YtsEnum); 212 212 vx=(double*)xmalloc(numberofnodes*sizeof(double)); 213 213 vy=(double*)xmalloc(numberofnodes*sizeof(double)); … … 252 252 if(!isstokes){ 253 253 if(ismacayealpattyn){ 254 fem_dh->parameters->FindParam(&numberofnodes, "numberofnodes");254 fem_dh->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 255 255 VecToMPISerial(&partition,fem_dh->partition->vector); 256 256 } 257 257 else{ 258 fem_dhu->parameters->FindParam(&numberofnodes, "numberofnodes");258 fem_dhu->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 259 259 VecToMPISerial(&partition,fem_dhu->partition->vector); 260 260 } 261 261 } 262 262 else{ 263 fem_ds->parameters->FindParam(&numberofnodes, "numberofnodes");263 fem_ds->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 264 264 VecToMPISerial(&partition,fem_ds->partition->vector); 265 265 } … … 285 285 result->GetField(&t_g); 286 286 VecToMPISerial(&t_g_serial,t_g); 287 fem_t->parameters->FindParam(&numberofnodes, "numberofnodes");287 fem_t->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 288 288 VecToMPISerial(&partition,fem_t->partition->vector); 289 289 … … 309 309 result->GetField(&grad_g); 310 310 VecToMPISerial(&grad_g_serial,grad_g); 311 fem_c->parameters->FindParam(&numberofnodes, "numberofnodes");311 fem_c->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 312 312 VecToMPISerial(&partition,fem_c->partition->vector); 313 313 … … 332 332 result->GetField(&m_g); 333 333 VecToMPISerial(&m_g_serial,m_g); 334 fem_t->parameters->FindParam(&numberofnodes, "numberofnodes");335 fem_t->parameters->FindParam(&yts, "yts");334 fem_t->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 335 fem_t->parameters->FindParam(&yts,YtsEnum); 336 336 VecToMPISerial(&partition,fem_t->partition->vector); 337 337 … … 356 356 result->GetField(&h_g); 357 357 VecToMPISerial(&h_g_serial,h_g); 358 fem_p->parameters->FindParam(&numberofvertices, "numberofvertices");358 fem_p->parameters->FindParam(&numberofvertices,NumberOfVerticesEnum); 359 359 VecToMPISerial(&partition,fem_p->partition->vector); 360 360 … … 379 379 result->GetField(&v_g); 380 380 VecToMPISerial(&v_g_serial,v_g); 381 fem_p->parameters->FindParam(&numberofnodes, "numberofnodes");381 fem_p->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 382 382 VecToMPISerial(&partition,fem_p->partition->vector); 383 383 … … 402 402 result->GetField(&s_g); 403 403 VecToMPISerial(&s_g_serial,s_g); 404 fem_p->parameters->FindParam(&numberofnodes, "numberofnodes");404 fem_p->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 405 405 VecToMPISerial(&partition,fem_p->partition->vector); 406 406 … … 425 425 result->GetField(&sx_g); 426 426 VecToMPISerial(&sx_g_serial,sx_g); 427 fem_sl->parameters->FindParam(&numberofnodes, "numberofnodes");427 fem_sl->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 428 428 VecToMPISerial(&partition,fem_sl->partition->vector); 429 429 … … 448 448 result->GetField(&sy_g); 449 449 VecToMPISerial(&sy_g_serial,sy_g); 450 fem_sl->parameters->FindParam(&numberofnodes, "numberofnodes");450 fem_sl->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 451 451 VecToMPISerial(&partition,fem_sl->partition->vector); 452 452 … … 471 471 result->GetField(&b_g); 472 472 VecToMPISerial(&b_g_serial,b_g); 473 fem_p->parameters->FindParam(&numberofnodes, "numberofnodes");473 fem_p->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 474 474 VecToMPISerial(&partition,fem_p->partition->vector); 475 475 … … 493 493 /*easy, param_g is of size numberofnodes, on 1 dof, just repartition: */ 494 494 result->GetField(¶m_g); 495 fem_dh->parameters->FindParam(&numberofnodes, "numberofnodes");495 fem_dh->parameters->FindParam(&numberofnodes,NumberOfNodesEnum); 496 496 VecToMPISerial(&partition,fem_dh->partition->vector); 497 497 … … 513 513 else if(strcmp(result->GetFieldName(),"riftproperties")==0){ 514 514 result->GetField(&riftproperties); 515 fem_dh->parameters->FindParam(&numrifts, "numrifts");515 fem_dh->parameters->FindParam(&numrifts,NumRiftsEnum); 516 516 VecToMPISerial(&riftproperties_serial,riftproperties); 517 517 … … 524 524 else if(strcmp(result->GetFieldName(),"sigma_zz")==0){ 525 525 /*easy, param_g is of size numberofelements, on 1 dof, just repartition: */ 526 fem_ds->parameters->FindParam(&numberofelements, "numberofelements");526 fem_ds->parameters->FindParam(&numberofelements,NumberOfElementsEnum); 527 527 result->GetField(&sigma_zz); 528 528 VecToMPISerial(&sigma_zz_serial,sigma_zz);
Note:
See TracChangeset
for help on using the changeset viewer.