double code(double v, double w, double r) {
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 v, double w, double r) {
return ((double) (3.0 + ((double) (((double) (2.0 * ((double) pow(r, -2.0)))) + ((double) (((double) (((double) (((double) cbrt(r)) * ((double) cbrt(r)))) * ((double) (((double) (r * w)) * ((double) (((double) (w * ((double) (((double) (0.125 / ((double) (1.0 - v)))) * ((double) (((double) (2.0 * v)) - 3.0)))))) * ((double) (((double) cbrt(((double) (((double) cbrt(r)) * ((double) cbrt(r)))))) * ((double) cbrt(((double) cbrt(r)))))))))))) - 4.5))))));
}



Bits error versus v



Bits error versus w



Bits error versus r
Results
Initial program 12.7
Simplified8.7
rmApplied associate-*r*2.5
rmApplied add-cube-cbrt2.7
Applied associate-*l*2.7
Simplified2.0
rmApplied add-cube-cbrt2.1
Applied cbrt-prod2.1
rmApplied div-inv2.1
Simplified2.0
Final simplification2.0
herbie shell --seed 2020191
(FPCore (v w r)
:name "Rosa's TurbineBenchmark"
:precision binary64
(- (- (+ 3.0 (/ 2.0 (* r r))) (/ (* (* 0.125 (- 3.0 (* 2.0 v))) (* (* (* w w) r) r)) (- 1.0 v))) 4.5))