6 May
2020
6 May
'20
7:18 p.m.
On Wed, 6 May 2020 at 19:12, Andrew Hicks <ahick17@lsu.edu> wrote:
Hi Matt,
Actually, the “hypre” pc-type worked just fine, it took far less time than originally, but still a bit slow.
I tried the “gamg” pc-type but got the following error from which I’ve cut out the middle part:
petsc4py.PETSc.Error: error code 83
…
[0] Require AIJ matrix type
I’m working with three spatial dimensions but the vector I’m using is 5D, so it’s the PDE consists of solving for a function from R^3 into R^5. Is that why this won’t work?
If you use a VectorFunctionSpace, by default we make a baij (block aij) matrix. Gamg needs aij instead. If you additionally add "mat_type": "aij" to the solve parameters it should work. Lawrence