[petsc-dev] VECCUSP needs to be zeroed?

[Jed Brown]

On Wed, Nov 9, 2011 at 13:25, Matthew Knepley <knepley at gmail.com> wrote:

> The simpler answer is yes, I always initialize the Mat or Vec to zero when
> I use it. Even if its initialized when created,
> I almost always get something from DMGetGlobalVector() etc. which could
> have been used before.
>

Since we keep track of access that could mutate, it would be easy to make
VecZeroEntries a no-op when the array is already zeroed. (In debug mode, it
could actually check that the array is still zero.)

I agree with Matt that you should call VecZeroEntries() when you want to
rely on the Vec being zeroed.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111109/44cf9169/attachment.html>