\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1.0}\frac{\frac{\left(\left(\alpha + \beta\right) + i\right) \cdot i}{\left(\alpha + \beta\right) + i \cdot 2}}{\left(\left(\alpha + \beta\right) + i \cdot 2\right) + \sqrt{1.0}} \cdot \frac{\sqrt{\frac{\left(\left(\alpha + \beta\right) + i\right) \cdot i + \beta \cdot \alpha}{\left(\alpha + \beta\right) + i \cdot 2}}}{\frac{\left(\left(\alpha + \beta\right) + i \cdot 2\right) - \sqrt{1.0}}{\sqrt{\frac{\left(\left(\alpha + \beta\right) + i\right) \cdot i + \beta \cdot \alpha}{\left(\alpha + \beta\right) + i \cdot 2}}}}double f(double alpha, double beta, double i) {
double r4360309 = i;
double r4360310 = alpha;
double r4360311 = beta;
double r4360312 = r4360310 + r4360311;
double r4360313 = r4360312 + r4360309;
double r4360314 = r4360309 * r4360313;
double r4360315 = r4360311 * r4360310;
double r4360316 = r4360315 + r4360314;
double r4360317 = r4360314 * r4360316;
double r4360318 = 2.0;
double r4360319 = r4360318 * r4360309;
double r4360320 = r4360312 + r4360319;
double r4360321 = r4360320 * r4360320;
double r4360322 = r4360317 / r4360321;
double r4360323 = 1.0;
double r4360324 = r4360321 - r4360323;
double r4360325 = r4360322 / r4360324;
return r4360325;
}
double f(double alpha, double beta, double i) {
double r4360326 = alpha;
double r4360327 = beta;
double r4360328 = r4360326 + r4360327;
double r4360329 = i;
double r4360330 = r4360328 + r4360329;
double r4360331 = r4360330 * r4360329;
double r4360332 = 2.0;
double r4360333 = r4360329 * r4360332;
double r4360334 = r4360328 + r4360333;
double r4360335 = r4360331 / r4360334;
double r4360336 = 1.0;
double r4360337 = sqrt(r4360336);
double r4360338 = r4360334 + r4360337;
double r4360339 = r4360335 / r4360338;
double r4360340 = r4360327 * r4360326;
double r4360341 = r4360331 + r4360340;
double r4360342 = r4360341 / r4360334;
double r4360343 = sqrt(r4360342);
double r4360344 = r4360334 - r4360337;
double r4360345 = r4360344 / r4360343;
double r4360346 = r4360343 / r4360345;
double r4360347 = r4360339 * r4360346;
return r4360347;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
Initial program 52.6
rmApplied add-sqr-sqrt52.6
Applied difference-of-squares52.6
Applied times-frac38.7
Applied times-frac36.7
rmApplied add-sqr-sqrt36.8
Applied associate-/l*36.8
Final simplification36.8
herbie shell --seed 2019152 +o rules:numerics
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/4"
:pre (and (> alpha -1) (> beta -1) (> i 1))
(/ (/ (* (* i (+ (+ alpha beta) i)) (+ (* beta alpha) (* i (+ (+ alpha beta) i)))) (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i)))) (- (* (+ (+ alpha beta) (* 2 i)) (+ (+ alpha beta) (* 2 i))) 1.0)))