[petsc-dev] PCASM with custom overlap/local matrices

Pierre Jolivet pierre.jolivet at enseeiht.fr
Sat Jan 4 05:30:03 CST 2020

I’d like to bypass the call to MatCreateSubMatrices during PCSetUp_PCASM because I’m using a custom (geometric) overlap (PCASMSetLocalSubdomains) and I have a direct access to the corresponding subdomain matrices with optimized boundary conditions.
Currently, I’m using PCSetModifySubMatrices which just overrides the local Dirichlet matrix with my own, so there is no real need for the call to MatCreateSubMatrices.
Do you see of a way to do this?


More information about the petsc-dev mailing list