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