Support for precise event times and adaptive solvers

This version is an overhaul of the solver runtime to use adaptive solvers and root finder routine for precise interpolation of event times.