[REF] Execute tests in float64

[f-dangel]

I repeatedly find myself having to bump tolerances in tests when I insert code that generates random numbers.
This is annoying and leads to long-term degradation of the tests because tolerances always increase, but never shrink.
I think this should be fixed by carrying out all tests in float64, which will also allow to use small tolerances.