[petsc-users] DMCreateSubDM() not available in petsc4py

[Justin Chang]

Hi all,

I'm writing a petsc4py routine to manually create nested fieldsplits using
index sets, and it looks like whenever I move onto the next level of splits
I need to rescale the IS's.

>From the PCFieldSplitSetDefault() routine, it looks like DMCreateSubDM()
does exactly this here
<https://bitbucket.org/petsc/petsc/src/048b3c40cbdab0ededa3294e9c569ef147398dbc/src/ksp/pc/impls/fieldsplit/fieldsplit.c?at=master&fileviewer=file-view-default#fieldsplit.c-433>.
However, I see no such Python wrapper for this function. While I could
write my own, I kind of have my hands tied behind my back by using the
petsc4py that's attached to FEniCS - I wonder if there's a general strategy
or workaround for writing your own recursive fieldsplitsetIS.

Thanks,
Justin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190321/e8a209b6/attachment.html>