\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}\begin{array}{l}
\mathbf{if}\;\beta \le 4.670506575327117809269966612039077622849 \cdot 10^{208}:\\
\;\;\;\;\frac{\frac{i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i} \cdot \frac{\frac{\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) - \sqrt{1}}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + \sqrt{1}}\\
\mathbf{else}:\\
\;\;\;\;0\\
\end{array}double f(double alpha, double beta, double i) {
double r122823 = i;
double r122824 = alpha;
double r122825 = beta;
double r122826 = r122824 + r122825;
double r122827 = r122826 + r122823;
double r122828 = r122823 * r122827;
double r122829 = r122825 * r122824;
double r122830 = r122829 + r122828;
double r122831 = r122828 * r122830;
double r122832 = 2.0;
double r122833 = r122832 * r122823;
double r122834 = r122826 + r122833;
double r122835 = r122834 * r122834;
double r122836 = r122831 / r122835;
double r122837 = 1.0;
double r122838 = r122835 - r122837;
double r122839 = r122836 / r122838;
return r122839;
}
double f(double alpha, double beta, double i) {
double r122840 = beta;
double r122841 = 4.670506575327118e+208;
bool r122842 = r122840 <= r122841;
double r122843 = i;
double r122844 = alpha;
double r122845 = r122844 + r122840;
double r122846 = r122845 + r122843;
double r122847 = r122843 * r122846;
double r122848 = 2.0;
double r122849 = r122848 * r122843;
double r122850 = r122845 + r122849;
double r122851 = r122847 / r122850;
double r122852 = r122840 * r122844;
double r122853 = r122852 + r122847;
double r122854 = r122853 / r122850;
double r122855 = 1.0;
double r122856 = sqrt(r122855);
double r122857 = r122850 - r122856;
double r122858 = r122854 / r122857;
double r122859 = r122851 * r122858;
double r122860 = r122850 + r122856;
double r122861 = r122859 / r122860;
double r122862 = 0.0;
double r122863 = r122842 ? r122861 : r122862;
return r122863;
}



Bits error versus alpha



Bits error versus beta



Bits error versus i
Results
if beta < 4.670506575327118e+208Initial program 52.9
rmApplied add-sqr-sqrt52.9
Applied difference-of-squares52.9
Applied times-frac38.1
Applied times-frac35.9
rmApplied associate-*l/35.9
if 4.670506575327118e+208 < beta Initial program 64.0
rmApplied add-sqr-sqrt64.0
Applied difference-of-squares64.0
Applied times-frac57.2
Applied times-frac56.2
rmApplied add-sqr-sqrt56.2
Applied associate-/l*56.2
Taylor expanded around inf 44.4
Final simplification36.9
herbie shell --seed 2019325
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/4"
:precision binary64
: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)))