
(FPCore (w l) :precision binary64 (* (exp (- w)) (pow l (exp w))))
double code(double w, double l) {
return exp(-w) * pow(l, exp(w));
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
code = exp(-w) * (l ** exp(w))
end function
public static double code(double w, double l) {
return Math.exp(-w) * Math.pow(l, Math.exp(w));
}
def code(w, l): return math.exp(-w) * math.pow(l, math.exp(w))
function code(w, l) return Float64(exp(Float64(-w)) * (l ^ exp(w))) end
function tmp = code(w, l) tmp = exp(-w) * (l ^ exp(w)); end
code[w_, l_] := N[(N[Exp[(-w)], $MachinePrecision] * N[Power[l, N[Exp[w], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{-w} \cdot {\ell}^{\left(e^{w}\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (w l) :precision binary64 (* (exp (- w)) (pow l (exp w))))
double code(double w, double l) {
return exp(-w) * pow(l, exp(w));
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
code = exp(-w) * (l ** exp(w))
end function
public static double code(double w, double l) {
return Math.exp(-w) * Math.pow(l, Math.exp(w));
}
def code(w, l): return math.exp(-w) * math.pow(l, math.exp(w))
function code(w, l) return Float64(exp(Float64(-w)) * (l ^ exp(w))) end
function tmp = code(w, l) tmp = exp(-w) * (l ^ exp(w)); end
code[w_, l_] := N[(N[Exp[(-w)], $MachinePrecision] * N[Power[l, N[Exp[w], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
e^{-w} \cdot {\ell}^{\left(e^{w}\right)}
\end{array}
(FPCore (w l) :precision binary64 (* (/ 1.0 (cbrt (exp (+ w w)))) (/ (pow l (exp w)) (cbrt (exp w)))))
double code(double w, double l) {
return (1.0 / cbrt(exp((w + w)))) * (pow(l, exp(w)) / cbrt(exp(w)));
}
public static double code(double w, double l) {
return (1.0 / Math.cbrt(Math.exp((w + w)))) * (Math.pow(l, Math.exp(w)) / Math.cbrt(Math.exp(w)));
}
function code(w, l) return Float64(Float64(1.0 / cbrt(exp(Float64(w + w)))) * Float64((l ^ exp(w)) / cbrt(exp(w)))) end
code[w_, l_] := N[(N[(1.0 / N[Power[N[Exp[N[(w + w), $MachinePrecision]], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] * N[(N[Power[l, N[Exp[w], $MachinePrecision]], $MachinePrecision] / N[Power[N[Exp[w], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\sqrt[3]{e^{w + w}}} \cdot \frac{{\ell}^{\left(e^{w}\right)}}{\sqrt[3]{e^{w}}}
\end{array}
Initial program 99.6%
exp-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
Simplified99.6%
*-un-lft-identity99.6%
add-cube-cbrt99.6%
times-frac99.6%
cbrt-unprod99.6%
prod-exp99.6%
Applied egg-rr99.6%
Final simplification99.6%
(FPCore (w l) :precision binary64 (let* ((t_0 (exp (* w -0.3333333333333333)))) (* (pow l (exp w)) (* t_0 (pow t_0 2.0)))))
double code(double w, double l) {
double t_0 = exp((w * -0.3333333333333333));
return pow(l, exp(w)) * (t_0 * pow(t_0, 2.0));
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: t_0
t_0 = exp((w * (-0.3333333333333333d0)))
code = (l ** exp(w)) * (t_0 * (t_0 ** 2.0d0))
end function
public static double code(double w, double l) {
double t_0 = Math.exp((w * -0.3333333333333333));
return Math.pow(l, Math.exp(w)) * (t_0 * Math.pow(t_0, 2.0));
}
def code(w, l): t_0 = math.exp((w * -0.3333333333333333)) return math.pow(l, math.exp(w)) * (t_0 * math.pow(t_0, 2.0))
function code(w, l) t_0 = exp(Float64(w * -0.3333333333333333)) return Float64((l ^ exp(w)) * Float64(t_0 * (t_0 ^ 2.0))) end
function tmp = code(w, l) t_0 = exp((w * -0.3333333333333333)); tmp = (l ^ exp(w)) * (t_0 * (t_0 ^ 2.0)); end
code[w_, l_] := Block[{t$95$0 = N[Exp[N[(w * -0.3333333333333333), $MachinePrecision]], $MachinePrecision]}, N[(N[Power[l, N[Exp[w], $MachinePrecision]], $MachinePrecision] * N[(t$95$0 * N[Power[t$95$0, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{w \cdot -0.3333333333333333}\\
{\ell}^{\left(e^{w}\right)} \cdot \left(t_0 \cdot {t_0}^{2}\right)
\end{array}
\end{array}
Initial program 99.6%
exp-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
Simplified99.6%
Taylor expanded in l around 0 94.3%
pow-to-exp99.6%
add-cbrt-cube99.6%
exp-sum99.6%
cbrt-prod99.6%
associate-/l/99.6%
div-inv99.6%
div-inv99.6%
associate-*l*99.6%
add-exp-log99.6%
rec-exp99.6%
pow1/399.6%
log-pow99.6%
add-log-exp99.6%
Applied egg-rr99.6%
Taylor expanded in w around -inf 99.6%
*-commutative99.6%
Simplified99.6%
Taylor expanded in w around -inf 99.6%
Final simplification99.6%
(FPCore (w l) :precision binary64 (* (pow l (exp w)) (cbrt (exp (* w (- 3.0))))))
double code(double w, double l) {
return pow(l, exp(w)) * cbrt(exp((w * -3.0)));
}
public static double code(double w, double l) {
return Math.pow(l, Math.exp(w)) * Math.cbrt(Math.exp((w * -3.0)));
}
function code(w, l) return Float64((l ^ exp(w)) * cbrt(exp(Float64(w * Float64(-3.0))))) end
code[w_, l_] := N[(N[Power[l, N[Exp[w], $MachinePrecision]], $MachinePrecision] * N[Power[N[Exp[N[(w * (-3.0)), $MachinePrecision]], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
{\ell}^{\left(e^{w}\right)} \cdot \sqrt[3]{e^{w \cdot \left(-3\right)}}
\end{array}
Initial program 99.6%
exp-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
Simplified99.6%
*-un-lft-identity99.6%
add-cube-cbrt99.6%
times-frac99.6%
cbrt-unprod99.6%
prod-exp99.6%
Applied egg-rr99.6%
Taylor expanded in w around inf 99.6%
metadata-eval99.6%
distribute-lft-neg-in99.6%
unpow1/399.6%
prod-exp99.6%
distribute-lft-neg-in99.6%
metadata-eval99.6%
rec-exp99.6%
distribute-lft1-in99.6%
metadata-eval99.6%
Simplified99.6%
Final simplification99.6%
(FPCore (w l) :precision binary64 (/ (pow l (exp w)) (exp w)))
double code(double w, double l) {
return pow(l, exp(w)) / exp(w);
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
code = (l ** exp(w)) / exp(w)
end function
public static double code(double w, double l) {
return Math.pow(l, Math.exp(w)) / Math.exp(w);
}
def code(w, l): return math.pow(l, math.exp(w)) / math.exp(w)
function code(w, l) return Float64((l ^ exp(w)) / exp(w)) end
function tmp = code(w, l) tmp = (l ^ exp(w)) / exp(w); end
code[w_, l_] := N[(N[Power[l, N[Exp[w], $MachinePrecision]], $MachinePrecision] / N[Exp[w], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{{\ell}^{\left(e^{w}\right)}}{e^{w}}
\end{array}
Initial program 99.6%
exp-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
Simplified99.6%
Final simplification99.6%
(FPCore (w l) :precision binary64 (/ l (exp w)))
double code(double w, double l) {
return l / exp(w);
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
code = l / exp(w)
end function
public static double code(double w, double l) {
return l / Math.exp(w);
}
def code(w, l): return l / math.exp(w)
function code(w, l) return Float64(l / exp(w)) end
function tmp = code(w, l) tmp = l / exp(w); end
code[w_, l_] := N[(l / N[Exp[w], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell}{e^{w}}
\end{array}
Initial program 99.6%
exp-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
Simplified99.6%
Taylor expanded in w around 0 97.3%
Final simplification97.3%
(FPCore (w l) :precision binary64 (- l (* w l)))
double code(double w, double l) {
return l - (w * l);
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
code = l - (w * l)
end function
public static double code(double w, double l) {
return l - (w * l);
}
def code(w, l): return l - (w * l)
function code(w, l) return Float64(l - Float64(w * l)) end
function tmp = code(w, l) tmp = l - (w * l); end
code[w_, l_] := N[(l - N[(w * l), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\ell - w \cdot \ell
\end{array}
Initial program 99.6%
exp-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
Simplified99.6%
Taylor expanded in w around 0 97.3%
Taylor expanded in w around 0 67.2%
+-commutative67.2%
mul-1-neg67.2%
unsub-neg67.2%
*-commutative67.2%
Simplified67.2%
Final simplification67.2%
(FPCore (w l) :precision binary64 l)
double code(double w, double l) {
return l;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
code = l
end function
public static double code(double w, double l) {
return l;
}
def code(w, l): return l
function code(w, l) return l end
function tmp = code(w, l) tmp = l; end
code[w_, l_] := l
\begin{array}{l}
\\
\ell
\end{array}
Initial program 99.6%
exp-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
Simplified99.6%
*-un-lft-identity99.6%
add-cube-cbrt99.6%
times-frac99.6%
cbrt-unprod99.6%
prod-exp99.6%
Applied egg-rr99.6%
Taylor expanded in w around 0 61.7%
Final simplification61.7%
herbie shell --seed 2023203
(FPCore (w l)
:name "exp-w (used to crash)"
:precision binary64
(* (exp (- w)) (pow l (exp w))))