\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{\frac{2}{r}}{r}\right) - \left(0.125 \cdot \left|w \cdot r\right|\right) \cdot \frac{3 - 2 \cdot v}{\frac{1 - v}{\left|w \cdot r\right|}}\right) - 4.5double f(double v, double w, double r) {
double r27739 = 3.0;
double r27740 = 2.0;
double r27741 = r;
double r27742 = r27741 * r27741;
double r27743 = r27740 / r27742;
double r27744 = r27739 + r27743;
double r27745 = 0.125;
double r27746 = v;
double r27747 = r27740 * r27746;
double r27748 = r27739 - r27747;
double r27749 = r27745 * r27748;
double r27750 = w;
double r27751 = r27750 * r27750;
double r27752 = r27751 * r27741;
double r27753 = r27752 * r27741;
double r27754 = r27749 * r27753;
double r27755 = 1.0;
double r27756 = r27755 - r27746;
double r27757 = r27754 / r27756;
double r27758 = r27744 - r27757;
double r27759 = 4.5;
double r27760 = r27758 - r27759;
return r27760;
}
double f(double v, double w, double r) {
double r27761 = 3.0;
double r27762 = 2.0;
double r27763 = r;
double r27764 = r27762 / r27763;
double r27765 = r27764 / r27763;
double r27766 = r27761 + r27765;
double r27767 = 0.125;
double r27768 = w;
double r27769 = r27768 * r27763;
double r27770 = fabs(r27769);
double r27771 = r27767 * r27770;
double r27772 = v;
double r27773 = r27762 * r27772;
double r27774 = r27761 - r27773;
double r27775 = 1.0;
double r27776 = r27775 - r27772;
double r27777 = r27776 / r27770;
double r27778 = r27774 / r27777;
double r27779 = r27771 * r27778;
double r27780 = r27766 - r27779;
double r27781 = 4.5;
double r27782 = r27780 - r27781;
return r27782;
}



Bits error versus v



Bits error versus w



Bits error versus r
Results
Initial program 12.8
rmApplied add-sqr-sqrt12.8
Simplified12.8
Simplified6.8
rmApplied associate-/l*0.4
Simplified0.4
rmApplied associate-/r*0.4
rmApplied add-sqr-sqrt0.6
Applied unpow-prod-down0.6
Applied *-un-lft-identity0.6
Applied times-frac0.5
Applied times-frac0.5
Simplified0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019323 +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))