[petsc-users] PCs and MATIS

Barry Smith bsmith at petsc.dev
Tue Dec 14 15:00:32 CST 2021

   Do the "overlapping meshes" "match" in the overlap regions or are you connecting completely different meshes discretizations by boundary conditions along the edges of all the sub meshes?  In other words, will your global linear system be defined by your overlapping meshes? 

   If they match but you want to use more general boundary conditions for the subproblems than PCASM supports by default you might be able to use PCSetModifySubMatrices() to allow you to modify the sub matrices before they are used in the preconditioned; you can for example modify the entries along the boundary of the domain to represent Robin's conditions. Or you can put whatever you want into the entire submatrix if modifying them is too tedious.

   But if the meshes don't match then you don't really need a new PC you need to even define what your nonlinear system is and you have a very big project to write a PDE solver for non-matching overlapping grids using PETSc.


> On Dec 14, 2021, at 3:37 PM, Eric Chamberland <Eric.Chamberland at giref.ulaval.ca> wrote:
> Hi,
> We want to use an Additive Schwarz preconditioner (like PCASM) combined with overlapping meshes *and* specific boundary conditions to local (MATIS) matrices.
> At first sight, MATIS is only supported for BDDC and FETI-DP and is not working with PCASM.
> Do we have to write a new PC from scratch to combine the use of mesh overlap, MATIS and customized local matrices?
> ...or is there any working example we should look at to start from? :)
> Thanks a lot!
> Eric
> -- 
> Eric Chamberland, ing., M. Ing
> Professionnel de recherche
> GIREF/Université Laval
> (418) 656-2131 poste 41 22 42

More information about the petsc-users mailing list