
(FPCore (alpha beta) :precision binary64 (/ (+ (/ (- beta 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;
}
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
public static double code(double alpha, double beta) {
return (((beta - alpha) / ((alpha + beta) + 2.0)) + 1.0) / 2.0;
}
def code(alpha, beta): return (((beta - alpha) / ((alpha + beta) + 2.0)) + 1.0) / 2.0
function code(alpha, beta) return Float64(Float64(Float64(Float64(beta - alpha) / Float64(Float64(alpha + beta) + 2.0)) + 1.0) / 2.0) end
function tmp = code(alpha, beta) tmp = (((beta - alpha) / ((alpha + beta) + 2.0)) + 1.0) / 2.0; 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]
\begin{array}{l}
\\
\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (alpha beta) :precision binary64 (/ (+ (/ (- beta 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;
}
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
public static double code(double alpha, double beta) {
return (((beta - alpha) / ((alpha + beta) + 2.0)) + 1.0) / 2.0;
}
def code(alpha, beta): return (((beta - alpha) / ((alpha + beta) + 2.0)) + 1.0) / 2.0
function code(alpha, beta) return Float64(Float64(Float64(Float64(beta - alpha) / Float64(Float64(alpha + beta) + 2.0)) + 1.0) / 2.0) end
function tmp = code(alpha, beta) tmp = (((beta - alpha) / ((alpha + beta) + 2.0)) + 1.0) / 2.0; 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]
\begin{array}{l}
\\
\frac{\frac{\beta - \alpha}{\left(\alpha + \beta\right) + 2} + 1}{2}
\end{array}
(FPCore (alpha beta) :precision binary64 (if (<= (/ (- beta alpha) (+ (+ beta alpha) 2.0)) -0.8) (/ (+ beta 1.0) alpha) (+ 0.5 (* (- alpha beta) (/ -0.5 (+ beta (+ alpha 2.0)))))))
double code(double alpha, double beta) {
double tmp;
if (((beta - alpha) / ((beta + alpha) + 2.0)) <= -0.8) {
tmp = (beta + 1.0) / alpha;
} else {
tmp = 0.5 + ((alpha - beta) * (-0.5 / (beta + (alpha + 2.0))));
}
return tmp;
}
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.8d0)) then
tmp = (beta + 1.0d0) / alpha
else
tmp = 0.5d0 + ((alpha - beta) * ((-0.5d0) / (beta + (alpha + 2.0d0))))
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (((beta - alpha) / ((beta + alpha) + 2.0)) <= -0.8) {
tmp = (beta + 1.0) / alpha;
} else {
tmp = 0.5 + ((alpha - beta) * (-0.5 / (beta + (alpha + 2.0))));
}
return tmp;
}
def code(alpha, beta): tmp = 0 if ((beta - alpha) / ((beta + alpha) + 2.0)) <= -0.8: tmp = (beta + 1.0) / alpha else: tmp = 0.5 + ((alpha - beta) * (-0.5 / (beta + (alpha + 2.0)))) return tmp
function code(alpha, beta) tmp = 0.0 if (Float64(Float64(beta - alpha) / Float64(Float64(beta + alpha) + 2.0)) <= -0.8) tmp = Float64(Float64(beta + 1.0) / alpha); else tmp = Float64(0.5 + Float64(Float64(alpha - beta) * Float64(-0.5 / Float64(beta + Float64(alpha + 2.0))))); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (((beta - alpha) / ((beta + alpha) + 2.0)) <= -0.8) tmp = (beta + 1.0) / alpha; else tmp = 0.5 + ((alpha - beta) * (-0.5 / (beta + (alpha + 2.0)))); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[N[(N[(beta - alpha), $MachinePrecision] / N[(N[(beta + alpha), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision], -0.8], N[(N[(beta + 1.0), $MachinePrecision] / alpha), $MachinePrecision], N[(0.5 + N[(N[(alpha - beta), $MachinePrecision] * N[(-0.5 / N[(beta + N[(alpha + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{\beta - \alpha}{\left(\beta + \alpha\right) + 2} \leq -0.8:\\
\;\;\;\;\frac{\beta + 1}{\alpha}\\
\mathbf{else}:\\
\;\;\;\;0.5 + \left(\alpha - \beta\right) \cdot \frac{-0.5}{\beta + \left(\alpha + 2\right)}\\
\end{array}
\end{array}
if (/.f64 (-.f64 beta alpha) (+.f64 (+.f64 alpha beta) #s(literal 2 binary64))) < -0.80000000000000004Initial program 6.7%
+-commutative6.7%
Simplified6.7%
Taylor expanded in alpha around inf 99.3%
*-commutative99.3%
Simplified99.3%
Taylor expanded in beta around 0 99.3%
+-commutative99.3%
*-lft-identity99.3%
associate-*l/99.2%
*-rgt-identity99.2%
distribute-lft-in99.2%
associate-*l/99.3%
*-lft-identity99.3%
Simplified99.3%
if -0.80000000000000004 < (/.f64 (-.f64 beta alpha) (+.f64 (+.f64 alpha beta) #s(literal 2 binary64))) Initial program 100.0%
+-commutative100.0%
sub-neg100.0%
+-commutative100.0%
neg-sub0100.0%
associate-+l-100.0%
sub0-neg100.0%
distribute-frac-neg100.0%
+-commutative100.0%
sub-neg100.0%
div-sub100.0%
sub-neg100.0%
metadata-eval100.0%
neg-mul-1100.0%
*-commutative100.0%
+-commutative100.0%
associate-/l/100.0%
associate-*l/100.0%
Simplified100.0%
Final simplification99.9%
(FPCore (alpha beta) :precision binary64 (if (<= alpha 2.9e+17) (+ 0.5 (* (- alpha beta) (/ 0.5 (- -2.0 beta)))) (/ (+ beta 1.0) alpha)))
double code(double alpha, double beta) {
double tmp;
if (alpha <= 2.9e+17) {
tmp = 0.5 + ((alpha - beta) * (0.5 / (-2.0 - beta)));
} else {
tmp = (beta + 1.0) / alpha;
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (alpha <= 2.9d+17) then
tmp = 0.5d0 + ((alpha - beta) * (0.5d0 / ((-2.0d0) - beta)))
else
tmp = (beta + 1.0d0) / alpha
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (alpha <= 2.9e+17) {
tmp = 0.5 + ((alpha - beta) * (0.5 / (-2.0 - beta)));
} else {
tmp = (beta + 1.0) / alpha;
}
return tmp;
}
def code(alpha, beta): tmp = 0 if alpha <= 2.9e+17: tmp = 0.5 + ((alpha - beta) * (0.5 / (-2.0 - beta))) else: tmp = (beta + 1.0) / alpha return tmp
function code(alpha, beta) tmp = 0.0 if (alpha <= 2.9e+17) tmp = Float64(0.5 + Float64(Float64(alpha - beta) * Float64(0.5 / Float64(-2.0 - beta)))); else tmp = Float64(Float64(beta + 1.0) / alpha); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (alpha <= 2.9e+17) tmp = 0.5 + ((alpha - beta) * (0.5 / (-2.0 - beta))); else tmp = (beta + 1.0) / alpha; end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[alpha, 2.9e+17], N[(0.5 + N[(N[(alpha - beta), $MachinePrecision] * N[(0.5 / N[(-2.0 - beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(beta + 1.0), $MachinePrecision] / alpha), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\alpha \leq 2.9 \cdot 10^{+17}:\\
\;\;\;\;0.5 + \left(\alpha - \beta\right) \cdot \frac{0.5}{-2 - \beta}\\
\mathbf{else}:\\
\;\;\;\;\frac{\beta + 1}{\alpha}\\
\end{array}
\end{array}
if alpha < 2.9e17Initial program 99.8%
+-commutative99.8%
sub-neg99.8%
+-commutative99.8%
neg-sub099.8%
associate-+l-99.8%
sub0-neg99.8%
distribute-frac-neg99.8%
+-commutative99.8%
sub-neg99.8%
div-sub99.8%
sub-neg99.8%
metadata-eval99.8%
neg-mul-199.8%
*-commutative99.8%
+-commutative99.8%
associate-/l/99.8%
associate-*l/99.8%
Simplified99.8%
Taylor expanded in alpha around 0 98.6%
metadata-eval98.6%
distribute-neg-frac98.6%
distribute-neg-frac298.6%
distribute-neg-in98.6%
metadata-eval98.6%
unsub-neg98.6%
Simplified98.6%
if 2.9e17 < alpha Initial program 20.7%
+-commutative20.7%
Simplified20.7%
Taylor expanded in alpha around inf 84.9%
*-commutative84.9%
Simplified84.9%
Taylor expanded in beta around 0 85.0%
+-commutative85.0%
*-lft-identity85.0%
associate-*l/84.9%
*-rgt-identity84.9%
distribute-lft-in84.9%
associate-*l/85.0%
*-lft-identity85.0%
Simplified85.0%
Final simplification95.3%
(FPCore (alpha beta) :precision binary64 (if (<= alpha 1.92e+62) (/ (+ 1.0 (/ beta (+ beta 2.0))) 2.0) (/ (+ beta 1.0) alpha)))
double code(double alpha, double beta) {
double tmp;
if (alpha <= 1.92e+62) {
tmp = (1.0 + (beta / (beta + 2.0))) / 2.0;
} else {
tmp = (beta + 1.0) / alpha;
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (alpha <= 1.92d+62) then
tmp = (1.0d0 + (beta / (beta + 2.0d0))) / 2.0d0
else
tmp = (beta + 1.0d0) / alpha
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (alpha <= 1.92e+62) {
tmp = (1.0 + (beta / (beta + 2.0))) / 2.0;
} else {
tmp = (beta + 1.0) / alpha;
}
return tmp;
}
def code(alpha, beta): tmp = 0 if alpha <= 1.92e+62: tmp = (1.0 + (beta / (beta + 2.0))) / 2.0 else: tmp = (beta + 1.0) / alpha return tmp
function code(alpha, beta) tmp = 0.0 if (alpha <= 1.92e+62) tmp = Float64(Float64(1.0 + Float64(beta / Float64(beta + 2.0))) / 2.0); else tmp = Float64(Float64(beta + 1.0) / alpha); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (alpha <= 1.92e+62) tmp = (1.0 + (beta / (beta + 2.0))) / 2.0; else tmp = (beta + 1.0) / alpha; end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[alpha, 1.92e+62], N[(N[(1.0 + N[(beta / N[(beta + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision], N[(N[(beta + 1.0), $MachinePrecision] / alpha), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\alpha \leq 1.92 \cdot 10^{+62}:\\
\;\;\;\;\frac{1 + \frac{\beta}{\beta + 2}}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\beta + 1}{\alpha}\\
\end{array}
\end{array}
if alpha < 1.9199999999999999e62Initial program 97.4%
+-commutative97.4%
Simplified97.4%
Taylor expanded in alpha around 0 95.9%
+-commutative95.9%
+-commutative95.9%
Simplified95.9%
if 1.9199999999999999e62 < alpha Initial program 14.8%
+-commutative14.8%
Simplified14.8%
Taylor expanded in alpha around inf 91.1%
*-commutative91.1%
Simplified91.1%
Taylor expanded in beta around 0 91.1%
+-commutative91.1%
*-lft-identity91.1%
associate-*l/91.0%
*-rgt-identity91.0%
distribute-lft-in91.0%
associate-*l/91.1%
*-lft-identity91.1%
Simplified91.1%
Final simplification94.9%
(FPCore (alpha beta) :precision binary64 (if (<= beta 2.0) (/ (+ 1.0 (* beta 0.5)) 2.0) 1.0))
double code(double alpha, double beta) {
double tmp;
if (beta <= 2.0) {
tmp = (1.0 + (beta * 0.5)) / 2.0;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 2.0d0) then
tmp = (1.0d0 + (beta * 0.5d0)) / 2.0d0
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 2.0) {
tmp = (1.0 + (beta * 0.5)) / 2.0;
} else {
tmp = 1.0;
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 2.0: tmp = (1.0 + (beta * 0.5)) / 2.0 else: tmp = 1.0 return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 2.0) tmp = Float64(Float64(1.0 + Float64(beta * 0.5)) / 2.0); else tmp = 1.0; end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 2.0) tmp = (1.0 + (beta * 0.5)) / 2.0; else tmp = 1.0; end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 2.0], N[(N[(1.0 + N[(beta * 0.5), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2:\\
\;\;\;\;\frac{1 + \beta \cdot 0.5}{2}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if beta < 2Initial program 76.3%
+-commutative76.3%
Simplified76.3%
Taylor expanded in alpha around 0 74.1%
+-commutative74.1%
+-commutative74.1%
Simplified74.1%
Taylor expanded in beta around 0 73.4%
if 2 < beta Initial program 89.2%
+-commutative89.2%
sub-neg89.2%
+-commutative89.2%
neg-sub089.2%
associate-+l-89.2%
sub0-neg89.2%
distribute-frac-neg89.2%
+-commutative89.2%
sub-neg89.2%
div-sub89.2%
sub-neg89.2%
metadata-eval89.2%
neg-mul-189.2%
*-commutative89.2%
+-commutative89.2%
associate-/l/89.2%
associate-*l/89.2%
Simplified89.2%
Taylor expanded in beta around inf 88.9%
Final simplification78.3%
(FPCore (alpha beta) :precision binary64 (if (<= beta 2.1) 0.5 1.0))
double code(double alpha, double beta) {
double tmp;
if (beta <= 2.1) {
tmp = 0.5;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 2.1d0) then
tmp = 0.5d0
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 2.1) {
tmp = 0.5;
} else {
tmp = 1.0;
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 2.1: tmp = 0.5 else: tmp = 1.0 return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 2.1) tmp = 0.5; else tmp = 1.0; end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 2.1) tmp = 0.5; else tmp = 1.0; end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 2.1], 0.5, 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2.1:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if beta < 2.10000000000000009Initial program 76.3%
+-commutative76.3%
Simplified76.3%
Taylor expanded in beta around 0 75.2%
Taylor expanded in alpha around 0 73.1%
if 2.10000000000000009 < beta Initial program 89.2%
+-commutative89.2%
sub-neg89.2%
+-commutative89.2%
neg-sub089.2%
associate-+l-89.2%
sub0-neg89.2%
distribute-frac-neg89.2%
+-commutative89.2%
sub-neg89.2%
div-sub89.2%
sub-neg89.2%
metadata-eval89.2%
neg-mul-189.2%
*-commutative89.2%
+-commutative89.2%
associate-/l/89.2%
associate-*l/89.2%
Simplified89.2%
Taylor expanded in beta around inf 88.9%
Final simplification78.0%
(FPCore (alpha beta) :precision binary64 0.5)
double code(double alpha, double beta) {
return 0.5;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
code = 0.5d0
end function
public static double code(double alpha, double beta) {
return 0.5;
}
def code(alpha, beta): return 0.5
function code(alpha, beta) return 0.5 end
function tmp = code(alpha, beta) tmp = 0.5; end
code[alpha_, beta_] := 0.5
\begin{array}{l}
\\
0.5
\end{array}
Initial program 80.3%
+-commutative80.3%
Simplified80.3%
Taylor expanded in beta around 0 56.3%
Taylor expanded in alpha around 0 55.7%
Final simplification55.7%
herbie shell --seed 2024076
(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))