FEM interface
It would be cool to have a working FEM (finite element method) discretization in the long run.
I am unsure, whether it's best to rely on some (likely more performant) third-party library (e.g. scikit-fem
?) or whether we could simply assemble the matrices ourselves.