Changeset 16218


Ignore:
Timestamp:
09/20/13 15:44:30 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: cleaning up updateelements

Location:
issm/trunk-jpl/src
Files:
19 edited

Legend:

Unmodified
Added
Removed

  • TabularUnified issm/trunk-jpl/src/c/analyses/stressbalance_core.cpp

    r16158 r16218  
    4242                InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuVyEnum,VyEnum);
    4343                InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuVzEnum,VzEnum);
    44                 InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuPressureEnum,PressureEnum);
     44                if(isFS) InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuPressureEnum,PressureEnum);
    4545        }
    4646

  • TabularUnified issm/trunk-jpl/src/c/analyses/transient_core.cpp

    r16215 r16218  
    2121        /*parameters: */
    2222        IssmDouble starttime,finaltime,dt,yts;
    23         bool   isstressbalance,ismasstransport,isthermal,isgroundingline,isenthalpy,isdelta18o,isgia;
     23        bool   isstressbalance,ismasstransport,isFS,isthermal,isgroundingline,isenthalpy,isdelta18o,isgia;
    2424        bool   save_results,dakota_analysis;
    2525        bool   time_adapt=false;
    … …  
    4848        femmodel->parameters->FindParam(&isgroundingline,TransientIsgroundinglineEnum);
    4949        femmodel->parameters->FindParam(&isenthalpy,ThermalIsenthalpyEnum);
     50        femmodel->parameters->FindParam(&isFS,FlowequationIsFSEnum);
    5051        if(isgroundingline) femmodel->parameters->FindParam(&groundingline_migration,GroundinglineMigrationEnum);
    5152        femmodel->parameters->FindParam(&numoutputs,TransientNumRequestedOutputsEnum);
    … …  
    6263                        InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuVxEnum,VxEnum);
    6364                        InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuVyEnum,VyEnum);
    64                         InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuVzEnum,VzEnum);
    65                         InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuPressureEnum,PressureEnum);
     65                        if(dim==3){
     66                                InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuVzEnum,VzEnum);
     67                                if(isFS)InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuPressureEnum,PressureEnum);
     68                        }
    6669                }
    67                 if(ismasstransport | isgroundingline){
     70                if(ismasstransport || isgroundingline){
    6871                        InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuThicknessEnum,ThicknessEnum);
    6972                        InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,QmuSurfaceEnum,SurfaceEnum);

  • TabularUnified issm/trunk-jpl/src/c/classes/Elements/Penta.cpp

    r16215 r16218  
    34983498
    34993499        /*Checks if debuging*/
    3500         /*{{{*/
    35013500        _assert_(iomodel->elements);
    3502         /*}}}*/
    35033501
    35043502        /*Recover element type*/
    … …  
    36563654                case StressbalanceAnalysisEnum:
    36573655
    3658                         /*default vx,vy and vz: either observation or 0 */
    3659                         if(!iomodel->Data(VxEnum)){
    3660                                 for(i=0;i<6;i++)nodeinputs[i]=0;
    3661                                 this->inputs->AddInput(new PentaInput(VxEnum,nodeinputs,P1Enum));
    3662                                 if(dakota_analysis) this->inputs->AddInput(new PentaInput(QmuVxEnum,nodeinputs,P1Enum));
    3663                         }
    3664                         if(!iomodel->Data(VyEnum)){
    3665                                 for(i=0;i<6;i++)nodeinputs[i]=0;
    3666                                 this->inputs->AddInput(new PentaInput(VyEnum,nodeinputs,P1Enum));
    3667                                 if(dakota_analysis) this->inputs->AddInput(new PentaInput(QmuVyEnum,nodeinputs,P1Enum));
    3668                         }
    3669                         if(!iomodel->Data(VzEnum)){
    3670                                 for(i=0;i<6;i++)nodeinputs[i]=0;
    3671                                 this->inputs->AddInput(new PentaInput(VzEnum,nodeinputs,P1Enum));
    3672                                 if(dakota_analysis) this->inputs->AddInput(new PentaInput(QmuVzEnum,nodeinputs,P1Enum));
    3673                         }
    3674                         if(!iomodel->Data(PressureEnum)){
    3675                                 for(i=0;i<6;i++)nodeinputs[i]=0;
    3676                                 if(dakota_analysis){
    3677                                         this->inputs->AddInput(new PentaInput(PressureEnum,nodeinputs,P1Enum));
    3678                                         this->inputs->AddInput(new PentaInput(QmuPressureEnum,nodeinputs,P1Enum));
    3679                                 }
    3680                                 if(isFS){
    3681                                         this->inputs->AddInput(new PentaInput(PressureEnum,nodeinputs,P1Enum));
    3682                                         this->inputs->AddInput(new PentaInput(PressurePicardEnum,nodeinputs,P1Enum));
    3683                                 }
    3684                         }
    36853656                        if(*(iomodel->Data(FlowequationElementEquationEnum)+index)==HOFSApproximationEnum){
    36863657                                /*Create VzHO and VzFS Enums*/
    … …  
    37133684                        break;
    37143685
    3715                 case ThermalAnalysisEnum:
    3716                         /*Initialize mesh velocity*/
    3717                         for(i=0;i<6;i++)nodeinputs[i]=0;
    3718                         this->inputs->AddInput(new PentaInput(VxMeshEnum,nodeinputs,P1Enum));
    3719                         this->inputs->AddInput(new PentaInput(VyMeshEnum,nodeinputs,P1Enum));
    3720                         this->inputs->AddInput(new PentaInput(VzMeshEnum,nodeinputs,P1Enum));
    3721                         if(dakota_analysis){
    3722                                 this->inputs->AddInput(new PentaInput(QmuVxMeshEnum,nodeinputs,P1Enum));
    3723                                 this->inputs->AddInput(new PentaInput(QmuVyMeshEnum,nodeinputs,P1Enum));
    3724                                 this->inputs->AddInput(new PentaInput(QmuVzMeshEnum,nodeinputs,P1Enum));
    3725                         }
    3726                         break;
    37273686
    37283687                case EnthalpyAnalysisEnum:
    37293688                        /*Initialize mesh velocity*/
    37303689                        for(i=0;i<6;i++)nodeinputs[i]=0;
    3731                         this->inputs->AddInput(new PentaInput(VxMeshEnum,nodeinputs,P1Enum));
    3732                         this->inputs->AddInput(new PentaInput(VyMeshEnum,nodeinputs,P1Enum));
    3733                         this->inputs->AddInput(new PentaInput(VzMeshEnum,nodeinputs,P1Enum));
    37343690                        if (iomodel->Data(TemperatureEnum) && iomodel->Data(WaterfractionEnum) && iomodel->Data(PressureEnum)) {
    37353691                                for(i=0;i<6;i++){

  • TabularUnified issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

    r16215 r16218  
    24932493        this->InputUpdateFromIoModel(index,iomodel);
    24942494
    2495         /*Defaults if not provided in iomodel*/
    2496         switch(analysis_type){
    2497 
    2498                 case StressbalanceAnalysisEnum:
    2499 
    2500                         /*default vx,vy and vz: either observation or 0 */
    2501                         if(!iomodel->Data(VxEnum)){
    2502                                 for(i=0;i<3;i++)nodeinputs[i]=0;
    2503                                 this->inputs->AddInput(new TriaInput(VxEnum,nodeinputs,P1Enum));
    2504                                 if(dakota_analysis)this->inputs->AddInput(new TriaInput(QmuVxEnum,nodeinputs,P1Enum));
    2505                         }
    2506                         if(!iomodel->Data(VyEnum)){
    2507                                 for(i=0;i<3;i++)nodeinputs[i]=0;
    2508                                 this->inputs->AddInput(new TriaInput(VyEnum,nodeinputs,P1Enum));
    2509                                 if(dakota_analysis)this->inputs->AddInput(new TriaInput(QmuVyEnum,nodeinputs,P1Enum));
    2510                         }
    2511                         if(!iomodel->Data(VzEnum)){
    2512                                 for(i=0;i<3;i++)nodeinputs[i]=0;
    2513                                 this->inputs->AddInput(new TriaInput(VzEnum,nodeinputs,P1Enum));
    2514                                 if(dakota_analysis)this->inputs->AddInput(new TriaInput(QmuVzEnum,nodeinputs,P1Enum));
    2515                         }
    2516                         if(!iomodel->Data(PressureEnum)){
    2517                                 for(i=0;i<3;i++)nodeinputs[i]=0;
    2518                                 if(dakota_analysis){
    2519                                         this->inputs->AddInput(new TriaInput(PressureEnum,nodeinputs,P1Enum));
    2520                                         this->inputs->AddInput(new TriaInput(QmuPressureEnum,nodeinputs,P1Enum));
    2521                                 }
    2522                         }
    2523                         break;
    2524 
    2525                 default:
    2526                         /*No update for other solution types*/
    2527                         break;
    2528 
    2529         }
    2530 
    2531         //this->parameters: we still can't point to it, it may not even exist. Configure will handle this.
    2532         this->parameters=NULL;
    25332495}
    25342496/*}}}*/
    … …  
    30793041        Input* vx_input=inputs->GetInput(VxEnum);                 _assert_(vx_input);
    30803042        Input* vy_input=inputs->GetInput(VyEnum);                 _assert_(vy_input);
    3081         Input* vz_input=inputs->GetInput(VzEnum);                 _assert_(vz_input);
    30823043        Input* gllevelset_input=NULL;
    30833044        parameters->FindParam(&analysis_type,AnalysisTypeEnum);
    … …  
    55515512
    55525513        /*Skip if water or ice shelf element*/
    5553         if(NoIceInElement() | IsFloating()) return NULL;
     5514        if(NoIceInElement() || IsFloating()) return NULL;
    55545515
    55555516        /*Fetch number of nodes and dof for this finite element*/
    … …  
    57765737
    57775738        /*Skip if water or ice shelf element*/
    5778         if(NoIceInElement() | IsFloating()) return NULL;
     5739        if(NoIceInElement() || IsFloating()) return NULL;
    57795740
    57805741        /*Fetch number of nodes and dof for this finite element*/

  • TabularUnified issm/trunk-jpl/src/c/modules/InputUpdateFromConstantx/InputUpdateFromConstantx.cpp

    r15849 r16218  
    7171
    7272}
     73void InputUpdateFromConstantx(Elements* elements,IssmDouble constant, int name){
     74
     75        int i;
     76        if(VerboseModule()) _printf0_("   Input updates from constant\n");
     77
     78        /*Elements and loads drive the update: */
     79        for(i=0;i<elements->Size();i++){
     80                Element* element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
     81                element->InputUpdateFromConstant(constant,name);
     82        }
     83}

  • TabularUnified issm/trunk-jpl/src/c/modules/InputUpdateFromConstantx/InputUpdateFromConstantx.h

    r15849 r16218  
    99
    1010/* local prototypes: */
    11 void InputUpdateFromConstantx(FemModel* femmode,bool   constant    ,int name);
    12 void InputUpdateFromConstantx(FemModel* femmode,int    constant    ,int name);
    13 void InputUpdateFromConstantx(FemModel* femmode,IssmDouble constant,int name);
     11void InputUpdateFromConstantx(FemModel* femmodel,bool       constant,int name);
     12void InputUpdateFromConstantx(FemModel* femmodel,int        constant,int name);
     13void InputUpdateFromConstantx(FemModel* femmodel,IssmDouble constant,int name);
     14void InputUpdateFromConstantx(Elements* elements,IssmDouble constant,int name);
    1415
    1516#endif  /* _UPDATEINPUTSFROMCONSTANTXX_H */

  • TabularUnified issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateFaces.cpp

    r15611 r16218  
    5656
    5757                        /*Go through all processed faces connected to v1 and check whether we have seen this edge yet*/
    58                         _assert_(v1>=0 & v1<iomodel->numberofvertices);
     58                        _assert_(v1>=0 && v1<iomodel->numberofvertices);
    5959                        for(int e=head_minv[v1]; e!=-1; e=next_edge[e]){
    6060                                if(facestemp[e*4+1]==v2){

  • TabularUnified issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/UpdateElementsEnthalpy.cpp

    r16026 r16218  
    77#include "../../../shared/shared.h"
    88#include "../../MeshPartitionx/MeshPartitionx.h"
     9#include "../../InputUpdateFromConstantx/InputUpdateFromConstantx.h"
    910#include "../ModelProcessorx.h"
    1011
    … …  
    2627                }
    2728        }
     29
     30        bool dakota_analysis;
     31        iomodel->Constant(&dakota_analysis,QmuIsdakotaEnum);
    2832
    2933        iomodel->FetchDataToInput(elements,ThicknessEnum);
    … …  
    4953        iomodel->FetchDataToInput(elements,VyEnum);
    5054        iomodel->FetchDataToInput(elements,VzEnum);
     55        InputUpdateFromConstantx(elements,0.,VxMeshEnum);
     56        InputUpdateFromConstantx(elements,0.,VyMeshEnum);
     57        InputUpdateFromConstantx(elements,0.,VzMeshEnum);
     58        if(dakota_analysis){
     59                elements->InputDuplicate(TemperatureEnum,QmuTemperatureEnum);
     60                elements->InputDuplicate(BasalforcingsMeltingRateEnum,QmuMeltingEnum);
     61                elements->InputDuplicate(VxMeshEnum,QmuVxMeshEnum);
     62                elements->InputDuplicate(VxMeshEnum,QmuVyMeshEnum);
     63                elements->InputDuplicate(VxMeshEnum,QmuVzMeshEnum);
     64        }
    5165
    5266        /*Free data: */

  • TabularUnified issm/trunk-jpl/src/c/modules/ModelProcessorx/FreeSurfaceBase/UpdateElementsFreeSurfaceBase.cpp

    r16211 r16218  
    2929        iomodel->FetchDataToInput(elements,SurfaceEnum);
    3030        iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum);
    31         iomodel->FetchDataToInput(elements,SurfaceforcingsMassBalanceEnum);
    3231        iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum);
    3332        iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateCorrectionEnum,0.);

  • TabularUnified issm/trunk-jpl/src/c/modules/ModelProcessorx/FreeSurfaceTop/UpdateElementsFreeSurfaceTop.cpp

    r16211 r16218  
    2828        iomodel->FetchDataToInput(elements,SurfaceEnum);
    2929        iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum);
    30         iomodel->FetchDataToInput(elements,SurfaceforcingsMassBalanceEnum);
    31         iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum);
    32         iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateCorrectionEnum,0.);
     30        iomodel->FetchDataToInput(elements,SurfaceforcingsMassBalanceEnum,0.);
    3331        iomodel->FetchDataToInput(elements,VxEnum);
    3432        iomodel->FetchDataToInput(elements,VyEnum);

  • TabularUnified issm/trunk-jpl/src/c/modules/ModelProcessorx/Masstransport/UpdateElementsMasstransport.cpp

    r16215 r16218  
    6767                iomodel->FetchDataToInput(elements,MeshElementonbedEnum);
    6868                iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum);
    69                 iomodel->FetchDataToInput(elements,TemperatureEnum);
    70                 iomodel->FetchDataToInput(elements,VzEnum);
    7169        }
     70
    7271        if(issmbgradients){
    7372                iomodel->FetchDataToInput(elements,SurfaceforcingsHrefEnum);
    … …  
    8887                }
    8988        }
    90         if(~ispdd & ~issmbgradients){
    91                 iomodel->FetchDataToInput(elements,SurfaceforcingsMassBalanceEnum);
     89        if(~ispdd && ~issmbgradients){
     90                iomodel->FetchDataToInput(elements,SurfaceforcingsMassBalanceEnum,0.);
    9291        }
    9392}

  • TabularUnified issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateConstraintsStressbalance.cpp

    r16166 r16218  
    5555
    5656        /*Now, is the flag macayaealHO on? otherwise, do nothing: */
    57         if(!isSSA & !isHO & !isFS & !isL1L2){
     57        if(!isSSA && !isHO && !isFS && !isL1L2){
    5858                *pconstraints=constraints;
    5959                return;

  • TabularUnified issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateLoadsStressbalance.cpp

    r16200 r16218  
    3030
    3131        /*Now, is the flag macayaealHO on? otherwise, do nothing: */
    32         if(!isSSA & !isHO & !isFS & !isL1L2) return;
     32        if(!isSSA && !isHO && !isFS && !isL1L2) return;
    3333
    3434        /*Initialize counter: */

  • TabularUnified issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp

    r16215 r16218  
    9797        iomodel->FetchDataToInput(elements,MaterialsRheologyNEnum);
    9898        iomodel->FetchDataToInput(elements,VxEnum);
     99        if(dakota_analysis)elements->InputDuplicate(VxEnum,QmuVxEnum);
    99100        iomodel->FetchDataToInput(elements,VyEnum);
     101        if(dakota_analysis)elements->InputDuplicate(VyEnum,QmuVyEnum);
    100102        iomodel->FetchDataToInput(elements,LoadingforceXEnum);
    101103        iomodel->FetchDataToInput(elements,LoadingforceYEnum);
    102104        iomodel->FetchDataToInput(elements,DamageDEnum);
    103         iomodel->FetchDataToInput(elements,PressureEnum);
    104         if(dakota_analysis)elements->InputDuplicate(PressureEnum,QmuPressureEnum);
    105105
    106106        if(iomodel->dim==3){
    107107                iomodel->FetchDataToInput(elements,MeshElementonbedEnum);
    108108                iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum);
    109                 iomodel->FetchDataToInput(elements,VzEnum);
    110                 iomodel->FetchDataToInput(elements,SurfaceforcingsMassBalanceEnum);
    111109                iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum);
    112110                iomodel->FetchDataToInput(elements,FlowequationBorderFSEnum);
    113111                iomodel->FetchDataToInput(elements,LoadingforceZEnum);
     112                iomodel->FetchDataToInput(elements,VzEnum);
     113                if(dakota_analysis)elements->InputDuplicate(VzEnum,QmuVzEnum);
     114                if(isFS){
     115                        iomodel->FetchDataToInput(elements,PressureEnum);
     116                        if(dakota_analysis)elements->InputDuplicate(PressureEnum,QmuPressureEnum);
     117                }
    114118        }
    115 
    116         if(control_analysis){
    117                 iomodel->FetchDataToInput(elements,InversionVxObsEnum);
    118                 iomodel->FetchDataToInput(elements,InversionVyObsEnum);
    119         }
    120 
    121         /*Post-process: */
    122         elements->InputDuplicate(VxEnum,VxPicardEnum);
    123         elements->InputDuplicate(VxEnum,InversionVxObsEnum);
    124         if(dakota_analysis)elements->InputDuplicate(VxEnum,QmuVxEnum);
    125119
    126120        #ifdef _HAVE_ANDROID_
    … …  
    128122        #endif
    129123
    130         elements->InputDuplicate(VyEnum,VyPicardEnum);
    131         elements->InputDuplicate(VyEnum,InversionVyObsEnum);
    132         if(dakota_analysis)elements->InputDuplicate(VyEnum,QmuVyEnum);
    133 
    134         if(iomodel->dim==3){
    135                 elements->InputDuplicate(VzEnum,VzPicardEnum);
    136                 elements->InputDuplicate(VzEnum,InversionVzObsEnum);
    137                 if(dakota_analysis)elements->InputDuplicate(VzEnum,QmuVzEnum);
    138         }
    139 
    140124        /*Free data: */
    141125        iomodel->DeleteData(1,FlowequationElementEquationEnum);

  • TabularUnified issm/trunk-jpl/src/c/modules/ModelProcessorx/StressbalanceVertical/UpdateElementsStressbalanceVertical.cpp

    r16210 r16218  
    3131        iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum);
    3232        iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum);
    33         iomodel->FetchDataToInput(elements,SurfaceforcingsMassBalanceEnum);
    3433        iomodel->FetchDataToInput(elements,VxEnum);
    3534        iomodel->FetchDataToInput(elements,VyEnum);

  • TabularUnified issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/UpdateElementsThermal.cpp

    r16215 r16218  
    77#include "../../../shared/shared.h"
    88#include "../../MeshPartitionx/MeshPartitionx.h"
     9#include "../../InputUpdateFromConstantx/InputUpdateFromConstantx.h"
    910#include "../ModelProcessorx.h"
    1011
    … …  
    2324                }
    2425        }
     26
     27        bool dakota_analysis;
     28        iomodel->Constant(&dakota_analysis,QmuIsdakotaEnum);
    2529
    2630        iomodel->FetchDataToInput(elements,ThicknessEnum);
    … …  
    4347        iomodel->FetchDataToInput(elements,VyEnum);
    4448        iomodel->FetchDataToInput(elements,VzEnum);
    45 
    46         bool dakota_analysis;
    47         iomodel->Constant(&dakota_analysis,QmuIsdakotaEnum);
     49        InputUpdateFromConstantx(elements,0.,VxMeshEnum);
     50        InputUpdateFromConstantx(elements,0.,VyMeshEnum);
     51        InputUpdateFromConstantx(elements,0.,VzMeshEnum);
    4852        if(dakota_analysis){
    4953                elements->InputDuplicate(TemperatureEnum,QmuTemperatureEnum);
    5054                elements->InputDuplicate(BasalforcingsMeltingRateEnum,QmuMeltingEnum);
     55                elements->InputDuplicate(VxMeshEnum,QmuVxMeshEnum);
     56                elements->InputDuplicate(VxMeshEnum,QmuVyMeshEnum);
     57                elements->InputDuplicate(VxMeshEnum,QmuVzMeshEnum);
    5158        }
    5259}

  • TabularUnified issm/trunk-jpl/src/c/modules/Scotchx/Scotchx.cpp

    r11401 r16218  
    260260    }
    261261  }
    262   if ((flagval & C_FLAGPART) != 0) {              /* If program run as the partitioner            */
     262  if ((flagval && C_FLAGPART) != 0) {              /* If program run as the partitioner            */
    263263    C_fileTab[3].name = C_fileTab[2].name;        /* Put provided file names at their right place */
    264264    C_fileTab[2].name = C_fileTab[1].name;
    … …  
    305305  clockStop (&runtime[0]);                        /* Get output time */
    306306
    307   if (flagval & C_FLAGVERBSTR) {
     307  if (flagval && C_FLAGVERBSTR) {
    308308    fprintf (C_filepntrlogout, "S\tStrat=");
    309309    SCOTCH_stratSave (&stradat, C_filepntrlogout);
    310310    putc ('\n', C_filepntrlogout);
    311311  }
    312   if (flagval & C_FLAGVERBTIM) {
     312  if (flagval && C_FLAGVERBTIM) {
    313313    fprintf (C_filepntrlogout, "T\tMapping\t\t%g\nT\tI/O\t\t%g\nT\tTotal\t\t%g\n",
    314314             (double) clockVal (&runtime[1]),
    … …  
    317317             (double) clockVal (&runtime[1]));
    318318  }
    319   if (flagval & C_FLAGVERBMAP)
     319  if (flagval && C_FLAGVERBMAP)
    320320    SCOTCH_graphMapView (&grafdat, &mapdat, C_filepntrlogout);
    321321

  • TabularUnified issm/trunk-jpl/src/c/shared/TriMesh/TriMeshUtils.cpp

    r14917 r16218  
    109109        int i;
    110110        for (i=0;i<nriftsegs;i++){
    111                 if ((*(riftsegments+4*i+0)==el) | (*(riftsegments+4*i+1)==el)){
     111                if ((*(riftsegments+4*i+0)==el) || (*(riftsegments+4*i+1)==el)){
    112112                        return 1;
    113113                }

  • TabularUnified issm/trunk-jpl/src/m/classes/damage.m

    r16198 r16218  
    66classdef damage
    77        properties (SetAccess=public) 
    8         % {{{
    98                %damage
    109                D                   = NaN;
    … …  
    2726                c4                  = NaN;
    2827                healing             = NaN;
    29 
    30         end % }}}
     28        end
    3129        methods
    3230                function obj = damage(varargin) % {{{
Note: See TracChangeset for help on using the changeset viewer.