<div class="gmail_extra">On Fri, May 4, 2012 at 2:27 PM, Yoo, Andy <span dir="ltr"><<a href="mailto:yoo2@llnl.gov" target="_blank">yoo2@llnl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div lang="EN-US" link="blue" vlink="purple"><p class="MsoNormal">Hi,<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">I have written a two-dimensional MatVec routine (<a href="http://dl.acm.org/citation.cfm?id=2063469" target="_blank">http://dl.acm.org/citation.cfm?id=2063469</a>) a while ago.</p>
</div></blockquote><div><br></div><div>Any chance of getting that into PETSc proper?</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="blue" vlink="purple">
<p class="MsoNormal"><u></u><u></u></p><p class="MsoNormal">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.<u></u><u></u></p><p class="MsoNormal">
I looked at the Petsc sources, but was not able to locate the default copy function for VecMPI.<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Can you please point me to the right direction?</p>
</div></blockquote><div><br></div><div>Its here:</div><div><br></div><div><a href="http://petsc.cs.iit.edu/petsc/petsc-dev/annotate/95406d5a1c14/src/vec/vec/impls/seq/bvec2.c#l190">http://petsc.cs.iit.edu/petsc/petsc-dev/annotate/95406d5a1c14/src/vec/vec/impls/seq/bvec2.c#l190</a></div>
<div><br></div><div>set here in the vtable:</div><div><br></div><div><a href="http://petsc.cs.iit.edu/petsc/petsc-dev/annotate/95406d5a1c14/src/vec/vec/impls/mpi/pbvec.c#l111">http://petsc.cs.iit.edu/petsc/petsc-dev/annotate/95406d5a1c14/src/vec/vec/impls/mpi/pbvec.c#l111</a></div>
<div><br></div><div>   Matt</div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="blue" vlink="purple"><p class="MsoNormal"> </p>
<p class="MsoNormal"><u></u></p><p class="MsoNormal">Thank you for your help in advance.<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Andy Yoo, Ph.D.<u></u><u></u></p><p class="MsoNormal">Center for Applied Scientific Computing<u></u><u></u></p><p class="MsoNormal">Lawrence Livermore National Laboratory<u></u><u></u></p><p class="MsoNormal">
<u></u> <u></u></p></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>
</div>