VecGhost trouble

Jed Brown jed at
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?)


-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 260 bytes
Desc: OpenPGP digital signature
URL: <>

More information about the petsc-dev mailing list