Changeset 27097


Ignore:
Timestamp:
06/23/22 17:26:19 (3 years ago)
Author:
Eric.Larour
Message:

CHG: fixing leaks.

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

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/analyses/SealevelchangeAnalysis.cpp

    r27061 r27097  
    99
    1010/*Model processing*/
     11SealevelchangeAnalysis::~SealevelchangeAnalysis(){/*{{{*/
     12
     13        BarystaticContributions* barycontrib=NULL;
     14        GenericParam<BarystaticContributions*>* barycontribparam=NULL;
     15
     16        barycontribparam = xDynamicCast<GenericParam<BarystaticContributions*>*>(femmodel->parameters->FindParamObject(BarystaticContributionsEnum));
     17        barycontrib=barycontribparam->GetParameterValue();
     18
     19        delete barycontrib;
     20}/*}}}*/
    1121void SealevelchangeAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/*{{{*/
    1222        /*No constraints*/
    … …  
    285295                        parameters->AddObject(new BoolParam(LoveIsTimeEnum,love_istime));
    286296
     297                        /*Free allocations:*/
     298                        xDelete<IssmDouble>(love_timefreq);
     299
    287300                        // AD performance is sensitive to calls to ensurecontiguous.
    288301                        // // Providing "t" will cause ensurecontiguous to be called.
    … …  
    532545                                nt=1; //in elastic, or if we run only selfattraction, we need only one step
    533546#ifdef _HAVE_AD_
    534                                 G_viscoelastic_interpolated=G_viscoelastic;
    535                                 U_viscoelastic_interpolated=U_viscoelastic;
    536                                 if(horiz)H_viscoelastic_interpolated=H_viscoelastic;
     547                                G_viscoelastic_interpolated=xNew<IssmDouble>(M,"t");
     548                                U_viscoelastic_interpolated=xNew<IssmDouble>(M,"t");
     549                                if(horiz) H_viscoelastic_interpolated=xNew<IssmDouble>(M,"t");
    537550#else
    538                                 G_viscoelastic_interpolated=G_viscoelastic;
    539                                 U_viscoelastic_interpolated=U_viscoelastic;
    540                                 if(horiz)H_viscoelastic_interpolated=H_viscoelastic;
     551                                G_viscoelastic_interpolated=xNew<IssmDouble>(M);
     552                                U_viscoelastic_interpolated=xNew<IssmDouble>(M);
     553                                if(horiz) H_viscoelastic_interpolated=xNew<IssmDouble>(M);
    541554#endif
     555
     556
     557                                xMemCpy<IssmDouble>(G_viscoelastic_interpolated,G_viscoelastic,M);
     558                                xMemCpy<IssmDouble>(U_viscoelastic_interpolated,U_viscoelastic,M);
     559                                if(horiz) xMemCpy<IssmDouble>(H_viscoelastic_interpolated,H_viscoelastic,M);
    542560
    543561                                if(rotation){ //if this cpu handles degree 2
    … …  
    558576#endif
    559577
    560                                         Pmtf_col_interpolated=love_pmtf_colinear;
    561                                         Pmtf_ortho_interpolated=love_pmtf_ortho;
     578                                        Pmtf_col_interpolated[0]=love_pmtf_colinear[0];
     579                                        Pmtf_ortho_interpolated[0]=love_pmtf_ortho[0];
    562580                                        Pmtf_z_interpolated[0]=1.0+love_k[2];
    563581                                        Love_tk2_interpolated[0]=love_tk[2];
    … …  
    595613                                xDelete<IssmDouble>(G_viscoelastic);
    596614                                xDelete<IssmDouble>(G_viscoelastic_local);
     615                                xDelete<IssmDouble>(G_viscoelastic_interpolated);
    597616                                xDelete<IssmDouble>(U_viscoelastic);
     617                                xDelete<IssmDouble>(U_viscoelastic_interpolated);
    598618                                xDelete<IssmDouble>(U_viscoelastic_local);
    599619                                if(horiz){
    600620                                        xDelete<IssmDouble>(H_viscoelastic);
     621                                        xDelete<IssmDouble>(H_viscoelastic_interpolated);
    601622                                        xDelete<IssmDouble>(H_viscoelastic_local);
    602623                                }
    … …  
    604625                                        xDelete<IssmDouble>(love_pmtf_colinear);
    605626                                        xDelete<IssmDouble>(love_pmtf_ortho);
     627                                       
     628                                        xDelete<IssmDouble>(Love_tk2_interpolated);
     629                                        xDelete<IssmDouble>(Love_th2_interpolated);
     630                                        xDelete<IssmDouble>(Pmtf_col_interpolated);
     631                                        xDelete<IssmDouble>(Pmtf_ortho_interpolated);
     632                                        xDelete<IssmDouble>(Pmtf_z_interpolated);
     633                                        if(horiz){
     634                                                xDelete<IssmDouble>(Love_tl2_interpolated);
     635                                        }
    606636
    607637                                }
    … …  
    632662        iomodel->DeleteData(&requestedoutputs,numoutputs,"md.solidearth.requested_outputs");
    633663        /*}}}*/
     664
    634665}/*}}}*/
    635666

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

    r27061 r27097  
    65836583                LoveRotRSL  = xNewZeroInit<IssmDouble>(nt);
    65846584                LoveRotU    = xNewZeroInit<IssmDouble>(nt);
    6585                 LoveRothoriz= xNewZeroInit<IssmDouble>(nt);
     6585                if(horiz)LoveRothoriz= xNewZeroInit<IssmDouble>(nt);
    65866586                Grot        = xNewZeroInit<IssmDouble>(3*3*nt); //3 polar motion components * 3 vertices * number of time steps
    65876587                GUrot       = xNewZeroInit<IssmDouble>(3*3*nt);
    … …  
    66696669                        }
    66706670                }
     6671                /*Free ressources:*/
     6672                xDelete<IssmDouble>(LoveRotRSL);
     6673                xDelete<IssmDouble>(LoveRotU);
     6674                if(horiz)xDelete<IssmDouble>(LoveRothoriz);
    66716675        }
    66726676        /*}}}*/
    … …  
    73927396        /*Keep track of barystatic contributions:*/
    73937397        barycontrib->Set(this->Sid(),bslcice,bslchydro,bslcbp);
     7398
     7399        /*Free ressources:*/
     7400        xDelete<IssmDouble>(areae);
    73947401
    73957402}
    … …  
    77087715                        /*Save viscous stack now that we updated the values:*/
    77097716                        this->inputs->SetArrayInput(viscousenum,this->lid,viscousfield,3*viscousnumsteps);
     7717
     7718                        /*Free ressources:*/
     7719                        xDelete<IssmDouble>(grdfieldinterp);
    77107720                }
    77117721
    … …  
    77267736                for(i=0;i<NUMVERTICES;i++) grdfieldout[i]=grdfield[i*nt+0];
    77277737        }
     7738        /*Free ressources:*/
     7739        xDelete<IssmDouble>(grdfield);
    77287740
    77297741} /*}}}*/

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

    r26800 r27097  
    1919GrdLoads::GrdLoads(int nel,SealevelGeometry* slgeom){ /*{{{*/
    2020
     21        /*allocate:*/
    2122        vloads=new Vector<IssmDouble>(nel);
    2223        for (int i=0;i<SLGEOM_NUMLOADS;i++) vsubloads[i]=new Vector<IssmDouble>(slgeom->nbar[i]);
    … …  
    2627
    2728        vsubsealevelloads=new Vector<IssmDouble>(slgeom->nbar[SLGEOM_OCEAN]);
     29
     30        /*make sure all pointers that are not allocated are NULL:*/
     31        loads=NULL;
     32        for (int i=0;i<SLGEOM_NUMLOADS;i++)subloads[i]=NULL;
     33        sealevelloads=NULL;
     34        subsealevelloads=NULL;
     35
    2836
    2937}; /*}}}*/
    … …  
    5058        }
    5159
     60        /*Avoid leaks:*/
     61        if(loads)xDelete<IssmDouble>(loads);
     62        for (int i=0;i<SLGEOM_NUMLOADS;i++){
     63                if(subloads[i])xDelete<IssmDouble>(subloads[i]);
     64        }
     65
     66        /*Serialize:*/
    5267        loads=vloads->ToMPISerial();
    5368        for (int i=0;i<SLGEOM_NUMLOADS;i++){
    … …  
    6479void GrdLoads::BroadcastSealevelLoads(void){ /*{{{*/
    6580
     81        /*Avoid leakds:*/
     82        if(sealevelloads)xDelete<IssmDouble>(sealevelloads);
     83        if(subsealevelloads)xDelete<IssmDouble>(subsealevelloads);
     84
     85        /*Serialize:*/
    6686        sealevelloads=vsealevelloads->ToMPISerial();
    6787        subsealevelloads=vsubsealevelloads->ToMPISerial();
    6888
    6989} /*}}}*/
    70 void GrdLoads::SHDegree2Coefficients(IssmDouble* deg2coeff, FemModel* femmodel, SealevelGeometry* slgeom){
     90void GrdLoads::SHDegree2Coefficients(IssmDouble* deg2coeff, FemModel* femmodel, SealevelGeometry* slgeom){ /*{{{*/
    7191
    7292        IssmDouble area,re, S;
    … …  
    119139        ISSM_MPI_Bcast(&deg2coeff[2],1,ISSM_MPI_DOUBLE,0,IssmComm::GetComm());
    120140
    121 }
     141} /*}}}*/

  • issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp

    r27061 r27097  
    8181        /*End profiler*/
    8282        femmodel->profiler->Stop(SLRCORE);
     83
     84        /*Free ressources:*/
     85        delete slgeom;
     86
    8387}
    8488/*}}}*/
    … …  
    364368
    365369                //convergence?
    366                 if(slcconvergence(loads->vsealevelloads,oldsealevelloads,eps_rel,eps_abs))break;
     370                if(slcconvergence(loads->vsealevelloads,oldsealevelloads,eps_rel,eps_abs)){
     371                        delete oldsealevelloads; break;
     372                }
    367373
    368374                //early return?
    369                 if(iterations>=max_nonlinear_iterations)break;
     375                if(iterations>=max_nonlinear_iterations){
     376                        delete oldsealevelloads; break;
     377                }
    370378                iterations++;
     379                delete oldsealevelloads;
    371380        }
    372381
    … …  
    388397        if(planethasocean){
    389398                if (!sealevelloading){ //we haven't done so before, so we need to compute the ocean average and area
    390                         loads->sealevelloads=NULL; //needed to trigger the calculation of areas
     399                        if(loads->sealevelloads)xDelete<IssmDouble>(loads->sealevelloads); loads->sealevelloads=NULL; //needed to trigger the calculation of areas
    391400                        /*retrieve sea level average  and ocean area:*/
    392401                        for(Object* & object : femmodel->elements->objects){
    … …  
    430439        }
    431440
     441        /*Free ressources:*/
     442        delete loads;
     443        delete subelementoceanareas;
     444        delete oceanareas;
     445        xDelete<IssmDouble>(sealevelpercpu);
     446
    432447}
    433448/*}}}*/
    … …  
    607622        #endif
    608623
     624        xDelete<IssmDouble>(xxe);
     625        xDelete<IssmDouble>(yye);
     626        xDelete<IssmDouble>(zze);
     627        xDelete<IssmDouble>(areae);
     628
    609629        return;
    610630
    … …  
    672692        femmodel->parameters->AddObject(new DoubleVecParam(ZzeEnum,zze,nel));
    673693        femmodel->parameters->AddObject(new DoubleVecParam(AreaeEnum,areae,nel));
     694
     695        /*Free ressources:*/
     696        xDelete<IssmDouble>(xxe);
     697        xDelete<IssmDouble>(yye);
     698        xDelete<IssmDouble>(zze);
     699        xDelete<IssmDouble>(areae);
    674700
    675701        return slgeom;
    … …  
    870896                        xDelete<IssmDouble>(m3interp);
    871897                }
    872         }
     898                xDelete<IssmDouble>(viscoustimes);
     899                xDelete<IssmDouble>(viscouspolarmotion);
     900        }
     901        xDelete<IssmDouble>(pmtf_col);
     902        xDelete<IssmDouble>(pmtf_ortho);
     903        xDelete<IssmDouble>(pmtf_z);
    873904
    874905} /*}}}*/

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

    r26327 r27097  
    223223                                                                                        {
    224224                                                                                                Tria* tria= xDynamicCast<Tria*>(element);
     225                                                                                                if(tria->hmaterial)delete tria->hmaterial;
    225226                                                                                                tria->hmaterial=NULL;
    226227                                                                                        }
Note: See TracChangeset for help on using the changeset viewer.