Re: [petsc-users] PCHPDDM and matrix type

2023-04-17 Thread Pierre Jolivet
1) PCHPDDM handles AIJ, BAIJ, SBAIJ, IS, NORMAL, NORMALHERMITIAN, SCHURCOMPLEMENT, HTOOL 2) This PC is based on domain decomposition, with no support yet for “over decomposition”. If you run with a single process, it’s like PCASM or PCBJACOBI, you’ll get the same behavior as if you were just

Re: [petsc-users] PCHPDDM and matrix type

2023-04-17 Thread Matthew Knepley
Yes, I cannot figure the error out. We will wait for Pierre to weigh in. Thanks, Matt On Mon, Apr 17, 2023 at 7:26 PM Alexander Lindsay wrote: > I don't really get much more of a stack trace out: > > [0]PETSC ERROR: [1]PETSC ERROR: - Error Message >

Re: [petsc-users] PCHPDDM and matrix type

2023-04-17 Thread Alexander Lindsay
I don't really get much more of a stack trace out: [0]PETSC ERROR: [1]PETSC ERROR: - Error Message -- [0]PETSC ERROR: Invalid argument [0]PETSC ERROR: - Error Message

Re: [petsc-users] PCHPDDM and matrix type

2023-04-17 Thread Matthew Knepley
I don't think so. Can you show the whole stack? THanks, Matt On Mon, Apr 17, 2023 at 6:24 PM Alexander Lindsay wrote: > 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

Re: [petsc-users] PCHPDDM and matrix type

2023-04-17 Thread 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 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

[petsc-users] PCHPDDM and matrix type

2023-04-17 Thread Alexander Lindsay
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