
(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 7 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.999999995) (+ (/ 1.0 alpha) (/ beta 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.999999995) {
tmp = (1.0 / alpha) + (beta / 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.999999995d0)) then
tmp = (1.0d0 / alpha) + (beta / 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.999999995) {
tmp = (1.0 / alpha) + (beta / 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.999999995: tmp = (1.0 / alpha) + (beta / 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.999999995) tmp = Float64(Float64(1.0 / alpha) + Float64(beta / 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.999999995) tmp = (1.0 / alpha) + (beta / 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.999999995], N[(N[(1.0 / alpha), $MachinePrecision] + N[(beta / alpha), $MachinePrecision]), $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.999999995:\\
\;\;\;\;\frac{1}{\alpha} + \frac{\beta}{\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.99999999500000003Initial program 6.1%
+-commutative6.1%
sub-neg6.1%
+-commutative6.1%
neg-sub06.1%
associate-+l-6.1%
sub0-neg6.1%
distribute-frac-neg6.1%
+-commutative6.1%
sub-neg6.1%
div-sub6.1%
sub-neg6.1%
metadata-eval6.1%
neg-mul-16.1%
*-commutative6.1%
+-commutative6.1%
associate-/l/6.1%
associate-*l/6.1%
Simplified6.1%
Taylor expanded in alpha around inf 99.5%
neg-mul-199.5%
associate--r+99.5%
sub-neg99.5%
distribute-neg-in99.5%
+-commutative99.5%
distribute-neg-in99.5%
metadata-eval99.5%
unsub-neg99.5%
Simplified99.5%
Taylor expanded in beta around 0 99.6%
if -0.99999999500000003 < (/.f64 (-.f64 beta alpha) (+.f64 (+.f64 alpha beta) #s(literal 2 binary64))) Initial program 99.7%
+-commutative99.7%
sub-neg99.7%
+-commutative99.7%
neg-sub099.7%
associate-+l-99.7%
sub0-neg99.7%
distribute-frac-neg99.7%
+-commutative99.7%
sub-neg99.7%
div-sub99.7%
sub-neg99.7%
metadata-eval99.7%
neg-mul-199.7%
*-commutative99.7%
+-commutative99.7%
associate-/l/99.7%
associate-*l/99.7%
Simplified99.7%
Final simplification99.7%
(FPCore (alpha beta) :precision binary64 (if (<= alpha 86000000.0) (/ (+ 1.0 (/ beta (+ beta 2.0))) 2.0) (+ (/ 1.0 alpha) (/ beta alpha))))
double code(double alpha, double beta) {
double tmp;
if (alpha <= 86000000.0) {
tmp = (1.0 + (beta / (beta + 2.0))) / 2.0;
} else {
tmp = (1.0 / alpha) + (beta / alpha);
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (alpha <= 86000000.0d0) then
tmp = (1.0d0 + (beta / (beta + 2.0d0))) / 2.0d0
else
tmp = (1.0d0 / alpha) + (beta / alpha)
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (alpha <= 86000000.0) {
tmp = (1.0 + (beta / (beta + 2.0))) / 2.0;
} else {
tmp = (1.0 / alpha) + (beta / alpha);
}
return tmp;
}
def code(alpha, beta): tmp = 0 if alpha <= 86000000.0: tmp = (1.0 + (beta / (beta + 2.0))) / 2.0 else: tmp = (1.0 / alpha) + (beta / alpha) return tmp
function code(alpha, beta) tmp = 0.0 if (alpha <= 86000000.0) tmp = Float64(Float64(1.0 + Float64(beta / Float64(beta + 2.0))) / 2.0); else tmp = Float64(Float64(1.0 / alpha) + Float64(beta / alpha)); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (alpha <= 86000000.0) tmp = (1.0 + (beta / (beta + 2.0))) / 2.0; else tmp = (1.0 / alpha) + (beta / alpha); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[alpha, 86000000.0], N[(N[(1.0 + N[(beta / N[(beta + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision], N[(N[(1.0 / alpha), $MachinePrecision] + N[(beta / alpha), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\alpha \leq 86000000:\\
\;\;\;\;\frac{1 + \frac{\beta}{\beta + 2}}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\alpha} + \frac{\beta}{\alpha}\\
\end{array}
\end{array}
if alpha < 8.6e7Initial program 100.0%
+-commutative100.0%
Simplified100.0%
Taylor expanded in alpha around 0 99.1%
+-commutative99.1%
Simplified99.1%
if 8.6e7 < alpha Initial program 25.6%
+-commutative25.6%
sub-neg25.6%
+-commutative25.6%
neg-sub025.6%
associate-+l-25.6%
sub0-neg25.6%
distribute-frac-neg25.6%
+-commutative25.6%
sub-neg25.6%
div-sub25.6%
sub-neg25.6%
metadata-eval25.6%
neg-mul-125.6%
*-commutative25.6%
+-commutative25.6%
associate-/l/25.6%
associate-*l/25.6%
Simplified25.7%
Taylor expanded in alpha around inf 80.5%
neg-mul-180.5%
associate--r+80.5%
sub-neg80.5%
distribute-neg-in80.5%
+-commutative80.5%
distribute-neg-in80.5%
metadata-eval80.5%
unsub-neg80.5%
Simplified80.5%
Taylor expanded in beta around 0 80.6%
Final simplification93.7%
(FPCore (alpha beta) :precision binary64 (if (<= beta 2.0) (+ 0.5 (* beta 0.25)) (+ 1.0 (/ -1.0 beta))))
double code(double alpha, double beta) {
double tmp;
if (beta <= 2.0) {
tmp = 0.5 + (beta * 0.25);
} else {
tmp = 1.0 + (-1.0 / beta);
}
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 = 0.5d0 + (beta * 0.25d0)
else
tmp = 1.0d0 + ((-1.0d0) / beta)
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 2.0) {
tmp = 0.5 + (beta * 0.25);
} else {
tmp = 1.0 + (-1.0 / beta);
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 2.0: tmp = 0.5 + (beta * 0.25) else: tmp = 1.0 + (-1.0 / beta) return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 2.0) tmp = Float64(0.5 + Float64(beta * 0.25)); else tmp = Float64(1.0 + Float64(-1.0 / beta)); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 2.0) tmp = 0.5 + (beta * 0.25); else tmp = 1.0 + (-1.0 / beta); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 2.0], N[(0.5 + N[(beta * 0.25), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-1.0 / beta), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2:\\
\;\;\;\;0.5 + \beta \cdot 0.25\\
\mathbf{else}:\\
\;\;\;\;1 + \frac{-1}{\beta}\\
\end{array}
\end{array}
if beta < 2Initial program 73.1%
+-commutative73.1%
Simplified73.1%
Taylor expanded in alpha around 0 71.6%
+-commutative71.6%
Simplified71.6%
Taylor expanded in beta around 0 71.6%
*-commutative71.6%
Simplified71.6%
if 2 < beta Initial program 87.8%
+-commutative87.8%
Simplified87.8%
Taylor expanded in alpha around 0 86.3%
+-commutative86.3%
Simplified86.3%
Taylor expanded in beta around inf 86.1%
Final simplification76.9%
(FPCore (alpha beta) :precision binary64 (if (<= beta 2.0) (+ 0.5 (* beta 0.25)) 1.0))
double code(double alpha, double beta) {
double tmp;
if (beta <= 2.0) {
tmp = 0.5 + (beta * 0.25);
} 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 = 0.5d0 + (beta * 0.25d0)
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 = 0.5 + (beta * 0.25);
} else {
tmp = 1.0;
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 2.0: tmp = 0.5 + (beta * 0.25) else: tmp = 1.0 return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 2.0) tmp = Float64(0.5 + Float64(beta * 0.25)); else tmp = 1.0; end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 2.0) tmp = 0.5 + (beta * 0.25); else tmp = 1.0; end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 2.0], N[(0.5 + N[(beta * 0.25), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2:\\
\;\;\;\;0.5 + \beta \cdot 0.25\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if beta < 2Initial program 73.1%
+-commutative73.1%
Simplified73.1%
Taylor expanded in alpha around 0 71.6%
+-commutative71.6%
Simplified71.6%
Taylor expanded in beta around 0 71.6%
*-commutative71.6%
Simplified71.6%
if 2 < beta Initial program 87.8%
+-commutative87.8%
Simplified87.8%
Taylor expanded in alpha around 0 86.3%
+-commutative86.3%
Simplified86.3%
Taylor expanded in beta around inf 85.0%
(FPCore (alpha beta) :precision binary64 (if (<= beta 2.0) 0.5 1.0))
double code(double alpha, double beta) {
double tmp;
if (beta <= 2.0) {
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.0d0) 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.0) {
tmp = 0.5;
} else {
tmp = 1.0;
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 2.0: tmp = 0.5 else: tmp = 1.0 return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 2.0) tmp = 0.5; else tmp = 1.0; end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 2.0) tmp = 0.5; else tmp = 1.0; end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 2.0], 0.5, 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if beta < 2Initial program 73.1%
+-commutative73.1%
Simplified73.1%
Taylor expanded in alpha around 0 71.6%
+-commutative71.6%
Simplified71.6%
Taylor expanded in beta around 0 71.6%
*-commutative71.6%
Simplified71.6%
Taylor expanded in beta around 0 71.4%
if 2 < beta Initial program 87.8%
+-commutative87.8%
Simplified87.8%
Taylor expanded in alpha around 0 86.3%
+-commutative86.3%
Simplified86.3%
Taylor expanded in beta around inf 85.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 78.5%
+-commutative78.5%
Simplified78.5%
Taylor expanded in alpha around 0 77.0%
+-commutative77.0%
Simplified77.0%
Taylor expanded in beta around 0 47.2%
*-commutative47.2%
Simplified47.2%
Taylor expanded in beta around 0 51.4%
(FPCore (alpha beta) :precision binary64 0.0)
double code(double alpha, double beta) {
return 0.0;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
code = 0.0d0
end function
public static double code(double alpha, double beta) {
return 0.0;
}
def code(alpha, beta): return 0.0
function code(alpha, beta) return 0.0 end
function tmp = code(alpha, beta) tmp = 0.0; end
code[alpha_, beta_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 78.5%
+-commutative78.5%
sub-neg78.5%
+-commutative78.5%
neg-sub078.5%
associate-+l-78.5%
sub0-neg78.5%
distribute-frac-neg78.5%
+-commutative78.5%
sub-neg78.5%
div-sub78.5%
sub-neg78.5%
metadata-eval78.5%
neg-mul-178.5%
*-commutative78.5%
+-commutative78.5%
associate-/l/78.5%
associate-*l/78.5%
Simplified78.5%
Taylor expanded in alpha around inf 3.5%
metadata-eval3.5%
Applied egg-rr3.5%
herbie shell --seed 2024160
(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))