[petsc-dev] I do not think this is the right solution

Jed Brown jedbrown at mcs.anl.gov
Sun Mar 18 19:10:35 CDT 2012


On Sun, Mar 18, 2012 at 18:55, Matthew Knepley <knepley at gmail.com> wrote:

> No off-diagonal.


So I think the MatGetArray() should be MatDenseGetArray() because it's not
directly useful otherwise.

For in-place modification (if that is actually useful), I'd do it with a
callback.

Note that for your specific purpose, you can do in-place modification with
MatMPIAIJGetAIJ() and MatGetArray() on each off-diagonal block.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120318/187c5e93/attachment.html>


More information about the petsc-dev mailing list