| Alternative 1 | |
|---|---|
| Error | 0.1 |
| Cost | 2116 |
(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.9998)
(/
(-
(/ (+ beta (+ beta 2.0)) alpha)
(* (/ (/ (+ 2.0 (+ beta beta)) alpha) alpha) (+ beta 2.0)))
2.0)
(/
(-
(/ beta (+ beta (+ alpha 2.0)))
(log (exp (+ -1.0 (/ alpha (+ alpha (+ beta 2.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.9998) {
tmp = (((beta + (beta + 2.0)) / alpha) - ((((2.0 + (beta + beta)) / alpha) / alpha) * (beta + 2.0))) / 2.0;
} else {
tmp = ((beta / (beta + (alpha + 2.0))) - log(exp((-1.0 + (alpha / (alpha + (beta + 2.0))))))) / 2.0;
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
code = (((beta - alpha) / ((alpha + beta) + 2.0d0)) + 1.0d0) / 2.0d0
end function
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (((beta - alpha) / ((beta + alpha) + 2.0d0)) <= (-0.9998d0)) then
tmp = (((beta + (beta + 2.0d0)) / alpha) - ((((2.0d0 + (beta + beta)) / alpha) / alpha) * (beta + 2.0d0))) / 2.0d0
else
tmp = ((beta / (beta + (alpha + 2.0d0))) - log(exp(((-1.0d0) + (alpha / (alpha + (beta + 2.0d0))))))) / 2.0d0
end if
code = tmp
end function
public static double code(double alpha, double beta) {
return (((beta - alpha) / ((alpha + beta) + 2.0)) + 1.0) / 2.0;
}
public static double code(double alpha, double beta) {
double tmp;
if (((beta - alpha) / ((beta + alpha) + 2.0)) <= -0.9998) {
tmp = (((beta + (beta + 2.0)) / alpha) - ((((2.0 + (beta + beta)) / alpha) / alpha) * (beta + 2.0))) / 2.0;
} else {
tmp = ((beta / (beta + (alpha + 2.0))) - Math.log(Math.exp((-1.0 + (alpha / (alpha + (beta + 2.0))))))) / 2.0;
}
return tmp;
}
def code(alpha, beta): return (((beta - alpha) / ((alpha + beta) + 2.0)) + 1.0) / 2.0
def code(alpha, beta): tmp = 0 if ((beta - alpha) / ((beta + alpha) + 2.0)) <= -0.9998: tmp = (((beta + (beta + 2.0)) / alpha) - ((((2.0 + (beta + beta)) / alpha) / alpha) * (beta + 2.0))) / 2.0 else: tmp = ((beta / (beta + (alpha + 2.0))) - math.log(math.exp((-1.0 + (alpha / (alpha + (beta + 2.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) tmp = 0.0 if (Float64(Float64(beta - alpha) / Float64(Float64(beta + alpha) + 2.0)) <= -0.9998) tmp = Float64(Float64(Float64(Float64(beta + Float64(beta + 2.0)) / alpha) - Float64(Float64(Float64(Float64(2.0 + Float64(beta + beta)) / alpha) / alpha) * Float64(beta + 2.0))) / 2.0); else tmp = Float64(Float64(Float64(beta / Float64(beta + Float64(alpha + 2.0))) - log(exp(Float64(-1.0 + Float64(alpha / Float64(alpha + Float64(beta + 2.0))))))) / 2.0); end return tmp end
function tmp = code(alpha, beta) tmp = (((beta - alpha) / ((alpha + beta) + 2.0)) + 1.0) / 2.0; end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (((beta - alpha) / ((beta + alpha) + 2.0)) <= -0.9998) tmp = (((beta + (beta + 2.0)) / alpha) - ((((2.0 + (beta + beta)) / alpha) / alpha) * (beta + 2.0))) / 2.0; else tmp = ((beta / (beta + (alpha + 2.0))) - log(exp((-1.0 + (alpha / (alpha + (beta + 2.0))))))) / 2.0; end tmp_2 = 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_] := If[LessEqual[N[(N[(beta - alpha), $MachinePrecision] / N[(N[(beta + alpha), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision], -0.9998], N[(N[(N[(N[(beta + N[(beta + 2.0), $MachinePrecision]), $MachinePrecision] / alpha), $MachinePrecision] - N[(N[(N[(N[(2.0 + N[(beta + beta), $MachinePrecision]), $MachinePrecision] / alpha), $MachinePrecision] / alpha), $MachinePrecision] * N[(beta + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision], N[(N[(N[(beta / N[(beta + N[(alpha + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[Log[N[Exp[N[(-1.0 + N[(alpha / N[(alpha + N[(beta + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]]
\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.9998:\\
\;\;\;\;\frac{\frac{\beta + \left(\beta + 2\right)}{\alpha} - \frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(\beta + 2\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\beta}{\beta + \left(\alpha + 2\right)} - \log \left(e^{-1 + \frac{\alpha}{\alpha + \left(\beta + 2\right)}}\right)}{2}\\
\end{array}
Results
if (/.f64 (-.f64 beta alpha) (+.f64 (+.f64 alpha beta) 2)) < -0.99980000000000002Initial program 59.4
Simplified59.4
[Start]59.4 | \[ \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}
\] |
|---|---|
+-commutative [=>]59.4 | \[ \frac{\frac{\beta - \alpha}{\color{blue}{\left(\beta + \alpha\right)} + 2} + 1}{2}
\] |
Taylor expanded in alpha around -inf 3.0
Simplified3.0
[Start]3.0 | \[ \frac{-1 \cdot \frac{-1 \cdot \beta - \left(\beta + 2\right)}{\alpha} + -1 \cdot \frac{{\left(\beta + 2\right)}^{2} + \beta \cdot \left(\beta + 2\right)}{{\alpha}^{2}}}{2}
\] |
|---|---|
distribute-lft-out [=>]3.0 | \[ \frac{\color{blue}{-1 \cdot \left(\frac{-1 \cdot \beta - \left(\beta + 2\right)}{\alpha} + \frac{{\left(\beta + 2\right)}^{2} + \beta \cdot \left(\beta + 2\right)}{{\alpha}^{2}}\right)}}{2}
\] |
mul-1-neg [=>]3.0 | \[ \frac{-1 \cdot \left(\frac{\color{blue}{\left(-\beta\right)} - \left(\beta + 2\right)}{\alpha} + \frac{{\left(\beta + 2\right)}^{2} + \beta \cdot \left(\beta + 2\right)}{{\alpha}^{2}}\right)}{2}
\] |
+-commutative [=>]3.0 | \[ \frac{-1 \cdot \left(\frac{\left(-\beta\right) - \color{blue}{\left(2 + \beta\right)}}{\alpha} + \frac{{\left(\beta + 2\right)}^{2} + \beta \cdot \left(\beta + 2\right)}{{\alpha}^{2}}\right)}{2}
\] |
+-commutative [=>]3.0 | \[ \frac{-1 \cdot \left(\frac{\left(-\beta\right) - \left(2 + \beta\right)}{\alpha} + \frac{{\color{blue}{\left(2 + \beta\right)}}^{2} + \beta \cdot \left(\beta + 2\right)}{{\alpha}^{2}}\right)}{2}
\] |
+-commutative [=>]3.0 | \[ \frac{-1 \cdot \left(\frac{\left(-\beta\right) - \left(2 + \beta\right)}{\alpha} + \frac{{\left(2 + \beta\right)}^{2} + \beta \cdot \color{blue}{\left(2 + \beta\right)}}{{\alpha}^{2}}\right)}{2}
\] |
unpow2 [=>]3.0 | \[ \frac{-1 \cdot \left(\frac{\left(-\beta\right) - \left(2 + \beta\right)}{\alpha} + \frac{{\left(2 + \beta\right)}^{2} + \beta \cdot \left(2 + \beta\right)}{\color{blue}{\alpha \cdot \alpha}}\right)}{2}
\] |
Applied egg-rr45.9
Taylor expanded in alpha around 0 3.0
Simplified0.2
[Start]3.0 | \[ \frac{-1 \cdot \left(-1 \cdot \frac{2 + 2 \cdot \beta}{\alpha} + \frac{\left(2 + 2 \cdot \beta\right) \cdot \left(\beta + 2\right)}{{\alpha}^{2}}\right)}{2}
\] |
|---|---|
+-commutative [=>]3.0 | \[ \frac{-1 \cdot \color{blue}{\left(\frac{\left(2 + 2 \cdot \beta\right) \cdot \left(\beta + 2\right)}{{\alpha}^{2}} + -1 \cdot \frac{2 + 2 \cdot \beta}{\alpha}\right)}}{2}
\] |
associate-/l* [=>]0.2 | \[ \frac{-1 \cdot \left(\color{blue}{\frac{2 + 2 \cdot \beta}{\frac{{\alpha}^{2}}{\beta + 2}}} + -1 \cdot \frac{2 + 2 \cdot \beta}{\alpha}\right)}{2}
\] |
count-2 [<=]0.2 | \[ \frac{-1 \cdot \left(\frac{2 + \color{blue}{\left(\beta + \beta\right)}}{\frac{{\alpha}^{2}}{\beta + 2}} + -1 \cdot \frac{2 + 2 \cdot \beta}{\alpha}\right)}{2}
\] |
unpow2 [=>]0.2 | \[ \frac{-1 \cdot \left(\frac{2 + \left(\beta + \beta\right)}{\frac{\color{blue}{\alpha \cdot \alpha}}{\beta + 2}} + -1 \cdot \frac{2 + 2 \cdot \beta}{\alpha}\right)}{2}
\] |
associate-/r/ [=>]0.2 | \[ \frac{-1 \cdot \left(\color{blue}{\frac{2 + \left(\beta + \beta\right)}{\alpha \cdot \alpha} \cdot \left(\beta + 2\right)} + -1 \cdot \frac{2 + 2 \cdot \beta}{\alpha}\right)}{2}
\] |
associate-/r* [=>]0.2 | \[ \frac{-1 \cdot \left(\color{blue}{\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha}} \cdot \left(\beta + 2\right) + -1 \cdot \frac{2 + 2 \cdot \beta}{\alpha}\right)}{2}
\] |
+-commutative [=>]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \color{blue}{\left(2 + \beta\right)} + -1 \cdot \frac{2 + 2 \cdot \beta}{\alpha}\right)}{2}
\] |
associate-*r/ [=>]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \color{blue}{\frac{-1 \cdot \left(2 + 2 \cdot \beta\right)}{\alpha}}\right)}{2}
\] |
count-2 [<=]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \frac{-1 \cdot \left(2 + \color{blue}{\left(\beta + \beta\right)}\right)}{\alpha}\right)}{2}
\] |
neg-mul-1 [<=]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \frac{\color{blue}{-\left(2 + \left(\beta + \beta\right)\right)}}{\alpha}\right)}{2}
\] |
neg-sub0 [=>]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \frac{\color{blue}{0 - \left(2 + \left(\beta + \beta\right)\right)}}{\alpha}\right)}{2}
\] |
associate-+r+ [=>]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \frac{0 - \color{blue}{\left(\left(2 + \beta\right) + \beta\right)}}{\alpha}\right)}{2}
\] |
+-commutative [<=]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \frac{0 - \left(\color{blue}{\left(\beta + 2\right)} + \beta\right)}{\alpha}\right)}{2}
\] |
associate--r+ [=>]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \frac{\color{blue}{\left(0 - \left(\beta + 2\right)\right) - \beta}}{\alpha}\right)}{2}
\] |
neg-sub0 [<=]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \frac{\color{blue}{\left(-\left(\beta + 2\right)\right)} - \beta}{\alpha}\right)}{2}
\] |
+-commutative [=>]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \frac{\left(-\color{blue}{\left(2 + \beta\right)}\right) - \beta}{\alpha}\right)}{2}
\] |
distribute-neg-in [=>]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \frac{\color{blue}{\left(\left(-2\right) + \left(-\beta\right)\right)} - \beta}{\alpha}\right)}{2}
\] |
metadata-eval [=>]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \frac{\left(\color{blue}{-2} + \left(-\beta\right)\right) - \beta}{\alpha}\right)}{2}
\] |
sub-neg [<=]0.2 | \[ \frac{-1 \cdot \left(\frac{\frac{2 + \left(\beta + \beta\right)}{\alpha}}{\alpha} \cdot \left(2 + \beta\right) + \frac{\color{blue}{\left(-2 - \beta\right)} - \beta}{\alpha}\right)}{2}
\] |
if -0.99980000000000002 < (/.f64 (-.f64 beta alpha) (+.f64 (+.f64 alpha beta) 2)) Initial program 0.0
Simplified0.0
[Start]0.0 | \[ \frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}
\] |
|---|---|
+-commutative [=>]0.0 | \[ \frac{\frac{\beta - \alpha}{\color{blue}{\left(\beta + \alpha\right)} + 2} + 1}{2}
\] |
Applied egg-rr0.0
Applied egg-rr0.0
Final simplification0.1
| Alternative 1 | |
|---|---|
| Error | 0.1 |
| Cost | 2116 |
| Alternative 2 | |
|---|---|
| Error | 0.1 |
| Cost | 1860 |
| Alternative 3 | |
|---|---|
| Error | 0.2 |
| Cost | 1604 |
| Alternative 4 | |
|---|---|
| Error | 0.1 |
| Cost | 1604 |
| Alternative 5 | |
|---|---|
| Error | 0.2 |
| Cost | 1476 |
| Alternative 6 | |
|---|---|
| Error | 4.2 |
| Cost | 964 |
| Alternative 7 | |
|---|---|
| Error | 7.5 |
| Cost | 708 |
| Alternative 8 | |
|---|---|
| Error | 4.5 |
| Cost | 708 |
| Alternative 9 | |
|---|---|
| Error | 18.2 |
| Cost | 580 |
| Alternative 10 | |
|---|---|
| Error | 18.0 |
| Cost | 580 |
| Alternative 11 | |
|---|---|
| Error | 18.4 |
| Cost | 196 |
| Alternative 12 | |
|---|---|
| Error | 32.7 |
| Cost | 64 |
herbie shell --seed 2023053
(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))