How to use pySDC with Fenics?

[zhyxue0]

Hello, everyone! I am a newer on pySDC. I notice that pySDC works with FEniCS, mpi4py-fft and PETSc (through petsc4py). Which versions of Fenics are compatible? Are there any examples of using it with Fenics? Is it suitable for firedrakes like Fenics?

Thanks a lot!

[pancetta]

Hi @zhyxue0, thanks for your inquiry. pySDC currently supports FEniCS 2019.x, not the new version which is currently being developed. I tested using FEniCSx a few months ago and for simple algorithms (like SDC) this seems to work. For multilevel algorithms like MLSDC or PFASST they were still missing interpolation operators and I could not figure out a way to get this to work. This may have changed since then, though.

To see how pySDC works with FEniCS 2019, see the corresponding tutorial.

Concerning firedrake: this does not work out of the box and probably also not inside the box.. I know that Jemma Shipton from Exeter and her team there played with this. She may know more about the status.

[pancetta]

Firedrake coupling now also works: https://github.com/Parallel-in-Time/pySDC/pull/513