\left(\left(3 - \frac{2}{r \cdot r}\right) - \frac{\left(0.125 \cdot \left(1 + 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 - v}\right) - 0.5\left(\left(3 - \frac{2}{r \cdot r}\right) - \frac{\left(0.125 \cdot \left(1 + 2 \cdot v\right)\right) \cdot \left(\left(\left(w \cdot w\right) \cdot r\right) \cdot r\right)}{1 - v}\right) - 0.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) (1.0 + ((double) (2.0 * v)))))) * ((double) (((double) (((double) (w * w)) * r)) * r)))) / ((double) (1.0 - v)))))) - 0.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) (1.0 + ((double) (2.0 * v)))))) * ((double) (((double) (((double) (w * w)) * r)) * r)))) / ((double) (1.0 - v)))))) - 0.5));
}



Bits error versus r



Bits error versus v



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