On Tue, Nov 8, 2011 at 10:37 AM, Jose E. Roman <jroman at dsic.upv.es> wrote:
> The following code prints a vector of uninitialized entries. > > ierr = VecCreate(comm,&v);CHKERRQ(ierr); > ierr = VecSetSizes(v,n,n);CHKERRQ(ierr); > ierr = VecSetType(v,VECSEQCUSP);CHKERRQ(ierr); > ierr = VecView(v,PETSC_NULL);CHKERRQ(ierr); > > In contrast, if the type is VECSEQ the vector is initialized to zero. So > it seems that VecSet(v,0.0) must always be called. Is this the expected > behaviour? There is no guarantee of initialization for any Vec type. Matt > > Jose > -- 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/20111108/68a810ad/attachment.html>
