\[\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}
\]
↓
\[\begin{array}{l}
t_0 := \frac{\beta + 2}{\alpha}\\
\mathbf{if}\;\frac{\beta - \alpha}{\left(\beta + \alpha\right) + 2} \leq -0.9996:\\
\;\;\;\;\frac{{t_0}^{3} + \left(\frac{\beta}{{\alpha}^{3}} \cdot {\left(\beta + 2\right)}^{2} + \mathsf{fma}\left(\frac{\left(-2 - \beta\right) - \beta}{\alpha}, t_0, \frac{\beta + \left(\beta + 2\right)}{\alpha}\right)\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\beta}{\beta + \left(\alpha + 2\right)} - \log \left(e^{\frac{\alpha}{\alpha + \left(\beta + 2\right)} + -1}\right)}{2}\\
\end{array}
\]
(FPCore (alpha beta)
:precision binary64
(/ (+ (/ (- beta alpha) (+ (+ alpha beta) 2.0)) 1.0) 2.0))
↓
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (/ (+ beta 2.0) alpha)))
(if (<= (/ (- beta alpha) (+ (+ beta alpha) 2.0)) -0.9996)
(/
(+
(pow t_0 3.0)
(+
(* (/ beta (pow alpha 3.0)) (pow (+ beta 2.0) 2.0))
(fma
(/ (- (- -2.0 beta) beta) alpha)
t_0
(/ (+ beta (+ beta 2.0)) alpha))))
2.0)
(/
(-
(/ beta (+ beta (+ alpha 2.0)))
(log (exp (+ (/ alpha (+ alpha (+ beta 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 t_0 = (beta + 2.0) / alpha;
double tmp;
if (((beta - alpha) / ((beta + alpha) + 2.0)) <= -0.9996) {
tmp = (pow(t_0, 3.0) + (((beta / pow(alpha, 3.0)) * pow((beta + 2.0), 2.0)) + fma((((-2.0 - beta) - beta) / alpha), t_0, ((beta + (beta + 2.0)) / alpha)))) / 2.0;
} else {
tmp = ((beta / (beta + (alpha + 2.0))) - log(exp(((alpha / (alpha + (beta + 2.0))) + -1.0)))) / 2.0;
}
return tmp;
}
function code(alpha, beta)
return Float64(Float64(Float64(Float64(beta - alpha) / Float64(Float64(alpha + beta) + 2.0)) + 1.0) / 2.0)
end
↓
function code(alpha, beta)
t_0 = Float64(Float64(beta + 2.0) / alpha)
tmp = 0.0
if (Float64(Float64(beta - alpha) / Float64(Float64(beta + alpha) + 2.0)) <= -0.9996)
tmp = Float64(Float64((t_0 ^ 3.0) + Float64(Float64(Float64(beta / (alpha ^ 3.0)) * (Float64(beta + 2.0) ^ 2.0)) + fma(Float64(Float64(Float64(-2.0 - beta) - beta) / alpha), t_0, Float64(Float64(beta + Float64(beta + 2.0)) / alpha)))) / 2.0);
else
tmp = Float64(Float64(Float64(beta / Float64(beta + Float64(alpha + 2.0))) - log(exp(Float64(Float64(alpha / Float64(alpha + Float64(beta + 2.0))) + -1.0)))) / 2.0);
end
return tmp
end
code[alpha_, beta_] := N[(N[(N[(N[(beta - alpha), $MachinePrecision] / N[(N[(alpha + beta), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / 2.0), $MachinePrecision]
↓
code[alpha_, beta_] := Block[{t$95$0 = N[(N[(beta + 2.0), $MachinePrecision] / alpha), $MachinePrecision]}, If[LessEqual[N[(N[(beta - alpha), $MachinePrecision] / N[(N[(beta + alpha), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision], -0.9996], N[(N[(N[Power[t$95$0, 3.0], $MachinePrecision] + N[(N[(N[(beta / N[Power[alpha, 3.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[(beta + 2.0), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] + N[(N[(N[(N[(-2.0 - beta), $MachinePrecision] - beta), $MachinePrecision] / alpha), $MachinePrecision] * t$95$0 + N[(N[(beta + N[(beta + 2.0), $MachinePrecision]), $MachinePrecision] / alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision], N[(N[(N[(beta / N[(beta + N[(alpha + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[Log[N[Exp[N[(N[(alpha / N[(alpha + N[(beta + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]]]
\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}
↓
\begin{array}{l}
t_0 := \frac{\beta + 2}{\alpha}\\
\mathbf{if}\;\frac{\beta - \alpha}{\left(\beta + \alpha\right) + 2} \leq -0.9996:\\
\;\;\;\;\frac{{t_0}^{3} + \left(\frac{\beta}{{\alpha}^{3}} \cdot {\left(\beta + 2\right)}^{2} + \mathsf{fma}\left(\frac{\left(-2 - \beta\right) - \beta}{\alpha}, t_0, \frac{\beta + \left(\beta + 2\right)}{\alpha}\right)\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\beta}{\beta + \left(\alpha + 2\right)} - \log \left(e^{\frac{\alpha}{\alpha + \left(\beta + 2\right)} + -1}\right)}{2}\\
\end{array}