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



Bits error versus x



Bits error versus v



Bits error versus dt



Bits error versus a



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