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