[petsc-users] Nested Fieldsplit for custom index sets

[Safin, Artur]

Barry, Lawrence,

> I think the SubKSPs (and therefore SubPCs) are not set up until you call KSPSetUp(ksp) which your code does not do explicitly and is therefore done in KSPSolve.

I added KSPSetUp(), but unfortunately the issue did not go away.



I have created a MWE that replicates the issue. The program tries to solve a tridiagonal system, where the first fieldsplit partitions the global matrix

[ P  x ]
[ x  T ],

and the nested fieldsplit partitions P into

[ A  x ]
[ x  B ].

Thanks for your help,

Artur

-------------- next part --------------
A non-text attachment was scrubbed...
Name: ex.c
Type: text/x-csrc
Size: 2395 bytes
Desc: ex.c
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160728/a019776b/attachment.c>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: run.sh
Type: application/x-shellscript
Size: 351 bytes
Desc: run.sh
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160728/a019776b/attachment.bin>