\left(\left(3 + \frac{2}{r \cdot r}\right) - \frac{\left(0.125 \cdot \left(3 - 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 - v}\right) - 4.5\left(\left(3 + \frac{2}{r \cdot r}\right) - \frac{\left(0.125 \cdot \left(3 - 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 - v}\right) - 4.5double code(double r, double v, double w) {
return ((double) (((double) (((double) (3.0 + ((double) (2.0 / ((double) (r * r)))))) - ((double) (((double) (((double) (0.125 * ((double) (3.0 - ((double) (2.0 * v)))))) * ((double) (((double) (((double) (w * w)) * r)) * r)))) / ((double) (1.0 - v)))))) - 4.5));
}
double code(double r, double v, double w) {
return ((double) (((double) (((double) (3.0 + ((double) (2.0 / ((double) (r * r)))))) - ((double) (((double) (((double) (0.125 * ((double) (3.0 - ((double) (2.0 * v)))))) * ((double) (((double) (((double) (w * w)) * r)) * r)))) / ((double) (1.0 - v)))))) - 4.5));
}



Bits error versus r



Bits error versus v



Bits error versus w
Results
Initial program 13.3
Final simplification13.3
herbie shell --seed 2020153
(FPCore (r v w)
:name "(- (- (+ 3 (/ 2 (* r r))) (/ (* (* 1/8 (- 3 (* 2 v))) (* (* (* w w) r) r)) (- 1 v))) 9/2)"
:precision binary64
(- (- (+ 3.0 (/ 2.0 (* r r))) (/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v))) 4.5))