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

Smith, Barry F. bsmith at mcs.anl.gov
Sat Jan 4 14:17:19 CST 2020

   Can you overload the MatCreateSubMatrices() to use your function instead of the default. Using MatSetOperation()?


> On Jan 4, 2020, at 5:30 AM, Pierre Jolivet <pierre.jolivet at enseeiht.fr> wrote:
> 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

More information about the petsc-dev mailing list