
(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 5 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 (let* ((t_0 (+ beta (+ alpha 2.0)))) (/ (- (/ 1.0 (* (/ 1.0 beta) t_0)) (+ (/ alpha t_0) -1.0)) 2.0)))
double code(double alpha, double beta) {
double t_0 = beta + (alpha + 2.0);
return ((1.0 / ((1.0 / beta) * t_0)) - ((alpha / t_0) + -1.0)) / 2.0;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: t_0
t_0 = beta + (alpha + 2.0d0)
code = ((1.0d0 / ((1.0d0 / beta) * t_0)) - ((alpha / t_0) + (-1.0d0))) / 2.0d0
end function
public static double code(double alpha, double beta) {
double t_0 = beta + (alpha + 2.0);
return ((1.0 / ((1.0 / beta) * t_0)) - ((alpha / t_0) + -1.0)) / 2.0;
}
def code(alpha, beta): t_0 = beta + (alpha + 2.0) return ((1.0 / ((1.0 / beta) * t_0)) - ((alpha / t_0) + -1.0)) / 2.0
function code(alpha, beta) t_0 = Float64(beta + Float64(alpha + 2.0)) return Float64(Float64(Float64(1.0 / Float64(Float64(1.0 / beta) * t_0)) - Float64(Float64(alpha / t_0) + -1.0)) / 2.0) end
function tmp = code(alpha, beta) t_0 = beta + (alpha + 2.0); tmp = ((1.0 / ((1.0 / beta) * t_0)) - ((alpha / t_0) + -1.0)) / 2.0; end
code[alpha_, beta_] := Block[{t$95$0 = N[(beta + N[(alpha + 2.0), $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(1.0 / N[(N[(1.0 / beta), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision] - N[(N[(alpha / t$95$0), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \beta + \left(\alpha + 2\right)\\
\frac{\frac{1}{\frac{1}{\beta} \cdot t\_0} - \left(\frac{\alpha}{t\_0} + -1\right)}{2}
\end{array}
\end{array}
Initial program 76.4%
div-subN/A
associate-+l-N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f6477.1
Applied egg-rr77.1%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6477.1
Applied egg-rr77.1%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6477.1
Applied egg-rr77.1%
(FPCore (alpha beta) :precision binary64 (let* ((t_0 (+ beta (+ alpha 2.0)))) (/ (fma (/ 1.0 t_0) beta (- 1.0 (/ alpha t_0))) 2.0)))
double code(double alpha, double beta) {
double t_0 = beta + (alpha + 2.0);
return fma((1.0 / t_0), beta, (1.0 - (alpha / t_0))) / 2.0;
}
function code(alpha, beta) t_0 = Float64(beta + Float64(alpha + 2.0)) return Float64(fma(Float64(1.0 / t_0), beta, Float64(1.0 - Float64(alpha / t_0))) / 2.0) end
code[alpha_, beta_] := Block[{t$95$0 = N[(beta + N[(alpha + 2.0), $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(1.0 / t$95$0), $MachinePrecision] * beta + N[(1.0 - N[(alpha / t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \beta + \left(\alpha + 2\right)\\
\frac{\mathsf{fma}\left(\frac{1}{t\_0}, \beta, 1 - \frac{\alpha}{t\_0}\right)}{2}
\end{array}
\end{array}
Initial program 76.4%
div-subN/A
associate-+l-N/A
sub-negN/A
div-invN/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f6477.1
Applied egg-rr77.1%
*-commutativeN/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
neg-sub0N/A
+-commutativeN/A
associate--r+N/A
metadata-evalN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6477.1
Applied egg-rr77.1%
(FPCore (alpha beta) :precision binary64 (let* ((t_0 (+ beta (+ alpha 2.0)))) (fma (/ beta t_0) 0.5 (* 0.5 (- (/ (- alpha) t_0) -1.0)))))
double code(double alpha, double beta) {
double t_0 = beta + (alpha + 2.0);
return fma((beta / t_0), 0.5, (0.5 * ((-alpha / t_0) - -1.0)));
}
function code(alpha, beta) t_0 = Float64(beta + Float64(alpha + 2.0)) return fma(Float64(beta / t_0), 0.5, Float64(0.5 * Float64(Float64(Float64(-alpha) / t_0) - -1.0))) end
code[alpha_, beta_] := Block[{t$95$0 = N[(beta + N[(alpha + 2.0), $MachinePrecision]), $MachinePrecision]}, N[(N[(beta / t$95$0), $MachinePrecision] * 0.5 + N[(0.5 * N[(N[((-alpha) / t$95$0), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \beta + \left(\alpha + 2\right)\\
\mathsf{fma}\left(\frac{\beta}{t\_0}, 0.5, 0.5 \cdot \left(\frac{-\alpha}{t\_0} - -1\right)\right)
\end{array}
\end{array}
Initial program 76.4%
div-subN/A
associate-+l-N/A
sub-negN/A
div-invN/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f6477.1
Applied egg-rr77.1%
div-invN/A
sub-negN/A
div-subN/A
sub-negN/A
div-invN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
div-invN/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr77.1%
Final simplification77.1%
(FPCore (alpha beta) :precision binary64 (fma (/ (- beta alpha) (+ beta (+ alpha 2.0))) 0.5 0.5))
double code(double alpha, double beta) {
return fma(((beta - alpha) / (beta + (alpha + 2.0))), 0.5, 0.5);
}
function code(alpha, beta) return fma(Float64(Float64(beta - alpha) / Float64(beta + Float64(alpha + 2.0))), 0.5, 0.5) end
code[alpha_, beta_] := N[(N[(N[(beta - alpha), $MachinePrecision] / N[(beta + N[(alpha + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\frac{\beta - \alpha}{\beta + \left(\alpha + 2\right)}, 0.5, 0.5\right)
\end{array}
Initial program 76.4%
clear-numN/A
associate-/r/N/A
distribute-rgt-inN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
metadata-eval76.4
Applied egg-rr76.4%
(FPCore (alpha beta) :precision binary64 (fma (- beta alpha) (/ 0.5 (+ beta (+ alpha 2.0))) 0.5))
double code(double alpha, double beta) {
return fma((beta - alpha), (0.5 / (beta + (alpha + 2.0))), 0.5);
}
function code(alpha, beta) return fma(Float64(beta - alpha), Float64(0.5 / Float64(beta + Float64(alpha + 2.0))), 0.5) end
code[alpha_, beta_] := N[(N[(beta - alpha), $MachinePrecision] * N[(0.5 / N[(beta + N[(alpha + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 0.5), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\beta - \alpha, \frac{0.5}{\beta + \left(\alpha + 2\right)}, 0.5\right)
\end{array}
Initial program 76.4%
div-subN/A
associate-+l-N/A
sub-negN/A
div-invN/A
accelerator-lowering-fma.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
neg-lowering-neg.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
+-lowering-+.f6477.1
Applied egg-rr77.1%
div-invN/A
div-invN/A
sub-negN/A
associate--r+N/A
div-subN/A
sub-negN/A
metadata-evalN/A
metadata-evalN/A
distribute-lft1-inN/A
associate-*l/N/A
associate-/l*N/A
accelerator-lowering-fma.f64N/A
Applied egg-rr76.3%
herbie shell --seed 2024208
(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))