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

[Pierre Jolivet]

Hello,
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?

Thanks,
Pierre