make integrators type stable

Created by: PhilipVinc

otherwise we upcast the type of parameters when integrating the RK equations, and since the tableaus are stored in double precision, this breaks the callback if the model has float32 parameters