[petsc-dev] DMDA_*PERIODIC and DMDA_XYZGHOSTED

Satish Balay balay at mcs.anl.gov
Thu Mar 10 17:49:44 CST 2011


patch now added to petsc-dev

satish

On Thu, 10 Mar 2011, Ethan Coon wrote:

> 
> 
> 
> > >
> > > Hmm, I didn't add the dd->idx, so I'm suprised it wasn't free'd already.
> > > I also wasn't sure when/if to log the memory, since it seemed to be in
> > > there for some stuff and not for others in the previous code.
> > >
> > 
> > You should just log it in the DM object, just uncomment the line in
> > the 2D and 3D setup code
> > 
> 
> Ok, done.
> 
> > > With the copy (your next email) to the IS, shouldn't that changeset leak
> > > memory?  I think it also needs a PetscFree(idx) then as well...
> > >
> > 
> > My two patches are cumulative... Take into account that DMDestroy_DA
> > applies to 1/2/3D, while the setup code is different. The free(idx) at
> > the common destoy() should be enough (provided you use COPY_VALUES for
> > the 1D code)
> > 
> 
> Ah, I see my confusion... previously the memory had been shared with the
> da->l2gmap, and so it was destroyed by that.  I missed the fact that
> your second change was in da1.  This passes valgrind's memcheck.  Thanks
> for clearing that up, Lisandro.
> 
> This patch does the change from DMDAPeriodicType -> DMDABoundaryType and
> deals with Lisandro's memory fixes.
> 
> The one thing I didn't change is the name of the DMDASetPeriodicity()
> method, which I wasn't sure what to go to.
> 
> Ethan
> 
> 
> 
> > >
> > > Will get these fixes in,
> > >
> > 
> > With both patches applied, petsc4py testsuite pass cleanly (and mem
> > leaks seems to go away). So I'm more or less confident on your (huge)
> > changes + my (tiny) fixes
> > 
> > 
> 
> 




More information about the petsc-dev mailing list