Changeset 16151


Ignore:
Timestamp:
09/17/13 14:11:06 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: more cppcheck errors

Location:
issm/trunk-jpl/src/c/classes
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/classes/Elements/Penta.cpp

    r16149 r16151  
    45444544
    45454545        /*Intermediaries*/
    4546         int    i,j,found=0;
    4547         int    friction_type,stabilization;
     4546        int        i,j,found=0;
     4547        int        stabilization;
    45484548        IssmDouble Jdet,phi,dt;
    45494549        IssmDouble rho_ice,heatcapacity;
     
    48134813
    48144814        /*Intermediaries*/
    4815         int    i,j,found=0;
    4816         int    friction_type,stabilization;
     4815        int        i,found=0;
     4816        int        stabilization;
    48174817        IssmDouble Jdet,phi,dt;
    48184818        IssmDouble rho_ice,heatcapacity;
     
    48444844        Input* vz_input=inputs->GetInput(VzEnum); _assert_(vz_input);
    48454845        Input* temperature_input=NULL;
    4846         if (reCast<bool,IssmDouble>(dt)) temperature_input=inputs->GetInput(TemperatureEnum); _assert_(inputs);
    4847         if (stabilization==2) diameter=MinEdgeLength(xyz_list);
     4846        if(reCast<bool,IssmDouble>(dt)) temperature_input=inputs->GetInput(TemperatureEnum); _assert_(inputs);
     4847        if(stabilization==2) diameter=MinEdgeLength(xyz_list);
    48484848
    48494849        /* Start  looping on the number of gaussian points: */
     
    55865586        IssmDouble epsilon[5]; /* epsilon=[exx,eyy,exy,exz,eyz];*/
    55875587        IssmDouble eps1[3],eps2[3];
    5588         IssmDouble phi[NUMVERTICES];
    55895588        IssmDouble dphi[3][NUMVERTICES];
    55905589        GaussPenta *gauss=NULL;
     
    61596158        /*dJ/dalpha = ∂L/∂alpha = ∂J/∂alpha + ∂/∂alpha(KU-F)*/
    61606159
    6161         int              i,approximation;
    6162         Tria*            tria=NULL;
     6160        int   approximation;
     6161        Tria* tria=NULL;
    61636162
    61646163        /*If on water, skip grad (=0): */
     
    62136212
    62146213        /*Now deal with ∂J/∂alpha*/
    6215         int        *responses = NULL;
    6216         int         num_responses,resp;
     6214        int *responses = NULL;
     6215        int num_responses,resp;
    62176216        this->parameters->FindParam(&num_responses,InversionNumCostFunctionsEnum);
    62186217        this->parameters->FindParam(&responses,NULL,NULL,StepResponsesEnum);
     
    80198018
    80208019        /*Intermediaries */
    8021         int         i,j;
    80228020        int         approximation;
    80238021        IssmDouble  xyz_list[NUMVERTICES][3];
     
    80648062
    80658063                D_scalar=2*newviscosity*gauss->weight*Jdet;
    8066                 for (i=0;i<5;i++) D[i*5+i]=D_scalar;
     8064                for(int i=0;i<5;i++) D[i*5+i]=D_scalar;
    80678065
    80688066                TripleMultiply(B,5,numdof,1,
     
    91279125
    91289126        /*Intermediaries*/
    9129         int         i,j;
    91309127        IssmDouble  Jdet;
    91319128        IssmDouble  slope[3]; //do not put 2! this goes into GetInputDerivativeValue, which addresses slope[3] also!
     
    91609157                driving_stress_baseline=matpar->GetRhoIce()*matpar->GetG();
    91619158
    9162                 for(i=0;i<numnodes;i++){
     9159                for(int i=0;i<numnodes;i++){
    91639160                        pe->values[i*NDOF2+0]+= -driving_stress_baseline*slope[0]*Jdet*gauss->weight*basis[i];
    91649161                        pe->values[i*NDOF2+1]+= -driving_stress_baseline*slope[1]*Jdet*gauss->weight*basis[i];
     
    93729369        /*Intermediaries*/
    93739370        int        i,j;
    9374         int        approximation;
    93759371        IssmDouble Jdet,gravity,rho_ice,B,D_scalar_stab,viscosity;
    93769372        IssmDouble forcex,forcey,forcez,diameter,FSreconditioning;
     
    93889384        IssmDouble dnodalbasis[6][6][3];
    93899385        IssmDouble SW[6][4][4];
    9390         int p,q,ii;
     9386        int p,ii;
    93919387        int c=3; //index of pressure
    93929388
     
    94639459
    94649460        /*Intermediaries*/
    9465         int        i,j;
     9461        int        i;
    94669462        int        approximation;
    94679463        IssmDouble Jdet,gravity,rho_ice;
     
    98669862
    98679863        /*Intermediaries */
    9868         int        i,j,approximation;
     9864        int        i,j;
    98699865        IssmDouble xyz_list[NUMVERTICES][3];
    98709866        IssmDouble Jdet;
  • issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

    r16149 r16151  
    18601860
    18611861        /*Intermediaries*/
    1862         int    i,j,t;
     1862        int    i,t;
    18631863        int    tria_vertex_ids[3];
    18641864        int    row;
     
    24182418
    24192419        /*Intermediaries*/
    2420         int        i,j;
     2420        int        i;
    24212421        int        tria_vertex_ids[3];
    24222422        IssmDouble nodeinputs[3];
     
    24292429
    24302430        /*Checks if debuging*/
    2431         /*{{{*/
    24322431        _assert_(iomodel->elements);
    2433         /*}}}*/
    24342432
    24352433        /*Fetch parameters: */
     
    29832981
    29842982        /*Intermediaries*/
    2985         int        i,j;
    29862983        IssmDouble xyz_list[NUMVERTICES][3];
    29872984        IssmDouble viscosity,newviscosity,oldviscosity;
     
    30273024                newviscosity=viscosity+viscosity_overshoot*(viscosity-oldviscosity);
    30283025                D_scalar=2.*newviscosity*thickness*gauss->weight*Jdet;
    3029                 for(i=0;i<3;i++) D[i*3+i]=D_scalar;
     3026                for(int i=0;i<3;i++) D[i*3+i]=D_scalar;
    30303027
    30313028                TripleMultiply(B,3,numdof,1,
     
    30513048        /*Intermediaries*/
    30523049        bool       mainlyfloating;
    3053         int        i,j;
    30543050        int        analysis_type,migration_style;
    30553051        int        point1;
     
    31143110                GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
    31153111                D_scalar=alpha2*gauss->weight*Jdet;
    3116                 for(i=0;i<2;i++) D[i*2+i]=D_scalar;
     3112                for(int i=0;i<2;i++) D[i*2+i]=D_scalar;
    31173113
    31183114                TripleMultiply(B,2,numdof,1,
     
    64706466        /*Intermediaries */
    64716467        int        stabilization;
    6472         int        dim;
    64736468        IssmDouble Jdet,D_scalar,dt,h;
    64746469        IssmDouble vel,vx,vy;
     
    65686563        /*Intermediaries */
    65696564        int        stabilization;
    6570         int        dim;
    65716565        IssmDouble Jdet,D_scalar,dt,h;
    65726566        IssmDouble vel,vx,vy;
     
    70157009
    70167010        /*Intermediaries */
    7017         int        stabilization;
    7018         int        i,j,dim;
     7011        int        stabilization,dim;
    70197012        IssmDouble Jdet,vx,vy,dvxdx,dvydy,vel,h;
    70207013        IssmDouble D_scalar;
  • issm/trunk-jpl/src/c/classes/Inputs/TransientInput.cpp

    r15737 r16151  
    267267void TransientInput::GetInputUpToCurrentTimeAverages(IssmDouble** pvalues, IssmDouble** ptimes, int* pnumtimes, IssmDouble currenttime){
    268268
    269         int i;
    270         IssmDouble* times=NULL;
    271         IssmDouble* values=NULL;
    272         int numsteps;
    273         int count;
    274         bool iscurrenttime_included=false;
     269        int         i;
     270        IssmDouble *times  = NULL;
     271        IssmDouble *values = NULL;
     272        int         numsteps;
     273        bool        iscurrenttime_included = false;
    275274
    276275        /*Figure out how many time steps we are going to return: */
  • issm/trunk-jpl/src/c/classes/Loads/Pengrid.cpp

    r16125 r16151  
    634634        int        reset_penalties = 0;
    635635        int        found           = 0;
    636         int        penalty_lock;
    637636        int        new_active;
    638637        IssmDouble pressure;
Note: See TracChangeset for help on using the changeset viewer.