\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}\begin{array}{l}
\mathbf{if}\;\frac{\beta - \alpha}{\left(\beta + \alpha\right) + 2} \leq -0.9999999999956046:\\
\;\;\;\;\left(\frac{\beta}{\alpha} + \frac{1}{\alpha}\right) - \left(3 \cdot \frac{\beta}{\alpha \cdot \alpha} + \left(\frac{\beta}{\alpha} \cdot \frac{\beta}{\alpha} + \frac{2}{\alpha \cdot \alpha}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;e^{\log \left(\frac{\frac{\beta}{\left(\beta + \alpha\right) + 2} - \left(\frac{\alpha}{\left(\beta + \alpha\right) + 2} - 1\right)}{2}\right)}\\
\end{array}(FPCore (alpha beta) :precision binary64 (/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0) 2.0))
(FPCore (alpha beta)
:precision binary64
(if (<= (/ (- beta alpha) (+ (+ beta alpha) 2.0)) -0.9999999999956046)
(-
(+ (/ beta alpha) (/ 1.0 alpha))
(+
(* 3.0 (/ beta (* alpha alpha)))
(+ (* (/ beta alpha) (/ beta alpha)) (/ 2.0 (* alpha alpha)))))
(exp
(log
(/
(-
(/ beta (+ (+ beta alpha) 2.0))
(- (/ alpha (+ (+ beta alpha) 2.0)) 1.0))
2.0)))))double code(double alpha, double beta) {
return (((beta - alpha) / ((alpha + beta) + 2.0)) + 1.0) / 2.0;
}
double code(double alpha, double beta) {
double tmp;
if (((beta - alpha) / ((beta + alpha) + 2.0)) <= -0.9999999999956046) {
tmp = ((beta / alpha) + (1.0 / alpha)) - ((3.0 * (beta / (alpha * alpha))) + (((beta / alpha) * (beta / alpha)) + (2.0 / (alpha * alpha))));
} else {
tmp = exp(log(((beta / ((beta + alpha) + 2.0)) - ((alpha / ((beta + alpha) + 2.0)) - 1.0)) / 2.0));
}
return tmp;
}












Bits error versus alpha












Bits error versus beta
Results
| Alternative 1 | |
|---|---|
| Error | 0.2 |
| Cost | 2689 |
| Alternative 2 | |
|---|---|
| Error | 0.2 |
| Cost | 2049 |
| Alternative 3 | |
|---|---|
| Error | 0.2 |
| Cost | 1793 |
| Alternative 4 | |
|---|---|
| Error | 0.2 |
| Cost | 1665 |
| Alternative 5 | |
|---|---|
| Error | 4.5 |
| Cost | 1025 |
| Alternative 6 | |
|---|---|
| Error | 4.5 |
| Cost | 897 |
| Alternative 7 | |
|---|---|
| Error | 27.0 |
| Cost | 641 |
| Alternative 8 | |
|---|---|
| Error | 40.0 |
| Cost | 64 |
| Alternative 9 | |
|---|---|
| Error | 61.6 |
| Cost | 64 |

if (/.f64 (-.f64 beta alpha) (+.f64 (+.f64 alpha beta) 2)) < -0.99999999999560463Initial program 60.4
rmApplied div-sub_binary64_247060.4
Applied associate-+l-_binary64_240058.4
Simplified58.4
rmApplied frac-2neg_binary64_247658.4
Simplified58.4
rmApplied neg-sub0_binary64_246058.4
Applied div-sub_binary64_247058.4
Applied associate--l-_binary64_240358.4
Simplified58.4
Taylor expanded around inf 3.2
Simplified0.0
Simplified0.0
if -0.99999999999560463 < (/.f64 (-.f64 beta alpha) (+.f64 (+.f64 alpha beta) 2)) Initial program 0.3
rmApplied div-sub_binary64_24700.3
Applied associate-+l-_binary64_24000.3
Simplified0.3
rmApplied add-exp-log_binary64_25030.3
Applied add-exp-log_binary64_25030.3
Applied div-exp_binary64_25160.3
Simplified0.3
Simplified0.3
Final simplification0.2
herbie shell --seed 2021014
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/1"
:precision binary64
:pre (and (> alpha -1.0) (> beta -1.0))
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0) 2.0))