Thorough check that ALL tests use epsilons for floating point comparison in boundary conditions and spatial params
This is the main reason for failing test. All tests should use a proper epsilon. * This epsilon usually is made of a base epsilon `static constexpr eps_ = 1.5e-7;` times a magnitude. In the case of grids this should be a characteristic cell size. * If the grid is in the scale of 1m there is no adjustment necessary. However a runtime check would make it flexible for using any grid * We could also introduce a preprocessor constant or global variable in a namespace for the epsilon * We can also use more of Dune::Float_Cmp
issue