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) {
double VAR;
if ((((double) (w * w)) <= 4.885314099265168e+300)) {
VAR = ((double) (3.0 + ((double) (((double) (2.0 / ((double) (r * r)))) + ((double) (((double) (r * ((double) (((double) (w * r)) * ((double) (w * ((double) (((double) (0.125 / ((double) (1.0 - v)))) * ((double) (((double) (2.0 * v)) - 3.0)))))))))) - 4.5))))));
} else {
VAR = ((double) (3.0 + ((double) (((double) (2.0 / ((double) (r * r)))) + ((double) (((double) (((double) (r * ((double) (w * ((double) (((double) cbrt(r)) * ((double) cbrt(r)))))))) * ((double) (((double) (w * ((double) (((double) (0.125 / ((double) (1.0 - v)))) * ((double) (((double) (2.0 * v)) - 3.0)))))) * ((double) cbrt(r)))))) - 4.5))))));
}
return VAR;
}



Bits error versus v



Bits error versus w



Bits error versus r
Results
if (* w w) < 4.88531409926516791e300Initial program 9.3
Simplified4.7
rmApplied associate-*r*0.3
rmApplied associate-*l*0.3
Simplified0.3
if 4.88531409926516791e300 < (* w w) Initial program 62.2
Simplified61.5
rmApplied associate-*r*30.1
rmApplied add-cube-cbrt30.2
Applied associate-*l*30.2
Simplified19.9
rmApplied associate-*r*1.0
Simplified1.0
Final simplification0.3
herbie shell --seed 2020190
(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))