Changeset 4765


Ignore:
Timestamp:
07/22/10 15:09:02 (15 years ago)
Author:
Mathieu Morlighem
Message:

moved Balancedthickness2 to Balancedthickness

Location:
issm/trunk/src/c
Files:
2 deleted
17 edited

Legend:

Unmodified
Added
Removed

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

    r4739 r4765  
    6262                case SurfaceSlopeXAnalysisEnum : return "SurfaceSlopeXAnalysis";
    6363                case SurfaceSlopeYAnalysisEnum : return "SurfaceSlopeYAnalysis";
    64                 case Balancedthickness2SolutionEnum : return "Balancedthickness2Solution";
    65                 case Balancedthickness2AnalysisEnum : return "Balancedthickness2Analysis";
    6664                case BalancedthicknessSolutionEnum : return "BalancedthicknessSolution";
    6765                case BalancedthicknessAnalysisEnum : return "BalancedthicknessAnalysis";
    … …  
    7169                case PrognosticAnalysisEnum : return "PrognosticAnalysis";
    7270                case MeltingAnalysisEnum : return "MeltingAnalysis";
    73                 case Mesh2gridAnalysisEnum : return "Mesh2gridAnalysis";
    7471                case ParametersSolutionEnum : return "ParametersSolution";
    7572                case ParametersAnalysisEnum : return "ParametersAnalysis";

  • issm/trunk/src/c/EnumDefinitions/EnumDefinitions.h

    r4739 r4765  
    6666        SurfaceSlopeYAnalysisEnum,
    6767        //prognostic
    68         Balancedthickness2SolutionEnum,
    69         Balancedthickness2AnalysisEnum,
    7068        BalancedthicknessSolutionEnum,
    7169        BalancedthicknessAnalysisEnum,
    … …  
    7674        //melting
    7775        MeltingAnalysisEnum,
    78         //mesh2grid
    79         Mesh2gridAnalysisEnum,
    8076        //parameters
    8177        ParametersSolutionEnum,

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

    r4739 r4765  
    6060        else if (strcmp(name,"SurfaceSlopeXAnalysis")==0) return SurfaceSlopeXAnalysisEnum;
    6161        else if (strcmp(name,"SurfaceSlopeYAnalysis")==0) return SurfaceSlopeYAnalysisEnum;
    62         else if (strcmp(name,"Balancedthickness2Solution")==0) return Balancedthickness2SolutionEnum;
    63         else if (strcmp(name,"Balancedthickness2Analysis")==0) return Balancedthickness2AnalysisEnum;
    6462        else if (strcmp(name,"BalancedthicknessSolution")==0) return BalancedthicknessSolutionEnum;
    6563        else if (strcmp(name,"BalancedthicknessAnalysis")==0) return BalancedthicknessAnalysisEnum;
    … …  
    6967        else if (strcmp(name,"PrognosticAnalysis")==0) return PrognosticAnalysisEnum;
    7068        else if (strcmp(name,"MeltingAnalysis")==0) return MeltingAnalysisEnum;
    71         else if (strcmp(name,"Mesh2gridAnalysis")==0) return Mesh2gridAnalysisEnum;
    7269        else if (strcmp(name,"ParametersSolution")==0) return ParametersSolutionEnum;
    7370        else if (strcmp(name,"ParametersAnalysis")==0) return ParametersAnalysisEnum;

  • issm/trunk/src/c/Makefile.am

    r4739 r4765  
    394394                                        ./modules/ModelProcessorx/Balancedthickness/CreateConstraintsBalancedthickness.cpp\
    395395                                        ./modules/ModelProcessorx/Balancedthickness/CreateLoadsBalancedthickness.cpp\
    396                                         ./modules/ModelProcessorx/Balancedthickness2/UpdateElementsBalancedthickness2.cpp\
    397                                         ./modules/ModelProcessorx/Balancedthickness2/CreateNodesBalancedthickness2.cpp\
    398                                         ./modules/ModelProcessorx/Balancedthickness2/CreateConstraintsBalancedthickness2.cpp\
    399                                         ./modules/ModelProcessorx/Balancedthickness2/CreateLoadsBalancedthickness2.cpp\
    400396                                        ./modules/ModelProcessorx/Balancedvelocities/UpdateElementsBalancedvelocities.cpp\
    401397                                        ./modules/ModelProcessorx/Balancedvelocities/CreateNodesBalancedvelocities.cpp\
    … …  
    950946                                        ./modules/ModelProcessorx/Balancedthickness/CreateConstraintsBalancedthickness.cpp\
    951947                                        ./modules/ModelProcessorx/Balancedthickness/CreateLoadsBalancedthickness.cpp\
    952                                         ./modules/ModelProcessorx/Balancedthickness2/UpdateElementsBalancedthickness2.cpp\
    953                                         ./modules/ModelProcessorx/Balancedthickness2/CreateNodesBalancedthickness2.cpp\
    954                                         ./modules/ModelProcessorx/Balancedthickness2/CreateConstraintsBalancedthickness2.cpp\
    955                                         ./modules/ModelProcessorx/Balancedthickness2/CreateLoadsBalancedthickness2.cpp\
    956948                                        ./modules/ModelProcessorx/Balancedvelocities/UpdateElementsBalancedvelocities.cpp\
    957949                                        ./modules/ModelProcessorx/Balancedvelocities/CreateNodesBalancedvelocities.cpp\
    … …  
    11291121                                        ./solutions/prognostic_core.cpp\
    11301122                                        ./solutions/balancedthickness_core.cpp\
    1131                                         ./solutions/balancedthickness2_core.cpp\
    11321123                                        ./solutions/balancedvelocities_core.cpp\
    11331124                                        ./solutions/surfaceslope_core.cpp\

  • issm/trunk/src/c/modules/ModelProcessorx/Balancedthickness/CreateConstraintsBalancedthickness.cpp

    r4236 r4765  
    2626        if(!constraints) constraints = new Constraints(ConstraintsEnum);
    2727
    28         /*Fetch data: */
    29         IoModelFetchData(&iomodel->spcthickness,NULL,NULL,iomodel_handle,"spcthickness");
     28        /*Do not add constraints in DG*/
     29        if(!iomodel->prognostic_DG){
    3030
    31         /*Initialize counter*/
    32         count=0;
     31                /*Fetch data: */
     32                IoModelFetchData(&iomodel->spcthickness,NULL,NULL,iomodel_handle,"spcthickness");
    3333
    34         /*Create spcs from x,y,z, as well as the spc values on those spcs: */
    35         for (i=0;i<iomodel->numberofvertices;i++){
    36                 if(iomodel->my_vertices[i]){
    37                        
    38                         if ((int)iomodel->spcthickness[2*i]){
    39                
    40                                 constraints->AddObject(new Spc(iomodel->constraintcounter+count+1,iomodel->nodecounter+i+1,1,
    41                                                                 *(iomodel->spcthickness+2*i+1),BalancedthicknessAnalysisEnum));//we enforce first translation degree of freedom, for temperature
    42                                 count++;
     34                /*Initialize counter*/
     35                count=0;
     36
     37                /*Create spcs from x,y,z, as well as the spc values on those spcs: */
     38                for (i=0;i<iomodel->numberofvertices;i++){
     39                        if(iomodel->my_vertices[i]){
     40
     41                                if ((int)iomodel->spcthickness[2*i]){
     42
     43                                        constraints->AddObject(new Spc(iomodel->constraintcounter+count+1,iomodel->nodecounter+i+1,1,
     44                                                                        *(iomodel->spcthickness+2*i+1),BalancedthicknessAnalysisEnum));//we enforce first translation degree of freedom, for temperature
     45                                        count++;
     46                                }
    4347                        }
    4448                }
     49
     50                /*Free data: */
     51                xfree((void**)&iomodel->spcthickness);
    4552        }
    46 
    47         /*Free data: */
    48         xfree((void**)&iomodel->spcthickness);
    4953
    5054        /*Assign output pointer: */

  • issm/trunk/src/c/modules/ModelProcessorx/Balancedthickness/CreateLoadsBalancedthickness.cpp

    r4236 r4765  
    1313void    CreateLoadsBalancedthickness(Loads** ploads, IoModel* iomodel,ConstDataHandle iomodel_handle){
    1414
    15         /*DataSet*/
     15        /*Intermediary*/
     16        int i;
     17        int element;
     18
     19        /*Output*/
    1620        Loads*    loads    = NULL;
    1721
    … …  
    2226        if(!loads) loads = new Loads(LoadsEnum);
    2327       
     28        /*Loads only in DG*/
     29        if (iomodel->prognostic_DG){
     30
     31                /*Get edges and elements*/
     32                IoModelFetchData(&iomodel->edges,&iomodel->numberofedges,NULL,iomodel_handle,"edges");
     33                IoModelFetchData(&iomodel->elements,NULL,NULL,iomodel_handle,"elements");
     34                IoModelFetchData(&iomodel->thickness,NULL,NULL,iomodel_handle,"thickness");
     35
     36                /*First load data:*/
     37                for (i=0;i<iomodel->numberofedges;i++){
     38
     39                        /*Get left and right elements*/
     40                        element=(int)iomodel->edges[4*i+2]-1; //edges are [node1 node2 elem1 elem2]
     41
     42                        /*Now, if this element is not in the partition, pass: */
     43                        if(!iomodel->my_elements[element]) continue;
     44
     45                        /* Add load */
     46                        loads->AddObject(new Numericalflux(iomodel->loadcounter+i+1,i,iomodel,BalancedthicknessAnalysisEnum));
     47                }
     48
     49                /*Free data: */
     50                xfree((void**)&iomodel->edges);
     51                xfree((void**)&iomodel->elements);
     52                xfree((void**)&iomodel->thickness);
     53        }
     54
    2455        /*Assign output pointer: */
    2556        *ploads=loads;

  • issm/trunk/src/c/modules/ModelProcessorx/Balancedthickness/CreateNodesBalancedthickness.cpp

    r4236 r4765  
    1616
    1717        /*Intermediary*/
    18         int i;
     18        int  i,j;
     19        int  node_id;
     20        int  vertex_id;
     21        int  io_index;
    1922        bool continuous_galerkin=true;
    2023
    … …  
    2932
    3033        /*Continuous Galerkin partition of nodes: */
     34        if(iomodel->prognostic_DG) continuous_galerkin=false;
    3135        NodesPartitioning(&iomodel->my_nodes,iomodel->my_elements, iomodel->my_vertices, iomodel->my_bordervertices, iomodel, iomodel_handle,continuous_galerkin);
     36
     37        /*Check in 3d*/
     38        if(iomodel->prognostic_DG && iomodel->dim==3) ISSMERROR("DG 3d not implemented yet");
    3239
    3340        /*First fetch data: */
    … …  
    3643                IoModelFetchData(&iomodel->uppernodes,NULL,NULL,iomodel_handle,"uppergrids");
    3744        }
     45        IoModelFetchData(&iomodel->elements,NULL,NULL,iomodel_handle,"elements");
    3846        IoModelFetchData(&iomodel->gridonbed,NULL,NULL,iomodel_handle,"gridonbed");
    3947        IoModelFetchData(&iomodel->gridonsurface,NULL,NULL,iomodel_handle,"gridonsurface");
    … …  
    4149        IoModelFetchData(&iomodel->gridoniceshelf,NULL,NULL,iomodel_handle,"gridoniceshelf");
    4250
    43         for (i=0;i<iomodel->numberofvertices;i++){
     51        if(continuous_galerkin){
    4452
    45                 if(iomodel->my_vertices[i]){
     53                /*Build Nodes dataset (Continuous Galerkin)*/
     54                for (i=0;i<iomodel->numberofvertices;i++){
    4655
    47                         /*Add node to nodes dataset: */
    48                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,i+1,i,iomodel,BalancedthicknessAnalysisEnum));
     56                        if(iomodel->my_vertices[i]){
    4957
     58                                /*Add node to nodes dataset: */
     59                                nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,i+1,i,iomodel,BalancedthicknessAnalysisEnum));
     60
     61                        }
     62                }
     63        }
     64        else{
     65
     66                /*Build Nodes dataset -> 3 for each element (Discontinuous Galerkin)*/
     67                for (i=0;i<iomodel->numberofelements;i++){
     68                        for (j=0;j<3;j++){
     69
     70                                if(iomodel->my_nodes[3*i+j]){
     71
     72                                        //Get index of the vertex on which the current node is located
     73                                        vertex_id=(int)*(iomodel->elements+3*i+j); //(Matlab indexing)
     74                                        io_index=vertex_id-1;                      //(C indexing)
     75                                        ISSMASSERT(vertex_id>0 && vertex_id<=iomodel->numberofvertices);
     76
     77                                        //Compute Node id
     78                                        node_id=iomodel->nodecounter+3*i+j+1;
     79
     80                                        /*Add node to nodes dataset: */
     81                                        nodes->AddObject(new Node(node_id,node_id-1,vertex_id,io_index,iomodel,BalancedthicknessAnalysisEnum));
     82
     83                                }
     84                        }
    5085                }
    5186        }
    5287
    5388        /*Clean fetched data: */
     89        xfree((void**)&iomodel->elements);
    5490        xfree((void**)&iomodel->deadgrids);
    5591        xfree((void**)&iomodel->gridonbed);

  • issm/trunk/src/c/modules/ModelProcessorx/CreateDataSets.cpp

    r4739 r4765  
    9999                        break;
    100100
    101                 case Balancedthickness2AnalysisEnum:
    102                         CreateNodesBalancedthickness2(pnodes, iomodel,iomodel_handle);
    103                         CreateConstraintsBalancedthickness2(pconstraints,iomodel,iomodel_handle);
    104                         CreateLoadsBalancedthickness2(ploads,iomodel,iomodel_handle);
    105                         UpdateElementsBalancedthickness2(elements,iomodel,iomodel_handle,analysis_counter,analysis_type);
    106                         break;
    107101                case BalancedvelocitiesAnalysisEnum:
    108102                        CreateNodesBalancedvelocities(pnodes, iomodel,iomodel_handle);

  • issm/trunk/src/c/modules/ModelProcessorx/ModelProcessorx.h

    r4739 r4765  
    8484void    UpdateElementsBalancedthickness(Elements* elements,IoModel* iomodel_handle,ConstDataHandle iomodel_handle,int analysis_counter,int analysis_type);
    8585
    86 void    CreateNodesBalancedthickness2(Nodes** pnodes,IoModel* iomodel_handle,ConstDataHandle iomodel_handle);
    87 void    CreateConstraintsBalancedthickness2(Constraints** pconstraints,IoModel* iomodel,ConstDataHandle iomodel_handle);
    88 void  CreateLoadsBalancedthickness2(Loads** ploads, IoModel* iomodel, ConstDataHandle iomodel_handle);
    89 void    UpdateElementsBalancedthickness2(Elements* elements,IoModel* iomodel_handle,ConstDataHandle iomodel_handle,int analysis_counter,int analysis_type);
    90 
    9186/*balancedvelocities:*/
    9287void    CreateNodesBalancedvelocities(Nodes** pnodes,IoModel* iomodel_handle,ConstDataHandle iomodel_handle);

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

    r4739 r4765  
    394394        else if (analysis_type==BalancedthicknessAnalysisEnum){
    395395                InputUpdateFromSolutionBalancedthickness( solution);
    396         }
    397         else if (analysis_type==Balancedthickness2AnalysisEnum){
    398                 InputUpdateFromSolutionBalancedthickness2( solution);
    399396        }
    400397        else if (analysis_type==BalancedvelocitiesAnalysisEnum){
    … …  
    26462643        /*Add thickness as inputs to the tria element: */
    26472644        this->inputs->AddInput(new PentaVertexInput(ThicknessEnum,values));
    2648 }
    2649 /*}}}*/
    2650 /*FUNCTION Penta::InputUpdateFromSolutionBalancedthickness2 {{{1*/
    2651 void  Penta::InputUpdateFromSolutionBalancedthickness2(double* solution){
    2652         ISSMERROR(" not supported yet!");
    26532645}
    26542646/*}}}*/

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

    r4739 r4765  
    175175                void    InputUpdateFromSolutionAdjointStokes( double* solutiong);
    176176                void    InputUpdateFromSolutionBalancedthickness( double* solutiong);
    177                 void    InputUpdateFromSolutionBalancedthickness2( double* solutiong);
    178177                void    InputUpdateFromSolutionBalancedvelocities( double* solutiong);
    179178                void    InputUpdateFromSolutionDiagnosticHoriz( double* solutiong);

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

    r4753 r4765  
    394394                InputUpdateFromSolutionBalancedthickness( solution);
    395395        }
    396         else if (analysis_type==Balancedthickness2AnalysisEnum){
    397                 InputUpdateFromSolutionBalancedthickness2( solution);
    398         }
    399396        else if (analysis_type==BalancedvelocitiesAnalysisEnum){
    400397                InputUpdateFromSolutionBalancedvelocities( solution);
    … …  
    688685        }
    689686        else if (analysis_type==BalancedthicknessAnalysisEnum){
    690                 CreateKMatrixBalancedthickness( Kgg);
    691         }
    692         else if (analysis_type==Balancedthickness2AnalysisEnum){
    693                 CreateKMatrixBalancedthickness2( Kgg);
     687                if (GetElementType()==P1Enum)
     688                 CreateKMatrixBalancedthickness_CG( Kgg);
     689                else if (GetElementType()==P1DGEnum)
     690                 CreateKMatrixBalancedthickness_DG( Kgg);
     691                else
     692                 ISSMERROR("Element type %s not supported yet",EnumAsString(GetElementType()));
    694693        }
    695694        else if (analysis_type==BalancedvelocitiesAnalysisEnum){
    … …  
    737736        }
    738737        else if (analysis_type==BalancedthicknessAnalysisEnum){
    739                 CreatePVectorBalancedthickness( pg);
    740         }
    741         else if (analysis_type==Balancedthickness2AnalysisEnum){
    742                 CreatePVectorBalancedthickness2( pg);
     738                if (GetElementType()==P1Enum)
     739                 CreatePVectorBalancedthickness_CG( pg);
     740                else if (GetElementType()==P1DGEnum)
     741                 CreatePVectorBalancedthickness_DG( pg);
     742                else
     743                 ISSMERROR("Element type %s not supported yet",EnumAsString(GetElementType()));
    743744        }
    744745        else if (analysis_type==BalancedvelocitiesAnalysisEnum){
    … …  
    22642265
    22652266        /*Recover element type*/
    2266         if ((analysis_type==PrognosticAnalysisEnum && iomodel->prognostic_DG) || analysis_type==Balancedthickness2AnalysisEnum){
     2267        if ((analysis_type==PrognosticAnalysisEnum || analysis_type==BalancedthicknessAnalysisEnum) && iomodel->prognostic_DG){
    22672268
    22682269                /*P1 Discontinuous Galerkin*/
    … …  
    24732474
    24742475/*Tria specific routines: */
    2475 /*FUNCTION Tria::CreateKMatrixBalancedthickness {{{1*/
    2476 void  Tria::CreateKMatrixBalancedthickness(Mat Kgg){
     2476/*FUNCTION Tria::CreateKMatrixBalancedthickness_CG {{{1*/
     2477void  Tria::CreateKMatrixBalancedthickness_CG(Mat Kgg){
    24772478
    24782479        /* local declarations */
    … …  
    26292630}
    26302631/*}}}*/
    2631 /*FUNCTION Tria::CreateKMatrixBalancedthickness2 {{{1*/
    2632 void  Tria::CreateKMatrixBalancedthickness2(Mat Kgg){
     2632/*FUNCTION Tria::CreateKMatrixBalancedthickness_DG {{{1*/
     2633void  Tria::CreateKMatrixBalancedthickness_DG(Mat Kgg){
    26332634
    26342635        /* local declarations */
    … …  
    38603861}
    38613862/*}}}*/
    3862 /*FUNCTION Tria::CreatePVectorBalancedthickness {{{1*/
    3863 void  Tria::CreatePVectorBalancedthickness(Vec pg ){
     3863/*FUNCTION Tria::CreatePVectorBalancedthickness_CG{{{1*/
     3864void  Tria::CreatePVectorBalancedthickness_CG(Vec pg ){
    38643865
    38653866
    … …  
    39343935}
    39353936/*}}}*/
    3936 /*FUNCTION Tria::CreatePVectorBalancedthickness2 {{{1*/
    3937 void  Tria::CreatePVectorBalancedthickness2(Vec pg){
     3937/*FUNCTION Tria::CreatePVectorBalancedthickness_DG {{{1*/
     3938void  Tria::CreatePVectorBalancedthickness_DG(Vec pg){
    39383939
    39393940
    … …  
    60286029        double       values[numdof];
    60296030        double       thickness[numvertices];
    6030 
    6031         int          dummy;
    6032 
    6033         /*Get dof list: */
    6034         GetDofList(&doflist[0],&dummy);
    6035 
    6036         /*Use the dof list to index into the solution vector: */
    6037         for(i=0;i<numdof;i++){
    6038                 values[i]=solution[doflist[i]];
    6039         }
    6040 
    6041         /*Add thickness as inputs to the tria element: */
    6042         this->inputs->AddInput(new TriaVertexInput(ThicknessEnum,values));
    6043 }
    6044 /*}}}*/
    6045 /*FUNCTION Tria::InputUpdateFromSolutionBalancedthickness2 {{{1*/
    6046 void  Tria::InputUpdateFromSolutionBalancedthickness2(double* solution){
    6047         int i;
    6048 
    6049         const int    numvertices=3;
    6050         const int    numdofpervertex=1;
    6051         const int    numdof=numdofpervertex*numvertices;
    6052 
    6053         int          doflist[numdof];
    6054         double       values[numdof];
    60556031
    60566032        int          dummy;

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

    r4753 r4765  
    111111                /*}}}*/
    112112                /*Tria specific routines:{{{1*/
    113                 void      CreateKMatrixBalancedthickness(Mat Kgg);
    114                 void      CreateKMatrixBalancedthickness2(Mat Kgg);
     113                void      CreateKMatrixBalancedthickness_DG(Mat Kgg);
     114                void      CreateKMatrixBalancedthickness_CG(Mat Kgg);
    115115                void      CreateKMatrixBalancedvelocities(Mat Kgg);
    116116                void      CreateKMatrixDiagnosticHoriz(Mat Kgg);
    … …  
    123123                void      CreateKMatrixSlope(Mat Kgg);
    124124                void      CreateKMatrixThermal(Mat Kgg);
    125                 void      CreatePVectorBalancedthickness(Vec pg);
    126                 void      CreatePVectorBalancedthickness2(Vec pg);
     125                void      CreatePVectorBalancedthickness_DG(Vec pg);
     126                void      CreatePVectorBalancedthickness_CG(Vec pg);
    127127                void      CreatePVectorBalancedvelocities(Vec pg);
    128128                void      CreatePVectorDiagnosticBaseVert(Vec pg);
    … …  
    169169                void      InputUpdateFromSolutionPrognostic( double* solution);
    170170                void      InputUpdateFromSolutionBalancedthickness( double* solution);
    171                 void      InputUpdateFromSolutionBalancedthickness2( double* solution);
    172171                void      InputUpdateFromSolutionBalancedvelocities( double* solution);
    173172                bool      IsInput(int name);

  • issm/trunk/src/c/objects/Loads/Numericalflux.cpp

    r4739 r4765  
    415415                parameters->FindParam(&dt,DtEnum);
    416416        }
    417         else if (analysis_type==Balancedthickness2AnalysisEnum){
     417        else if (analysis_type==BalancedthicknessAnalysisEnum){
    418418                /*No transient term is involved*/
    419419                dt=1;
    … …  
    447447               
    448448                UdotN=vx*normal[0]+vy*normal[1];
    449                 if (fabs(UdotN)<1.0e-9 && analysis_type==Balancedthickness2AnalysisEnum) printf("Edge number %i has a flux very small (u.n = %g ), which could lead to unaccurate results\n",id,UdotN);
     449                if (fabs(UdotN)<1.0e-9 && analysis_type==BalancedthicknessAnalysisEnum) printf("Edge number %i has a flux very small (u.n = %g ), which could lead to unaccurate results\n",id,UdotN);
    450450
    451451                /*Get L and B: */
    … …  
    532532                parameters->FindParam(&dt,DtEnum);
    533533        }
    534         else if (analysis_type==Balancedthickness2AnalysisEnum){
     534        else if (analysis_type==BalancedthicknessAnalysisEnum){
    535535                /*No transient term is involved*/
    536536                dt=1;
    … …  
    659659                parameters->FindParam(&dt,DtEnum);
    660660        }
    661         else if (analysis_type==Balancedthickness2AnalysisEnum){
     661        else if (analysis_type==BalancedthicknessAnalysisEnum){
    662662                /*No transient term is involved*/
    663663                dt=1;

  • issm/trunk/src/c/shared/Dofs/DistributeNumDofs.cpp

    r4739 r4765  
    4040                numdofs=1;
    4141        }
    42         else if (analysis_type==Balancedthickness2AnalysisEnum){
    43                 numdofs=1;
    44         }
    4542        else if (analysis_type==BalancedvelocitiesAnalysisEnum){
    4643                numdofs=1;

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

    r4739 r4765  
    7676                        break;
    7777
    78                 case Balancedthickness2SolutionEnum:
    79                         numanalyses=1;
    80                         analyses=(int*)xmalloc(numanalyses*sizeof(int));
    81                         analyses[0]=Balancedthickness2AnalysisEnum;
    82                         solutioncore=&balancedthickness2_core;
    83                         break;
    84 
    8578                case BalancedvelocitiesSolutionEnum:
    8679                        numanalyses=1;

  • issm/trunk/src/c/solutions/solutions.h

    r4739 r4765  
    2323void prognostic_core(FemModel* femmodel);
    2424void balancedthickness_core(FemModel* femmodel);
    25 void balancedthickness2_core(FemModel* femmodel);
    2625void balancedvelocities_core(FemModel* femmodel);
    2726void slopecompute_core(FemModel* femmodel);
Note: See TracChangeset for help on using the changeset viewer.