Implement difference method for directional restriction

This returns `f(t)-f(0)`. This is potentially more stable to evaluate
than conputing the difference manually.
4 jobs for master in 10 minutes and 25 seconds (queued for 2 seconds)
latest