
(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 13 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 (/ (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%
remove-double-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
remove-double-neg99.6%
Simplified99.6%
(FPCore (w l) :precision binary64 (if (<= w -220.0) (/ l (exp w)) (/ (pow l (exp w)) (+ 1.0 (* w (+ 1.0 (* w 0.5)))))))
double code(double w, double l) {
double tmp;
if (w <= -220.0) {
tmp = l / exp(w);
} else {
tmp = pow(l, exp(w)) / (1.0 + (w * (1.0 + (w * 0.5))));
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= (-220.0d0)) then
tmp = l / exp(w)
else
tmp = (l ** exp(w)) / (1.0d0 + (w * (1.0d0 + (w * 0.5d0))))
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= -220.0) {
tmp = l / Math.exp(w);
} else {
tmp = Math.pow(l, Math.exp(w)) / (1.0 + (w * (1.0 + (w * 0.5))));
}
return tmp;
}
def code(w, l): tmp = 0 if w <= -220.0: tmp = l / math.exp(w) else: tmp = math.pow(l, math.exp(w)) / (1.0 + (w * (1.0 + (w * 0.5)))) return tmp
function code(w, l) tmp = 0.0 if (w <= -220.0) tmp = Float64(l / exp(w)); else tmp = Float64((l ^ exp(w)) / Float64(1.0 + Float64(w * Float64(1.0 + Float64(w * 0.5))))); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= -220.0) tmp = l / exp(w); else tmp = (l ^ exp(w)) / (1.0 + (w * (1.0 + (w * 0.5)))); end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, -220.0], N[(l / N[Exp[w], $MachinePrecision]), $MachinePrecision], N[(N[Power[l, N[Exp[w], $MachinePrecision]], $MachinePrecision] / N[(1.0 + N[(w * N[(1.0 + N[(w * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq -220:\\
\;\;\;\;\frac{\ell}{e^{w}}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\ell}^{\left(e^{w}\right)}}{1 + w \cdot \left(1 + w \cdot 0.5\right)}\\
\end{array}
\end{array}
if w < -220Initial program 100.0%
exp-neg100.0%
remove-double-neg100.0%
associate-*l/100.0%
*-lft-identity100.0%
remove-double-neg100.0%
Simplified100.0%
add-sqr-sqrt0.0%
sqrt-unprod50.0%
sqr-neg50.0%
sqrt-unprod50.0%
add-sqr-sqrt50.0%
add-sqr-sqrt50.0%
sqrt-unprod50.0%
add-sqr-sqrt50.0%
sqrt-unprod50.0%
sqr-neg50.0%
sqrt-unprod0.0%
add-sqr-sqrt0.0%
pow10.0%
exp-neg0.0%
inv-pow0.0%
pow-prod-up100.0%
metadata-eval100.0%
metadata-eval100.0%
metadata-eval100.0%
Applied egg-rr100.0%
Taylor expanded in l around 0 100.0%
if -220 < w Initial program 99.5%
exp-neg99.5%
remove-double-neg99.5%
associate-*l/99.5%
*-lft-identity99.5%
remove-double-neg99.5%
Simplified99.5%
Taylor expanded in w around 0 99.4%
*-commutative99.4%
Simplified99.4%
(FPCore (w l) :precision binary64 (if (<= w -1.0) (/ l (exp w)) (/ (pow l (+ 1.0 (* w (+ 1.0 (* w 0.5))))) (+ w 1.0))))
double code(double w, double l) {
double tmp;
if (w <= -1.0) {
tmp = l / exp(w);
} else {
tmp = pow(l, (1.0 + (w * (1.0 + (w * 0.5))))) / (w + 1.0);
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= (-1.0d0)) then
tmp = l / exp(w)
else
tmp = (l ** (1.0d0 + (w * (1.0d0 + (w * 0.5d0))))) / (w + 1.0d0)
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= -1.0) {
tmp = l / Math.exp(w);
} else {
tmp = Math.pow(l, (1.0 + (w * (1.0 + (w * 0.5))))) / (w + 1.0);
}
return tmp;
}
def code(w, l): tmp = 0 if w <= -1.0: tmp = l / math.exp(w) else: tmp = math.pow(l, (1.0 + (w * (1.0 + (w * 0.5))))) / (w + 1.0) return tmp
function code(w, l) tmp = 0.0 if (w <= -1.0) tmp = Float64(l / exp(w)); else tmp = Float64((l ^ Float64(1.0 + Float64(w * Float64(1.0 + Float64(w * 0.5))))) / Float64(w + 1.0)); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= -1.0) tmp = l / exp(w); else tmp = (l ^ (1.0 + (w * (1.0 + (w * 0.5))))) / (w + 1.0); end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, -1.0], N[(l / N[Exp[w], $MachinePrecision]), $MachinePrecision], N[(N[Power[l, N[(1.0 + N[(w * N[(1.0 + N[(w * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(w + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq -1:\\
\;\;\;\;\frac{\ell}{e^{w}}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\ell}^{\left(1 + w \cdot \left(1 + w \cdot 0.5\right)\right)}}{w + 1}\\
\end{array}
\end{array}
if w < -1Initial program 100.0%
exp-neg100.0%
remove-double-neg100.0%
associate-*l/100.0%
*-lft-identity100.0%
remove-double-neg100.0%
Simplified100.0%
add-sqr-sqrt0.0%
sqrt-unprod50.0%
sqr-neg50.0%
sqrt-unprod50.0%
add-sqr-sqrt50.0%
add-sqr-sqrt50.0%
sqrt-unprod50.0%
add-sqr-sqrt50.0%
sqrt-unprod50.0%
sqr-neg50.0%
sqrt-unprod0.0%
add-sqr-sqrt0.0%
pow10.0%
exp-neg0.0%
inv-pow0.0%
pow-prod-up100.0%
metadata-eval100.0%
metadata-eval100.0%
metadata-eval100.0%
Applied egg-rr100.0%
Taylor expanded in l around 0 100.0%
if -1 < w Initial program 99.5%
exp-neg99.5%
remove-double-neg99.5%
associate-*l/99.5%
*-lft-identity99.5%
remove-double-neg99.5%
Simplified99.5%
Taylor expanded in w around 0 99.2%
+-commutative99.2%
Simplified99.2%
Taylor expanded in w around 0 99.2%
*-commutative99.4%
Simplified99.2%
(FPCore (w l) :precision binary64 (* l (/ 1.0 (exp w))))
double code(double w, double l) {
return l * (1.0 / exp(w));
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
code = l * (1.0d0 / exp(w))
end function
public static double code(double w, double l) {
return l * (1.0 / Math.exp(w));
}
def code(w, l): return l * (1.0 / math.exp(w))
function code(w, l) return Float64(l * Float64(1.0 / exp(w))) end
function tmp = code(w, l) tmp = l * (1.0 / exp(w)); end
code[w_, l_] := N[(l * N[(1.0 / N[Exp[w], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\ell \cdot \frac{1}{e^{w}}
\end{array}
Initial program 99.6%
exp-neg99.6%
remove-double-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
remove-double-neg99.6%
Simplified99.6%
add-sqr-sqrt38.9%
sqrt-unprod81.9%
sqr-neg81.9%
sqrt-unprod42.9%
add-sqr-sqrt80.7%
add-sqr-sqrt80.7%
sqrt-unprod80.7%
add-sqr-sqrt42.9%
sqrt-unprod67.4%
sqr-neg67.4%
sqrt-unprod24.5%
add-sqr-sqrt50.3%
pow150.3%
exp-neg50.3%
inv-pow50.3%
pow-prod-up97.9%
metadata-eval97.9%
metadata-eval97.9%
metadata-eval97.9%
Applied egg-rr97.9%
frac-2neg97.9%
div-inv97.9%
*-rgt-identity97.9%
Applied egg-rr97.9%
Final simplification97.9%
(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%
remove-double-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
remove-double-neg99.6%
Simplified99.6%
add-sqr-sqrt38.9%
sqrt-unprod81.9%
sqr-neg81.9%
sqrt-unprod42.9%
add-sqr-sqrt80.7%
add-sqr-sqrt80.7%
sqrt-unprod80.7%
add-sqr-sqrt42.9%
sqrt-unprod67.4%
sqr-neg67.4%
sqrt-unprod24.5%
add-sqr-sqrt50.3%
pow150.3%
exp-neg50.3%
inv-pow50.3%
pow-prod-up97.9%
metadata-eval97.9%
metadata-eval97.9%
metadata-eval97.9%
Applied egg-rr97.9%
Taylor expanded in l around 0 97.9%
(FPCore (w l) :precision binary64 (if (<= w 0.098) (+ l (* l (* w (+ -1.0 (* w (+ 0.5 (* w -0.16666666666666666))))))) (* l 0.0)))
double code(double w, double l) {
double tmp;
if (w <= 0.098) {
tmp = l + (l * (w * (-1.0 + (w * (0.5 + (w * -0.16666666666666666))))));
} else {
tmp = l * 0.0;
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= 0.098d0) then
tmp = l + (l * (w * ((-1.0d0) + (w * (0.5d0 + (w * (-0.16666666666666666d0)))))))
else
tmp = l * 0.0d0
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= 0.098) {
tmp = l + (l * (w * (-1.0 + (w * (0.5 + (w * -0.16666666666666666))))));
} else {
tmp = l * 0.0;
}
return tmp;
}
def code(w, l): tmp = 0 if w <= 0.098: tmp = l + (l * (w * (-1.0 + (w * (0.5 + (w * -0.16666666666666666)))))) else: tmp = l * 0.0 return tmp
function code(w, l) tmp = 0.0 if (w <= 0.098) tmp = Float64(l + Float64(l * Float64(w * Float64(-1.0 + Float64(w * Float64(0.5 + Float64(w * -0.16666666666666666))))))); else tmp = Float64(l * 0.0); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= 0.098) tmp = l + (l * (w * (-1.0 + (w * (0.5 + (w * -0.16666666666666666)))))); else tmp = l * 0.0; end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, 0.098], N[(l + N[(l * N[(w * N[(-1.0 + N[(w * N[(0.5 + N[(w * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l * 0.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq 0.098:\\
\;\;\;\;\ell + \ell \cdot \left(w \cdot \left(-1 + w \cdot \left(0.5 + w \cdot -0.16666666666666666\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot 0\\
\end{array}
\end{array}
if w < 0.098000000000000004Initial program 99.6%
exp-neg99.6%
remove-double-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
remove-double-neg99.6%
Simplified99.6%
add-sqr-sqrt28.9%
sqrt-unprod78.9%
sqr-neg78.9%
sqrt-unprod50.0%
add-sqr-sqrt78.4%
add-sqr-sqrt78.4%
sqrt-unprod78.4%
add-sqr-sqrt50.0%
sqrt-unprod78.4%
sqr-neg78.4%
sqrt-unprod28.4%
add-sqr-sqrt58.4%
pow158.4%
exp-neg58.4%
inv-pow58.4%
pow-prod-up98.4%
metadata-eval98.4%
metadata-eval98.4%
metadata-eval98.4%
Applied egg-rr98.4%
Taylor expanded in w around 0 85.5%
Taylor expanded in l around 0 88.1%
if 0.098000000000000004 < w Initial program 100.0%
exp-neg100.0%
remove-double-neg100.0%
associate-*l/100.0%
*-lft-identity100.0%
remove-double-neg100.0%
Simplified100.0%
add-sqr-sqrt100.0%
sqrt-unprod100.0%
sqr-neg100.0%
sqrt-unprod0.0%
add-sqr-sqrt94.7%
add-sqr-sqrt94.7%
sqrt-unprod94.7%
add-sqr-sqrt0.0%
sqrt-unprod0.4%
sqr-neg0.4%
sqrt-unprod0.4%
add-sqr-sqrt0.4%
pow10.4%
exp-neg0.4%
inv-pow0.4%
pow-prod-up94.8%
metadata-eval94.8%
metadata-eval94.8%
metadata-eval94.8%
Applied egg-rr94.8%
Taylor expanded in w around 0 3.3%
mul-1-neg3.3%
unsub-neg3.3%
Simplified3.3%
Taylor expanded in w around inf 3.3%
neg-mul-13.3%
distribute-rgt-neg-in3.3%
Simplified3.3%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod3.4%
add-sqr-sqrt3.4%
add-log-exp1.9%
pow11.9%
pow11.9%
add-sqr-sqrt1.9%
sqrt-unprod1.9%
sqr-neg1.9%
sqrt-unprod0.0%
add-sqr-sqrt1.7%
add-sqr-sqrt1.7%
sqrt-unprod1.7%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod2.7%
add-sqr-sqrt2.7%
pow12.7%
exp-neg2.7%
inv-pow2.7%
pow-prod-up97.2%
metadata-eval97.2%
metadata-eval97.2%
Applied egg-rr97.2%
Final simplification89.3%
(FPCore (w l) :precision binary64 (if (<= w 0.13) (+ l (* w (- (* w (* w (* l -0.16666666666666666))) l))) (* l 0.0)))
double code(double w, double l) {
double tmp;
if (w <= 0.13) {
tmp = l + (w * ((w * (w * (l * -0.16666666666666666))) - l));
} else {
tmp = l * 0.0;
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= 0.13d0) then
tmp = l + (w * ((w * (w * (l * (-0.16666666666666666d0)))) - l))
else
tmp = l * 0.0d0
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= 0.13) {
tmp = l + (w * ((w * (w * (l * -0.16666666666666666))) - l));
} else {
tmp = l * 0.0;
}
return tmp;
}
def code(w, l): tmp = 0 if w <= 0.13: tmp = l + (w * ((w * (w * (l * -0.16666666666666666))) - l)) else: tmp = l * 0.0 return tmp
function code(w, l) tmp = 0.0 if (w <= 0.13) tmp = Float64(l + Float64(w * Float64(Float64(w * Float64(w * Float64(l * -0.16666666666666666))) - l))); else tmp = Float64(l * 0.0); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= 0.13) tmp = l + (w * ((w * (w * (l * -0.16666666666666666))) - l)); else tmp = l * 0.0; end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, 0.13], N[(l + N[(w * N[(N[(w * N[(w * N[(l * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l * 0.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq 0.13:\\
\;\;\;\;\ell + w \cdot \left(w \cdot \left(w \cdot \left(\ell \cdot -0.16666666666666666\right)\right) - \ell\right)\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot 0\\
\end{array}
\end{array}
if w < 0.13Initial program 99.6%
exp-neg99.6%
remove-double-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
remove-double-neg99.6%
Simplified99.6%
add-sqr-sqrt28.9%
sqrt-unprod78.9%
sqr-neg78.9%
sqrt-unprod50.0%
add-sqr-sqrt78.4%
add-sqr-sqrt78.4%
sqrt-unprod78.4%
add-sqr-sqrt50.0%
sqrt-unprod78.4%
sqr-neg78.4%
sqrt-unprod28.4%
add-sqr-sqrt58.4%
pow158.4%
exp-neg58.4%
inv-pow58.4%
pow-prod-up98.4%
metadata-eval98.4%
metadata-eval98.4%
metadata-eval98.4%
Applied egg-rr98.4%
Taylor expanded in w around 0 85.5%
Taylor expanded in w around inf 85.5%
distribute-rgt-out85.5%
metadata-eval85.5%
+-commutative85.5%
mul-1-neg85.5%
distribute-rgt-out85.5%
metadata-eval85.5%
fma-undefine85.5%
neg-mul-185.5%
distribute-rgt-neg-in85.5%
fma-undefine85.5%
metadata-eval85.5%
distribute-rgt-out85.5%
mul-1-neg85.5%
+-commutative85.5%
metadata-eval85.5%
Simplified85.5%
if 0.13 < w Initial program 100.0%
exp-neg100.0%
remove-double-neg100.0%
associate-*l/100.0%
*-lft-identity100.0%
remove-double-neg100.0%
Simplified100.0%
add-sqr-sqrt100.0%
sqrt-unprod100.0%
sqr-neg100.0%
sqrt-unprod0.0%
add-sqr-sqrt94.7%
add-sqr-sqrt94.7%
sqrt-unprod94.7%
add-sqr-sqrt0.0%
sqrt-unprod0.4%
sqr-neg0.4%
sqrt-unprod0.4%
add-sqr-sqrt0.4%
pow10.4%
exp-neg0.4%
inv-pow0.4%
pow-prod-up94.8%
metadata-eval94.8%
metadata-eval94.8%
metadata-eval94.8%
Applied egg-rr94.8%
Taylor expanded in w around 0 3.3%
mul-1-neg3.3%
unsub-neg3.3%
Simplified3.3%
Taylor expanded in w around inf 3.3%
neg-mul-13.3%
distribute-rgt-neg-in3.3%
Simplified3.3%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod3.4%
add-sqr-sqrt3.4%
add-log-exp1.9%
pow11.9%
pow11.9%
add-sqr-sqrt1.9%
sqrt-unprod1.9%
sqr-neg1.9%
sqrt-unprod0.0%
add-sqr-sqrt1.7%
add-sqr-sqrt1.7%
sqrt-unprod1.7%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod2.7%
add-sqr-sqrt2.7%
pow12.7%
exp-neg2.7%
inv-pow2.7%
pow-prod-up97.2%
metadata-eval97.2%
metadata-eval97.2%
Applied egg-rr97.2%
(FPCore (w l) :precision binary64 (if (<= w 0.114) (* l (- (* w (- -1.0 (* w -0.5))) -1.0)) (* l 0.0)))
double code(double w, double l) {
double tmp;
if (w <= 0.114) {
tmp = l * ((w * (-1.0 - (w * -0.5))) - -1.0);
} else {
tmp = l * 0.0;
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= 0.114d0) then
tmp = l * ((w * ((-1.0d0) - (w * (-0.5d0)))) - (-1.0d0))
else
tmp = l * 0.0d0
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= 0.114) {
tmp = l * ((w * (-1.0 - (w * -0.5))) - -1.0);
} else {
tmp = l * 0.0;
}
return tmp;
}
def code(w, l): tmp = 0 if w <= 0.114: tmp = l * ((w * (-1.0 - (w * -0.5))) - -1.0) else: tmp = l * 0.0 return tmp
function code(w, l) tmp = 0.0 if (w <= 0.114) tmp = Float64(l * Float64(Float64(w * Float64(-1.0 - Float64(w * -0.5))) - -1.0)); else tmp = Float64(l * 0.0); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= 0.114) tmp = l * ((w * (-1.0 - (w * -0.5))) - -1.0); else tmp = l * 0.0; end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, 0.114], N[(l * N[(N[(w * N[(-1.0 - N[(w * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - -1.0), $MachinePrecision]), $MachinePrecision], N[(l * 0.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq 0.114:\\
\;\;\;\;\ell \cdot \left(w \cdot \left(-1 - w \cdot -0.5\right) - -1\right)\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot 0\\
\end{array}
\end{array}
if w < 0.114000000000000004Initial program 99.6%
exp-neg99.6%
remove-double-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
remove-double-neg99.6%
Simplified99.6%
add-sqr-sqrt28.9%
sqrt-unprod78.9%
sqr-neg78.9%
sqrt-unprod50.0%
add-sqr-sqrt78.4%
add-sqr-sqrt78.4%
sqrt-unprod78.4%
add-sqr-sqrt50.0%
sqrt-unprod78.4%
sqr-neg78.4%
sqrt-unprod28.4%
add-sqr-sqrt58.4%
pow158.4%
exp-neg58.4%
inv-pow58.4%
pow-prod-up98.4%
metadata-eval98.4%
metadata-eval98.4%
metadata-eval98.4%
Applied egg-rr98.4%
frac-2neg98.4%
div-inv98.4%
*-rgt-identity98.4%
Applied egg-rr98.4%
Taylor expanded in w around 0 84.6%
if 0.114000000000000004 < w Initial program 100.0%
exp-neg100.0%
remove-double-neg100.0%
associate-*l/100.0%
*-lft-identity100.0%
remove-double-neg100.0%
Simplified100.0%
add-sqr-sqrt100.0%
sqrt-unprod100.0%
sqr-neg100.0%
sqrt-unprod0.0%
add-sqr-sqrt94.7%
add-sqr-sqrt94.7%
sqrt-unprod94.7%
add-sqr-sqrt0.0%
sqrt-unprod0.4%
sqr-neg0.4%
sqrt-unprod0.4%
add-sqr-sqrt0.4%
pow10.4%
exp-neg0.4%
inv-pow0.4%
pow-prod-up94.8%
metadata-eval94.8%
metadata-eval94.8%
metadata-eval94.8%
Applied egg-rr94.8%
Taylor expanded in w around 0 3.3%
mul-1-neg3.3%
unsub-neg3.3%
Simplified3.3%
Taylor expanded in w around inf 3.3%
neg-mul-13.3%
distribute-rgt-neg-in3.3%
Simplified3.3%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod3.4%
add-sqr-sqrt3.4%
add-log-exp1.9%
pow11.9%
pow11.9%
add-sqr-sqrt1.9%
sqrt-unprod1.9%
sqr-neg1.9%
sqrt-unprod0.0%
add-sqr-sqrt1.7%
add-sqr-sqrt1.7%
sqrt-unprod1.7%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod2.7%
add-sqr-sqrt2.7%
pow12.7%
exp-neg2.7%
inv-pow2.7%
pow-prod-up97.2%
metadata-eval97.2%
metadata-eval97.2%
Applied egg-rr97.2%
Final simplification86.3%
(FPCore (w l) :precision binary64 (if (<= w 0.122) (+ l (* w (- (* w (* l 0.5)) l))) (* l 0.0)))
double code(double w, double l) {
double tmp;
if (w <= 0.122) {
tmp = l + (w * ((w * (l * 0.5)) - l));
} else {
tmp = l * 0.0;
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= 0.122d0) then
tmp = l + (w * ((w * (l * 0.5d0)) - l))
else
tmp = l * 0.0d0
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= 0.122) {
tmp = l + (w * ((w * (l * 0.5)) - l));
} else {
tmp = l * 0.0;
}
return tmp;
}
def code(w, l): tmp = 0 if w <= 0.122: tmp = l + (w * ((w * (l * 0.5)) - l)) else: tmp = l * 0.0 return tmp
function code(w, l) tmp = 0.0 if (w <= 0.122) tmp = Float64(l + Float64(w * Float64(Float64(w * Float64(l * 0.5)) - l))); else tmp = Float64(l * 0.0); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= 0.122) tmp = l + (w * ((w * (l * 0.5)) - l)); else tmp = l * 0.0; end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, 0.122], N[(l + N[(w * N[(N[(w * N[(l * 0.5), $MachinePrecision]), $MachinePrecision] - l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l * 0.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq 0.122:\\
\;\;\;\;\ell + w \cdot \left(w \cdot \left(\ell \cdot 0.5\right) - \ell\right)\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot 0\\
\end{array}
\end{array}
if w < 0.122Initial program 99.6%
exp-neg99.6%
remove-double-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
remove-double-neg99.6%
Simplified99.6%
add-sqr-sqrt28.9%
sqrt-unprod78.9%
sqr-neg78.9%
sqrt-unprod50.0%
add-sqr-sqrt78.4%
add-sqr-sqrt78.4%
sqrt-unprod78.4%
add-sqr-sqrt50.0%
sqrt-unprod78.4%
sqr-neg78.4%
sqrt-unprod28.4%
add-sqr-sqrt58.4%
pow158.4%
exp-neg58.4%
inv-pow58.4%
pow-prod-up98.4%
metadata-eval98.4%
metadata-eval98.4%
metadata-eval98.4%
Applied egg-rr98.4%
frac-2neg98.4%
div-inv98.4%
*-rgt-identity98.4%
Applied egg-rr98.4%
Taylor expanded in w around 0 80.3%
mul-1-neg80.3%
distribute-rgt-out80.3%
metadata-eval80.3%
distribute-rgt-neg-out80.3%
distribute-rgt-neg-in80.3%
metadata-eval80.3%
Simplified80.3%
if 0.122 < w Initial program 100.0%
exp-neg100.0%
remove-double-neg100.0%
associate-*l/100.0%
*-lft-identity100.0%
remove-double-neg100.0%
Simplified100.0%
add-sqr-sqrt100.0%
sqrt-unprod100.0%
sqr-neg100.0%
sqrt-unprod0.0%
add-sqr-sqrt94.7%
add-sqr-sqrt94.7%
sqrt-unprod94.7%
add-sqr-sqrt0.0%
sqrt-unprod0.4%
sqr-neg0.4%
sqrt-unprod0.4%
add-sqr-sqrt0.4%
pow10.4%
exp-neg0.4%
inv-pow0.4%
pow-prod-up94.8%
metadata-eval94.8%
metadata-eval94.8%
metadata-eval94.8%
Applied egg-rr94.8%
Taylor expanded in w around 0 3.3%
mul-1-neg3.3%
unsub-neg3.3%
Simplified3.3%
Taylor expanded in w around inf 3.3%
neg-mul-13.3%
distribute-rgt-neg-in3.3%
Simplified3.3%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod3.4%
add-sqr-sqrt3.4%
add-log-exp1.9%
pow11.9%
pow11.9%
add-sqr-sqrt1.9%
sqrt-unprod1.9%
sqr-neg1.9%
sqrt-unprod0.0%
add-sqr-sqrt1.7%
add-sqr-sqrt1.7%
sqrt-unprod1.7%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod2.7%
add-sqr-sqrt2.7%
pow12.7%
exp-neg2.7%
inv-pow2.7%
pow-prod-up97.2%
metadata-eval97.2%
metadata-eval97.2%
Applied egg-rr97.2%
(FPCore (w l) :precision binary64 (if (<= w -202.0) (* l (- w)) (if (<= w 0.135) l (* l 0.0))))
double code(double w, double l) {
double tmp;
if (w <= -202.0) {
tmp = l * -w;
} else if (w <= 0.135) {
tmp = l;
} else {
tmp = l * 0.0;
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= (-202.0d0)) then
tmp = l * -w
else if (w <= 0.135d0) then
tmp = l
else
tmp = l * 0.0d0
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= -202.0) {
tmp = l * -w;
} else if (w <= 0.135) {
tmp = l;
} else {
tmp = l * 0.0;
}
return tmp;
}
def code(w, l): tmp = 0 if w <= -202.0: tmp = l * -w elif w <= 0.135: tmp = l else: tmp = l * 0.0 return tmp
function code(w, l) tmp = 0.0 if (w <= -202.0) tmp = Float64(l * Float64(-w)); elseif (w <= 0.135) tmp = l; else tmp = Float64(l * 0.0); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= -202.0) tmp = l * -w; elseif (w <= 0.135) tmp = l; else tmp = l * 0.0; end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, -202.0], N[(l * (-w)), $MachinePrecision], If[LessEqual[w, 0.135], l, N[(l * 0.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq -202:\\
\;\;\;\;\ell \cdot \left(-w\right)\\
\mathbf{elif}\;w \leq 0.135:\\
\;\;\;\;\ell\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot 0\\
\end{array}
\end{array}
if w < -202Initial program 100.0%
exp-neg100.0%
remove-double-neg100.0%
associate-*l/100.0%
*-lft-identity100.0%
remove-double-neg100.0%
Simplified100.0%
add-sqr-sqrt0.0%
sqrt-unprod50.0%
sqr-neg50.0%
sqrt-unprod50.0%
add-sqr-sqrt50.0%
add-sqr-sqrt50.0%
sqrt-unprod50.0%
add-sqr-sqrt50.0%
sqrt-unprod50.0%
sqr-neg50.0%
sqrt-unprod0.0%
add-sqr-sqrt0.0%
pow10.0%
exp-neg0.0%
inv-pow0.0%
pow-prod-up100.0%
metadata-eval100.0%
metadata-eval100.0%
metadata-eval100.0%
Applied egg-rr100.0%
Taylor expanded in w around 0 28.9%
mul-1-neg28.9%
unsub-neg28.9%
Simplified28.9%
Taylor expanded in w around inf 28.9%
neg-mul-128.9%
distribute-rgt-neg-in28.9%
Simplified28.9%
if -202 < w < 0.13500000000000001Initial program 99.3%
Taylor expanded in w around 0 97.4%
if 0.13500000000000001 < w Initial program 100.0%
exp-neg100.0%
remove-double-neg100.0%
associate-*l/100.0%
*-lft-identity100.0%
remove-double-neg100.0%
Simplified100.0%
add-sqr-sqrt100.0%
sqrt-unprod100.0%
sqr-neg100.0%
sqrt-unprod0.0%
add-sqr-sqrt94.7%
add-sqr-sqrt94.7%
sqrt-unprod94.7%
add-sqr-sqrt0.0%
sqrt-unprod0.4%
sqr-neg0.4%
sqrt-unprod0.4%
add-sqr-sqrt0.4%
pow10.4%
exp-neg0.4%
inv-pow0.4%
pow-prod-up94.8%
metadata-eval94.8%
metadata-eval94.8%
metadata-eval94.8%
Applied egg-rr94.8%
Taylor expanded in w around 0 3.3%
mul-1-neg3.3%
unsub-neg3.3%
Simplified3.3%
Taylor expanded in w around inf 3.3%
neg-mul-13.3%
distribute-rgt-neg-in3.3%
Simplified3.3%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod3.4%
add-sqr-sqrt3.4%
add-log-exp1.9%
pow11.9%
pow11.9%
add-sqr-sqrt1.9%
sqrt-unprod1.9%
sqr-neg1.9%
sqrt-unprod0.0%
add-sqr-sqrt1.7%
add-sqr-sqrt1.7%
sqrt-unprod1.7%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod2.7%
add-sqr-sqrt2.7%
pow12.7%
exp-neg2.7%
inv-pow2.7%
pow-prod-up97.2%
metadata-eval97.2%
metadata-eval97.2%
Applied egg-rr97.2%
(FPCore (w l) :precision binary64 (if (<= w 0.135) (- l (* l w)) (* l 0.0)))
double code(double w, double l) {
double tmp;
if (w <= 0.135) {
tmp = l - (l * w);
} else {
tmp = l * 0.0;
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= 0.135d0) then
tmp = l - (l * w)
else
tmp = l * 0.0d0
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= 0.135) {
tmp = l - (l * w);
} else {
tmp = l * 0.0;
}
return tmp;
}
def code(w, l): tmp = 0 if w <= 0.135: tmp = l - (l * w) else: tmp = l * 0.0 return tmp
function code(w, l) tmp = 0.0 if (w <= 0.135) tmp = Float64(l - Float64(l * w)); else tmp = Float64(l * 0.0); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= 0.135) tmp = l - (l * w); else tmp = l * 0.0; end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, 0.135], N[(l - N[(l * w), $MachinePrecision]), $MachinePrecision], N[(l * 0.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq 0.135:\\
\;\;\;\;\ell - \ell \cdot w\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot 0\\
\end{array}
\end{array}
if w < 0.13500000000000001Initial program 99.6%
exp-neg99.6%
remove-double-neg99.6%
associate-*l/99.6%
*-lft-identity99.6%
remove-double-neg99.6%
Simplified99.6%
add-sqr-sqrt28.9%
sqrt-unprod78.9%
sqr-neg78.9%
sqrt-unprod50.0%
add-sqr-sqrt78.4%
add-sqr-sqrt78.4%
sqrt-unprod78.4%
add-sqr-sqrt50.0%
sqrt-unprod78.4%
sqr-neg78.4%
sqrt-unprod28.4%
add-sqr-sqrt58.4%
pow158.4%
exp-neg58.4%
inv-pow58.4%
pow-prod-up98.4%
metadata-eval98.4%
metadata-eval98.4%
metadata-eval98.4%
Applied egg-rr98.4%
Taylor expanded in w around 0 70.0%
mul-1-neg70.0%
unsub-neg70.0%
Simplified70.0%
if 0.13500000000000001 < w Initial program 100.0%
exp-neg100.0%
remove-double-neg100.0%
associate-*l/100.0%
*-lft-identity100.0%
remove-double-neg100.0%
Simplified100.0%
add-sqr-sqrt100.0%
sqrt-unprod100.0%
sqr-neg100.0%
sqrt-unprod0.0%
add-sqr-sqrt94.7%
add-sqr-sqrt94.7%
sqrt-unprod94.7%
add-sqr-sqrt0.0%
sqrt-unprod0.4%
sqr-neg0.4%
sqrt-unprod0.4%
add-sqr-sqrt0.4%
pow10.4%
exp-neg0.4%
inv-pow0.4%
pow-prod-up94.8%
metadata-eval94.8%
metadata-eval94.8%
metadata-eval94.8%
Applied egg-rr94.8%
Taylor expanded in w around 0 3.3%
mul-1-neg3.3%
unsub-neg3.3%
Simplified3.3%
Taylor expanded in w around inf 3.3%
neg-mul-13.3%
distribute-rgt-neg-in3.3%
Simplified3.3%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod3.4%
add-sqr-sqrt3.4%
add-log-exp1.9%
pow11.9%
pow11.9%
add-sqr-sqrt1.9%
sqrt-unprod1.9%
sqr-neg1.9%
sqrt-unprod0.0%
add-sqr-sqrt1.7%
add-sqr-sqrt1.7%
sqrt-unprod1.7%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod2.7%
add-sqr-sqrt2.7%
pow12.7%
exp-neg2.7%
inv-pow2.7%
pow-prod-up97.2%
metadata-eval97.2%
metadata-eval97.2%
Applied egg-rr97.2%
(FPCore (w l) :precision binary64 (if (<= w 0.135) l (* l 0.0)))
double code(double w, double l) {
double tmp;
if (w <= 0.135) {
tmp = l;
} else {
tmp = l * 0.0;
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= 0.135d0) then
tmp = l
else
tmp = l * 0.0d0
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= 0.135) {
tmp = l;
} else {
tmp = l * 0.0;
}
return tmp;
}
def code(w, l): tmp = 0 if w <= 0.135: tmp = l else: tmp = l * 0.0 return tmp
function code(w, l) tmp = 0.0 if (w <= 0.135) tmp = l; else tmp = Float64(l * 0.0); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= 0.135) tmp = l; else tmp = l * 0.0; end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, 0.135], l, N[(l * 0.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq 0.135:\\
\;\;\;\;\ell\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot 0\\
\end{array}
\end{array}
if w < 0.13500000000000001Initial program 99.6%
Taylor expanded in w around 0 60.1%
if 0.13500000000000001 < w Initial program 100.0%
exp-neg100.0%
remove-double-neg100.0%
associate-*l/100.0%
*-lft-identity100.0%
remove-double-neg100.0%
Simplified100.0%
add-sqr-sqrt100.0%
sqrt-unprod100.0%
sqr-neg100.0%
sqrt-unprod0.0%
add-sqr-sqrt94.7%
add-sqr-sqrt94.7%
sqrt-unprod94.7%
add-sqr-sqrt0.0%
sqrt-unprod0.4%
sqr-neg0.4%
sqrt-unprod0.4%
add-sqr-sqrt0.4%
pow10.4%
exp-neg0.4%
inv-pow0.4%
pow-prod-up94.8%
metadata-eval94.8%
metadata-eval94.8%
metadata-eval94.8%
Applied egg-rr94.8%
Taylor expanded in w around 0 3.3%
mul-1-neg3.3%
unsub-neg3.3%
Simplified3.3%
Taylor expanded in w around inf 3.3%
neg-mul-13.3%
distribute-rgt-neg-in3.3%
Simplified3.3%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod3.4%
add-sqr-sqrt3.4%
add-log-exp1.9%
pow11.9%
pow11.9%
add-sqr-sqrt1.9%
sqrt-unprod1.9%
sqr-neg1.9%
sqrt-unprod0.0%
add-sqr-sqrt1.7%
add-sqr-sqrt1.7%
sqrt-unprod1.7%
add-sqr-sqrt0.0%
sqrt-unprod2.7%
sqr-neg2.7%
sqrt-unprod2.7%
add-sqr-sqrt2.7%
pow12.7%
exp-neg2.7%
inv-pow2.7%
pow-prod-up97.2%
metadata-eval97.2%
metadata-eval97.2%
Applied egg-rr97.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%
Taylor expanded in w around 0 52.4%
herbie shell --seed 2024136
(FPCore (w l)
:name "exp-w (used to crash)"
:precision binary64
(* (exp (- w)) (pow l (exp w))))