Ignore:
Timestamp:
08/25/09 13:04:05 (16 years ago)
Author:
Eric.Larour
Message:

Model and FemModel are now classes in their own right.
This changes the cores quite a bit.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/parallel/ControlInitialization.cpp

    r1854 r1881  
    5252        int dof3[1]={3};
    5353
     54        /*first recover parameters common to all solutions:*/
     55        model->FindParam(&debug,"debug");
     56        model->FindParam(&dim,"dim");
     57        model->FindParam(&ishutter,"ishutter");
     58        model->FindParam(&ismacayealpattyn,"ismacayealpattyn");
     59        model->FindParam(&numberofnodes,"numberofnodes");
     60        model->FindParam(&isstokes,"isstokes");
     61        model->FindParam(&stokesreconditioning,"stokesreconditioning");
     62
    5463        /*recover fem models: */
    55         fem_dh=model->DiagnosticHorizontal();
    56         fem_dv=model->DiagnosticVertical();
    57         fem_ds=model->DiagnosticStokes();
    58         fem_dhu=model->DiagnosticHutter();
    59         fem_sl=model->Slope();
    60 
    61         //first recover parameters common to all solutions
    62         fem_dh->parameters->FindParam((void*)&debug,"debug");
    63         fem_dh->parameters->FindParam((void*)&dim,"dim");
    64         fem_dhu->parameters->FindParam((void*)&ishutter,"ishutter");
    65         fem_dh->parameters->FindParam((void*)&ismacayealpattyn,"ismacayealpattyn");
    66         fem_dh->parameters->FindParam((void*)&numberofnodes,"numberofnodes");
    67         fem_ds->parameters->FindParam((void*)&isstokes,"isstokes");
    68         fem_ds->parameters->FindParam((void*)&stokesreconditioning,"stokesreconditioning");
     64        fem_dh=model->GetFormulation(DiagnosticAnalysisEnum(),HorizAnalysisEnum());
     65        fem_dv=model->GetFormulation(DiagnosticAnalysisEnum(),VertAnalysisEnum());
     66        fem_ds=model->GetFormulation(DiagnosticAnalysisEnum(),StokesAnalysisEnum());
     67        fem_dhu=model->GetFormulation(DiagnosticAnalysisEnum(),HutterAnalysisEnum());
     68        fem_sl=model->GetFormulation(SlopeComputeAnalysisEnum());
    6969
    7070        //specific parameters for specific models
    71         fem_dh->parameters->FindParam((void*)&numberofdofspernode_dh,"numberofdofspernode");
    72         fem_sl->parameters->FindParam((void*)&numberofdofspernode_sl,"numberofdofspernode");
    73         fem_ds->parameters->FindParam((void*)&numberofdofspernode_ds,"numberofdofspernode");
     71        fem_dh->FindParam((void*)&numberofdofspernode_dh,"numberofdofspernode");
     72        fem_sl->FindParam((void*)&numberofdofspernode_sl,"numberofdofspernode");
     73        fem_ds->FindParam((void*)&numberofdofspernode_ds,"numberofdofspernode");
    7474
    7575        /*if no Stokes, assign output and return*/
    7676        if (!isstokes){
    77                 model->SetActive(fem_dh);
     77                model->SetActiveFormulation(fem_dh);
    7878                return;
    7979        }
     
    146146
    147147        /*Assign output*/
    148         model->SetActive(fem_ds);
     148        model->SetActiveFormulation(fem_ds);
    149149}
Note: See TracChangeset for help on using the changeset viewer.