Add option to output predicted values

(but code in MicromagEquations not yet implemented, so doesn't work yet)
Remove support for old IMR implementation from ArrayInterpolator

* Should be faster

* Simplifies code ready for adding varying time step support
Add interpolation of history values

In ArrayInterpolator