Changeset 8287


Ignore:
Timestamp:
05/16/11 08:19:00 (14 years ago)
Author:
Mathieu Morlighem
Message:

moving Balanced* to Balance* TBC

Location:
issm/trunk/src
Files:
5 added
5 deleted
32 edited
4 moved

Legend:

Unmodified
Added
Removed

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

    r8224 r8287  
    3838        AdjointSolutionEnum,
    3939        AdjointHorizAnalysisEnum,
    40         AdjointBalancedthicknessAnalysisEnum,
     40        AdjointBalancethicknessAnalysisEnum,
    4141        ControlAnalysisEnum,
    4242        GradientAnalysisEnum,
    … …  
    6666        SurfaceSlopeYAnalysisEnum,
    6767        //prognostic
    68         BalancedthicknessSolutionEnum,
    69         BalancedthicknessAnalysisEnum,
    70         BalancedvelocitiesSolutionEnum,
    71         BalancedvelocitiesAnalysisEnum,
     68        BalancethicknessSolutionEnum,
     69        BalancethicknessAnalysisEnum,
     70        BalancevelocitiesSolutionEnum,
     71        BalancevelocitiesAnalysisEnum,
    7272        PrognosticSolutionEnum,
    7373        PrognosticAnalysisEnum,

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

    r8285 r8287  
    459459                                        ./modules/ModelProcessorx/Prognostic/CreateConstraintsPrognostic.cpp\
    460460                                        ./modules/ModelProcessorx/Prognostic/CreateLoadsPrognostic.cpp\
    461                                         ./modules/ModelProcessorx/Balancedthickness/UpdateElementsBalancedthickness.cpp\
    462                                         ./modules/ModelProcessorx/Balancedthickness/CreateNodesBalancedthickness.cpp\
    463                                         ./modules/ModelProcessorx/Balancedthickness/CreateConstraintsBalancedthickness.cpp\
    464                                         ./modules/ModelProcessorx/Balancedthickness/CreateLoadsBalancedthickness.cpp\
    465                                         ./modules/ModelProcessorx/Balancedvelocities/UpdateElementsBalancedvelocities.cpp\
    466                                         ./modules/ModelProcessorx/Balancedvelocities/CreateNodesBalancedvelocities.cpp\
    467                                         ./modules/ModelProcessorx/Balancedvelocities/CreateConstraintsBalancedvelocities.cpp\
    468                                         ./modules/ModelProcessorx/Balancedvelocities/CreateLoadsBalancedvelocities.cpp\
     461                                        ./modules/ModelProcessorx/Balancethickness/UpdateElementsBalancethickness.cpp\
     462                                        ./modules/ModelProcessorx/Balancethickness/CreateNodesBalancethickness.cpp\
     463                                        ./modules/ModelProcessorx/Balancethickness/CreateConstraintsBalancethickness.cpp\
     464                                        ./modules/ModelProcessorx/Balancethickness/CreateLoadsBalancethickness.cpp\
     465                                        ./modules/ModelProcessorx/Balancevelocities/UpdateElementsBalancevelocities.cpp\
     466                                        ./modules/ModelProcessorx/Balancevelocities/CreateNodesBalancevelocities.cpp\
     467                                        ./modules/ModelProcessorx/Balancevelocities/CreateConstraintsBalancevelocities.cpp\
     468                                        ./modules/ModelProcessorx/Balancevelocities/CreateLoadsBalancevelocities.cpp\
    469469                                        ./modules/ModelProcessorx/Qmu/CreateParametersQmu.cpp\
    470470                                        ./modules/NodesDofx/NodesDofx.h\
    … …  
    10971097                                        ./modules/ModelProcessorx/Prognostic/CreateConstraintsPrognostic.cpp\
    10981098                                        ./modules/ModelProcessorx/Prognostic/CreateLoadsPrognostic.cpp\
    1099                                         ./modules/ModelProcessorx/Balancedthickness/UpdateElementsBalancedthickness.cpp\
    1100                                         ./modules/ModelProcessorx/Balancedthickness/CreateNodesBalancedthickness.cpp\
    1101                                         ./modules/ModelProcessorx/Balancedthickness/CreateConstraintsBalancedthickness.cpp\
    1102                                         ./modules/ModelProcessorx/Balancedthickness/CreateLoadsBalancedthickness.cpp\
    1103                                         ./modules/ModelProcessorx/Balancedvelocities/UpdateElementsBalancedvelocities.cpp\
    1104                                         ./modules/ModelProcessorx/Balancedvelocities/CreateNodesBalancedvelocities.cpp\
    1105                                         ./modules/ModelProcessorx/Balancedvelocities/CreateConstraintsBalancedvelocities.cpp\
    1106                                         ./modules/ModelProcessorx/Balancedvelocities/CreateLoadsBalancedvelocities.cpp\
     1099                                        ./modules/ModelProcessorx/Balancethickness/UpdateElementsBalancethickness.cpp\
     1100                                        ./modules/ModelProcessorx/Balancethickness/CreateNodesBalancethickness.cpp\
     1101                                        ./modules/ModelProcessorx/Balancethickness/CreateConstraintsBalancethickness.cpp\
     1102                                        ./modules/ModelProcessorx/Balancethickness/CreateLoadsBalancethickness.cpp\
     1103                                        ./modules/ModelProcessorx/Balancevelocities/UpdateElementsBalancevelocities.cpp\
     1104                                        ./modules/ModelProcessorx/Balancevelocities/CreateNodesBalancevelocities.cpp\
     1105                                        ./modules/ModelProcessorx/Balancevelocities/CreateConstraintsBalancevelocities.cpp\
     1106                                        ./modules/ModelProcessorx/Balancevelocities/CreateLoadsBalancevelocities.cpp\
    11071107                                        ./modules/ModelProcessorx/Qmu/CreateParametersQmu.cpp\
    11081108                                        ./modules/NodesDofx/NodesDofx.h\

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

    r8224 r8287  
    3636                case AdjointSolutionEnum : return "AdjointSolution";
    3737                case AdjointHorizAnalysisEnum : return "AdjointHorizAnalysis";
    38                 case AdjointBalancedthicknessAnalysisEnum : return "AdjointBalancedthicknessAnalysis";
     38                case AdjointBalancethicknessAnalysisEnum : return "AdjointBalancethicknessAnalysis";
    3939                case ControlAnalysisEnum : return "ControlAnalysis";
    4040                case GradientAnalysisEnum : return "GradientAnalysis";
    … …  
    5959                case SurfaceSlopeXAnalysisEnum : return "SurfaceSlopeXAnalysis";
    6060                case SurfaceSlopeYAnalysisEnum : return "SurfaceSlopeYAnalysis";
    61                 case BalancedthicknessSolutionEnum : return "BalancedthicknessSolution";
    62                 case BalancedthicknessAnalysisEnum : return "BalancedthicknessAnalysis";
    63                 case BalancedvelocitiesSolutionEnum : return "BalancedvelocitiesSolution";
    64                 case BalancedvelocitiesAnalysisEnum : return "BalancedvelocitiesAnalysis";
     61                case BalancethicknessSolutionEnum : return "BalancethicknessSolution";
     62                case BalancethicknessAnalysisEnum : return "BalancethicknessAnalysis";
     63                case BalancevelocitiesSolutionEnum : return "BalancevelocitiesSolution";
     64                case BalancevelocitiesAnalysisEnum : return "BalancevelocitiesAnalysis";
    6565                case PrognosticSolutionEnum : return "PrognosticSolution";
    6666                case PrognosticAnalysisEnum : return "PrognosticAnalysis";

  • issm/trunk/src/c/modules/ModelProcessorx/Balancevelocities/CreateConstraintsBalancedvelocities.cpp

    r6237 r8287  
    11/*
    2  * CreateConstraintsBalancedvelocities.c:
     2 * CreateConstraintsBalancevelocities.c:
    33 */
    44
    … …  
    1111#include "../ModelProcessorx.h"
    1212
    13 void    CreateConstraintsBalancedvelocities(Constraints** pconstraints, IoModel* iomodel,ConstDataHandle iomodel_handle){
     13void    CreateConstraintsBalancevelocities(Constraints** pconstraints, IoModel* iomodel,ConstDataHandle iomodel_handle){
    1414
    1515        /*Intermediary*/
    … …  
    4040                                /*This grid needs to be spc'd: */
    4141                                constraints->AddObject(new Spc(iomodel->constraintcounter+count+1,iomodel->nodecounter+i+1,1,
    42                                                                 pow( pow(*(iomodel->spcvelocity+6*i+4),2.0) + pow(*(iomodel->spcvelocity+6*i+5),2.0) ,0.5),BalancedvelocitiesAnalysisEnum));
     42                                                                pow( pow(*(iomodel->spcvelocity+6*i+4),2.0) + pow(*(iomodel->spcvelocity+6*i+5),2.0) ,0.5),BalancevelocitiesAnalysisEnum));
    4343                                count++;
    4444                        }

  • issm/trunk/src/c/modules/ModelProcessorx/Balancevelocities/CreateLoadsBalancedvelocities.cpp

    r4236 r8287  
    1 /*! \file CreateLoadsBalancedvelocities.c:
     1/*! \file CreateLoadsBalancevelocities.c:
    22 */
    33
    … …  
    1111#include "../ModelProcessorx.h"
    1212
    13 void    CreateLoadsBalancedvelocities(Loads** ploads, IoModel* iomodel,ConstDataHandle iomodel_handle){
     13void    CreateLoadsBalancevelocities(Loads** ploads, IoModel* iomodel,ConstDataHandle iomodel_handle){
    1414
    1515        /*DataSet*/

  • issm/trunk/src/c/modules/ModelProcessorx/Balancevelocities/CreateNodesBalancedvelocities.cpp

    r7516 r8287  
    11/*
    2  * CreateNodesBalancedvelocities.c:
     2 * CreateNodesBalancevelocities.c:
    33 */
    44
    … …  
    1313#include "../ModelProcessorx.h"
    1414
    15 void    CreateNodesBalancedvelocities(Nodes** pnodes, IoModel* iomodel,ConstDataHandle iomodel_handle){
     15void    CreateNodesBalancevelocities(Nodes** pnodes, IoModel* iomodel,ConstDataHandle iomodel_handle){
    1616
    1717        /*Intermediary*/
    … …  
    4343                       
    4444                        /*Add node to nodes dataset: */
    45                         nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,i+1,i,iomodel,BalancedvelocitiesAnalysisEnum));
     45                        nodes->AddObject(new Node(iomodel->nodecounter+i+1,i,i+1,i,iomodel,BalancevelocitiesAnalysisEnum));
    4646
    4747                }

  • issm/trunk/src/c/modules/ModelProcessorx/Balancevelocities/UpdateElementsBalancedvelocities.cpp

    r4441 r8287  
    11/*
    2  * UpdateElementsBalancedvelocities:
     2 * UpdateElementsBalancevelocities:
    33 */
    44
    … …  
    1313#include "../ModelProcessorx.h"
    1414
    15 void    UpdateElementsBalancedvelocities(Elements* elements, IoModel* iomodel,ConstDataHandle iomodel_handle,int analysis_counter,int analysis_type){
     15void    UpdateElementsBalancevelocities(Elements* elements, IoModel* iomodel,ConstDataHandle iomodel_handle,int analysis_counter,int analysis_type){
    1616
    1717        /*Intermediary*/

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

    r8224 r8287  
    9494                        break;
    9595
    96                 case BalancedthicknessAnalysisEnum:
    97                         CreateNodesBalancedthickness(pnodes, iomodel,iomodel_handle);
    98                         CreateConstraintsBalancedthickness(pconstraints,iomodel,iomodel_handle);
    99                         CreateLoadsBalancedthickness(ploads,iomodel,iomodel_handle);
    100                         UpdateElementsBalancedthickness(elements,iomodel,iomodel_handle,analysis_counter,analysis_type);
     96                case BalancethicknessAnalysisEnum:
     97                        CreateNodesBalancethickness(pnodes, iomodel,iomodel_handle);
     98                        CreateConstraintsBalancethickness(pconstraints,iomodel,iomodel_handle);
     99                        CreateLoadsBalancethickness(ploads,iomodel,iomodel_handle);
     100                        UpdateElementsBalancethickness(elements,iomodel,iomodel_handle,analysis_counter,analysis_type);
    101101                        break;
    102102
    103                 case BalancedvelocitiesAnalysisEnum:
    104                         CreateNodesBalancedvelocities(pnodes, iomodel,iomodel_handle);
    105                         CreateConstraintsBalancedvelocities(pconstraints,iomodel,iomodel_handle);
    106                         CreateLoadsBalancedvelocities(ploads,iomodel,iomodel_handle);
    107                         UpdateElementsBalancedvelocities(elements,iomodel,iomodel_handle,analysis_counter,analysis_type);
     103                case BalancevelocitiesAnalysisEnum:
     104                        CreateNodesBalancevelocities(pnodes, iomodel,iomodel_handle);
     105                        CreateConstraintsBalancevelocities(pconstraints,iomodel,iomodel_handle);
     106                        CreateLoadsBalancevelocities(ploads,iomodel,iomodel_handle);
     107                        UpdateElementsBalancevelocities(elements,iomodel,iomodel_handle,analysis_counter,analysis_type);
    108108                        break;
    109109                default:

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

    r8224 r8287  
    8484                numdofs=1;
    8585        }
    86         else if (analysis_type==BalancedthicknessAnalysisEnum){
     86        else if (analysis_type==BalancethicknessAnalysisEnum){
    8787                numdofs=1;
    8888        }
    89         else if (analysis_type==BalancedvelocitiesAnalysisEnum){
     89        else if (analysis_type==BalancevelocitiesAnalysisEnum){
    9090                numdofs=1;
    9191        }

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

    r7640 r8287  
    8181
    8282/*balancedthickness:*/
    83 void    CreateNodesBalancedthickness(Nodes** pnodes,IoModel* iomodel,ConstDataHandle iomodel_handle);
    84 void    CreateConstraintsBalancedthickness(Constraints** pconstraints,IoModel* iomodel,ConstDataHandle iomodel_handle);
    85 void  CreateLoadsBalancedthickness(Loads** ploads, IoModel* iomodel, ConstDataHandle iomodel_handle);
    86 void    UpdateElementsBalancedthickness(Elements* elements,IoModel* iomodel,ConstDataHandle iomodel_handle,int analysis_counter,int analysis_type);
     83void    CreateNodesBalancethickness(Nodes** pnodes,IoModel* iomodel,ConstDataHandle iomodel_handle);
     84void    CreateConstraintsBalancethickness(Constraints** pconstraints,IoModel* iomodel,ConstDataHandle iomodel_handle);
     85void  CreateLoadsBalancethickness(Loads** ploads, IoModel* iomodel, ConstDataHandle iomodel_handle);
     86void    UpdateElementsBalancethickness(Elements* elements,IoModel* iomodel,ConstDataHandle iomodel_handle,int analysis_counter,int analysis_type);
    8787
    8888/*balancedvelocities:*/
    89 void    CreateNodesBalancedvelocities(Nodes** pnodes,IoModel* iomodel,ConstDataHandle iomodel_handle);
    90 void    CreateConstraintsBalancedvelocities(Constraints** pconstraints,IoModel* iomodel,ConstDataHandle iomodel_handle);
    91 void  CreateLoadsBalancedvelocities(Loads** ploads, IoModel* iomodel, ConstDataHandle iomodel_handle);
    92 void    UpdateElementsBalancedvelocities(Elements* elements,IoModel* iomodel,ConstDataHandle iomodel_handle,int analysis_counter,int analysis_type);
     89void    CreateNodesBalancevelocities(Nodes** pnodes,IoModel* iomodel,ConstDataHandle iomodel_handle);
     90void    CreateConstraintsBalancevelocities(Constraints** pconstraints,IoModel* iomodel,ConstDataHandle iomodel_handle);
     91void  CreateLoadsBalancevelocities(Loads** ploads, IoModel* iomodel, ConstDataHandle iomodel_handle);
     92void    UpdateElementsBalancevelocities(Elements* elements,IoModel* iomodel,ConstDataHandle iomodel_handle,int analysis_counter,int analysis_type);
    9393
    9494/*partitioning: */

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

    r8224 r8287  
    3434        else if (strcmp(name,"AdjointSolution")==0) return AdjointSolutionEnum;
    3535        else if (strcmp(name,"AdjointHorizAnalysis")==0) return AdjointHorizAnalysisEnum;
    36         else if (strcmp(name,"AdjointBalancedthicknessAnalysis")==0) return AdjointBalancedthicknessAnalysisEnum;
     36        else if (strcmp(name,"AdjointBalancethicknessAnalysis")==0) return AdjointBalancethicknessAnalysisEnum;
    3737        else if (strcmp(name,"ControlAnalysis")==0) return ControlAnalysisEnum;
    3838        else if (strcmp(name,"GradientAnalysis")==0) return GradientAnalysisEnum;
    … …  
    5757        else if (strcmp(name,"SurfaceSlopeXAnalysis")==0) return SurfaceSlopeXAnalysisEnum;
    5858        else if (strcmp(name,"SurfaceSlopeYAnalysis")==0) return SurfaceSlopeYAnalysisEnum;
    59         else if (strcmp(name,"BalancedthicknessSolution")==0) return BalancedthicknessSolutionEnum;
    60         else if (strcmp(name,"BalancedthicknessAnalysis")==0) return BalancedthicknessAnalysisEnum;
    61         else if (strcmp(name,"BalancedvelocitiesSolution")==0) return BalancedvelocitiesSolutionEnum;
    62         else if (strcmp(name,"BalancedvelocitiesAnalysis")==0) return BalancedvelocitiesAnalysisEnum;
     59        else if (strcmp(name,"BalancethicknessSolution")==0) return BalancethicknessSolutionEnum;
     60        else if (strcmp(name,"BalancethicknessAnalysis")==0) return BalancethicknessAnalysisEnum;
     61        else if (strcmp(name,"BalancevelocitiesSolution")==0) return BalancevelocitiesSolutionEnum;
     62        else if (strcmp(name,"BalancevelocitiesAnalysis")==0) return BalancevelocitiesAnalysisEnum;
    6363        else if (strcmp(name,"PrognosticSolution")==0) return PrognosticSolutionEnum;
    6464        else if (strcmp(name,"PrognosticAnalysis")==0) return PrognosticAnalysisEnum;

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

    r8224 r8287  
    561561                        Ke=CreateKMatrixPrognostic();
    562562                        break;
    563                 case BalancedthicknessAnalysisEnum:
    564                         Ke=CreateKMatrixBalancedthickness();
     563                case BalancethicknessAnalysisEnum:
     564                        Ke=CreateKMatrixBalancethickness();
    565565                        break;
    566                 case BalancedvelocitiesAnalysisEnum:
    567                         Ke=CreateKMatrixBalancedvelocities();
     566                case BalancevelocitiesAnalysisEnum:
     567                        Ke=CreateKMatrixBalancevelocities();
    568568                        break;
    569569                case ThermalAnalysisEnum:
    … …  
    589589}
    590590/*}}}*/
    591 /*FUNCTION Penta::CreateKMatrixBalancedthickness {{{1*/
    592 ElementMatrix* Penta::CreateKMatrixBalancedthickness(void){
     591/*FUNCTION Penta::CreateKMatrixBalancethickness {{{1*/
     592ElementMatrix* Penta::CreateKMatrixBalancethickness(void){
    593593
    594594        /*Figure out if this penta is collapsed. If so, then bailout, except if it is at the
    … …  
    603603        /*Spawn Tria element from the base of the Penta: */
    604604        Tria* tria=(Tria*)SpawnTria(0,1,2); //grids 0, 1 and 2 make the new tria.
    605         ElementMatrix* Ke=tria->CreateKMatrixBalancedthickness();
     605        ElementMatrix* Ke=tria->CreateKMatrixBalancethickness();
    606606        delete tria->matice; delete tria;
    607607
    … …  
    614614}
    615615/*}}}*/
    616 /*FUNCTION Penta::CreateKMatrixBalancedvelocities {{{1*/
    617 ElementMatrix* Penta::CreateKMatrixBalancedvelocities(void){
     616/*FUNCTION Penta::CreateKMatrixBalancevelocities {{{1*/
     617ElementMatrix* Penta::CreateKMatrixBalancevelocities(void){
    618618
    619619        /*Figure out if this penta is collapsed. If so, then bailout, except if it is at the
    … …  
    628628        /*Spawn Tria element from the base of the Penta: */
    629629        Tria* tria=(Tria*)SpawnTria(0,1,2); //grids 0, 1 and 2 make the new tria.
    630         ElementMatrix* Ke=tria->CreateKMatrixBalancedvelocities();
     630        ElementMatrix* Ke=tria->CreateKMatrixBalancevelocities();
    631631        delete tria->matice; delete tria;
    632632
    … …  
    19081908                        pe=CreatePVectorPrognostic();
    19091909                        break;
    1910                 case BalancedthicknessAnalysisEnum:
    1911                         pe=CreatePVectorBalancedthickness();
     1910                case BalancethicknessAnalysisEnum:
     1911                        pe=CreatePVectorBalancethickness();
    19121912                        break;
    1913                 case BalancedvelocitiesAnalysisEnum:
    1914                         pe=CreatePVectorBalancedvelocities();
     1913                case BalancevelocitiesAnalysisEnum:
     1914                        pe=CreatePVectorBalancevelocities();
    19151915                        break;
    19161916                case ThermalAnalysisEnum:
    … …  
    19791979}
    19801980/*}}}*/
    1981 /*FUNCTION Penta::CreatePVectorBalancedthickness {{{1*/
    1982 ElementVector* Penta::CreatePVectorBalancedthickness(void){
     1981/*FUNCTION Penta::CreatePVectorBalancethickness {{{1*/
     1982ElementVector* Penta::CreatePVectorBalancethickness(void){
    19831983
    19841984        if (!IsOnBed()) return NULL;
    … …  
    19901990        /*Call Tria function*/
    19911991        Tria* tria=(Tria*)SpawnTria(0,1,2); //grids 0, 1 and 2 make the new tria.
    1992         ElementVector* pe=tria->CreatePVectorBalancedthickness();
     1992        ElementVector* pe=tria->CreatePVectorBalancethickness();
    19931993        delete tria->matice; delete tria;
    19941994
    … …  
    20012001}
    20022002/*}}}*/
    2003 /*FUNCTION Penta::CreatePVectorBalancedvelocities {{{1*/
    2004 ElementVector* Penta::CreatePVectorBalancedvelocities(void){
     2003/*FUNCTION Penta::CreatePVectorBalancevelocities {{{1*/
     2004ElementVector* Penta::CreatePVectorBalancevelocities(void){
    20052005
    20062006        if (!IsOnBed()) return NULL;
    … …  
    20122012        /*Call Tria function*/
    20132013        Tria* tria=(Tria*)SpawnTria(0,1,2); //grids 0, 1 and 2 make the new tria.
    2014         ElementVector* pe=tria->CreatePVectorBalancedvelocities();
     2014        ElementVector* pe=tria->CreatePVectorBalancevelocities();
    20152015        delete tria->matice; delete tria;
    20162016
    … …  
    41834183                InputUpdateFromSolutionPrognostic(solution);
    41844184        }
    4185         else if (analysis_type==BalancedthicknessAnalysisEnum){
     4185        else if (analysis_type==BalancethicknessAnalysisEnum){
    41864186                InputUpdateFromSolutionOneDofCollapsed(solution,ThicknessEnum);
    41874187        }
    4188         else if (analysis_type==BalancedvelocitiesAnalysisEnum){
     4188        else if (analysis_type==BalancevelocitiesAnalysisEnum){
    41894189                InputUpdateFromSolutionOneDofCollapsed(solution,VelEnum);
    41904190        }
    … …  
    61556155
    61566156        /*Recover element type*/
    6157         if ((analysis_type==PrognosticAnalysisEnum || analysis_type==BalancedthicknessAnalysisEnum) && iomodel->prognostic_DG){
     6157        if ((analysis_type==PrognosticAnalysisEnum || analysis_type==BalancethicknessAnalysisEnum) && iomodel->prognostic_DG){
    61586158                /*P1 Discontinuous Galerkin*/
    61596159                penta_type=P1DGEnum;

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

    r7391 r8287  
    137137                /*Penta specific routines:{{{1*/
    138138                void      BedNormal(double* bed_normal, double xyz_list[3][3]);
    139                 ElementMatrix* CreateKMatrixBalancedthickness(void);
    140                 ElementMatrix* CreateKMatrixBalancedvelocities(void);
     139                ElementMatrix* CreateKMatrixBalancethickness(void);
     140                ElementMatrix* CreateKMatrixBalancevelocities(void);
    141141                ElementMatrix* CreateKMatrixCouplingMacAyealPattyn(void);
    142142                ElementMatrix* CreateKMatrixCouplingMacAyealPattynViscous(void);
    … …  
    172172                ElementMatrix* CreateKMatrixThermalVolume(void);
    173173                ElementMatrix* CreateKMatrixThermalShelf(void);
    174                 ElementVector* CreatePVectorBalancedthickness(void);
    175                 ElementVector* CreatePVectorBalancedvelocities(void);
     174                ElementVector* CreatePVectorBalancethickness(void);
     175                ElementVector* CreatePVectorBalancevelocities(void);
    176176                ElementVector* CreatePVectorAdjointHoriz(void);
    177177                ElementVector* CreatePVectorAdjointMacAyeal(void);

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

    r8224 r8287  
    408408                        Ke=CreateKMatrixHydrology();
    409409                        break;
    410                 case BalancedthicknessAnalysisEnum:
    411                         Ke=CreateKMatrixBalancedthickness();
     410                case BalancethicknessAnalysisEnum:
     411                        Ke=CreateKMatrixBalancethickness();
    412412                        break;
    413                 case AdjointBalancedthicknessAnalysisEnum:
    414                         Ke=CreateKMatrixAdjointBalancedthickness();
     413                case AdjointBalancethicknessAnalysisEnum:
     414                        Ke=CreateKMatrixAdjointBalancethickness();
    415415                        break;
    416                 case BalancedvelocitiesAnalysisEnum:
    417                         Ke=CreateKMatrixBalancedvelocities();
     416                case BalancevelocitiesAnalysisEnum:
     417                        Ke=CreateKMatrixBalancevelocities();
    418418                        break;
    419419                default:
    … …  
    428428}
    429429/*}}}*/
    430 /*FUNCTION Tria::CreateKMatrixAdjointBalancedthickness {{{1*/
    431 ElementMatrix* Tria::CreateKMatrixAdjointBalancedthickness(void){
     430/*FUNCTION Tria::CreateKMatrixAdjointBalancethickness {{{1*/
     431ElementMatrix* Tria::CreateKMatrixAdjointBalancethickness(void){
    432432
    433433        ElementMatrix* Ke=NULL;
    … …  
    436436        switch(GetElementType()){
    437437                case P1Enum:
    438                         Ke=CreateKMatrixBalancedthickness_CG();
     438                        Ke=CreateKMatrixBalancethickness_CG();
    439439                        break;
    440440                case P1DGEnum:
    441                         Ke=CreateKMatrixBalancedthickness_DG();
     441                        Ke=CreateKMatrixBalancethickness_DG();
    442442                        break;
    443443                default:
    … …  
    450450}
    451451/*}}}*/
    452 /*FUNCTION Tria::CreateKMatrixBalancedthickness {{{1*/
    453 ElementMatrix* Tria::CreateKMatrixBalancedthickness(void){
     452/*FUNCTION Tria::CreateKMatrixBalancethickness {{{1*/
     453ElementMatrix* Tria::CreateKMatrixBalancethickness(void){
    454454
    455455        switch(GetElementType()){
    456456                case P1Enum:
    457                         return CreateKMatrixBalancedthickness_CG();
     457                        return CreateKMatrixBalancethickness_CG();
    458458                case P1DGEnum:
    459                         return CreateKMatrixBalancedthickness_DG();
     459                        return CreateKMatrixBalancethickness_DG();
    460460                default:
    461461                        _error_("Element type %s not supported yet",EnumToStringx(GetElementType()));
    … …  
    464464}
    465465/*}}}*/
    466 /*FUNCTION Tria::CreateKMatrixBalancedthickness_CG {{{1*/
    467 ElementMatrix* Tria::CreateKMatrixBalancedthickness_CG(void){
     466/*FUNCTION Tria::CreateKMatrixBalancethickness_CG {{{1*/
     467ElementMatrix* Tria::CreateKMatrixBalancethickness_CG(void){
    468468
    469469        /*Constants*/
    … …  
    576576}
    577577/*}}}*/
    578 /*FUNCTION Tria::CreateKMatrixBalancedthickness_DG {{{1*/
    579 ElementMatrix* Tria::CreateKMatrixBalancedthickness_DG(void){
     578/*FUNCTION Tria::CreateKMatrixBalancethickness_DG {{{1*/
     579ElementMatrix* Tria::CreateKMatrixBalancethickness_DG(void){
    580580
    581581        /*Constants*/
    … …  
    633633}
    634634/*}}}*/
    635 /*FUNCTION Tria::CreateKMatrixBalancedvelocities {{{1*/
    636 ElementMatrix* Tria::CreateKMatrixBalancedvelocities(void){
     635/*FUNCTION Tria::CreateKMatrixBalancevelocities {{{1*/
     636ElementMatrix* Tria::CreateKMatrixBalancevelocities(void){
    637637
    638638        /*Constants*/
    … …  
    16181618                        pe=CreatePVectorHydrology();
    16191619                        break;
    1620                 case BalancedthicknessAnalysisEnum:
    1621                         pe=CreatePVectorBalancedthickness();
     1620                case BalancethicknessAnalysisEnum:
     1621                        pe=CreatePVectorBalancethickness();
    16221622                        break;
    1623                 case AdjointBalancedthicknessAnalysisEnum:
    1624                         pe=CreatePVectorAdjointBalancedthickness();
     1623                case AdjointBalancethicknessAnalysisEnum:
     1624                        pe=CreatePVectorAdjointBalancethickness();
    16251625                        break;
    1626                 case BalancedvelocitiesAnalysisEnum:
    1627                         pe=CreatePVectorBalancedvelocities();
     1626                case BalancevelocitiesAnalysisEnum:
     1627                        pe=CreatePVectorBalancevelocities();
    16281628                        break;
    16291629                default:
    … …  
    16381638}
    16391639/*}}}*/
    1640 /*FUNCTION Tria::CreatePVectorBalancedthickness{{{1*/
    1641 ElementVector* Tria::CreatePVectorBalancedthickness(void){
     1640/*FUNCTION Tria::CreatePVectorBalancethickness{{{1*/
     1641ElementVector* Tria::CreatePVectorBalancethickness(void){
    16421642
    16431643        switch(GetElementType()){
    16441644                case P1Enum:
    1645                         return CreatePVectorBalancedthickness_CG();
     1645                        return CreatePVectorBalancethickness_CG();
    16461646                        break;
    16471647                case P1DGEnum:
    1648                         return CreatePVectorBalancedthickness_DG();
     1648                        return CreatePVectorBalancethickness_DG();
    16491649                default:
    16501650                        _error_("Element type %s not supported yet",EnumToStringx(GetElementType()));
    … …  
    16521652}
    16531653/*}}}*/
    1654 /*FUNCTION Tria::CreatePVectorBalancedthickness_CG{{{1*/
    1655 ElementVector* Tria::CreatePVectorBalancedthickness_CG(void){
     1654/*FUNCTION Tria::CreatePVectorBalancethickness_CG{{{1*/
     1655ElementVector* Tria::CreatePVectorBalancethickness_CG(void){
    16561656
    16571657        /*Constants*/
    … …  
    16951695}
    16961696/*}}}*/
    1697 /*FUNCTION Tria::CreatePVectorBalancedthickness_DG {{{1*/
    1698 ElementVector* Tria::CreatePVectorBalancedthickness_DG(void){
     1697/*FUNCTION Tria::CreatePVectorBalancethickness_DG {{{1*/
     1698ElementVector* Tria::CreatePVectorBalancethickness_DG(void){
    16991699
    17001700        /*Constants*/
    … …  
    17381738}
    17391739/*}}}*/
    1740 /*FUNCTION Tria::CreatePVectorBalancedvelocities {{{1*/
    1741 ElementVector* Tria::CreatePVectorBalancedvelocities(void){
     1740/*FUNCTION Tria::CreatePVectorBalancevelocities {{{1*/
     1741ElementVector* Tria::CreatePVectorBalancevelocities(void){
    17421742
    17431743        /*Constants*/
    … …  
    19051905}
    19061906/*}}}*/
    1907 /*FUNCTION Tria::CreatePVectorAdjointBalancedthickness{{{1*/
    1908 ElementVector* Tria::CreatePVectorAdjointBalancedthickness(void){
     1907/*FUNCTION Tria::CreatePVectorAdjointBalancethickness{{{1*/
     1908ElementVector* Tria::CreatePVectorAdjointBalancethickness(void){
    19091909
    19101910        /*Constants*/
    … …  
    19461946
    19471947                for(i=0;i<numdof;i++) pe->values[i]+=(thicknessobs-thickness)*weight*Jdet*gauss->weight*l1l2l3[i];
    1948                 /*Regularization of the constraint: 20000000*/
    1949                 //for(i=0;i<numdof;i++) pe->values[i]+= - 1*2000000*dH[0]*dbasis[0][i]*Jdet*gauss->weight;
    1950                 //for(i=0;i<numdof;i++) pe->values[i]+= - 1*2000000*dH[1]*dbasis[1][i]*weight*Jdet*gauss->weight;
     1948                /*Regularization of the constraint: 2000000 79 N*/
     1949                //for(i=0;i<numdof;i++) pe->values[i]+= - 1*100000*dH[0]*dbasis[0][i]*Jdet*gauss->weight;
     1950                //for(i=0;i<numdof;i++) pe->values[i]+= - 1*100000*dH[1]*dbasis[1][i]*weight*Jdet*gauss->weight;
    19511951        }
    19521952
    … …  
    40224022                        InputUpdateFromSolutionHydrology(solution);
    40234023                        break;
    4024                 case BalancedthicknessAnalysisEnum:
     4024                case BalancethicknessAnalysisEnum:
    40254025                        InputUpdateFromSolutionOneDof(solution,ThicknessEnum);
    40264026                        break;
    4027                 case AdjointBalancedthicknessAnalysisEnum:
    4028                         InputUpdateFromSolutionAdjointBalancedthickness( solution);
     4027                case AdjointBalancethicknessAnalysisEnum:
     4028                        InputUpdateFromSolutionAdjointBalancethickness( solution);
    40294029                        break;
    4030                 case BalancedvelocitiesAnalysisEnum:
     4030                case BalancevelocitiesAnalysisEnum:
    40314031                        InputUpdateFromSolutionOneDof(solution,VelEnum);
    40324032                        break;
    … …  
    40364036}
    40374037/*}}}*/
    4038 /*FUNCTION Tria::InputUpdateFromSolutionAdjointBalancedthickness {{{1*/
    4039 void  Tria::InputUpdateFromSolutionAdjointBalancedthickness(double* solution){
     4038/*FUNCTION Tria::InputUpdateFromSolutionAdjointBalancethickness {{{1*/
     4039void  Tria::InputUpdateFromSolutionAdjointBalancethickness(double* solution){
    40404040
    40414041        const int numdof=NDOF1*NUMVERTICES;
    … …  
    57825782
    57835783        /*Recover element type*/
    5784         if ((analysis_type==PrognosticAnalysisEnum || analysis_type==BalancedthicknessAnalysisEnum) && iomodel->prognostic_DG){
     5784        if ((analysis_type==PrognosticAnalysisEnum || analysis_type==BalancethicknessAnalysisEnum) && iomodel->prognostic_DG){
    57855785                /*P1 Discontinuous Galerkin*/
    57865786                tria_type=P1DGEnum;

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

    r7640 r8287  
    141141                /*}}}*/
    142142                /*Tria specific routines:{{{1*/
    143                 ElementMatrix* CreateKMatrixAdjointBalancedthickness(void);
    144                 ElementMatrix* CreateKMatrixBalancedthickness(void);
    145                 ElementMatrix* CreateKMatrixBalancedthickness_DG(void);
    146                 ElementMatrix* CreateKMatrixBalancedthickness_CG(void);
    147                 ElementMatrix* CreateKMatrixBalancedvelocities(void);
     143                ElementMatrix* CreateKMatrixAdjointBalancethickness(void);
     144                ElementMatrix* CreateKMatrixBalancethickness(void);
     145                ElementMatrix* CreateKMatrixBalancethickness_DG(void);
     146                ElementMatrix* CreateKMatrixBalancethickness_CG(void);
     147                ElementMatrix* CreateKMatrixBalancevelocities(void);
    148148                ElementMatrix* CreateKMatrixDiagnosticMacAyeal(void);
    149149                ElementMatrix* CreateKMatrixDiagnosticMacAyealViscous(void);
    … …  
    160160                ElementMatrix* CreateKMatrixSlope(void);
    161161                ElementMatrix* CreateKMatrixThermal(void);
    162                 ElementVector* CreatePVectorBalancedthickness(void);
    163                 ElementVector* CreatePVectorBalancedthickness_DG(void);
    164                 ElementVector* CreatePVectorBalancedthickness_CG(void);
    165                 ElementVector* CreatePVectorBalancedvelocities(void);
     162                ElementVector* CreatePVectorBalancethickness(void);
     163                ElementVector* CreatePVectorBalancethickness_DG(void);
     164                ElementVector* CreatePVectorBalancethickness_CG(void);
     165                ElementVector* CreatePVectorBalancevelocities(void);
    166166                ElementVector* CreatePVectorDiagnosticBaseVert(void);
    167167                ElementVector* CreatePVectorDiagnosticMacAyeal(void);
    168168                ElementVector* CreatePVectorAdjointHoriz(void);
    169169                ElementVector* CreatePVectorAdjointStokes(void);
    170                 ElementVector* CreatePVectorAdjointBalancedthickness(void);
     170                ElementVector* CreatePVectorAdjointBalancethickness(void);
    171171                ElementVector* CreatePVectorDiagnosticHutter(void);
    172172                ElementVector* CreatePVectorHydrology(void);
    … …  
    190190                void    GetStrainRate2d(double* epsilon,double* xyz_list, GaussTria* gauss, Input* vx_input, Input* vy_input);
    191191                void      GradjDragStokes(Vec gradient);
    192                 void      InputUpdateFromSolutionAdjointBalancedthickness( double* solution);
     192                void      InputUpdateFromSolutionAdjointBalancethickness( double* solution);
    193193                void      InputUpdateFromSolutionAdjointHoriz( double* solution);
    194194                void      InputUpdateFromSolutionDiagnosticHoriz( double* solution);

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

    r8224 r8287  
    339339                        Ke=CreateKMatrixPrognostic();
    340340                        break;
    341                 case BalancedthicknessAnalysisEnum:
    342                         Ke=CreateKMatrixBalancedthickness();
     341                case BalancethicknessAnalysisEnum:
     342                        Ke=CreateKMatrixBalancethickness();
    343343                        break;
    344                 case AdjointBalancedthicknessAnalysisEnum:
    345                         Ke=CreateKMatrixAdjointBalancedthickness();
     344                case AdjointBalancethicknessAnalysisEnum:
     345                        Ke=CreateKMatrixAdjointBalancethickness();
    346346                        break;
    347347                default:
    … …  
    369369                        pe=CreatePVectorPrognostic();
    370370                        break;
    371                 case BalancedthicknessAnalysisEnum:
    372                         pe=CreatePVectorBalancedthickness();
     371                case BalancethicknessAnalysisEnum:
     372                        pe=CreatePVectorBalancethickness();
    373373                        break;
    374                 case AdjointBalancedthicknessAnalysisEnum:
    375                         pe=CreatePVectorAdjointBalancedthickness();
     374                case AdjointBalancethicknessAnalysisEnum:
     375                        pe=CreatePVectorAdjointBalancethickness();
    376376                        break;
    377377                default:
    … …  
    564564}
    565565/*}}}*/
    566 /*FUNCTION Numericalflux::CreateKMatrixBalancedthickness{{{1*/
    567 ElementMatrix* Numericalflux::CreateKMatrixBalancedthickness(void){
     566/*FUNCTION Numericalflux::CreateKMatrixBalancethickness{{{1*/
     567ElementMatrix* Numericalflux::CreateKMatrixBalancethickness(void){
    568568
    569569        int type;
    … …  
    572572        switch(type){
    573573                case InternalEnum:
    574                         return CreateKMatrixBalancedthicknessInternal();
     574                        return CreateKMatrixBalancethicknessInternal();
    575575                case BoundaryEnum:
    576                         return CreateKMatrixBalancedthicknessBoundary();
     576                        return CreateKMatrixBalancethicknessBoundary();
    577577                default:
    578578                        _error_("type not supported yet");
    … …  
    580580}
    581581/*}}}*/
    582 /*FUNCTION Numericalflux::CreateKMatrixBalancedthicknessInternal {{{1*/
    583 ElementMatrix* Numericalflux::CreateKMatrixBalancedthicknessInternal(void){
     582/*FUNCTION Numericalflux::CreateKMatrixBalancethicknessInternal {{{1*/
     583ElementMatrix* Numericalflux::CreateKMatrixBalancethicknessInternal(void){
    584584
    585585        /* constants*/
    … …  
    644644}
    645645/*}}}*/
    646 /*FUNCTION Numericalflux::CreateKMatrixBalancedthicknessBoundary {{{1*/
    647 ElementMatrix* Numericalflux::CreateKMatrixBalancedthicknessBoundary(void){
     646/*FUNCTION Numericalflux::CreateKMatrixBalancethicknessBoundary {{{1*/
     647ElementMatrix* Numericalflux::CreateKMatrixBalancethicknessBoundary(void){
    648648
    649649        /* constants*/
    … …  
    715715}
    716716/*}}}*/
    717 /*FUNCTION Numericalflux::CreateKMatrixAdjointBalancedthickness{{{1*/
    718 ElementMatrix* Numericalflux::CreateKMatrixAdjointBalancedthickness(void){
     717/*FUNCTION Numericalflux::CreateKMatrixAdjointBalancethickness{{{1*/
     718ElementMatrix* Numericalflux::CreateKMatrixAdjointBalancethickness(void){
    719719
    720720        int type;
    … …  
    723723        switch(type){
    724724                case InternalEnum:
    725                         return CreateKMatrixAdjointBalancedthicknessInternal();
     725                        return CreateKMatrixAdjointBalancethicknessInternal();
    726726                case BoundaryEnum:
    727                         return CreateKMatrixAdjointBalancedthicknessBoundary();
     727                        return CreateKMatrixAdjointBalancethicknessBoundary();
    728728                default:
    729729                        _error_("type not supported yet");
    … …  
    731731}
    732732/*}}}*/
    733 /*FUNCTION Numericalflux::CreateKMatrixAdjointBalancedthicknessInternal {{{1*/
    734 ElementMatrix* Numericalflux::CreateKMatrixAdjointBalancedthicknessInternal(void){
    735 
    736         ElementMatrix* Ke=CreateKMatrixBalancedthicknessInternal();
     733/*FUNCTION Numericalflux::CreateKMatrixAdjointBalancethicknessInternal {{{1*/
     734ElementMatrix* Numericalflux::CreateKMatrixAdjointBalancethicknessInternal(void){
     735
     736        ElementMatrix* Ke=CreateKMatrixBalancethicknessInternal();
    737737        if (Ke) Ke->Transpose();
    738738        return Ke;
    739739}
    740740/*}}}*/
    741 /*FUNCTION Numericalflux::CreateKMatrixAdjointBalancedthicknessBoundary {{{1*/
    742 ElementMatrix* Numericalflux::CreateKMatrixAdjointBalancedthicknessBoundary(void){
    743 
    744         ElementMatrix* Ke=CreateKMatrixBalancedthicknessBoundary();
     741/*FUNCTION Numericalflux::CreateKMatrixAdjointBalancethicknessBoundary {{{1*/
     742ElementMatrix* Numericalflux::CreateKMatrixAdjointBalancethicknessBoundary(void){
     743
     744        ElementMatrix* Ke=CreateKMatrixBalancethicknessBoundary();
    745745        if(Ke) Ke->Transpose();
    746746        return Ke;
    … …  
    839839}
    840840/*}}}*/
    841 /*FUNCTION Numericalflux::CreatePVectorBalancedthickness{{{1*/
    842 ElementVector* Numericalflux::CreatePVectorBalancedthickness(void){
     841/*FUNCTION Numericalflux::CreatePVectorBalancethickness{{{1*/
     842ElementVector* Numericalflux::CreatePVectorBalancethickness(void){
    843843
    844844        int type;
    … …  
    847847        switch(type){
    848848                case InternalEnum:
    849                         return CreatePVectorBalancedthicknessInternal();
     849                        return CreatePVectorBalancethicknessInternal();
    850850                case BoundaryEnum:
    851                         return CreatePVectorBalancedthicknessBoundary();
     851                        return CreatePVectorBalancethicknessBoundary();
    852852                default:
    853853                        _error_("type not supported yet");
    … …  
    855855}
    856856/*}}}*/
    857 /*FUNCTION Numericalflux::CreatePVectorBalancedthicknessInternal{{{1*/
    858 ElementVector* Numericalflux::CreatePVectorBalancedthicknessInternal(void){
     857/*FUNCTION Numericalflux::CreatePVectorBalancethicknessInternal{{{1*/
     858ElementVector* Numericalflux::CreatePVectorBalancethicknessInternal(void){
    859859
    860860        /*Nothing added to PVector*/
    … …  
    863863}
    864864/*}}}*/
    865 /*FUNCTION Numericalflux::CreatePVectorBalancedthicknessBoundary{{{1*/
    866 ElementVector* Numericalflux::CreatePVectorBalancedthicknessBoundary(void){
     865/*FUNCTION Numericalflux::CreatePVectorBalancethicknessBoundary{{{1*/
     866ElementVector* Numericalflux::CreatePVectorBalancethicknessBoundary(void){
    867867
    868868        /* constants*/
    … …  
    929929}
    930930/*}}}*/
    931 /*FUNCTION Numericalflux::CreatePVectorAdjointBalancedthickness{{{1*/
    932 ElementVector* Numericalflux::CreatePVectorAdjointBalancedthickness(void){
     931/*FUNCTION Numericalflux::CreatePVectorAdjointBalancethickness{{{1*/
     932ElementVector* Numericalflux::CreatePVectorAdjointBalancethickness(void){
    933933
    934934        /*No PVector for the Adjoint*/

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

    r6412 r8287  
    7878                ElementMatrix* CreateKMatrixPrognosticInternal(void);
    7979                ElementMatrix* CreateKMatrixPrognosticBoundary(void);
    80                 ElementMatrix* CreateKMatrixBalancedthickness(void);
    81                 ElementMatrix* CreateKMatrixBalancedthicknessInternal(void);
    82                 ElementMatrix* CreateKMatrixBalancedthicknessBoundary(void);
    83                 ElementMatrix* CreateKMatrixAdjointBalancedthickness(void);
    84                 ElementMatrix* CreateKMatrixAdjointBalancedthicknessInternal(void);
    85                 ElementMatrix* CreateKMatrixAdjointBalancedthicknessBoundary(void);
     80                ElementMatrix* CreateKMatrixBalancethickness(void);
     81                ElementMatrix* CreateKMatrixBalancethicknessInternal(void);
     82                ElementMatrix* CreateKMatrixBalancethicknessBoundary(void);
     83                ElementMatrix* CreateKMatrixAdjointBalancethickness(void);
     84                ElementMatrix* CreateKMatrixAdjointBalancethicknessInternal(void);
     85                ElementMatrix* CreateKMatrixAdjointBalancethicknessBoundary(void);
    8686                ElementVector* CreatePVectorPrognostic(void);
    8787                ElementVector* CreatePVectorPrognosticInternal(void);
    8888                ElementVector* CreatePVectorPrognosticBoundary(void);
    89                 ElementVector* CreatePVectorBalancedthickness(void);
    90                 ElementVector* CreatePVectorBalancedthicknessInternal(void);
    91                 ElementVector* CreatePVectorBalancedthicknessBoundary(void);
    92                 ElementVector* CreatePVectorAdjointBalancedthickness(void);
     89                ElementVector* CreatePVectorBalancethickness(void);
     90                ElementVector* CreatePVectorBalancethicknessInternal(void);
     91                ElementVector* CreatePVectorBalancethicknessBoundary(void);
     92                ElementVector* CreatePVectorAdjointBalancethickness(void);
    9393                /*}}}*/
    9494

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

    r8260 r8287  
    108108        }
    109109
    110         /*Prognostic/ Melting/ Slopecompute/ Balancedvelocities/ Balancedthickness*/
     110        /*Prognostic/ Melting/ Slopecompute/ Balancevelocities/ Balancethickness*/
    111111        if (
    112112                                analysis_type==PrognosticAnalysisEnum ||
    … …  
    114114                                analysis_type==BedSlopeAnalysisEnum ||
    115115                                analysis_type==SurfaceSlopeAnalysisEnum ||
    116                                 analysis_type==BalancedvelocitiesAnalysisEnum ||
    117                                 analysis_type==BalancedthicknessAnalysisEnum
     116                                analysis_type==BalancevelocitiesAnalysisEnum ||
     117                                analysis_type==BalancethicknessAnalysisEnum
    118118                                ){
    119119                if (iomodel->dim==3){

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

    r8224 r8287  
    3131                        adjointcore=&adjointdiagnostic_core;
    3232                        break;
    33                 case BalancedthicknessSolutionEnum:
     33                case BalancethicknessSolutionEnum:
    3434                        adjointcore=&adjointbalancedthickness_core;
    3535                        break;

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

    r8224 r8287  
    3737                        solutioncore=&prognostic_core;
    3838                        break;
    39                 case BalancedthicknessSolutionEnum:
     39                case BalancethicknessSolutionEnum:
    4040                        solutioncore=&balancedthickness_core;
    4141                        break;
    42                 case BalancedvelocitiesSolutionEnum:
     42                case BalancevelocitiesSolutionEnum:
    4343                        solutioncore=&balancedvelocities_core;
    4444                        break;

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

    r8224 r8287  
    7575                        break;
    7676
    77                 case BalancedthicknessSolutionEnum:
     77                case BalancethicknessSolutionEnum:
    7878                        numanalyses=1;
    7979                        analyses=(int*)xmalloc(numanalyses*sizeof(int));
    80                         analyses[0]=BalancedthicknessAnalysisEnum;
     80                        analyses[0]=BalancethicknessAnalysisEnum;
    8181                        break;
    8282
    83                 case BalancedvelocitiesSolutionEnum:
     83                case BalancevelocitiesSolutionEnum:
    8484                        numanalyses=1;
    8585                        analyses=(int*)xmalloc(numanalyses*sizeof(int));
    86                         analyses[0]=BalancedvelocitiesAnalysisEnum;
     86                        analyses[0]=BalancevelocitiesAnalysisEnum;
    8787                        break;
    8888

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

    r6412 r8287  
    2424        /*compute thickness */
    2525        _printf_(VerboseSolution(),"%s\n","   computing thickness");
    26         femmodel->SetCurrentConfiguration(BalancedthicknessAnalysisEnum);
     26        femmodel->SetCurrentConfiguration(BalancethicknessAnalysisEnum);
    2727        solver_linear(femmodel);
    2828
    … …  
    3232        /*compute adjoint*/
    3333        _printf_(VerboseSolution(),"%s\n","   computing adjoint");
    34         femmodel->SetCurrentConfiguration(BalancedthicknessAnalysisEnum,AdjointBalancedthicknessAnalysisEnum);
     34        femmodel->SetCurrentConfiguration(BalancethicknessAnalysisEnum,AdjointBalancethicknessAnalysisEnum);
    3535        solver_adjoint_linear(femmodel);
    3636       

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

    r8285 r8287  
    2020
    2121        /*activate formulation: */
    22         femmodel->SetCurrentConfiguration(BalancedthicknessAnalysisEnum);
     22        femmodel->SetCurrentConfiguration(BalancethicknessAnalysisEnum);
    2323       
    2424        /*recover parameters: */
    … …  
    3030        solver_linear(femmodel);
    3131
    32         if(solution_type==BalancedthicknessSolutionEnum && !control_analysis){
     32        if(solution_type==BalancethicknessSolutionEnum && !control_analysis){
    3333                _printf_(VerboseSolution(),"   saving results\n");
    3434                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,ThicknessEnum);

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

    r8285 r8287  
    1919
    2020        /*activate formulation: */
    21         femmodel->SetCurrentConfiguration(BalancedvelocitiesAnalysisEnum);
     21        femmodel->SetCurrentConfiguration(BalancevelocitiesAnalysisEnum);
    2222       
    2323        /*recover parameters: */
    … …  
    2929        solver_linear(femmodel);
    3030
    31         if(solution_type==BalancedvelocitiesSolutionEnum && !control_analysis){
     31        if(solution_type==BalancevelocitiesSolutionEnum && !control_analysis){
    3232                _printf_(VerboseSolution(),"   saving results\n");
    3333                InputToResultx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,VelEnum);

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

    r8224 r8287  
    5151                femmodel->SetCurrentConfiguration(DiagnosticHorizAnalysisEnum);
    5252        }
    53         else if (solution_type==BalancedthicknessSolutionEnum){
    54                 femmodel->SetCurrentConfiguration(BalancedthicknessAnalysisEnum);
     53        else if (solution_type==BalancethicknessSolutionEnum){
     54                femmodel->SetCurrentConfiguration(BalancethicknessAnalysisEnum);
    5555        }
    5656        else{
    … …  
    6868                solver_nonlinear(femmodel,conserve_loads);
    6969        }
    70         else if (solution_type==BalancedthicknessSolutionEnum){
     70        else if (solution_type==BalancethicknessSolutionEnum){
    7171                solver_linear(femmodel);
    7272        }

  • issm/trunk/src/m/enum/EnumToModelField.m

    r8248 r8287  
    1111switch enum,
    1212
    13                 case ThicknessEnum(), string='thickness'; return
    14                 case DragCoefficientEnum(), string='drag_coefficient'; return
    15                 case RheologyBEnum(), string='rheology_B'; return
    16                 case RheologyBbarEnum(), string='rheology_B'; return
    17                 case DhDtEnum(), string='dhdt'; return
    18                 case VxEnum(), string='vx'; return
    19                 case VxObsEnum(), string='vx_obs'; return
    20                 case VyEnum(), string='vy'; return
    21                 case VyObsEnum(), string='vy_obs'; return
    22                 case GroundingLineMigrationEnum(), string='gl_migration'; return
    23                 otherwise, error(['Enum ' num2str(enum)  ' not found associated to any model field']);
    24 
    25 end

  • issm/trunk/src/m/model/ismodelselfconsistent.m

    r8249 r8287  
    199199
    200200        %OBSERVED VELOCITIES
    201         if md.solution_type==BalancedthicknessSolutionEnum
     201        if md.solution_type==BalancethicknessSolutionEnum
    202202                fields={'thickness_obs'};
    203203                checksize(md,fields,[md.numberofgrids 1]);
    … …  
    484484                        % No checks for now
    485485                        %}}}
    486                 case BalancedthicknessAnalysisEnum,
     486                case BalancethicknessAnalysisEnum,
    487487                        % {{{2
    488488                        %VELOCITIES MELTING AND ACCUMULATION
    … …  
    503503                        end
    504504                        %}}}
    505                 case BalancedvelocitiesAnalysisEnum,
     505                case BalancevelocitiesAnalysisEnum,
    506506                        % {{{2
    507507                        %VELOCITIES MELTING AND ACCUMULATION

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

    r5377 r8287  
    99        case DiagnosticSolutionEnum,         adjointcore='adjointdiagnostic_core';
    1010        case SteadystateSolutionEnum,        adjointcore='adjointdiagnostic_core';
    11         case BalancedthicknessSolutionEnum,  adjointcore='adjointbalancedthickness_core';
     11        case BalancethicknessSolutionEnum,   adjointcore='adjointbalancethickness_core';
    1212        otherwise error('%s%s%s',' solution type: ',EnumToString(solutiontype),' not supported yet!');
    1313

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

    r7687 r8287  
    1111        case ThermalSolutionEnum,                  solutioncore='thermal_core';
    1212        case PrognosticSolutionEnum,               solutioncore='prognostic_core';
    13         case BalancedthicknessSolutionEnum,        solutioncore='balancedthickness_core';
    14         case BalancedvelocitiesSolutionEnum,       solutioncore='balancedvelocities_core';
     13        case BalancethicknessSolutionEnum,         solutioncore='balancethickness_core';
     14        case BalancevelocitiesSolutionEnum,        solutioncore='balancevelocities_core';
    1515        case SurfaceSlopeSolutionEnum,             solutioncore='surfaceslope_core';
    1616        case BedSlopeSolutionEnum,                 solutioncore='bedslope_core';

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

    r7687 r8287  
    2626                analyses=[PrognosticAnalysisEnum];
    2727
    28         case BalancedthicknessSolutionEnum,
     28        case BalancethicknessSolutionEnum,
    2929                numanalyses=1;
    30                 analyses=[BalancedthicknessAnalysisEnum];
     30                analyses=[BalancethicknessAnalysisEnum];
    3131
    32         case BalancedvelocitiesSolutionEnum,
     32        case BalancevelocitiesSolutionEnum,
    3333                numanalyses=1;
    34                 analyses=[BalancedvelocitiesAnalysisEnum];
     34                analyses=[BalancevelocitiesAnalysisEnum];
    3535
    3636        case SurfaceSlopeSolutionEnum,

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

    r6323 r8287  
    1 function femmodel=adjointbalancedthickness_core(femmodel),
     1function femmodel=adjointbalancethickness_core(femmodel),
    22%ADJOINT_CORE - compute inverse method adjoint state
    33%
    44%   Usage:
    5 %      femmodel=adjointbalancedthickness_core(femmodel)
     5%      femmodel=adjointbalancethickness_core(femmodel)
    66
    77        %recover parameters common to all solutions
    … …  
    1111        %set analysis type to compute velocity:
    1212        issmprintf(VerboseSolution,'   computing thickness');
    13         femmodel=SetCurrentConfiguration(femmodel,BalancedthicknessAnalysisEnum);
     13        femmodel=SetCurrentConfiguration(femmodel,BalancethicknessAnalysisEnum);
    1414        femmodel=solver_linear(femmodel);
    1515
    … …  
    1919        %compute adjoint
    2020        issmprintf(VerboseSolution,'   computing adjoint');
    21         femmodel=SetCurrentConfiguration(femmodel,BalancedthicknessAnalysisEnum,AdjointBalancedthicknessAnalysisEnum);
     21        femmodel=SetCurrentConfiguration(femmodel,BalancethicknessAnalysisEnum,AdjointBalancethicknessAnalysisEnum);
    2222        femmodel=solver_adjoint_linear(femmodel);
    2323

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

    r6323 r8287  
    1 function femmodel=balancedthickness_core(femmodel)
     1function femmodel=balancethickness_core(femmodel)
    22%BALANCEDTHICKNESS_CORE - linear solution sequence
    33%
    44%   Usage:
    5 %      femmodel=balancedthickness_core(femmode)
     5%      femmodel=balancethickness_core(femmode)
    66
    77        %recover parameters common to all solutions
    … …  
    1111
    1212        %Activate formulation
    13         femmodel=SetCurrentConfiguration(femmodel,BalancedthicknessAnalysisEnum);
     13        femmodel=SetCurrentConfiguration(femmodel,BalancethicknessAnalysisEnum);
    1414
    1515        issmprintf(VerboseSolution,'   call computational core');
    1616        femmodel=solver_linear(femmodel);
    1717
    18         if (solution_type==BalancedthicknessSolutionEnum & ~control_analysis),
     18        if (solution_type==BalancethicknessSolutionEnum & ~control_analysis),
    1919                issmprintf(VerboseSolution,'   saving results');
    2020                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);

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

    r6323 r8287  
    1 function femmodel=balancedvelocities_core(femmdoel)
     1function femmodel=balancevelocities_core(femmdoel)
    22%BALANCEDVELOCITIES_CORE - linear solution sequence
    33%
    44%   Usage:
    5 %      femmodel=balancedvelocities_core(femmodel)
     5%      femmodel=balancevelocities_core(femmodel)
    66
    77        %recover parameters common to all solutions
    … …  
    1111
    1212        %Activate formulation
    13         femmodel=SetCurrentConfiguration(femmodel,BalancedvelocitiesAnalysisEnum);
     13        femmodel=SetCurrentConfiguration(femmodel,BalancevelocitiesAnalysisEnum);
    1414
    1515        issmprintf(VerboseSolution,'   call computational core');
    1616        femmodel=solver_linear(femmodel);
    1717       
    18         if (solution_type==BalancedvelocitiesSolutionEnum & ~control_analysis),
     18        if (solution_type==BalancevelocitiesSolutionEnum & ~control_analysis),
    1919                issmprintf(VerboseSolution,'   saving results');
    2020                femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VelEnum);

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

    r6226 r8287  
    1212if (solution_type==DiagnosticSolutionEnum | solution_type==SteadystateSolutionEnum)
    1313        femmodel=SetCurrentConfiguration(femmodel,DiagnosticHorizAnalysisEnum);
    14 elseif (solution_type==BalancedthicknessSolutionEnum)
    15         femmodel=SetCurrentConfiguration(femmodel,BalancedthicknessAnalysisEnum);
     14elseif (solution_type==BalancethicknessSolutionEnum)
     15        femmodel=SetCurrentConfiguration(femmodel,BalancethicknessAnalysisEnum);
    1616else
    1717        error('solution not implemented yet');
    … …  
    2626elseif (solution_type==SteadystateSolutionEnum)
    2727        femmodel=diagnostic_core(femmodel);  %We need a 3D velocity!! (vz is required for the next thermal run)
    28 elseif (solution_type==BalancedthicknessSolutionEnum)
     28elseif (solution_type==BalancethicknessSolutionEnum)
    2929        femmodel=solver_linear(femmodel);
    3030else
Note: See TracChangeset for help on using the changeset viewer.