Is your feature request related to a problem? Please describe.
We want internal calculations for initial conditions to be done at compiler precision, instead of casting them at the end.
Describe the solution you'd like
See discussion regarding precision in #50