[petsc-users] Can a MatShell that operates on ghosted vectors work in place?

Jed Brown jedbrown at mcs.anl.gov
Mon Aug 20 10:15:31 CDT 2012

On Mon, Aug 20, 2012 at 10:12 AM, Matthew Knepley <knepley at gmail.com> wrote:

> Personally, I think this is quite cumbersome compared to the Local/Global
> vector strategy we use in DM, and
> does not have any real advantages.

For example, storage for the ghost values is allocated for every vector
duplicated from a VecGhost (e.g. the entire Krylov basis). This isn't an
issue if you have large subdomains, but for smallish subdomains using
methods with a fat stencil, that can double the storage needed for the
Krylov basis.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120820/9afc032d/attachment.html>

More information about the petsc-users mailing list