[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