
(FPCore (x) :precision binary64 (exp (- (- 1.0 (* x x)))))
double code(double x) {
return exp(-(1.0 - (x * x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = exp(-(1.0d0 - (x * x)))
end function
public static double code(double x) {
return Math.exp(-(1.0 - (x * x)));
}
def code(x): return math.exp(-(1.0 - (x * x)))
function code(x) return exp(Float64(-Float64(1.0 - Float64(x * x)))) end
function tmp = code(x) tmp = exp(-(1.0 - (x * x))); end
code[x_] := N[Exp[(-N[(1.0 - N[(x * x), $MachinePrecision]), $MachinePrecision])], $MachinePrecision]
\begin{array}{l}
\\
e^{-\left(1 - x \cdot x\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 11 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (exp (- (- 1.0 (* x x)))))
double code(double x) {
return exp(-(1.0 - (x * x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = exp(-(1.0d0 - (x * x)))
end function
public static double code(double x) {
return Math.exp(-(1.0 - (x * x)));
}
def code(x): return math.exp(-(1.0 - (x * x)))
function code(x) return exp(Float64(-Float64(1.0 - Float64(x * x)))) end
function tmp = code(x) tmp = exp(-(1.0 - (x * x))); end
code[x_] := N[Exp[(-N[(1.0 - N[(x * x), $MachinePrecision]), $MachinePrecision])], $MachinePrecision]
\begin{array}{l}
\\
e^{-\left(1 - x \cdot x\right)}
\end{array}
(FPCore (x) :precision binary64 (/ (pow (exp (- x)) (- x)) (E)))
\begin{array}{l}
\\
\frac{{\left(e^{-x}\right)}^{\left(-x\right)}}{\mathsf{E}\left(\right)}
\end{array}
Initial program 99.9%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
/-rgt-identityN/A
clear-numN/A
lift-pow.f64N/A
pow-flipN/A
pow-flipN/A
exp-to-powN/A
lift-exp.f64N/A
rem-log-expN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
rem-log-expN/A
lift-exp.f64N/A
neg-logN/A
exp-to-powN/A
lower-pow.f64N/A
inv-powN/A
lower-pow.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
lift-pow.f64N/A
unpow-1N/A
lift-exp.f64N/A
rec-expN/A
lift-neg.f64N/A
lower-exp.f64100.0
Applied rewrites100.0%
(FPCore (x) :precision binary64 (/ (pow (exp x) x) (E)))
\begin{array}{l}
\\
\frac{{\left(e^{x}\right)}^{x}}{\mathsf{E}\left(\right)}
\end{array}
Initial program 99.9%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
(FPCore (x) :precision binary64 (pow (exp (- x -1.0)) (- x 1.0)))
double code(double x) {
return pow(exp((x - -1.0)), (x - 1.0));
}
real(8) function code(x)
real(8), intent (in) :: x
code = exp((x - (-1.0d0))) ** (x - 1.0d0)
end function
public static double code(double x) {
return Math.pow(Math.exp((x - -1.0)), (x - 1.0));
}
def code(x): return math.pow(math.exp((x - -1.0)), (x - 1.0))
function code(x) return exp(Float64(x - -1.0)) ^ Float64(x - 1.0) end
function tmp = code(x) tmp = exp((x - -1.0)) ^ (x - 1.0); end
code[x_] := N[Power[N[Exp[N[(x - -1.0), $MachinePrecision]], $MachinePrecision], N[(x - 1.0), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
{\left(e^{x - -1}\right)}^{\left(x - 1\right)}
\end{array}
Initial program 99.9%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
/-rgt-identityN/A
clear-numN/A
lift-pow.f64N/A
pow-flipN/A
pow-flipN/A
exp-to-powN/A
lift-exp.f64N/A
rem-log-expN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
rem-log-expN/A
lift-exp.f64N/A
neg-logN/A
exp-to-powN/A
lower-pow.f64N/A
inv-powN/A
lower-pow.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
lift-/.f64N/A
lift-pow.f64N/A
lift-pow.f64N/A
pow-powN/A
lift-neg.f64N/A
neg-mul-1N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
lift-exp.f64N/A
*-rgt-identityN/A
pow-expN/A
lift-*.f64N/A
lift-E.f64N/A
e-exp-1N/A
div-expN/A
lift-*.f64N/A
difference-of-sqr-1N/A
lift-+.f64N/A
lift--.f64N/A
exp-prodN/A
lower-pow.f64N/A
Applied rewrites100.0%
(FPCore (x) :precision binary64 (/ (pow (sqrt (E)) (* 2.0 (* x x))) (E)))
\begin{array}{l}
\\
\frac{{\left(\sqrt{\mathsf{E}\left(\right)}\right)}^{\left(2 \cdot \left(x \cdot x\right)\right)}}{\mathsf{E}\left(\right)}
\end{array}
Initial program 99.9%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
lift-pow.f64N/A
lift-exp.f64N/A
pow-expN/A
*-lft-identityN/A
exp-prodN/A
e-exp-1N/A
add-sqr-sqrtN/A
pow2N/A
pow-powN/A
lower-pow.f64N/A
lift-E.f64N/A
lower-sqrt.f64N/A
lower-*.f64N/A
lower-*.f6499.9
Applied rewrites99.9%
(FPCore (x) :precision binary64 (/ (exp (* x x)) (E)))
\begin{array}{l}
\\
\frac{e^{x \cdot x}}{\mathsf{E}\left(\right)}
\end{array}
Initial program 99.9%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
lift-pow.f64N/A
lift-exp.f64N/A
pow-expN/A
lower-exp.f64N/A
lower-*.f6499.9
Applied rewrites99.9%
(FPCore (x) :precision binary64 (exp (fma x x -1.0)))
double code(double x) {
return exp(fma(x, x, -1.0));
}
function code(x) return exp(fma(x, x, -1.0)) end
code[x_] := N[Exp[N[(x * x + -1.0), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
e^{\mathsf{fma}\left(x, x, -1\right)}
\end{array}
Initial program 99.9%
lift-neg.f64N/A
neg-sub0N/A
lift--.f64N/A
associate--r-N/A
metadata-evalN/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6499.9
Applied rewrites99.9%
(FPCore (x) :precision binary64 (pow (E) (- x 1.0)))
\begin{array}{l}
\\
{\mathsf{E}\left(\right)}^{\left(x - 1\right)}
\end{array}
Initial program 99.9%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
/-rgt-identityN/A
clear-numN/A
lift-pow.f64N/A
pow-flipN/A
pow-flipN/A
exp-to-powN/A
lift-exp.f64N/A
rem-log-expN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
rem-log-expN/A
lift-exp.f64N/A
neg-logN/A
exp-to-powN/A
lower-pow.f64N/A
inv-powN/A
lower-pow.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
lift-/.f64N/A
lift-pow.f64N/A
lift-pow.f64N/A
pow-powN/A
lift-neg.f64N/A
neg-mul-1N/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
lift-exp.f64N/A
*-rgt-identityN/A
pow-expN/A
lift-*.f64N/A
lift-E.f64N/A
e-exp-1N/A
div-expN/A
lift-*.f64N/A
difference-of-sqr-1N/A
lift-+.f64N/A
lift--.f64N/A
exp-prodN/A
lower-pow.f64N/A
Applied rewrites100.0%
Taylor expanded in x around 0
exp-1-eN/A
lower-E.f6472.7
Applied rewrites72.7%
(FPCore (x) :precision binary64 (if (<= (* x x) 1e+304) (/ (/ (* (fma x x -1.0) (fma x x 1.0)) (fma x x -1.0)) (E)) (/ (* x x) (E))))
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \cdot x \leq 10^{+304}:\\
\;\;\;\;\frac{\frac{\mathsf{fma}\left(x, x, -1\right) \cdot \mathsf{fma}\left(x, x, 1\right)}{\mathsf{fma}\left(x, x, -1\right)}}{\mathsf{E}\left(\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot x}{\mathsf{E}\left(\right)}\\
\end{array}
\end{array}
if (*.f64 x x) < 9.9999999999999994e303Initial program 99.9%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
/-rgt-identityN/A
clear-numN/A
lift-pow.f64N/A
pow-flipN/A
pow-flipN/A
exp-to-powN/A
lift-exp.f64N/A
rem-log-expN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
rem-log-expN/A
lift-exp.f64N/A
neg-logN/A
exp-to-powN/A
lower-pow.f64N/A
inv-powN/A
lower-pow.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in x around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6464.3
Applied rewrites64.3%
Applied rewrites82.4%
if 9.9999999999999994e303 < (*.f64 x x) Initial program 100.0%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
/-rgt-identityN/A
clear-numN/A
lift-pow.f64N/A
pow-flipN/A
pow-flipN/A
exp-to-powN/A
lift-exp.f64N/A
rem-log-expN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
rem-log-expN/A
lift-exp.f64N/A
neg-logN/A
exp-to-powN/A
lower-pow.f64N/A
inv-powN/A
lower-pow.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in x around 0
+-commutativeN/A
unpow2N/A
lower-fma.f64100.0
Applied rewrites100.0%
Taylor expanded in x around inf
Applied rewrites100.0%
(FPCore (x) :precision binary64 (if (<= (* x x) 1.0) (/ 1.0 (E)) (/ (* x x) (E))))
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \cdot x \leq 1:\\
\;\;\;\;\frac{1}{\mathsf{E}\left(\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot x}{\mathsf{E}\left(\right)}\\
\end{array}
\end{array}
if (*.f64 x x) < 1Initial program 100.0%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
Taylor expanded in x around 0
Applied rewrites98.7%
if 1 < (*.f64 x x) Initial program 99.9%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
/-rgt-identityN/A
clear-numN/A
lift-pow.f64N/A
pow-flipN/A
pow-flipN/A
exp-to-powN/A
lift-exp.f64N/A
rem-log-expN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
rem-log-expN/A
lift-exp.f64N/A
neg-logN/A
exp-to-powN/A
lower-pow.f64N/A
inv-powN/A
lower-pow.f64N/A
lower-neg.f64100.0
Applied rewrites100.0%
Taylor expanded in x around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6445.0
Applied rewrites45.0%
Taylor expanded in x around inf
Applied rewrites45.0%
(FPCore (x) :precision binary64 (/ (fma x x 1.0) (E)))
\begin{array}{l}
\\
\frac{\mathsf{fma}\left(x, x, 1\right)}{\mathsf{E}\left(\right)}
\end{array}
Initial program 99.9%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
Taylor expanded in x around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6471.9
Applied rewrites71.9%
(FPCore (x) :precision binary64 (/ 1.0 (E)))
\begin{array}{l}
\\
\frac{1}{\mathsf{E}\left(\right)}
\end{array}
Initial program 99.9%
lift-exp.f64N/A
lift-neg.f64N/A
exp-negN/A
lift--.f64N/A
exp-diffN/A
clear-numN/A
lower-/.f64N/A
lift-*.f64N/A
exp-prodN/A
lower-pow.f64N/A
lower-exp.f64N/A
exp-1-eN/A
lower-E.f64100.0
Applied rewrites100.0%
Taylor expanded in x around 0
Applied rewrites50.6%
herbie shell --seed 2024313
(FPCore (x)
:name "exp neg sub"
:precision binary64
(exp (- (- 1.0 (* x x)))))