
(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 8 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)}^{x}}{\mathsf{E}\left(\right)}
\end{array}
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%
(FPCore (x) :precision binary64 (if (<= (exp (+ (* x x) -1.0)) 0.5) (/ 1.0 (E)) (/ (* x x) (E))))
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;e^{x \cdot x + -1} \leq 0.5:\\
\;\;\;\;\frac{1}{\mathsf{E}\left(\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot x}{\mathsf{E}\left(\right)}\\
\end{array}
\end{array}
if (exp.f64 (neg.f64 (-.f64 #s(literal 1 binary64) (*.f64 x x)))) < 0.5Initial 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 rewrites96.9%
if 0.5 < (exp.f64 (neg.f64 (-.f64 #s(literal 1 binary64) (*.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%
Taylor expanded in x around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6454.8
Applied rewrites54.8%
Taylor expanded in x around inf
Applied rewrites54.8%
Final simplification78.0%
(FPCore (x) :precision binary64 (if (<= (* x x) 0.2) (/ 1.0 (/ (E) (fma x x 1.0))) (exp (* x x))))
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \cdot x \leq 0.2:\\
\;\;\;\;\frac{1}{\frac{\mathsf{E}\left(\right)}{\mathsf{fma}\left(x, x, 1\right)}}\\
\mathbf{else}:\\
\;\;\;\;e^{x \cdot x}\\
\end{array}
\end{array}
if (*.f64 x x) < 0.20000000000000001Initial 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
+-commutativeN/A
unpow2N/A
lower-fma.f6497.9
Applied rewrites97.9%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6497.9
Applied rewrites97.9%
if 0.20000000000000001 < (*.f64 x x) Initial program 100.0%
Taylor expanded in x around inf
unpow2N/A
lower-*.f64100.0
Applied rewrites100.0%
(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 100.0%
lift-neg.f64N/A
neg-sub0N/A
lift--.f64N/A
associate--r-N/A
metadata-evalN/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f64100.0
Applied rewrites100.0%
(FPCore (x) :precision binary64 (pow (E) (- x 1.0)))
\begin{array}{l}
\\
{\mathsf{E}\left(\right)}^{\left(x - 1\right)}
\end{array}
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%
lift-/.f64N/A
div-invN/A
lift-pow.f64N/A
pow-to-expN/A
lift-E.f64N/A
e-exp-1N/A
rec-expN/A
metadata-evalN/A
prod-expN/A
lift-exp.f64N/A
rem-log-expN/A
difference-of-sqr--1N/A
exp-prodN/A
lower-pow.f64N/A
+-commutativeN/A
exp-sumN/A
e-exp-1N/A
lift-E.f64N/A
lift-exp.f64N/A
lower-*.f64N/A
lower--.f64100.0
Applied rewrites100.0%
Taylor expanded in x around 0
lower-E.f6475.7
Applied rewrites75.7%
(FPCore (x) :precision binary64 (/ 1.0 (/ (E) (fma x x 1.0))))
\begin{array}{l}
\\
\frac{1}{\frac{\mathsf{E}\left(\right)}{\mathsf{fma}\left(x, x, 1\right)}}
\end{array}
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%
Taylor expanded in x around 0
+-commutativeN/A
unpow2N/A
lower-fma.f6478.5
Applied rewrites78.5%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6478.5
Applied rewrites78.5%
(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 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
+-commutativeN/A
unpow2N/A
lower-fma.f6478.5
Applied rewrites78.5%
(FPCore (x) :precision binary64 (/ 1.0 (E)))
\begin{array}{l}
\\
\frac{1}{\mathsf{E}\left(\right)}
\end{array}
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%
Taylor expanded in x around 0
Applied rewrites54.8%
herbie shell --seed 2024255
(FPCore (x)
:name "exp neg sub"
:precision binary64
(exp (- (- 1.0 (* x x)))))