\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(3 + \frac{\frac{2}{r}}{r}\right) - \left(\frac{0.125 \cdot \left(3 - 2 \cdot v\right)}{\frac{1 - v}{{\left(\left|w \cdot r\right|\right)}^{2}}} + 4.5\right)double f(double v, double w, double r) {
double r15162 = 3.0;
double r15163 = 2.0;
double r15164 = r;
double r15165 = r15164 * r15164;
double r15166 = r15163 / r15165;
double r15167 = r15162 + r15166;
double r15168 = 0.125;
double r15169 = v;
double r15170 = r15163 * r15169;
double r15171 = r15162 - r15170;
double r15172 = r15168 * r15171;
double r15173 = w;
double r15174 = r15173 * r15173;
double r15175 = r15174 * r15164;
double r15176 = r15175 * r15164;
double r15177 = r15172 * r15176;
double r15178 = 1.0;
double r15179 = r15178 - r15169;
double r15180 = r15177 / r15179;
double r15181 = r15167 - r15180;
double r15182 = 4.5;
double r15183 = r15181 - r15182;
return r15183;
}
double f(double v, double w, double r) {
double r15184 = 3.0;
double r15185 = 2.0;
double r15186 = r;
double r15187 = r15185 / r15186;
double r15188 = r15187 / r15186;
double r15189 = r15184 + r15188;
double r15190 = 0.125;
double r15191 = v;
double r15192 = r15185 * r15191;
double r15193 = r15184 - r15192;
double r15194 = r15190 * r15193;
double r15195 = 1.0;
double r15196 = r15195 - r15191;
double r15197 = w;
double r15198 = r15197 * r15186;
double r15199 = fabs(r15198);
double r15200 = 2.0;
double r15201 = pow(r15199, r15200);
double r15202 = r15196 / r15201;
double r15203 = r15194 / r15202;
double r15204 = 4.5;
double r15205 = r15203 + r15204;
double r15206 = r15189 - r15205;
return r15206;
}



Bits error versus v



Bits error versus w



Bits error versus r
Results
Initial program 12.8
Simplified8.4
rmApplied associate-*l*2.6
Simplified2.6
rmApplied add-sqr-sqrt2.6
Simplified2.6
Simplified0.3
rmApplied associate-/r*0.3
rmApplied fma-udef0.3
Simplified0.4
Final simplification0.4
herbie shell --seed 2020047 +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))