But the main changes of this update are bug fixes to saveat and improved Rosenbrock methods. There was a bug in saveat such that for large changes in dt there could be oscillations in the output. Also, errors would occur if values of saveat matched tstops. Both of these bugs have been fixed, and saveat should work properly with all combinations of tstops, save_timeseries, and dense, letting you mix and match the behavior.
In addition, the Rosenbrock23 and Rosenbrock32 methods have received some updates. There have been some minor performance improvements, but more importantly the Rosenbrock methods now have a “stiff-safe” interpolant. Before it had used the standard 3rd order Hermite interpolation, which “in theory” since it’s order 3 was “good enough”, but in reality for really stiff equations this showed oscillations. A new interpolation has been placed on this method which utilizes the internal steps to not oscillate (and have less error) on stiff equations, making these Rosenbrock methods now a very good general method for stiff equations with dense (continuous) output.