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



Bits error versus v



Bits error versus w



Bits error versus r
Results
if (* w w) < 2.0062619874181415e217Initial program 8.8
Simplified4.6
rmApplied associate-*r*4.6
rmApplied associate-*r*0.3
Simplified0.3
if 2.0062619874181415e217 < (* w w) Initial program 43.2
Simplified37.3
rmApplied associate-*r*37.3
rmApplied associate-*r*19.2
Simplified19.2
rmApplied associate-*r*0.3
Final simplification0.3
herbie shell --seed 2020181
(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))