Changeset 4548


Ignore:
Timestamp:
07/13/10 09:21:07 (15 years ago)
Author:
Eric.Larour
Message:

New processing of units in funcions such as MinVel, MaxVx, etc ...
issm.cpp: added extern declarations for my_rank and num_procs, extremely important!

Location:
issm/trunk
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/issm.cpp

    r4491 r4548  
    66
    77int main(int argc,char* *argv){
    8        
     8
     9        /*MPI: */
     10        extern int my_rank;
     11        extern int num_procs;
     12
    913        /*I/O: */
    1014        FILE     *fid              = NULL;

  • issm/trunk/src/c/modules/Qmux/SpawnCoreParallel.cpp

    r4454 r4548  
    6565        /*Run the core solution sequence: */
    6666        if(verbose)_printf_("%s%s%s\n","Starting ",EnumAsString(solution_type)," core:");
    67         SolutionConfiguration(NULL,NULL,&solutioncore,solution_type); solutioncore(femmodel);
     67        SolutionConfiguration(NULL,NULL,&solutioncore,solution_type);
     68        solutioncore(femmodel);
    6869
    6970        /*compute responses: */

  • issm/trunk/src/c/objects/Elements/Beam.cpp

    r4546 r4548  
    481481        inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numgrids,VxEnum);
    482482
     483        /*process units if requested: */
     484        if(process_units)NodalValuesUnitConversion(&vx_values[0],numgrids,VxEnum,this->parameters);
     485
    483486        /*now, compute maximum:*/
    484487        maxabsvx=fabs(vx_values[0]);
    … …  
    507510        inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numgrids,VyEnum);
    508511
     512        /*process units if requested: */
     513        if(process_units)NodalValuesUnitConversion(&vy_values[0],numgrids,VyEnum,this->parameters);
     514
    509515        /*now, compute maximum:*/
    510516        maxabsvy=fabs(vy_values[0]);
    … …  
    533539        inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numgrids,VzEnum);
    534540
     541        /*process units if requested: */
     542        if(process_units)NodalValuesUnitConversion(&vz_values[0],numgrids,VzEnum,this->parameters);
     543
    535544        /*now, compute maximum:*/
    536545        maxabsvz=fabs(vz_values[0]);
    … …  
    559568        inputs->GetParameterValues(&vel_values[0],&gaussgrids[0][0],numgrids,VelEnum);
    560569
     570        /*process units if requested: */
     571        if(process_units)NodalValuesUnitConversion(&vel_values[0],numgrids,VelEnum,this->parameters);
     572
    561573        /*now, compute maximum:*/
    562574        maxvel=vel_values[0];
    … …  
    586598        inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numgrids,VxEnum);
    587599
     600        /*process units if requested: */
     601        if(process_units)NodalValuesUnitConversion(&vx_values[0],numgrids,VxEnum,this->parameters);
     602
    588603        /*now, compute maximum:*/
    589604        maxvx=vx_values[0];
    … …  
    613628        inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numgrids,VyEnum);
    614629
     630        /*process units if requested: */
     631        if(process_units)NodalValuesUnitConversion(&vy_values[0],numgrids,VyEnum,this->parameters);
     632
    615633        /*now, compute maximum:*/
    616634        maxvy=vy_values[0];
    … …  
    640658        inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numgrids,VzEnum);
    641659
     660        /*process units if requested: */
     661        if(process_units)NodalValuesUnitConversion(&vz_values[0],numgrids,VzEnum,this->parameters);
     662
    642663        /*now, compute maximum:*/
    643664        maxvz=vz_values[0];
    … …  
    667688        inputs->GetParameterValues(&vel_values[0],&gaussgrids[0][0],numgrids,VelEnum);
    668689
     690        /*process units if requested: */
     691        if(process_units)NodalValuesUnitConversion(&vel_values[0],numgrids,VelEnum,this->parameters);
     692
    669693        /*now, compute minimum:*/
    670694        minvel=vel_values[0];
    … …  
    694718        inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numgrids,VxEnum);
    695719
     720        /*process units if requested: */
     721        if(process_units)NodalValuesUnitConversion(&vx_values[0],numgrids,VxEnum,this->parameters);
     722
    696723        /*now, compute minimum:*/
    697724        minvx=vx_values[0];
    … …  
    721748        inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numgrids,VyEnum);
    722749
     750        /*process units if requested: */
     751        if(process_units)NodalValuesUnitConversion(&vy_values[0],numgrids,VyEnum,this->parameters);
     752
    723753        /*now, compute minimum:*/
    724754        minvy=vy_values[0];
    … …  
    747777        /*retrive velocity values at nodes */
    748778        inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numgrids,VzEnum);
     779
     780        /*process units if requested: */
     781        if(process_units)NodalValuesUnitConversion(&vz_values[0],numgrids,VzEnum,this->parameters);
    749782
    750783        /*now, compute minimum:*/

  • issm/trunk/src/c/objects/Elements/Penta.cpp

    r4546 r4548  
    12401240        inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numgrids,VxEnum);
    12411241
     1242        /*process units if requested: */
     1243        if(process_units)NodalValuesUnitConversion(&vx_values[0],numgrids,VxEnum,this->parameters);
     1244
    12421245        /*now, compute maximum:*/
    12431246        maxabsvx=fabs(vx_values[0]);
    … …  
    12661269        inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numgrids,VyEnum);
    12671270
     1271        /*process units if requested: */
     1272        if(process_units)NodalValuesUnitConversion(&vy_values[0],numgrids,VyEnum,this->parameters);
     1273
    12681274        /*now, compute maximum:*/
    12691275        maxabsvy=fabs(vy_values[0]);
    … …  
    12921298        inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numgrids,VzEnum);
    12931299
     1300        /*process units if requested: */
     1301        if(process_units)NodalValuesUnitConversion(&vz_values[0],numgrids,VzEnum,this->parameters);
     1302
    12941303        /*now, compute maximum:*/
    12951304        maxabsvz=fabs(vz_values[0]);
    … …  
    13181327        inputs->GetParameterValues(&vel_values[0],&gaussgrids[0][0],numgrids,VelEnum);
    13191328
     1329        /*process units if requested: */
     1330        if(process_units)NodalValuesUnitConversion(&vel_values[0],numgrids,VelEnum,this->parameters);
     1331
    13201332        /*now, compute maximum:*/
    13211333        maxvel=vel_values[0];
    … …  
    13451357        inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numgrids,VxEnum);
    13461358
     1359        /*process units if requested: */
     1360        if(process_units)NodalValuesUnitConversion(&vx_values[0],numgrids,VxEnum,this->parameters);
     1361
    13471362        /*now, compute maximum:*/
    13481363        maxvx=vx_values[0];
    … …  
    13721387        inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numgrids,VyEnum);
    13731388
     1389        /*process units if requested: */
     1390        if(process_units)NodalValuesUnitConversion(&vy_values[0],numgrids,VyEnum,this->parameters);
     1391
    13741392        /*now, compute maximum:*/
    13751393        maxvy=vy_values[0];
    … …  
    13991417        inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numgrids,VzEnum);
    14001418
     1419        /*process units if requested: */
     1420        if(process_units)NodalValuesUnitConversion(&vz_values[0],numgrids,VzEnum,this->parameters);
     1421
    14011422        /*now, compute maximum:*/
    14021423        maxvz=vz_values[0];
    … …  
    14261447        inputs->GetParameterValues(&vel_values[0],&gaussgrids[0][0],numgrids,VelEnum);
    14271448
     1449        /*process units if requested: */
     1450        if(process_units)NodalValuesUnitConversion(&vel_values[0],numgrids,VelEnum,this->parameters);
     1451
    14281452        /*now, compute minimum:*/
    14291453        minvel=vel_values[0];
    … …  
    14531477        inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numgrids,VxEnum);
    14541478
     1479        /*process units if requested: */
     1480        if(process_units)NodalValuesUnitConversion(&vx_values[0],numgrids,VxEnum,this->parameters);
     1481
    14551482        /*now, compute minimum:*/
    14561483        minvx=vx_values[0];
    … …  
    14801507        inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numgrids,VyEnum);
    14811508
     1509        /*process units if requested: */
     1510        if(process_units)NodalValuesUnitConversion(&vy_values[0],numgrids,VyEnum,this->parameters);
     1511
    14821512        /*now, compute minimum:*/
    14831513        minvy=vy_values[0];
    … …  
    15061536        /*retrive velocity values at nodes */
    15071537        inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numgrids,VzEnum);
     1538
     1539        /*process units if requested: */
     1540        if(process_units)NodalValuesUnitConversion(&vz_values[0],numgrids,VzEnum,this->parameters);
    15081541
    15091542        /*now, compute minimum:*/

  • issm/trunk/src/c/objects/Elements/Sing.cpp

    r4546 r4548  
    445445        /*retrive velocity values at nodes */
    446446        inputs->GetParameterValue(&maxabsvx,VxEnum);
     447
     448        /*process units if requested: */
     449        if(process_units)NodalValuesUnitConversion(&maxabsvx,1,VxEnum,this->parameters);
     450
    447451        maxabsvx=fabs(maxabsvx);
    448452
    … …  
    462466        /*retrive velocity values at nodes */
    463467        inputs->GetParameterValue(&maxabsvy,VyEnum);
     468
     469        /*process units if requested: */
     470        if(process_units)NodalValuesUnitConversion(&maxabsvy,1,VyEnum,this->parameters);
     471       
    464472        maxabsvy=fabs(maxabsvy);
    465473
    … …  
    479487        /*retrive velocity values at nodes */
    480488        inputs->GetParameterValue(&maxabsvz,VzEnum);
     489
     490        /*process units if requested: */
     491        if(process_units)NodalValuesUnitConversion(&maxabsvz,1,VzEnum,this->parameters);
     492
    481493        maxabsvz=fabs(maxabsvz);
    482494
    … …  
    497509        inputs->GetParameterValue(&maxvel,VelEnum);
    498510
     511        /*process units if requested: */
     512        if(process_units)NodalValuesUnitConversion(&maxvel,1,VelEnum,this->parameters);
     513
    499514        /*Assign output pointers:*/
    500515        *pmaxvel=maxvel;
    … …  
    514529        inputs->GetParameterValue(&maxvx,VxEnum);
    515530
     531        /*process units if requested: */
     532        if(process_units)NodalValuesUnitConversion(&maxvx,1,VxEnum,this->parameters);
     533
    516534        /*Assign output pointers:*/
    517535        *pmaxvx=maxvx;
    … …  
    531549        inputs->GetParameterValue(&maxvy,VyEnum);
    532550
     551        /*process units if requested: */
     552        if(process_units)NodalValuesUnitConversion(&maxvy,1,VyEnum,this->parameters);
     553
    533554        /*Assign output pointers:*/
    534555        *pmaxvy=maxvy;
    … …  
    548569        inputs->GetParameterValue(&maxvz,VzEnum);
    549570
     571        /*process units if requested: */
     572        if(process_units)NodalValuesUnitConversion(&maxvz,1,VzEnum,this->parameters);
     573
    550574        /*Assign output pointers:*/
    551575        *pmaxvz=maxvz;
    … …  
    565589        inputs->GetParameterValue(&minvel,VelEnum);
    566590
     591        /*process units if requested: */
     592        if(process_units)NodalValuesUnitConversion(&minvel,1,VelEnum,this->parameters);
     593
    567594        /*Assign output pointers:*/
    568595        *pminvel=minvel;
    … …  
    582609        inputs->GetParameterValue(&minvx,VxEnum);
    583610
     611        /*process units if requested: */
     612        if(process_units)NodalValuesUnitConversion(&minvx,1,VxEnum,this->parameters);
     613
    584614        /*Assign output pointers:*/
    585615        *pminvx=minvx;
    … …  
    599629        inputs->GetParameterValue(&minvy,VyEnum);
    600630
     631        /*process units if requested: */
     632        if(process_units)NodalValuesUnitConversion(&minvy,1,VyEnum,this->parameters);
     633
    601634        /*Assign output pointers:*/
    602635        *pminvy=minvy;
    … …  
    615648        /*retrive velocity values at nodes */
    616649        inputs->GetParameterValue(&minvz,VzEnum);
     650
     651        /*process units if requested: */
     652        if(process_units)NodalValuesUnitConversion(&minvz,1,VzEnum,this->parameters);
    617653
    618654        /*Assign output pointers:*/

  • issm/trunk/src/c/objects/Elements/Tria.cpp

    r4546 r4548  
    15901590        inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numgrids,VxEnum);
    15911591
     1592        /*process units if requested: */
     1593        if(process_units)NodalValuesUnitConversion(&vx_values[0],numgrids,VxEnum,this->parameters);
     1594
    15921595        /*now, compute maximum:*/
    15931596        maxabsvx=fabs(vx_values[0]);
    … …  
    16161619        inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numgrids,VyEnum);
    16171620
     1621        /*process units if requested: */
     1622        if(process_units)NodalValuesUnitConversion(&vy_values[0],numgrids,VyEnum,this->parameters);
     1623
    16181624        /*now, compute maximum:*/
    16191625        maxabsvy=fabs(vy_values[0]);
    … …  
    16421648        inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numgrids,VzEnum);
    16431649
     1650        /*process units if requested: */
     1651        if(process_units)NodalValuesUnitConversion(&vz_values[0],numgrids,VzEnum,this->parameters);
     1652
    16441653        /*now, compute maximum:*/
    16451654        maxabsvz=fabs(vz_values[0]);
    … …  
    16681677        inputs->GetParameterValues(&vel_values[0],&gaussgrids[0][0],numgrids,VelEnum);
    16691678
     1679        /*process units if requested: */
     1680        if(process_units)NodalValuesUnitConversion(&vel_values[0],numgrids,VelEnum,this->parameters);
     1681
    16701682        /*now, compute maximum:*/
    16711683        maxvel=vel_values[0];
    … …  
    16951707        inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numgrids,VxEnum);
    16961708
     1709        /*process units if requested: */
     1710        if(process_units)NodalValuesUnitConversion(&vx_values[0],numgrids,VxEnum,this->parameters);
     1711
    16971712        /*now, compute maximum:*/
    16981713        maxvx=vx_values[0];
    … …  
    17221737        inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numgrids,VyEnum);
    17231738
     1739        /*process units if requested: */
     1740        if(process_units)NodalValuesUnitConversion(&vy_values[0],numgrids,VyEnum,this->parameters);
     1741
    17241742        /*now, compute maximum:*/
    17251743        maxvy=vy_values[0];
    … …  
    17491767        inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numgrids,VzEnum);
    17501768
     1769        /*process units if requested: */
     1770        if(process_units)NodalValuesUnitConversion(&vz_values[0],numgrids,VzEnum,this->parameters);
     1771
    17511772        /*now, compute maximum:*/
    17521773        maxvz=vz_values[0];
    … …  
    17761797        inputs->GetParameterValues(&vel_values[0],&gaussgrids[0][0],numgrids,VelEnum);
    17771798
     1799        /*process units if requested: */
     1800        if(process_units)NodalValuesUnitConversion(&vel_values[0],numgrids,VelEnum,this->parameters);
     1801
    17781802        /*now, compute minimum:*/
    17791803        minvel=vel_values[0];
    … …  
    18031827        inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numgrids,VxEnum);
    18041828
     1829        /*process units if requested: */
     1830        if(process_units)NodalValuesUnitConversion(&vx_values[0],numgrids,VxEnum,this->parameters);
     1831
    18051832        /*now, compute minimum:*/
    18061833        minvx=vx_values[0];
    … …  
    18301857        inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numgrids,VyEnum);
    18311858
     1859        /*process units if requested: */
     1860        if(process_units)NodalValuesUnitConversion(&vy_values[0],numgrids,VyEnum,this->parameters);
     1861
    18321862        /*now, compute minimum:*/
    18331863        minvy=vy_values[0];
    … …  
    18561886        /*retrive velocity values at nodes */
    18571887        inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numgrids,VzEnum);
     1888
     1889        /*process units if requested: */
     1890        if(process_units)NodalValuesUnitConversion(&vz_values[0],numgrids,VzEnum,this->parameters);
    18581891
    18591892        /*now, compute minimum:*/

  • issm/trunk/src/m/solutions/SpawnCore.m

    r4532 r4548  
    1616%first update the inputs to the femmodel using the variables provided to us by dakota.
    1717[femmodel.elements femmodel.loads femmodel.materials]=InputUpdateFromDakota(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,femmodel.part,variables,variabledescriptors);
    18 Echo(femmodel.materials)
    1918
    2019%now run the core solution
    … …  
    2625
    2726format long
    28 variables
    2927responses

  • issm/trunk/test/Verification/test25_IceSheetIceFrontM2dDakota/Square.par

    r4532 r4548  
    3939
    4040
    41 
    4241%Dakota options
    4342md.variables.nuv=normal_uncertain.empty();
Note: See TracChangeset for help on using the changeset viewer.