[petsc-dev] Petsc VecCopy code.

[Yoo, Andy]

Thanks, Mathew.

From: petsc-dev-bounces at mcs.anl.gov [mailto:petsc-dev-bounces at mcs.anl.gov] On Behalf Of Matthew Knepley
Sent: Friday, May 04, 2012 11:41 AM
To: For users of the development version of PETSc
Subject: Re: [petsc-dev] Petsc VecCopy code.

On Fri, May 4, 2012 at 2:27 PM, Yoo, Andy <yoo2 at llnl.gov<mailto:yoo2 at llnl.gov>> wrote:
Hi,

I have written a two-dimensional MatVec routine (http://dl.acm.org/citation.cfm?id=2063469) a while ago.

Any chance of getting that into PETSc proper?

Now, I am noticing there is a small memory leak that appears to come from VecCopy. It is small, as it only appears like after 300000 MatVecs.
I looked at the Petsc sources, but was not able to locate the default copy function for VecMPI.

Can you please point me to the right direction?

Its here:

http://petsc.cs.iit.edu/petsc/petsc-dev/annotate/95406d5a1c14/src/vec/vec/impls/seq/bvec2.c#l190

set here in the vtable:

http://petsc.cs.iit.edu/petsc/petsc-dev/annotate/95406d5a1c14/src/vec/vec/impls/mpi/pbvec.c#l111

   Matt


Thank you for your help in advance.



Andy Yoo, Ph.D.
Center for Applied Scientific Computing
Lawrence Livermore National Laboratory




--
What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120504/8633d9fb/attachment.html>