On Wed, May 6, 2020 at 11:42 AM Lawrence Mitchell <wence@gmx.li> wrote:
Hi Andrew,

> On 6 May 2020, at 16:39, Andrew Hicks <ahick17@lsu.edu> wrote:
>
> Hello all,

> I have a simple question. I’m trying to solve an elliptic PDE in 3D, but the regular solver in firedrake is just too slow. What is the best 3D solver in firedrake for this? Perhaps there is an algebraic multigrid solver or something along those lines? Any suggestions would be highly appreciated!

There are a number of multigrid options available (via PETSc).

Try with "pc_type": "hypre" in your solver_parameters dictionary for the solver.

It may need some tuning for best performance for 3D problems, but if you have not been setting any parameters at all up to now, it should be much better.

Hypre is the best for 2D. In 3D, agglomeration multigrid becomes competitive, -pc_type gamg.

  Thanks,

     Matt
 
Cheers,

Lawrence
_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake


--
What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
-- Norbert Wiener