VecGhost trouble

Jed Brown jed at 59A2.org
Wed May 20 11:49:29 CDT 2009


I have some heavy use of VecGhost with MatShell (so the ghosts are
required to define MatMult).  While exercising this framework, I have
noticed some operations that do not work correctly, in particular

* MatComputeExplicitOperator uses VecCreateMPI rather than MatGetVecs

* same pattern in MatNullSpaceTest (if there are no vectors to
duplicate in the null space, i.e. just the constant)

I think these cases can just use MatGetVecs which I can use to produce
the correct vectors.  Any problems with this?  (Any reason the current
code calls VecCreateMPI directly?)

Jed

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 260 bytes
Desc: OpenPGP digital signature
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20090520/64249ead/attachment.sig>


More information about the petsc-dev mailing list