\left(x + v \cdot dt\right) + \left(\left(a \cdot \frac{2}{3} - aprev \cdot \frac{1}{6}\right) \cdot dt\right) \cdot dtdt \cdot \left(v + \left(a \cdot \frac{2}{3} - aprev \cdot \frac{1}{6}\right) \cdot dt\right) + xdouble code(double x, double v, double dt, double a, double aprev) {
return ((double) (((double) (x + ((double) (v * dt)))) + ((double) (((double) (((double) (((double) (a * ((double) (2.0 / 3.0)))) - ((double) (aprev * ((double) (1.0 / 6.0)))))) * dt)) * dt))));
}
double code(double x, double v, double dt, double a, double aprev) {
return ((double) (((double) (dt * ((double) (v + ((double) (((double) (((double) (a * ((double) (2.0 / 3.0)))) - ((double) (aprev * ((double) (1.0 / 6.0)))))) * dt)))))) + x));
}



Bits error versus x



Bits error versus v



Bits error versus dt



Bits error versus a



Bits error versus aprev
Results
Initial program 0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2020153
(FPCore (x v dt a aprev)
:name "(+ (+ x (* v dt)) (* (* (- (* a (/ 2 3)) (* aprev (/ 1 6))) dt) dt))"
:precision binary64
(+ (+ x (* v dt)) (* (* (- (* a (/ 2.0 3.0)) (* aprev (/ 1.0 6.0))) dt) dt)))