[petsc-dev] [petsc-users] Is xout always already zero'ed when being called back on PCShellSetApply ?

Matthew Knepley knepley at gmail.com
Sun May 21 11:41:10 CDT 2017


On Sun, May 21, 2017 at 11:11 AM, Franck Houssen <franck.houssen at inria.fr>
wrote:

> When using http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/
> PCShellSetApply.html, do I have the guaranty that xout from "
> PetscErrorCode
> <http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Sys/PetscErrorCode.html#PetscErrorCode>
> apply (PC
> <http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PC.html#PC>
> pc,Vec
> <http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/Vec.html#Vec>
> xin,Vec
> <http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/Vec.html#Vec>
> xout)" has always been previously filled with zeros ?
> I may have to fill xout by blocks (I would += several possibly overlapping
> blocks: I need to make sure xout is first filled with zero to get the
> correct result).
>

No, we do not initialize the output vector. Yo ucan call VecSet(xout, 0.0);

  Thanks,

    Matt


> Franck
>



-- 
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

http://www.caam.rice.edu/~mk51/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20170521/bf822ced/attachment.html>


More information about the petsc-dev mailing list