Summary
Instead of first creating a global basis and then passing this together with a grid to a ProblemStat, the latter can be constructed directly from the underlying pre-basis factory. This allows to write code like
Grid grid = ...;
ProblemStat prob{"name", grid, lagrange<1>()};