Geometric Multigrid solver changes b/f on lower levels
The geometric multigrid solver currently changes the right hand side (f/b) on lower levels. Not the one that is passed to solve(...)
This might lead to unexpected results in the application.
Otherwise another temporary function must be created or passed to the solver.