[petsc-users] DMComposite and Matrix construction

[Jed Brown]

Gideon Simpson <gideon.simpson at gmail.com> writes:
> I suspect this is because DMCreateMatrix is picking a sparsity pattern which is not consistent with what I need.

If the issue is in the diagonal blocks, then make sure the component DMs
are correct.  If in the off-diagonal blocks, you can use
DMCompositeSetCoupling to add nonzeros for the coupling terms.  This
interface is a bit ugly and not generic, but hopefully it will work for
you.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150729/09a8bf8c/attachment.pgp>