Mesh adaption
This will probably come only after we generalized the discretization schemes (#11 (closed)), but in the long run FEM and FD could support adaptive mesh refinement. That could give a huge performance boost.
Until then, any discussion can happen here.