Gideon Simpson <[email protected]> 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.
signature.asc
Description: PGP signature
