
(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 14 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.7%
exp-neg99.7%
remove-double-neg99.7%
associate-*l/99.7%
*-lft-identity99.7%
remove-double-neg99.7%
Simplified99.7%
(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(Float64(-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}
\\
\ell \cdot e^{-w}
\end{array}
Initial program 99.7%
add-sqr-sqrt41.3%
sqrt-unprod85.0%
sqr-neg85.0%
sqrt-unprod43.7%
add-sqr-sqrt84.7%
add-sqr-sqrt84.7%
sqrt-unprod84.7%
add-sqr-sqrt43.7%
sqrt-unprod68.7%
sqr-neg68.7%
sqrt-unprod25.0%
add-sqr-sqrt52.0%
pow152.0%
exp-neg52.0%
inv-pow52.0%
pow-prod-up98.4%
metadata-eval98.4%
metadata-eval98.4%
metadata-eval98.4%
Applied egg-rr98.4%
Final simplification98.4%
(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.7%
add-sqr-sqrt41.3%
sqrt-unprod85.0%
sqr-neg85.0%
sqrt-unprod43.7%
add-sqr-sqrt84.7%
add-sqr-sqrt84.7%
sqrt-unprod84.7%
add-sqr-sqrt43.7%
sqrt-unprod68.7%
sqr-neg68.7%
sqrt-unprod25.0%
add-sqr-sqrt52.0%
pow152.0%
exp-neg52.0%
inv-pow52.0%
pow-prod-up98.4%
metadata-eval98.4%
metadata-eval98.4%
metadata-eval98.4%
Applied egg-rr98.4%
Taylor expanded in w around inf 98.4%
rem-exp-log94.0%
prod-exp94.0%
sub-neg94.0%
div-exp94.0%
rem-exp-log98.4%
Simplified98.4%
(FPCore (w l)
:precision binary64
(let* ((t_0 (* w (* l (- (* w -0.16666666666666666) 0.5))))
(t_1 (+ 0.5 (* w -0.16666666666666666))))
(if (<= w -2.95e+54)
(+ l (* l (* w (+ (* w t_1) -1.0))))
(if (<= w -0.24)
(+ l (* w (- (/ (* (* w (* l t_1)) t_0) t_0) l)))
(/ l (+ 1.0 (* w (+ 1.0 (* w (+ 0.5 (* w 0.16666666666666666)))))))))))
double code(double w, double l) {
double t_0 = w * (l * ((w * -0.16666666666666666) - 0.5));
double t_1 = 0.5 + (w * -0.16666666666666666);
double tmp;
if (w <= -2.95e+54) {
tmp = l + (l * (w * ((w * t_1) + -1.0)));
} else if (w <= -0.24) {
tmp = l + (w * ((((w * (l * t_1)) * t_0) / t_0) - l));
} else {
tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666))))));
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = w * (l * ((w * (-0.16666666666666666d0)) - 0.5d0))
t_1 = 0.5d0 + (w * (-0.16666666666666666d0))
if (w <= (-2.95d+54)) then
tmp = l + (l * (w * ((w * t_1) + (-1.0d0))))
else if (w <= (-0.24d0)) then
tmp = l + (w * ((((w * (l * t_1)) * t_0) / t_0) - l))
else
tmp = l / (1.0d0 + (w * (1.0d0 + (w * (0.5d0 + (w * 0.16666666666666666d0))))))
end if
code = tmp
end function
public static double code(double w, double l) {
double t_0 = w * (l * ((w * -0.16666666666666666) - 0.5));
double t_1 = 0.5 + (w * -0.16666666666666666);
double tmp;
if (w <= -2.95e+54) {
tmp = l + (l * (w * ((w * t_1) + -1.0)));
} else if (w <= -0.24) {
tmp = l + (w * ((((w * (l * t_1)) * t_0) / t_0) - l));
} else {
tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666))))));
}
return tmp;
}
def code(w, l): t_0 = w * (l * ((w * -0.16666666666666666) - 0.5)) t_1 = 0.5 + (w * -0.16666666666666666) tmp = 0 if w <= -2.95e+54: tmp = l + (l * (w * ((w * t_1) + -1.0))) elif w <= -0.24: tmp = l + (w * ((((w * (l * t_1)) * t_0) / t_0) - l)) else: tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666)))))) return tmp
function code(w, l) t_0 = Float64(w * Float64(l * Float64(Float64(w * -0.16666666666666666) - 0.5))) t_1 = Float64(0.5 + Float64(w * -0.16666666666666666)) tmp = 0.0 if (w <= -2.95e+54) tmp = Float64(l + Float64(l * Float64(w * Float64(Float64(w * t_1) + -1.0)))); elseif (w <= -0.24) tmp = Float64(l + Float64(w * Float64(Float64(Float64(Float64(w * Float64(l * t_1)) * t_0) / t_0) - l))); else tmp = Float64(l / Float64(1.0 + Float64(w * Float64(1.0 + Float64(w * Float64(0.5 + Float64(w * 0.16666666666666666))))))); end return tmp end
function tmp_2 = code(w, l) t_0 = w * (l * ((w * -0.16666666666666666) - 0.5)); t_1 = 0.5 + (w * -0.16666666666666666); tmp = 0.0; if (w <= -2.95e+54) tmp = l + (l * (w * ((w * t_1) + -1.0))); elseif (w <= -0.24) tmp = l + (w * ((((w * (l * t_1)) * t_0) / t_0) - l)); else tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666)))))); end tmp_2 = tmp; end
code[w_, l_] := Block[{t$95$0 = N[(w * N[(l * N[(N[(w * -0.16666666666666666), $MachinePrecision] - 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(0.5 + N[(w * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[w, -2.95e+54], N[(l + N[(l * N[(w * N[(N[(w * t$95$1), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[w, -0.24], N[(l + N[(w * N[(N[(N[(N[(w * N[(l * t$95$1), $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision] / t$95$0), $MachinePrecision] - l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l / N[(1.0 + N[(w * N[(1.0 + N[(w * N[(0.5 + N[(w * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := w \cdot \left(\ell \cdot \left(w \cdot -0.16666666666666666 - 0.5\right)\right)\\
t_1 := 0.5 + w \cdot -0.16666666666666666\\
\mathbf{if}\;w \leq -2.95 \cdot 10^{+54}:\\
\;\;\;\;\ell + \ell \cdot \left(w \cdot \left(w \cdot t\_1 + -1\right)\right)\\
\mathbf{elif}\;w \leq -0.24:\\
\;\;\;\;\ell + w \cdot \left(\frac{\left(w \cdot \left(\ell \cdot t\_1\right)\right) \cdot t\_0}{t\_0} - \ell\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{1 + w \cdot \left(1 + w \cdot \left(0.5 + w \cdot 0.16666666666666666\right)\right)}\\
\end{array}
\end{array}
if w < -2.9499999999999999e54Initial program 100.0%
add-sqr-sqrt0.0%
sqrt-unprod55.4%
sqr-neg55.4%
sqrt-unprod55.4%
add-sqr-sqrt55.4%
add-sqr-sqrt55.4%
sqrt-unprod55.4%
add-sqr-sqrt55.4%
sqrt-unprod55.4%
sqr-neg55.4%
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 81.1%
Taylor expanded in l around 0 86.9%
if -2.9499999999999999e54 < w < -0.23999999999999999Initial program 99.1%
add-sqr-sqrt0.0%
sqrt-unprod47.4%
sqr-neg47.4%
sqrt-unprod47.4%
add-sqr-sqrt47.4%
add-sqr-sqrt47.4%
sqrt-unprod47.4%
add-sqr-sqrt47.4%
sqrt-unprod47.4%
sqr-neg47.4%
sqrt-unprod0.0%
add-sqr-sqrt1.0%
pow11.0%
exp-neg1.0%
inv-pow1.0%
pow-prod-up87.7%
metadata-eval87.7%
metadata-eval87.7%
metadata-eval87.7%
Applied egg-rr87.7%
Taylor expanded in w around 0 11.7%
distribute-lft-in11.7%
flip-+15.9%
*-commutative15.9%
associate-*l*15.9%
*-commutative15.9%
associate-*l*15.9%
*-commutative15.9%
*-commutative15.9%
*-commutative15.9%
associate-*l*15.9%
*-commutative15.9%
Applied egg-rr15.9%
difference-of-squares42.5%
distribute-lft-out42.5%
+-commutative42.5%
distribute-lft-in42.5%
distribute-lft-out--42.5%
distribute-lft-out--42.5%
distribute-lft-out--42.5%
distribute-lft-out--42.5%
Simplified42.5%
if -0.23999999999999999 < w Initial program 99.6%
add-sqr-sqrt60.0%
sqrt-unprod99.1%
sqr-neg99.1%
sqrt-unprod39.0%
add-sqr-sqrt98.7%
add-sqr-sqrt98.7%
sqrt-unprod98.7%
add-sqr-sqrt39.0%
sqrt-unprod75.5%
sqr-neg75.5%
sqrt-unprod36.4%
add-sqr-sqrt75.5%
pow175.5%
exp-neg75.5%
inv-pow75.5%
pow-prod-up98.8%
metadata-eval98.8%
metadata-eval98.8%
metadata-eval98.8%
Applied egg-rr98.8%
Taylor expanded in w around inf 98.8%
rem-exp-log92.4%
prod-exp92.4%
sub-neg92.4%
div-exp92.4%
rem-exp-log98.8%
Simplified98.8%
Taylor expanded in w around 0 95.4%
*-commutative95.4%
Simplified95.4%
Final simplification90.1%
(FPCore (w l) :precision binary64 (if (<= w 2e-7) (* l (+ 1.0 (* w (+ (* w (+ 0.5 (* w -0.16666666666666666))) -1.0)))) (/ l (+ 1.0 (* w (+ 1.0 (* w (+ 0.5 (* w 0.16666666666666666)))))))))
double code(double w, double l) {
double tmp;
if (w <= 2e-7) {
tmp = l * (1.0 + (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0)));
} else {
tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666))))));
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= 2d-7) then
tmp = l * (1.0d0 + (w * ((w * (0.5d0 + (w * (-0.16666666666666666d0)))) + (-1.0d0))))
else
tmp = l / (1.0d0 + (w * (1.0d0 + (w * (0.5d0 + (w * 0.16666666666666666d0))))))
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= 2e-7) {
tmp = l * (1.0 + (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0)));
} else {
tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666))))));
}
return tmp;
}
def code(w, l): tmp = 0 if w <= 2e-7: tmp = l * (1.0 + (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0))) else: tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666)))))) return tmp
function code(w, l) tmp = 0.0 if (w <= 2e-7) tmp = Float64(l * Float64(1.0 + Float64(w * Float64(Float64(w * Float64(0.5 + Float64(w * -0.16666666666666666))) + -1.0)))); else tmp = Float64(l / Float64(1.0 + Float64(w * Float64(1.0 + Float64(w * Float64(0.5 + Float64(w * 0.16666666666666666))))))); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= 2e-7) tmp = l * (1.0 + (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0))); else tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666)))))); end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, 2e-7], N[(l * N[(1.0 + N[(w * N[(N[(w * N[(0.5 + N[(w * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l / N[(1.0 + N[(w * N[(1.0 + N[(w * N[(0.5 + N[(w * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq 2 \cdot 10^{-7}:\\
\;\;\;\;\ell \cdot \left(1 + w \cdot \left(w \cdot \left(0.5 + w \cdot -0.16666666666666666\right) + -1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{1 + w \cdot \left(1 + w \cdot \left(0.5 + w \cdot 0.16666666666666666\right)\right)}\\
\end{array}
\end{array}
if w < 1.9999999999999999e-7Initial program 99.6%
add-sqr-sqrt30.1%
sqrt-unprod82.1%
sqr-neg82.1%
sqrt-unprod52.0%
add-sqr-sqrt81.8%
add-sqr-sqrt81.8%
sqrt-unprod81.8%
add-sqr-sqrt52.0%
sqrt-unprod81.8%
sqr-neg81.8%
sqrt-unprod29.8%
add-sqr-sqrt61.9%
pow161.9%
exp-neg61.9%
inv-pow61.9%
pow-prod-up98.1%
metadata-eval98.1%
metadata-eval98.1%
metadata-eval98.1%
Applied egg-rr98.1%
Taylor expanded in w around 0 88.9%
if 1.9999999999999999e-7 < w Initial program 100.0%
add-sqr-sqrt100.0%
sqrt-unprod100.0%
sqr-neg100.0%
sqrt-unprod0.0%
add-sqr-sqrt100.0%
add-sqr-sqrt100.0%
sqrt-unprod100.0%
add-sqr-sqrt0.0%
sqrt-unprod0.0%
sqr-neg0.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 inf 100.0%
rem-exp-log100.0%
prod-exp100.0%
sub-neg100.0%
div-exp100.0%
rem-exp-log100.0%
Simplified100.0%
Taylor expanded in w around 0 85.5%
*-commutative85.5%
Simplified85.5%
Final simplification88.3%
(FPCore (w l) :precision binary64 (if (<= w 2e-7) (+ l (* l (* w (+ (* w (+ 0.5 (* w -0.16666666666666666))) -1.0)))) (/ l (+ 1.0 (* w (+ 1.0 (* w (+ 0.5 (* w 0.16666666666666666)))))))))
double code(double w, double l) {
double tmp;
if (w <= 2e-7) {
tmp = l + (l * (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0)));
} else {
tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666))))));
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= 2d-7) then
tmp = l + (l * (w * ((w * (0.5d0 + (w * (-0.16666666666666666d0)))) + (-1.0d0))))
else
tmp = l / (1.0d0 + (w * (1.0d0 + (w * (0.5d0 + (w * 0.16666666666666666d0))))))
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= 2e-7) {
tmp = l + (l * (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0)));
} else {
tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666))))));
}
return tmp;
}
def code(w, l): tmp = 0 if w <= 2e-7: tmp = l + (l * (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0))) else: tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666)))))) return tmp
function code(w, l) tmp = 0.0 if (w <= 2e-7) tmp = Float64(l + Float64(l * Float64(w * Float64(Float64(w * Float64(0.5 + Float64(w * -0.16666666666666666))) + -1.0)))); else tmp = Float64(l / Float64(1.0 + Float64(w * Float64(1.0 + Float64(w * Float64(0.5 + Float64(w * 0.16666666666666666))))))); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= 2e-7) tmp = l + (l * (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0))); else tmp = l / (1.0 + (w * (1.0 + (w * (0.5 + (w * 0.16666666666666666)))))); end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, 2e-7], N[(l + N[(l * N[(w * N[(N[(w * N[(0.5 + N[(w * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l / N[(1.0 + N[(w * N[(1.0 + N[(w * N[(0.5 + N[(w * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq 2 \cdot 10^{-7}:\\
\;\;\;\;\ell + \ell \cdot \left(w \cdot \left(w \cdot \left(0.5 + w \cdot -0.16666666666666666\right) + -1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{1 + w \cdot \left(1 + w \cdot \left(0.5 + w \cdot 0.16666666666666666\right)\right)}\\
\end{array}
\end{array}
if w < 1.9999999999999999e-7Initial program 99.6%
add-sqr-sqrt30.1%
sqrt-unprod82.1%
sqr-neg82.1%
sqrt-unprod52.0%
add-sqr-sqrt81.8%
add-sqr-sqrt81.8%
sqrt-unprod81.8%
add-sqr-sqrt52.0%
sqrt-unprod81.8%
sqr-neg81.8%
sqrt-unprod29.8%
add-sqr-sqrt61.9%
pow161.9%
exp-neg61.9%
inv-pow61.9%
pow-prod-up98.1%
metadata-eval98.1%
metadata-eval98.1%
metadata-eval98.1%
Applied egg-rr98.1%
Taylor expanded in w around 0 87.1%
Taylor expanded in l around 0 88.9%
if 1.9999999999999999e-7 < w Initial program 100.0%
add-sqr-sqrt100.0%
sqrt-unprod100.0%
sqr-neg100.0%
sqrt-unprod0.0%
add-sqr-sqrt100.0%
add-sqr-sqrt100.0%
sqrt-unprod100.0%
add-sqr-sqrt0.0%
sqrt-unprod0.0%
sqr-neg0.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 inf 100.0%
rem-exp-log100.0%
prod-exp100.0%
sub-neg100.0%
div-exp100.0%
rem-exp-log100.0%
Simplified100.0%
Taylor expanded in w around 0 85.5%
*-commutative85.5%
Simplified85.5%
Final simplification88.3%
(FPCore (w l) :precision binary64 (if (<= w -12.2) (+ l (* l (* w (+ (* w (+ 0.5 (* w -0.16666666666666666))) -1.0)))) (/ l (+ 1.0 (* w (+ 1.0 (* w 0.5)))))))
double code(double w, double l) {
double tmp;
if (w <= -12.2) {
tmp = l + (l * (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0)));
} else {
tmp = l / (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 <= (-12.2d0)) then
tmp = l + (l * (w * ((w * (0.5d0 + (w * (-0.16666666666666666d0)))) + (-1.0d0))))
else
tmp = l / (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 <= -12.2) {
tmp = l + (l * (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0)));
} else {
tmp = l / (1.0 + (w * (1.0 + (w * 0.5))));
}
return tmp;
}
def code(w, l): tmp = 0 if w <= -12.2: tmp = l + (l * (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0))) else: tmp = l / (1.0 + (w * (1.0 + (w * 0.5)))) return tmp
function code(w, l) tmp = 0.0 if (w <= -12.2) tmp = Float64(l + Float64(l * Float64(w * Float64(Float64(w * Float64(0.5 + Float64(w * -0.16666666666666666))) + -1.0)))); else tmp = Float64(l / 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 <= -12.2) tmp = l + (l * (w * ((w * (0.5 + (w * -0.16666666666666666))) + -1.0))); else tmp = l / (1.0 + (w * (1.0 + (w * 0.5)))); end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, -12.2], N[(l + N[(l * N[(w * N[(N[(w * N[(0.5 + N[(w * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l / 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 -12.2:\\
\;\;\;\;\ell + \ell \cdot \left(w \cdot \left(w \cdot \left(0.5 + w \cdot -0.16666666666666666\right) + -1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{1 + w \cdot \left(1 + w \cdot 0.5\right)}\\
\end{array}
\end{array}
if w < -12.199999999999999Initial program 100.0%
add-sqr-sqrt0.0%
sqrt-unprod55.1%
sqr-neg55.1%
sqrt-unprod55.1%
add-sqr-sqrt55.1%
add-sqr-sqrt55.1%
sqrt-unprod55.1%
add-sqr-sqrt55.1%
sqrt-unprod55.1%
sqr-neg55.1%
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 69.6%
Taylor expanded in l around 0 74.4%
if -12.199999999999999 < w Initial program 99.6%
add-sqr-sqrt59.4%
sqrt-unprod98.0%
sqr-neg98.0%
sqrt-unprod38.7%
add-sqr-sqrt97.7%
add-sqr-sqrt97.7%
sqrt-unprod97.7%
add-sqr-sqrt38.7%
sqrt-unprod74.7%
sqr-neg74.7%
sqrt-unprod36.0%
add-sqr-sqrt74.7%
pow174.7%
exp-neg74.7%
inv-pow74.7%
pow-prod-up97.8%
metadata-eval97.8%
metadata-eval97.8%
metadata-eval97.8%
Applied egg-rr97.8%
Taylor expanded in w around inf 97.8%
rem-exp-log91.4%
prod-exp91.4%
sub-neg91.4%
div-exp91.4%
rem-exp-log97.8%
Simplified97.8%
Taylor expanded in w around 0 93.5%
*-commutative93.5%
Simplified93.5%
Final simplification87.7%
(FPCore (w l) :precision binary64 (if (<= w -1.4) (* l (+ 1.0 (* w (+ -1.0 (* w 0.5))))) (/ l (+ 1.0 (* w (+ 1.0 (* w 0.5)))))))
double code(double w, double l) {
double tmp;
if (w <= -1.4) {
tmp = l * (1.0 + (w * (-1.0 + (w * 0.5))));
} else {
tmp = l / (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 <= (-1.4d0)) then
tmp = l * (1.0d0 + (w * ((-1.0d0) + (w * 0.5d0))))
else
tmp = l / (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 <= -1.4) {
tmp = l * (1.0 + (w * (-1.0 + (w * 0.5))));
} else {
tmp = l / (1.0 + (w * (1.0 + (w * 0.5))));
}
return tmp;
}
def code(w, l): tmp = 0 if w <= -1.4: tmp = l * (1.0 + (w * (-1.0 + (w * 0.5)))) else: tmp = l / (1.0 + (w * (1.0 + (w * 0.5)))) return tmp
function code(w, l) tmp = 0.0 if (w <= -1.4) tmp = Float64(l * Float64(1.0 + Float64(w * Float64(-1.0 + Float64(w * 0.5))))); else tmp = Float64(l / 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 <= -1.4) tmp = l * (1.0 + (w * (-1.0 + (w * 0.5)))); else tmp = l / (1.0 + (w * (1.0 + (w * 0.5)))); end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, -1.4], N[(l * N[(1.0 + N[(w * N[(-1.0 + N[(w * 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l / 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 -1.4:\\
\;\;\;\;\ell \cdot \left(1 + w \cdot \left(-1 + w \cdot 0.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{1 + w \cdot \left(1 + w \cdot 0.5\right)}\\
\end{array}
\end{array}
if w < -1.3999999999999999Initial program 100.0%
add-sqr-sqrt0.0%
sqrt-unprod55.1%
sqr-neg55.1%
sqrt-unprod55.1%
add-sqr-sqrt55.1%
add-sqr-sqrt55.1%
sqrt-unprod55.1%
add-sqr-sqrt55.1%
sqrt-unprod55.1%
sqr-neg55.1%
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 64.5%
if -1.3999999999999999 < w Initial program 99.6%
add-sqr-sqrt59.4%
sqrt-unprod98.0%
sqr-neg98.0%
sqrt-unprod38.7%
add-sqr-sqrt97.7%
add-sqr-sqrt97.7%
sqrt-unprod97.7%
add-sqr-sqrt38.7%
sqrt-unprod74.7%
sqr-neg74.7%
sqrt-unprod36.0%
add-sqr-sqrt74.7%
pow174.7%
exp-neg74.7%
inv-pow74.7%
pow-prod-up97.8%
metadata-eval97.8%
metadata-eval97.8%
metadata-eval97.8%
Applied egg-rr97.8%
Taylor expanded in w around inf 97.8%
rem-exp-log91.4%
prod-exp91.4%
sub-neg91.4%
div-exp91.4%
rem-exp-log97.8%
Simplified97.8%
Taylor expanded in w around 0 93.5%
*-commutative93.5%
Simplified93.5%
Final simplification84.7%
(FPCore (w l) :precision binary64 (if (<= w -7.3) (- l (* w (+ l (* w (* l -0.5))))) (/ l (+ 1.0 (* w (+ 1.0 (* w 0.5)))))))
double code(double w, double l) {
double tmp;
if (w <= -7.3) {
tmp = l - (w * (l + (w * (l * -0.5))));
} else {
tmp = l / (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 <= (-7.3d0)) then
tmp = l - (w * (l + (w * (l * (-0.5d0)))))
else
tmp = l / (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 <= -7.3) {
tmp = l - (w * (l + (w * (l * -0.5))));
} else {
tmp = l / (1.0 + (w * (1.0 + (w * 0.5))));
}
return tmp;
}
def code(w, l): tmp = 0 if w <= -7.3: tmp = l - (w * (l + (w * (l * -0.5)))) else: tmp = l / (1.0 + (w * (1.0 + (w * 0.5)))) return tmp
function code(w, l) tmp = 0.0 if (w <= -7.3) tmp = Float64(l - Float64(w * Float64(l + Float64(w * Float64(l * -0.5))))); else tmp = Float64(l / 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 <= -7.3) tmp = l - (w * (l + (w * (l * -0.5)))); else tmp = l / (1.0 + (w * (1.0 + (w * 0.5)))); end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, -7.3], N[(l - N[(w * N[(l + N[(w * N[(l * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(l / 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 -7.3:\\
\;\;\;\;\ell - w \cdot \left(\ell + w \cdot \left(\ell \cdot -0.5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{1 + w \cdot \left(1 + w \cdot 0.5\right)}\\
\end{array}
\end{array}
if w < -7.29999999999999982Initial program 100.0%
add-sqr-sqrt0.0%
sqrt-unprod55.1%
sqr-neg55.1%
sqrt-unprod55.1%
add-sqr-sqrt55.1%
add-sqr-sqrt55.1%
sqrt-unprod55.1%
add-sqr-sqrt55.1%
sqrt-unprod55.1%
sqr-neg55.1%
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 inf 100.0%
rem-exp-log100.0%
prod-exp100.0%
sub-neg100.0%
div-exp100.0%
rem-exp-log100.0%
Simplified100.0%
Taylor expanded in w around 0 56.1%
distribute-rgt-out56.1%
metadata-eval56.1%
metadata-eval56.1%
distribute-rgt-neg-in56.1%
associate-*r*56.1%
neg-mul-156.1%
distribute-rgt-neg-in56.1%
metadata-eval56.1%
Simplified56.1%
if -7.29999999999999982 < w Initial program 99.6%
add-sqr-sqrt59.4%
sqrt-unprod98.0%
sqr-neg98.0%
sqrt-unprod38.7%
add-sqr-sqrt97.7%
add-sqr-sqrt97.7%
sqrt-unprod97.7%
add-sqr-sqrt38.7%
sqrt-unprod74.7%
sqr-neg74.7%
sqrt-unprod36.0%
add-sqr-sqrt74.7%
pow174.7%
exp-neg74.7%
inv-pow74.7%
pow-prod-up97.8%
metadata-eval97.8%
metadata-eval97.8%
metadata-eval97.8%
Applied egg-rr97.8%
Taylor expanded in w around inf 97.8%
rem-exp-log91.4%
prod-exp91.4%
sub-neg91.4%
div-exp91.4%
rem-exp-log97.8%
Simplified97.8%
Taylor expanded in w around 0 93.5%
*-commutative93.5%
Simplified93.5%
Final simplification82.1%
(FPCore (w l) :precision binary64 (if (<= w -0.122) (* l (- w)) (/ l (+ 1.0 (* w (+ 1.0 (* w 0.5)))))))
double code(double w, double l) {
double tmp;
if (w <= -0.122) {
tmp = l * -w;
} else {
tmp = l / (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 <= (-0.122d0)) then
tmp = l * -w
else
tmp = l / (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 <= -0.122) {
tmp = l * -w;
} else {
tmp = l / (1.0 + (w * (1.0 + (w * 0.5))));
}
return tmp;
}
def code(w, l): tmp = 0 if w <= -0.122: tmp = l * -w else: tmp = l / (1.0 + (w * (1.0 + (w * 0.5)))) return tmp
function code(w, l) tmp = 0.0 if (w <= -0.122) tmp = Float64(l * Float64(-w)); else tmp = Float64(l / 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 <= -0.122) tmp = l * -w; else tmp = l / (1.0 + (w * (1.0 + (w * 0.5)))); end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, -0.122], N[(l * (-w)), $MachinePrecision], N[(l / 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 -0.122:\\
\;\;\;\;\ell \cdot \left(-w\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{1 + w \cdot \left(1 + w \cdot 0.5\right)}\\
\end{array}
\end{array}
if w < -0.122Initial program 99.8%
add-sqr-sqrt0.0%
sqrt-unprod53.9%
sqr-neg53.9%
sqrt-unprod53.9%
add-sqr-sqrt53.9%
add-sqr-sqrt53.9%
sqrt-unprod53.9%
add-sqr-sqrt53.9%
sqrt-unprod53.9%
sqr-neg53.9%
sqrt-unprod0.0%
add-sqr-sqrt0.2%
pow10.2%
exp-neg0.2%
inv-pow0.2%
pow-prod-up97.7%
metadata-eval97.7%
metadata-eval97.7%
metadata-eval97.7%
Applied egg-rr97.7%
Taylor expanded in w around inf 97.7%
rem-exp-log97.7%
prod-exp97.7%
sub-neg97.7%
div-exp97.7%
rem-exp-log97.7%
Simplified97.7%
Taylor expanded in w around 0 32.1%
*-rgt-identity32.1%
neg-mul-132.1%
distribute-rgt-neg-in32.1%
distribute-lft-out32.1%
sub-neg32.1%
Simplified32.1%
Taylor expanded in w around inf 32.1%
associate-*r*32.1%
mul-1-neg32.1%
Simplified32.1%
if -0.122 < w Initial program 99.6%
add-sqr-sqrt60.0%
sqrt-unprod99.1%
sqr-neg99.1%
sqrt-unprod39.0%
add-sqr-sqrt98.7%
add-sqr-sqrt98.7%
sqrt-unprod98.7%
add-sqr-sqrt39.0%
sqrt-unprod75.5%
sqr-neg75.5%
sqrt-unprod36.4%
add-sqr-sqrt75.5%
pow175.5%
exp-neg75.5%
inv-pow75.5%
pow-prod-up98.8%
metadata-eval98.8%
metadata-eval98.8%
metadata-eval98.8%
Applied egg-rr98.8%
Taylor expanded in w around inf 98.8%
rem-exp-log92.4%
prod-exp92.4%
sub-neg92.4%
div-exp92.4%
rem-exp-log98.8%
Simplified98.8%
Taylor expanded in w around 0 94.5%
*-commutative94.5%
Simplified94.5%
Final simplification75.0%
(FPCore (w l) :precision binary64 (if (<= w -0.04) (* l (- w)) (/ l (+ w 1.0))))
double code(double w, double l) {
double tmp;
if (w <= -0.04) {
tmp = l * -w;
} else {
tmp = l / (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 <= (-0.04d0)) then
tmp = l * -w
else
tmp = l / (w + 1.0d0)
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= -0.04) {
tmp = l * -w;
} else {
tmp = l / (w + 1.0);
}
return tmp;
}
def code(w, l): tmp = 0 if w <= -0.04: tmp = l * -w else: tmp = l / (w + 1.0) return tmp
function code(w, l) tmp = 0.0 if (w <= -0.04) tmp = Float64(l * Float64(-w)); else tmp = Float64(l / Float64(w + 1.0)); end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= -0.04) tmp = l * -w; else tmp = l / (w + 1.0); end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, -0.04], N[(l * (-w)), $MachinePrecision], N[(l / N[(w + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq -0.04:\\
\;\;\;\;\ell \cdot \left(-w\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{w + 1}\\
\end{array}
\end{array}
if w < -0.0400000000000000008Initial program 99.8%
add-sqr-sqrt0.0%
sqrt-unprod53.9%
sqr-neg53.9%
sqrt-unprod53.9%
add-sqr-sqrt53.9%
add-sqr-sqrt53.9%
sqrt-unprod53.9%
add-sqr-sqrt53.9%
sqrt-unprod53.9%
sqr-neg53.9%
sqrt-unprod0.0%
add-sqr-sqrt0.2%
pow10.2%
exp-neg0.2%
inv-pow0.2%
pow-prod-up97.7%
metadata-eval97.7%
metadata-eval97.7%
metadata-eval97.7%
Applied egg-rr97.7%
Taylor expanded in w around inf 97.7%
rem-exp-log97.7%
prod-exp97.7%
sub-neg97.7%
div-exp97.7%
rem-exp-log97.7%
Simplified97.7%
Taylor expanded in w around 0 32.1%
*-rgt-identity32.1%
neg-mul-132.1%
distribute-rgt-neg-in32.1%
distribute-lft-out32.1%
sub-neg32.1%
Simplified32.1%
Taylor expanded in w around inf 32.1%
associate-*r*32.1%
mul-1-neg32.1%
Simplified32.1%
if -0.0400000000000000008 < w Initial program 99.6%
add-sqr-sqrt60.0%
sqrt-unprod99.1%
sqr-neg99.1%
sqrt-unprod39.0%
add-sqr-sqrt98.7%
add-sqr-sqrt98.7%
sqrt-unprod98.7%
add-sqr-sqrt39.0%
sqrt-unprod75.5%
sqr-neg75.5%
sqrt-unprod36.4%
add-sqr-sqrt75.5%
pow175.5%
exp-neg75.5%
inv-pow75.5%
pow-prod-up98.8%
metadata-eval98.8%
metadata-eval98.8%
metadata-eval98.8%
Applied egg-rr98.8%
Taylor expanded in w around inf 98.8%
rem-exp-log92.4%
prod-exp92.4%
sub-neg92.4%
div-exp92.4%
rem-exp-log98.8%
Simplified98.8%
Taylor expanded in w around 0 89.7%
+-commutative89.7%
Simplified89.7%
Final simplification71.7%
(FPCore (w l) :precision binary64 (if (<= w -0.011) (* l (- w)) l))
double code(double w, double l) {
double tmp;
if (w <= -0.011) {
tmp = l * -w;
} else {
tmp = l;
}
return tmp;
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
real(8) :: tmp
if (w <= (-0.011d0)) then
tmp = l * -w
else
tmp = l
end if
code = tmp
end function
public static double code(double w, double l) {
double tmp;
if (w <= -0.011) {
tmp = l * -w;
} else {
tmp = l;
}
return tmp;
}
def code(w, l): tmp = 0 if w <= -0.011: tmp = l * -w else: tmp = l return tmp
function code(w, l) tmp = 0.0 if (w <= -0.011) tmp = Float64(l * Float64(-w)); else tmp = l; end return tmp end
function tmp_2 = code(w, l) tmp = 0.0; if (w <= -0.011) tmp = l * -w; else tmp = l; end tmp_2 = tmp; end
code[w_, l_] := If[LessEqual[w, -0.011], N[(l * (-w)), $MachinePrecision], l]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;w \leq -0.011:\\
\;\;\;\;\ell \cdot \left(-w\right)\\
\mathbf{else}:\\
\;\;\;\;\ell\\
\end{array}
\end{array}
if w < -0.010999999999999999Initial program 99.8%
add-sqr-sqrt0.0%
sqrt-unprod53.9%
sqr-neg53.9%
sqrt-unprod53.9%
add-sqr-sqrt53.9%
add-sqr-sqrt53.9%
sqrt-unprod53.9%
add-sqr-sqrt53.9%
sqrt-unprod53.9%
sqr-neg53.9%
sqrt-unprod0.0%
add-sqr-sqrt0.2%
pow10.2%
exp-neg0.2%
inv-pow0.2%
pow-prod-up97.7%
metadata-eval97.7%
metadata-eval97.7%
metadata-eval97.7%
Applied egg-rr97.7%
Taylor expanded in w around inf 97.7%
rem-exp-log97.7%
prod-exp97.7%
sub-neg97.7%
div-exp97.7%
rem-exp-log97.7%
Simplified97.7%
Taylor expanded in w around 0 32.1%
*-rgt-identity32.1%
neg-mul-132.1%
distribute-rgt-neg-in32.1%
distribute-lft-out32.1%
sub-neg32.1%
Simplified32.1%
Taylor expanded in w around inf 32.1%
associate-*r*32.1%
mul-1-neg32.1%
Simplified32.1%
if -0.010999999999999999 < w Initial program 99.6%
Taylor expanded in w around 0 76.8%
Final simplification62.9%
(FPCore (w l) :precision binary64 (* l (- 1.0 w)))
double code(double w, double l) {
return l * (1.0 - w);
}
real(8) function code(w, l)
real(8), intent (in) :: w
real(8), intent (in) :: l
code = l * (1.0d0 - w)
end function
public static double code(double w, double l) {
return l * (1.0 - w);
}
def code(w, l): return l * (1.0 - w)
function code(w, l) return Float64(l * Float64(1.0 - w)) end
function tmp = code(w, l) tmp = l * (1.0 - w); end
code[w_, l_] := N[(l * N[(1.0 - w), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\ell \cdot \left(1 - w\right)
\end{array}
Initial program 99.7%
add-sqr-sqrt41.3%
sqrt-unprod85.0%
sqr-neg85.0%
sqrt-unprod43.7%
add-sqr-sqrt84.7%
add-sqr-sqrt84.7%
sqrt-unprod84.7%
add-sqr-sqrt43.7%
sqrt-unprod68.7%
sqr-neg68.7%
sqrt-unprod25.0%
add-sqr-sqrt52.0%
pow152.0%
exp-neg52.0%
inv-pow52.0%
pow-prod-up98.4%
metadata-eval98.4%
metadata-eval98.4%
metadata-eval98.4%
Applied egg-rr98.4%
Taylor expanded in w around inf 98.4%
rem-exp-log94.0%
prod-exp94.0%
sub-neg94.0%
div-exp94.0%
rem-exp-log98.4%
Simplified98.4%
Taylor expanded in w around 0 62.5%
*-rgt-identity62.5%
neg-mul-162.5%
distribute-rgt-neg-in62.5%
distribute-lft-out62.5%
sub-neg62.5%
Simplified62.5%
(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.7%
Taylor expanded in w around 0 54.1%
herbie shell --seed 2024140
(FPCore (w l)
:name "exp-w (used to crash)"
:precision binary64
(* (exp (- w)) (pow l (exp w))))