\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}\;\beta \le 1.57967512348392177 \cdot 10^{214}:\\
\;\;\;\;\frac{\frac{\frac{1}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot 1}} \cdot \frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{\sqrt{\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{\frac{\frac{1}{\sqrt{\left(\alpha + \beta\right) + 2 \cdot 1}} \cdot \left(\sqrt{0.5} \cdot \left(\beta + 0.75 \cdot \alpha\right) + \left(1 \cdot \sqrt{0.5} - 0.125 \cdot \frac{\beta}{\sqrt{0.5}}\right)\right)}{\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 ((double) (((double) (((double) (((double) (((double) (((double) (alpha + beta)) + ((double) (beta * alpha)))) + 1.0)) / ((double) (((double) (alpha + beta)) + ((double) (2.0 * 1.0)))))) / ((double) (((double) (alpha + beta)) + ((double) (2.0 * 1.0)))))) / ((double) (((double) (((double) (alpha + beta)) + ((double) (2.0 * 1.0)))) + 1.0))));
}
double code(double alpha, double beta) {
double VAR;
if ((beta <= 1.5796751234839218e+214)) {
VAR = ((double) (((double) (((double) (((double) (1.0 / ((double) sqrt(((double) (((double) (alpha + beta)) + ((double) (2.0 * 1.0)))))))) * ((double) (((double) (((double) (((double) (alpha + beta)) + ((double) (beta * alpha)))) + 1.0)) / ((double) sqrt(((double) (((double) (alpha + beta)) + ((double) (2.0 * 1.0)))))))))) / ((double) (((double) (alpha + beta)) + ((double) (2.0 * 1.0)))))) / ((double) (((double) (((double) (alpha + beta)) + ((double) (2.0 * 1.0)))) + 1.0))));
} else {
VAR = ((double) (((double) (((double) (((double) (1.0 / ((double) sqrt(((double) (((double) (alpha + beta)) + ((double) (2.0 * 1.0)))))))) * ((double) (((double) (((double) sqrt(0.5)) * ((double) (beta + ((double) (0.75 * alpha)))))) + ((double) (((double) (1.0 * ((double) sqrt(0.5)))) - ((double) (0.125 * ((double) (beta / ((double) sqrt(0.5)))))))))))) / ((double) (((double) (alpha + beta)) + ((double) (2.0 * 1.0)))))) / ((double) (((double) (((double) (alpha + beta)) + ((double) (2.0 * 1.0)))) + 1.0))));
}
return VAR;
}



Bits error versus alpha



Bits error versus beta
Results
if beta < 1.5796751234839218e+214Initial program 1.8
rmApplied add-sqr-sqrt2.3
Applied *-un-lft-identity2.3
Applied times-frac2.3
if 1.5796751234839218e+214 < beta Initial program 20.5
rmApplied add-sqr-sqrt20.5
Applied *-un-lft-identity20.5
Applied times-frac20.5
Taylor expanded around 0 5.9
Simplified5.9
Final simplification2.7
herbie shell --seed 2020148
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/3"
:precision binary64
:pre (and (> alpha -1.0) (> beta -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)))