Changeset 13813


Ignore:
Timestamp:
10/23/12 13:52:13 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: more cppcheck issues

Location:
issm/trunk-jpl/src/c
Files:
16 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/Container/Elements.cpp

    r13790 r13813  
    169169        int num_procs;
    170170
    171         Patch  *patch        = NULL;
    172         int    *resultsenums = NULL;
    173         int    *resultssizes  = NULL;
    174         int    *resultssteps  = NULL;
    175         IssmDouble *resultstimes = NULL;
    176         IssmDouble *vector_serial= NULL;
    177         Vector<IssmDouble>*     vector      = NULL;
    178         bool   io_gather;
    179         bool   results_as_patches;
    180         int    numberofvertices,numberofelements;
    181         int    numberofresults,vectorsize;
    182         int    rank;
    183         int    minrank;
     171        Patch               *patch              = NULL;
     172        int                 *resultsenums      = NULL;
     173        int                 *resultssizes       = NULL;
     174        int                 *resultssteps       = NULL;
     175        IssmDouble          *resultstimes      = NULL;
     176        IssmDouble          *vector_serial      = NULL;
     177        Vector<IssmDouble> *vector = NULL;
     178        bool                io_gather;
     179        bool                results_as_patches;
     180        int                 numberofvertices     ,numberofelements;
     181        int                 numberofresults      ,vectorsize;
     182        int                 rank;
     183        int                 minrank;
    184184
    185185        /*recover my_rank:*/

  • issm/trunk-jpl/src/c/Container/Inputs.cpp

    r13798 r13813  
    363363void  Inputs::DuplicateInput(int original_enum,int new_enum){
    364364
    365         Input* original=NULL;
    366         Input* copy=NULL;
    367 
    368365        /*Make a copy of the original input: */
    369         original=dynamic_cast<Input*>(this->GetInput(original_enum));
     366        Input* original=dynamic_cast<Input*>(this->GetInput(original_enum));
    370367        if(!original)_error_("could not find input with enum: " << EnumToStringx(original_enum));
    371         copy=dynamic_cast<Input*>(original->copy());
     368        Input* copy=dynamic_cast<Input*>(original->copy());
    372369
    373370        /*Change copy enum to reinitialized_enum: */
    … …  
    407404void  Inputs::AXPY(int MeshYEnum, IssmDouble scalar, int MeshXEnum){
    408405
    409         Input* xinput=NULL;
    410         Input* yinput=NULL;
    411 
    412406        /*Find x and y inputs: */
    413         xinput=dynamic_cast<Input*>(this->GetInput(MeshXEnum));
    414         yinput=dynamic_cast<Input*>(this->GetInput(MeshYEnum));
     407        Input* xinput=dynamic_cast<Input*>(this->GetInput(MeshXEnum));
     408        Input* yinput=dynamic_cast<Input*>(this->GetInput(MeshYEnum));
    415409
    416410        /*some checks: */

  • issm/trunk-jpl/src/c/Container/Nodes.cpp

    r13798 r13813  
    153153        int i;
    154154        int num_procs;
    155 
    156         int  numnodes;
     155        int numnodes;
    157156
    158157        /*recover num_procs: */

  • issm/trunk-jpl/src/c/classes/FemModel.cpp

    r13798 r13813  
    827827
    828828        /*Intermediary*/
    829         int      i;
    830829        int      num_responses;
    831830        Element *element       = NULL;

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

    r13787 r13813  
    45654565ElementMatrix* Penta::CreateKMatrixAdjointPattyn(void){
    45664566
    4567         /*Constants*/
    4568         const int    numdof=NDOF2*NUMVERTICES;
    4569 
    45704567        /*Intermediaries */
    45714568        int        i,j;
    … …  
    67816778
    67826779        /*Intermediaries */
    6783         int        i,j,approximation;
     6780        int        i,approximation;
    67846781        IssmDouble Jdet,viscosity,stokesreconditioning;
    67856782        IssmDouble xyz_list[NUMVERTICES][3];
    … …  
    69326929
    69336930        /*Intermediaries */
    6934         int         i,j;
    69356931        IssmDouble  Jdet;
    69366932        IssmDouble  xyz_list[NUMVERTICES][3];
    … …  
    71887184
    71897185        /*Intermediaries */
    7190         int         i,j;
     7186        int         i;
    71917187        int         approximation;
    71927188        IssmDouble  viscosity,Jdet;
    … …  
    73977393ElementVector* Penta::CreatePVectorDiagnosticHutter(void){
    73987394
    7399         /*Constants*/
    7400         const int numdofs=NDOF2*NUMVERTICES;
    7401 
    74027395        /*Intermediaries*/
    74037396        int          i,j;

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

    r13787 r13813  
    405405
    406406        /*Intermediaries */
    407         int        i,j,dim;
     407        int        dim;
    408408        IssmDouble xyz_list[NUMVERTICES][3];
    409409        IssmDouble Jdettria,dt,vx,vy;
    … …  
    593593
    594594        /*Intermediaries */
    595         int        i,j;
    596595        IssmDouble Jdettria,dt;
    597596        IssmDouble surface_mass_balance_g,basal_melting_g,basal_melting_correction_g,thickness_g;
    … …  
    629628                 basal_melting_correction_g=0.;
    630629
    631                 for(i=0;i<numdof;i++) pe->values[i]+=Jdettria*gauss->weight*(thickness_g+dt*(surface_mass_balance_g-basal_melting_g-basal_melting_correction_g))*L[i];
     630                for(int i=0;i<numdof;i++) pe->values[i]+=Jdettria*gauss->weight*(thickness_g+dt*(surface_mass_balance_g-basal_melting_g-basal_melting_correction_g))*L[i];
    632631        }
    633632
    … …  
    644643
    645644        /*Intermediaries */
    646         int        i,j;
    647645        IssmDouble Jdettria,dt;
    648646        IssmDouble surface_mass_balance_g,basal_melting_g,thickness_g;
    … …  
    674672                thickness_input->GetInputValue(&thickness_g,gauss);
    675673
    676                 for(i=0;i<numdof;i++) pe->values[i]+=Jdettria*gauss->weight*(thickness_g+dt*(surface_mass_balance_g-basal_melting_g))*L[i];
     674                for(int i=0;i<numdof;i++) pe->values[i]+=Jdettria*gauss->weight*(thickness_g+dt*(surface_mass_balance_g-basal_melting_g))*L[i];
    677675        }
    678676
    … …  
    689687
    690688        /*Intermediaries */
    691         int        i,j;
     689        int        i;
    692690        int        analysis_type;
    693691        IssmDouble Jdet;
    … …  
    13751373void  Tria::InputToResult(int enum_type,int step,IssmDouble time){
    13761374
    1377         int    i;
    13781375        Input *input = NULL;   
    13791376
    13801377        /*Go through all the input objects, and find the one corresponding to enum_type, if it exists: */
    1381         if (enum_type==MaterialsRheologyBbarEnum || enum_type==MaterialsRheologyZbarEnum) input=this->material->inputs->GetInput(enum_type);
    1382         else input=this->inputs->GetInput(enum_type);
     1378        if (enum_type==MaterialsRheologyBbarEnum || enum_type==MaterialsRheologyZbarEnum)
     1379         input=this->material->inputs->GetInput(enum_type);
     1380        else
     1381         input=this->inputs->GetInput(enum_type);
    13831382        //if (!input) _error_("Input " << EnumToStringx(enum_type) << " not found in tria->inputs");
    13841383        if(!input)return;
    … …  
    25642563IssmDouble Tria::MassFlux( IssmDouble* segment,bool process_units){
    25652564
    2566         const int    numdofs=2;
    2567 
    2568         int        i,dim;
    2569         IssmDouble     mass_flux=0;
    2570         IssmDouble     xyz_list[NUMVERTICES][3];
    2571         IssmDouble     normal[2];
    2572         IssmDouble     length,rho_ice;
    2573         IssmDouble     x1,y1,x2,y2,h1,h2;
    2574         IssmDouble     vx1,vx2,vy1,vy2;
     2565        const int  numdofs=2;
     2566        int        dim;
     2567        IssmDouble mass_flux=0;
     2568        IssmDouble xyz_list[NUMVERTICES][3];
     2569        IssmDouble normal[2];
     2570        IssmDouble length,rho_ice;
     2571        IssmDouble x1,y1,x2,y2,h1,h2;
     2572        IssmDouble vx1,vx2,vy1,vy2;
    25752573        GaussTria* gauss_1=NULL;
    25762574        GaussTria* gauss_2=NULL;
    … …  
    39873985
    39883986        /* Intermediaries */
    3989         int        ig;
    3990         IssmDouble     Jelem = 0;
    3991         IssmDouble     weight;
    3992         IssmDouble     Jdet;
    3993         IssmDouble     xyz_list[NUMVERTICES][3];
    3994         IssmDouble     dp[NDOF2];
     3987        IssmDouble Jelem = 0;
     3988        IssmDouble weight;
     3989        IssmDouble Jdet;
     3990        IssmDouble xyz_list[NUMVERTICES][3];
     3991        IssmDouble dp[NDOF2];
    39953992        GaussTria *gauss = NULL;
    39963993
    … …  
    40304027IssmDouble Tria::SurfaceAverageVelMisfit(bool process_units,int weight_index){
    40314028
    4032         const int    numdof=2*NUMVERTICES;
    4033 
    4034         int        i;
    40354029        IssmDouble Jelem=0,S,Jdet;
    40364030        IssmDouble misfit;
    … …  
    40914085IssmDouble Tria::SurfaceLogVelMisfit(bool process_units,int weight_index){
    40924086
    4093         const int    numdof=NDOF2*NUMVERTICES;
    4094 
    4095         int        i;
    40964087        IssmDouble Jelem=0;
    40974088        IssmDouble misfit,Jdet;
    … …  
    41564147IssmDouble Tria::SurfaceLogVxVyMisfit(bool process_units,int weight_index){
    41574148
    4158         const int    numdof=NDOF2*NUMVERTICES;
    4159 
    4160         int        i;
    41614149        int        fit=-1;
    41624150        IssmDouble Jelem=0, S=0;
    … …  
    42224210IssmDouble Tria::SurfaceAbsVelMisfit(bool process_units,int weight_index){
    42234211
    4224         const int    numdof=NDOF2*NUMVERTICES;
    4225 
    42264212        IssmDouble Jelem=0;
    42274213        IssmDouble misfit,Jdet;
    … …  
    42814267/*FUNCTION Tria::SurfaceRelVelMisfit {{{*/
    42824268IssmDouble Tria::SurfaceRelVelMisfit(bool process_units,int weight_index){
    4283         const int  numdof=2*NUMVERTICES;
    42844269
    42854270        IssmDouble  Jelem=0;
    … …  
    43464331
    43474332        /* Intermediaries */
    4348         int        ig;
    4349         IssmDouble     Jelem = 0;
    4350         IssmDouble     weight;
    4351         IssmDouble     Jdet;
    4352         IssmDouble     xyz_list[NUMVERTICES][3];
    4353         IssmDouble     dp[NDOF2];
     4333        IssmDouble Jelem = 0;
     4334        IssmDouble weight;
     4335        IssmDouble Jdet;
     4336        IssmDouble xyz_list[NUMVERTICES][3];
     4337        IssmDouble dp[NDOF2];
    43544338        GaussTria *gauss = NULL;
    43554339
    … …  
    43904374
    43914375        /* Intermediaries */
    4392         int         ig;
    43934376        IssmDouble  Jelem = 0;
    43944377        IssmDouble  weight;
    … …  
    44424425
    44434426        /* Intermediaries */
    4444         int         ig;
    44454427        IssmDouble  Jelem = 0;
    44464428        IssmDouble  weight;
    … …  
    44944476
    44954477        /*Intermediaries*/
    4496         int        i;
    44974478        IssmDouble thickness,thicknessobs,weight;
    44984479        IssmDouble Jdet;
    … …  
    49944975
    49954976        /* Intermediaries */
    4996         int        ig;
    4997         IssmDouble     Jelem = 0;
    4998         IssmDouble     weight;
    4999         IssmDouble     Jdet;
    5000         IssmDouble     xyz_list[NUMVERTICES][3];
    5001         IssmDouble     dp[NDOF2];
     4977        IssmDouble Jelem = 0;
     4978        IssmDouble weight;
     4979        IssmDouble Jdet;
     4980        IssmDouble xyz_list[NUMVERTICES][3];
     4981        IssmDouble dp[NDOF2];
    50024982        GaussTria *gauss = NULL;
    50034983
    … …  
    53165296        /*Intermediaries */
    53175297        IssmDouble diffusivity;
    5318         int        i,j;
    53195298        IssmDouble Jdettria,DL_scalar,dt,h;
    53205299        IssmDouble vx,vy,vel,dvxdx,dvydy;

  • issm/trunk-jpl/src/c/classes/objects/Inputs/ControlInput.cpp

    r13622 r13813  
    100100        output->control_id=this->control_id;
    101101
    102         if(values)      output->values=(Input*)this->values->copy();
    103         if(savedvalues) output->savedvalues=(Input*)this->savedvalues->copy();
    104         if(minvalues)   output->minvalues=(Input*)this->minvalues->copy();
    105         if(maxvalues)   output->maxvalues=(Input*)this->maxvalues->copy();
    106         if(gradient)    output->gradient=(Input*)this->gradient->copy();
     102        if(values)      output->values      = dynamic_cast<Input*>(this->values->copy());
     103        if(savedvalues) output->savedvalues = dynamic_cast<Input*>(this->savedvalues->copy());
     104        if(minvalues)   output->minvalues   = dynamic_cast<Input*>(this->minvalues->copy());
     105        if(maxvalues)   output->maxvalues   = dynamic_cast<Input*>(this->maxvalues->copy());
     106        if(gradient)    output->gradient    = dynamic_cast<Input*>(this->gradient->copy());
    107107
    108108        return output;
    … …  
    255255
    256256        if(savedvalues) delete this->savedvalues;
    257         this->savedvalues=(Input*)this->values->copy();
     257        this->savedvalues=dynamic_cast<Input*>(this->values->copy());
    258258}/*}}}*/
    259259/*FUNCTION ControlInput::UpdateValue{{{*/
    … …  
    263263
    264264        if(values) delete this->values;
    265         this->values=(Input*)this->savedvalues->copy();
     265        this->values=dynamic_cast<Input*>(this->savedvalues->copy());
    266266        this->values->AXPY(gradient,scalar);
    267267}/*}}}*/

  • issm/trunk-jpl/src/c/classes/objects/Inputs/DatasetInput.cpp

    r13622 r13813  
    8383        outinput=new DatasetInput();
    8484        outinput->enum_type=this->enum_type;
    85         outinput->inputs=(Inputs*)this->inputs->SpawnTriaInputs(indices);
     85        outinput->inputs=dynamic_cast<Inputs*>(this->inputs->SpawnTriaInputs(indices));
    8686
    8787        /*Assign output*/
    … …  
    110110        /*Get requested input within dataset*/
    111111        if(index<0 || index > inputs->Size()-1) _error_("index requested (" << index << ") exceeds dataset size (" << inputs->Size() << ")");
    112         Input* input=(Input*)this->inputs->GetObjectByOffset(index);
     112        Input* input=dynamic_cast<Input*>(this->inputs->GetObjectByOffset(index));
    113113
    114114        input->GetInputValue(pvalue,gauss);

  • issm/trunk-jpl/src/c/classes/objects/Loads/Icefront.cpp

    r13761 r13813  
    472472ElementVector* Icefront::CreatePVectorDiagnosticMacAyeal3d(void){
    473473
    474         Icefront* icefront=NULL;
    475         Penta*    penta=NULL;
    476         Tria*     tria=NULL;
    477         bool      onbed;
     474        Icefront *icefront = NULL;
     475        Penta    *penta    = NULL;
     476        Tria     *tria     = NULL;
    478477
    479478        /*Cast element onto Penta*/

  • issm/trunk-jpl/src/c/modules/ContourToNodesx/ContourToNodesx.cpp

    r13787 r13813  
    1313        double*  xc=NULL;
    1414        double*  yc=NULL;
    15         double   value;
    1615
    1716        /*output: */
    … …  
    3332        /*Assign output pointers: */
    3433        *pflags=flags;
    35 
    3634        return 1;
    3735}

  • issm/trunk-jpl/src/c/modules/ModelProcessorx/Dakota/CreateParametersDakota.cpp

    r13622 r13813  
    1616
    1717        /*variable declarations: {{{*/
    18         int i,j,k;
    19 
     18        int i;
    2019        Parameters* parameters = NULL;
    21         int      second_count;
    2220
    2321        int*     part=NULL;
    … …  
    3937        char tag[50];
    4038
    41         int      M;
    42         int      m,n;
    4339        bool     dakota_analysis=false;
    4440        char*    name=NULL;

  • issm/trunk-jpl/src/c/modules/Solverx/SolverxPetsc.cpp

    r13762 r13813  
    5858
    5959        /*Display message*/
     60        PetscSynchronizedPrintf(MPI_COMM_WORLD,"-------------- file: SolverxPetsc.cpp line: %i\n",__LINE__);
     61        PetscSynchronizedFlush(MPI_COMM_WORLD);
    6062        if(VerboseModule()) _pprintLine_("   Solving");
    6163        #if _PETSC_MAJOR_ < 3 || (_PETSC_MAJOR_ == 3 && _PETSC_MINOR_ < 2)
    … …  
    7173                *puf=NewVec(0,IssmComm::GetComm()); return;
    7274        }
     75        PetscSynchronizedPrintf(MPI_COMM_WORLD,"-------------- file: SolverxPetsc.cpp line: %i\n",__LINE__);
     76        PetscSynchronizedFlush(MPI_COMM_WORLD);
    7377
    7478        /*Initial guess */
    … …  
    8589        }
    8690        else{
     91                PetscSynchronizedPrintf(MPI_COMM_WORLD,"-------------- file: SolverxPetsc.cpp line: %i\n",__LINE__);
     92                PetscSynchronizedFlush(MPI_COMM_WORLD);
    8793                MatGetLocalSize(Kff,&local_m,&local_n);uf=NewVec(local_n,IssmComm::GetComm(),fromlocalsize);
    8894        }
    … …  
    104110        KSPSetOperators(ksp,Kff,Kff,DIFFERENT_NONZERO_PATTERN);
    105111        KSPSetFromOptions(ksp);
     112        PetscSynchronizedPrintf(MPI_COMM_WORLD,"-------------- file: SolverxPetsc.cpp line: %i\n",__LINE__);
     113        PetscSynchronizedFlush(MPI_COMM_WORLD);
    106114
    107115        #if _PETSC_MAJOR_==3
    … …  
    115123                #endif
    116124        }
     125        PetscSynchronizedPrintf(MPI_COMM_WORLD,"-------------- file: SolverxPetsc.cpp line: %i\n",__LINE__);
     126        PetscSynchronizedFlush(MPI_COMM_WORLD);
    117127
    118128        /*Stokes: */
    … …  
    146156        /*Solve: */
    147157        if(VerboseSolver())KSPView(ksp,PETSC_VIEWER_STDOUT_WORLD);
     158        PetscSynchronizedPrintf(MPI_COMM_WORLD,"-------------- file: SolverxPetsc.cpp line: %i\n",__LINE__);
     159        PetscSynchronizedFlush(MPI_COMM_WORLD);
    148160        KSPSolve(ksp,pf,uf);
     161        PetscSynchronizedPrintf(MPI_COMM_WORLD,"-------------- file: SolverxPetsc.cpp line: %i\n",__LINE__);
     162        PetscSynchronizedFlush(MPI_COMM_WORLD);
    149163
    150164        /*Check convergence*/
    … …  
    153167
    154168        /*Free resources:*/
     169        PetscSynchronizedPrintf(MPI_COMM_WORLD,"-------------- file: SolverxPetsc.cpp line: %i\n",__LINE__);
     170        PetscSynchronizedFlush(MPI_COMM_WORLD);
    155171        KSPFree(&ksp);
    156172

  • issm/trunk-jpl/src/c/modules/Solverx/SolverxSeq.cpp

    r13530 r13813  
    2424#ifdef _HAVE_GSL_
    2525        /*Intermediary: */
    26         int M,N,N2,s;
     26        int M,N,N2;
    2727        SeqVec<IssmDouble> *uf = NULL;
    2828

  • issm/trunk-jpl/src/c/shared/Exp/exp.h

    r13758 r13813  
    170170        doubletype         **pprofx        = NULL;
    171171        doubletype         **pprofy        = NULL;
    172         Contour<doubletype> *contour       = NULL;
    173172
    174173        /*output: */

  • issm/trunk-jpl/src/c/solutions/DakotaSpawnCore.cpp

    r13799 r13813  
    4343int DakotaSpawnCore(double* d_responses, int d_numresponses, double* d_variables, char** d_variables_descriptors,int d_numvariables, void* void_femmodel,int counter){
    4444
    45         int i;
    46 
    47         char   **responses_descriptors     = NULL; //these are our!  there are only numresponsedescriptors of them, not d_numresponses!!!
    48         int      numresponsedescriptors;
    49         char    *string                    = NULL;
    50         int      string_length;
    51         int      solution_type;
    52         bool     control_analysis          = false;
    53         void (*solutioncore)(FemModel*)    = NULL;
    54         FemModel* femmodel                 = NULL;
    55         bool      nodakotacore             = true;
     45        char     **responses_descriptors    = NULL;      //these are our! there are only numresponsedescriptors of them, not d_numresponses!!!
     46        int        numresponsedescriptors;
     47        char      *string                   = NULL;
     48        int        solution_type;
     49        bool       control_analysis         = false;
     50        void     (*solutioncore)(FemModel*) = NULL;
     51        FemModel  *femmodel                 = NULL;
     52        bool       nodakotacore             = true;
    5653
    5754        /*If counter==-1 on cpu0, it means that the dakota runs are done. In which case, bail out and return 0: */

  • issm/trunk-jpl/src/c/solutions/issm.cpp

    r13621 r13813  
    2121        /*Output results: */
    2222        femmodel->OutputResults();
     23        _error_("STOP HERE");
    2324
    2425        /*Wrap up: */
    … …  
    2627
    2728        /*Finalize environment:*/
    28         EnvironmentFinalize();
     29        //EnvironmentFinalize();
    2930
    3031        /*Finalize exception trapping: */
Note: See TracChangeset for help on using the changeset viewer.