\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\begin{array}{l}
\mathbf{if}\;\alpha \le 1.49557542553151027 \cdot 10^{178}:\\
\;\;\;\;\frac{\frac{\frac{1 \cdot \left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\alpha + \beta\right) + 2 \cdot 1}}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\\
\mathbf{else}:\\
\;\;\;\;\frac{0.25 \cdot \alpha + \left(0.5 + 0.25 \cdot \beta\right)}{\left(\alpha + \beta\right) \cdot \left(\alpha + \beta\right) - \left(2 \cdot 1\right) \cdot \left(2 \cdot 1\right)} \cdot \frac{\left(\alpha + \beta\right) - 2 \cdot 1}{\left(\left(\alpha + \beta\right) + 2 \cdot 1\right) + 1}\\
\end{array}double code(double alpha, double beta) {
return ((((((alpha + beta) + (beta * alpha)) + 1.0) / ((alpha + beta) + (2.0 * 1.0))) / ((alpha + beta) + (2.0 * 1.0))) / (((alpha + beta) + (2.0 * 1.0)) + 1.0));
}
double code(double alpha, double beta) {
double temp;
if ((alpha <= 1.4955754255315103e+178)) {
temp = (((((1.0 * ((alpha + beta) + (beta * alpha))) + 1.0) / ((alpha + beta) + (2.0 * 1.0))) / ((alpha + beta) + (2.0 * 1.0))) / (((alpha + beta) + (2.0 * 1.0)) + 1.0));
} else {
temp = ((((0.25 * alpha) + (0.5 + (0.25 * beta))) / (((alpha + beta) * (alpha + beta)) - ((2.0 * 1.0) * (2.0 * 1.0)))) * (((alpha + beta) - (2.0 * 1.0)) / (((alpha + beta) + (2.0 * 1.0)) + 1.0)));
}
return temp;
}



Bits error versus alpha



Bits error versus beta
Results
if alpha < 1.4955754255315103e+178Initial program 1.5
rmApplied *-un-lft-identity1.5
if 1.4955754255315103e+178 < alpha Initial program 16.9
rmApplied *-un-lft-identity16.9
Applied flip-+17.5
Applied associate-/r/17.5
Applied times-frac17.5
Simplified17.5
Taylor expanded around 0 7.2
Final simplification2.3
herbie shell --seed 2020049
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/3"
:precision binary64
:pre (and (> alpha -1) (> beta -1))
(/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1) (+ (+ alpha beta) (* 2 1))) (+ (+ alpha beta) (* 2 1))) (+ (+ (+ alpha beta) (* 2 1)) 1)))