<div dir="ltr"><div dir="ltr"><div dir="ltr">Hi all,<div><br></div><div>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. </div><div><br></div><div>From the PCFieldSplitSetDefault() routine, it looks like DMCreateSubDM() does exactly this <a href="https://bitbucket.org/petsc/petsc/src/048b3c40cbdab0ededa3294e9c569ef147398dbc/src/ksp/pc/impls/fieldsplit/fieldsplit.c?at=master&fileviewer=file-view-default#fieldsplit.c-433">here</a>. 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.</div><div><br></div><div>Thanks,</div><div>Justin</div></div></div></div>