\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) - \left(\sqrt{\frac{0.125 \cdot \left(3 - 2 \cdot v\right)}{1 - v}} \cdot \left|\left|w \cdot r\right|\right|\right) \cdot \sqrt{\frac{0.125 \cdot \left(3 - 2 \cdot v\right)}{1 - v} \cdot {\left(\left|w \cdot r\right|\right)}^{2}}\right) - 4.5double f(double v, double w, double r) {
double r27329 = 3.0;
double r27330 = 2.0;
double r27331 = r;
double r27332 = r27331 * r27331;
double r27333 = r27330 / r27332;
double r27334 = r27329 + r27333;
double r27335 = 0.125;
double r27336 = v;
double r27337 = r27330 * r27336;
double r27338 = r27329 - r27337;
double r27339 = r27335 * r27338;
double r27340 = w;
double r27341 = r27340 * r27340;
double r27342 = r27341 * r27331;
double r27343 = r27342 * r27331;
double r27344 = r27339 * r27343;
double r27345 = 1.0;
double r27346 = r27345 - r27336;
double r27347 = r27344 / r27346;
double r27348 = r27334 - r27347;
double r27349 = 4.5;
double r27350 = r27348 - r27349;
return r27350;
}
double f(double v, double w, double r) {
double r27351 = 3.0;
double r27352 = 2.0;
double r27353 = r;
double r27354 = r27353 * r27353;
double r27355 = r27352 / r27354;
double r27356 = r27351 + r27355;
double r27357 = 0.125;
double r27358 = v;
double r27359 = r27352 * r27358;
double r27360 = r27351 - r27359;
double r27361 = r27357 * r27360;
double r27362 = 1.0;
double r27363 = r27362 - r27358;
double r27364 = r27361 / r27363;
double r27365 = sqrt(r27364);
double r27366 = w;
double r27367 = r27366 * r27353;
double r27368 = fabs(r27367);
double r27369 = fabs(r27368);
double r27370 = r27365 * r27369;
double r27371 = 2.0;
double r27372 = pow(r27368, r27371);
double r27373 = r27364 * r27372;
double r27374 = sqrt(r27373);
double r27375 = r27370 * r27374;
double r27376 = r27356 - r27375;
double r27377 = 4.5;
double r27378 = r27376 - r27377;
return r27378;
}



Bits error versus v



Bits error versus w



Bits error versus r
Results
Initial program 12.8
Simplified8.5
rmApplied add-sqr-sqrt8.6
Simplified8.5
Simplified0.4
rmApplied fma-udef0.4
Applied associate--r+0.4
Simplified0.4
rmApplied add-sqr-sqrt0.4
rmApplied sqrt-prod0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019326 +o rules:numerics
(FPCore (v w r)
:name "Rosa's TurbineBenchmark"
:precision binary64
(- (- (+ 3 (/ 2 (* r r))) (/ (* (* 0.125 (- 3 (* 2 v))) (* (* (* w w) r) r)) (- 1 v))) 4.5))