Ignore:
Timestamp:
04/08/14 20:39:20 (11 years ago)
Author:
Mathieu Morlighem
Message:

CHG: following Chris' suggestions, we adopt the following convention
md.mesh.dimension -> DomainDimensionEnum
md.mesh.domaintype ->DomainTypeEnum
md.mesh.elementtype->MeshElementtypeEnum

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/analyses/L2ProjectionEPLAnalysis.cpp

    r17674 r17686  
    66
    77/*Model processing*/
    8 int  L2ProjectionEPLAnalysis::DofsPerNode(int** doflist,int meshxdim,int approximation){/*{{{*/
     8int  L2ProjectionEPLAnalysis::DofsPerNode(int** doflist,int domaintype,int approximation){/*{{{*/
    99        return 1;
    1010}/*}}}*/
    … …  
    3636        iomodel->FetchDataToInput(elements,EplHeadEnum);
    3737        iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum);
    38         if(iomodel->meshxdim==Mesh3DEnum){
     38        if(iomodel->domaintype==Mesh3DEnum){
    3939                iomodel->FetchDataToInput(elements,MeshElementonbaseEnum);
    4040        }
    … …  
    5151        if(!isefficientlayer) return;
    5252
    53         if(iomodel->meshxdim==Mesh3DEnum){
     53        if(iomodel->domaintype==Mesh3DEnum){
    5454                iomodel->FetchData(1,MeshVertexonbaseEnum);
    5555        }
    56         else if(iomodel->meshxdim==Mesh2DverticalEnum){
     56        else if(iomodel->domaintype==Mesh2DverticalEnum){
    5757                iomodel->FetchData(1,MeshVertexonbaseEnum);
    5858        }
    … …  
    8383
    8484        /*Intermediaries*/
    85         int      meshxdim;
     85        int      domaintype;
    8686        bool     active_element;
    8787        Element* basalelement;
    8888
    8989        /*Get basal element*/
    90         element->FindParam(&meshxdim,MeshXDimEnum);
    91         switch(meshxdim){
     90        element->FindParam(&domaintype,DomainTypeEnum);
     91        switch(domaintype){
    9292                case Mesh2DhorizontalEnum:
    9393                        basalelement = element;
    … …  
    101101                        basalelement = element->SpawnBasalElement();
    102102                        break;
    103                 default: _error_("mesh "<<EnumToStringx(meshxdim)<<" not supported yet");
     103                default: _error_("mesh "<<EnumToStringx(domaintype)<<" not supported yet");
    104104        }
    105105
    … …  
    109109        /* Check that all nodes are active, else return empty matrix */
    110110        if(!active_element){
    111                 if(meshxdim!=Mesh2DhorizontalEnum){
     111                if(domaintype!=Mesh2DhorizontalEnum){
    112112                        basalelement->DeleteMaterials();
    113113                        delete basalelement;
    … …  
    149149        xDelete<IssmDouble>(basis);
    150150        delete gauss;
    151         if(meshxdim!=Mesh2DhorizontalEnum){basalelement->DeleteMaterials(); delete basalelement;};
     151        if(domaintype!=Mesh2DhorizontalEnum){basalelement->DeleteMaterials(); delete basalelement;};
    152152        return Ke;
    153153}/*}}}*/
    … …  
    155155
    156156        /*Intermediaries*/
    157         int      meshxdim;
     157        int      domaintype;
    158158        bool     active_element;
    159159        Element* basalelement;
    160160
    161161        /*Get basal element*/
    162         element->FindParam(&meshxdim,MeshXDimEnum);
    163         switch(meshxdim){
     162        element->FindParam(&domaintype,DomainTypeEnum);
     163        switch(domaintype){
    164164                case Mesh2DhorizontalEnum:
    165165                        basalelement = element;
    … …  
    169169                        basalelement = element->SpawnBasalElement();
    170170                        break;
    171                 default: _error_("mesh "<<EnumToStringx(meshxdim)<<" not supported yet");
     171                default: _error_("mesh "<<EnumToStringx(domaintype)<<" not supported yet");
    172172        }
    173173
    … …  
    177177        /*Check that all nodes are active, else return empty matrix*/
    178178        if(!active_element) {
    179                 if(meshxdim!=Mesh2DhorizontalEnum){
     179                if(domaintype!=Mesh2DhorizontalEnum){
    180180                        basalelement->DeleteMaterials();
    181181                        delete basalelement;
    … …  
    222222        xDelete<IssmDouble>(basis);
    223223        delete gauss;
    224         if(meshxdim!=Mesh2DhorizontalEnum){basalelement->DeleteMaterials(); delete basalelement;};
     224        if(domaintype!=Mesh2DhorizontalEnum){basalelement->DeleteMaterials(); delete basalelement;};
    225225        return pe;
    226226}/*}}}*/
    … …  
    229229}/*}}}*/
    230230void L2ProjectionEPLAnalysis::InputUpdateFromSolution(IssmDouble* solution,Element* element){/*{{{*/
    231         int inputenum,meshxdim;
     231        int inputenum,domaintype;
    232232
    233233        element->FindParam(&inputenum,InputToL2ProjectEnum);
    234         element->FindParam(&meshxdim,MeshXDimEnum);
    235         switch(meshxdim){
     234        element->FindParam(&domaintype,DomainTypeEnum);
     235        switch(domaintype){
    236236                case Mesh2DhorizontalEnum:
    237237                        element->InputUpdateFromSolutionOneDof(solution,inputenum);
    … …  
    243243                        element->InputUpdateFromSolutionOneDofCollapsed(solution,inputenum);
    244244                        break;
    245                 default: _error_("mesh "<<EnumToStringx(meshxdim)<<" not supported yet");
     245                default: _error_("mesh "<<EnumToStringx(domaintype)<<" not supported yet");
    246246        }
    247247}/*}}}*/
Note: See TracChangeset for help on using the changeset viewer.