double f(double alpha, double beta) {
double r6923281 = beta;
double r6923282 = alpha;
double r6923283 = r6923281 - r6923282;
double r6923284 = r6923282 + r6923281;
double r6923285 = 2.0;
double r6923286 = r6923284 + r6923285;
double r6923287 = r6923283 / r6923286;
double r6923288 = 1.0;
double r6923289 = r6923287 + r6923288;
double r6923290 = r6923289 / r6923285;
return r6923290;
}
double f(double alpha, double beta) {
double r6923291 = alpha;
double r6923292 = 1.2004402010900637e+17;
bool r6923293 = r6923291 <= r6923292;
double r6923294 = beta;
double r6923295 = r6923294 + r6923291;
double r6923296 = 2.0;
double r6923297 = r6923295 + r6923296;
double r6923298 = r6923294 / r6923297;
double r6923299 = r6923298 * r6923298;
double r6923300 = r6923298 * r6923299;
double r6923301 = r6923291 / r6923297;
double r6923302 = 1.0;
double r6923303 = r6923301 - r6923302;
double r6923304 = r6923303 * r6923303;
double r6923305 = r6923304 * r6923303;
double r6923306 = r6923300 - r6923305;
double r6923307 = log(r6923306);
double r6923308 = r6923298 * r6923303;
double r6923309 = r6923304 + r6923308;
double r6923310 = r6923299 + r6923309;
double r6923311 = log(r6923310);
double r6923312 = r6923307 - r6923311;
double r6923313 = exp(r6923312);
double r6923314 = r6923313 / r6923296;
double r6923315 = 4.0;
double r6923316 = r6923291 * r6923291;
double r6923317 = r6923315 / r6923316;
double r6923318 = r6923296 / r6923291;
double r6923319 = r6923317 - r6923318;
double r6923320 = 8.0;
double r6923321 = r6923320 / r6923291;
double r6923322 = r6923321 / r6923316;
double r6923323 = r6923319 - r6923322;
double r6923324 = r6923298 - r6923323;
double r6923325 = r6923324 / r6923296;
double r6923326 = r6923293 ? r6923314 : r6923325;
return r6923326;
}
\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2.0} + 1.0}{2.0}\begin{array}{l}
\mathbf{if}\;\alpha \le 1.2004402010900637 \cdot 10^{+17}:\\
\;\;\;\;\frac{e^{\log \left(\frac{\beta}{\left(\beta + \alpha\right) + 2.0} \cdot \left(\frac{\beta}{\left(\beta + \alpha\right) + 2.0} \cdot \frac{\beta}{\left(\beta + \alpha\right) + 2.0}\right) - \left(\left(\frac{\alpha}{\left(\beta + \alpha\right) + 2.0} - 1.0\right) \cdot \left(\frac{\alpha}{\left(\beta + \alpha\right) + 2.0} - 1.0\right)\right) \cdot \left(\frac{\alpha}{\left(\beta + \alpha\right) + 2.0} - 1.0\right)\right) - \log \left(\frac{\beta}{\left(\beta + \alpha\right) + 2.0} \cdot \frac{\beta}{\left(\beta + \alpha\right) + 2.0} + \left(\left(\frac{\alpha}{\left(\beta + \alpha\right) + 2.0} - 1.0\right) \cdot \left(\frac{\alpha}{\left(\beta + \alpha\right) + 2.0} - 1.0\right) + \frac{\beta}{\left(\beta + \alpha\right) + 2.0} \cdot \left(\frac{\alpha}{\left(\beta + \alpha\right) + 2.0} - 1.0\right)\right)\right)}}{2.0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\beta}{\left(\beta + \alpha\right) + 2.0} - \left(\left(\frac{4.0}{\alpha \cdot \alpha} - \frac{2.0}{\alpha}\right) - \frac{\frac{8.0}{\alpha}}{\alpha \cdot \alpha}\right)}{2.0}\\
\end{array}


Bits error versus alpha



Bits error versus beta
if alpha < 1.2004402010900637e+17Initial program 0.5
rmApplied div-sub0.5
Applied associate-+l-0.5
rmApplied add-exp-log0.5
rmApplied flip3--0.5
Applied log-div0.5
Simplified0.5
if 1.2004402010900637e+17 < alpha Initial program 50.7
rmApplied div-sub50.7
Applied associate-+l-49.1
Taylor expanded around -inf 17.7
Simplified17.7
Final simplification6.1
herbie shell --seed 2019101
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/1"
:pre (and (> alpha -1) (> beta -1))
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0) 2.0))