[petsc-users] PCHPDDM and matrix type

[Alexander Lindsay]

If it helps: if I use those exact same options in serial, then no errors
and the linear solve is beautiful :-)

On Mon, Apr 17, 2023 at 4:22 PM Alexander Lindsay <alexlindsay239 at gmail.com>
wrote:

> I'm likely revealing a lot of ignorance, but in order to use HPDDM as a
> preconditioner does my system matrix (I am using the same matrix for A and
> P) need to be block type, e.g. baij or sbaij ? In MOOSE our default is aij
> and I am currently getting
>
> [1]PETSC ERROR: #1 buildTwo() at
> /raid/lindad/moose/petsc/arch-moose/include/HPDDM_schwarz.hpp:1012
>
> with options:
>
> -pc_type hpddm -pc_hpddm_block_splitting -pc_hpddm_coarse_mat_type baij
> -pc_hpddm_coarse_pc_type lu -pc_hpddm_define_subdomains
> -pc_hpddm_levels_1_eps_gen_non_hermitian -pc_hpddm_levels_1_eps_nev 50
> -pc_hpddm_levels_1_eps_threshold 0.1 -pc_hpddm_levels_1_st_matstructure
> SAME -pc_hpddm_levels_1_st_share_sub_ksp
> -pc_hpddm_levels_1_sub_pc_factor_mat_solver_type mumps
> -pc_hpddm_levels_1_sub_pc_type lu
>
> Alex
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230417/3efe26bd/attachment.html>