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

[Smith, Barry F.]

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

   Barry


> 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