
(FPCore (J l K U) :precision binary64 (+ (* (* J (- (exp l) (exp (- l)))) (cos (/ K 2.0))) U))
double code(double J, double l, double K, double U) {
return ((J * (exp(l) - exp(-l))) * cos((K / 2.0))) + U;
}
real(8) function code(j, l, k, u)
real(8), intent (in) :: j
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8), intent (in) :: u
code = ((j * (exp(l) - exp(-l))) * cos((k / 2.0d0))) + u
end function
public static double code(double J, double l, double K, double U) {
return ((J * (Math.exp(l) - Math.exp(-l))) * Math.cos((K / 2.0))) + U;
}
def code(J, l, K, U): return ((J * (math.exp(l) - math.exp(-l))) * math.cos((K / 2.0))) + U
function code(J, l, K, U) return Float64(Float64(Float64(J * Float64(exp(l) - exp(Float64(-l)))) * cos(Float64(K / 2.0))) + U) end
function tmp = code(J, l, K, U) tmp = ((J * (exp(l) - exp(-l))) * cos((K / 2.0))) + U; end
code[J_, l_, K_, U_] := N[(N[(N[(J * N[(N[Exp[l], $MachinePrecision] - N[Exp[(-l)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision]
\begin{array}{l}
\\
\left(J \cdot \left(e^{\ell} - e^{-\ell}\right)\right) \cdot \cos \left(\frac{K}{2}\right) + U
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 24 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (J l K U) :precision binary64 (+ (* (* J (- (exp l) (exp (- l)))) (cos (/ K 2.0))) U))
double code(double J, double l, double K, double U) {
return ((J * (exp(l) - exp(-l))) * cos((K / 2.0))) + U;
}
real(8) function code(j, l, k, u)
real(8), intent (in) :: j
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8), intent (in) :: u
code = ((j * (exp(l) - exp(-l))) * cos((k / 2.0d0))) + u
end function
public static double code(double J, double l, double K, double U) {
return ((J * (Math.exp(l) - Math.exp(-l))) * Math.cos((K / 2.0))) + U;
}
def code(J, l, K, U): return ((J * (math.exp(l) - math.exp(-l))) * math.cos((K / 2.0))) + U
function code(J, l, K, U) return Float64(Float64(Float64(J * Float64(exp(l) - exp(Float64(-l)))) * cos(Float64(K / 2.0))) + U) end
function tmp = code(J, l, K, U) tmp = ((J * (exp(l) - exp(-l))) * cos((K / 2.0))) + U; end
code[J_, l_, K_, U_] := N[(N[(N[(J * N[(N[Exp[l], $MachinePrecision] - N[Exp[(-l)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision]
\begin{array}{l}
\\
\left(J \cdot \left(e^{\ell} - e^{-\ell}\right)\right) \cdot \cos \left(\frac{K}{2}\right) + U
\end{array}
(FPCore (J l K U) :precision binary64 (fma (* (* (cos (/ K -2.0)) J) 2.0) (sinh l) U))
double code(double J, double l, double K, double U) {
return fma(((cos((K / -2.0)) * J) * 2.0), sinh(l), U);
}
function code(J, l, K, U) return fma(Float64(Float64(cos(Float64(K / -2.0)) * J) * 2.0), sinh(l), U) end
code[J_, l_, K_, U_] := N[(N[(N[(N[Cos[N[(K / -2.0), $MachinePrecision]], $MachinePrecision] * J), $MachinePrecision] * 2.0), $MachinePrecision] * N[Sinh[l], $MachinePrecision] + U), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\left(\cos \left(\frac{K}{-2}\right) \cdot J\right) \cdot 2, \sinh \ell, U\right)
\end{array}
Initial program 85.6%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift--.f64N/A
lift-exp.f64N/A
lift-exp.f64N/A
lift-neg.f64N/A
sinh-undefN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites99.9%
(FPCore (J l K U)
:precision binary64
(let* ((t_0 (+ (* (* J (- (exp l) (exp (- l)))) (cos (/ K 2.0))) U)))
(if (<= t_0 (- INFINITY))
(fma (* (* l l) (* l 2.0)) J U)
(if (<= t_0 1e+291) (fma (+ l l) J U) (fma (* l l) 4.0 U)))))
double code(double J, double l, double K, double U) {
double t_0 = ((J * (exp(l) - exp(-l))) * cos((K / 2.0))) + U;
double tmp;
if (t_0 <= -((double) INFINITY)) {
tmp = fma(((l * l) * (l * 2.0)), J, U);
} else if (t_0 <= 1e+291) {
tmp = fma((l + l), J, U);
} else {
tmp = fma((l * l), 4.0, U);
}
return tmp;
}
function code(J, l, K, U) t_0 = Float64(Float64(Float64(J * Float64(exp(l) - exp(Float64(-l)))) * cos(Float64(K / 2.0))) + U) tmp = 0.0 if (t_0 <= Float64(-Inf)) tmp = fma(Float64(Float64(l * l) * Float64(l * 2.0)), J, U); elseif (t_0 <= 1e+291) tmp = fma(Float64(l + l), J, U); else tmp = fma(Float64(l * l), 4.0, U); end return tmp end
code[J_, l_, K_, U_] := Block[{t$95$0 = N[(N[(N[(J * N[(N[Exp[l], $MachinePrecision] - N[Exp[(-l)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision]}, If[LessEqual[t$95$0, (-Infinity)], N[(N[(N[(l * l), $MachinePrecision] * N[(l * 2.0), $MachinePrecision]), $MachinePrecision] * J + U), $MachinePrecision], If[LessEqual[t$95$0, 1e+291], N[(N[(l + l), $MachinePrecision] * J + U), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * 4.0 + U), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(J \cdot \left(e^{\ell} - e^{-\ell}\right)\right) \cdot \cos \left(\frac{K}{2}\right) + U\\
\mathbf{if}\;t\_0 \leq -\infty:\\
\;\;\;\;\mathsf{fma}\left(\left(\ell \cdot \ell\right) \cdot \left(\ell \cdot 2\right), J, U\right)\\
\mathbf{elif}\;t\_0 \leq 10^{+291}:\\
\;\;\;\;\mathsf{fma}\left(\ell + \ell, J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\ell \cdot \ell, 4, U\right)\\
\end{array}
\end{array}
if (+.f64 (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) U) < -inf.0Initial program 100.0%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval36.6
Applied rewrites36.6%
Taylor expanded in K around 0
Applied rewrites28.9%
Applied rewrites28.9%
Applied rewrites53.5%
if -inf.0 < (+.f64 (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) U) < 9.9999999999999996e290Initial program 72.0%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval98.3
Applied rewrites98.3%
Taylor expanded in K around 0
Applied rewrites83.9%
Applied rewrites83.9%
if 9.9999999999999996e290 < (+.f64 (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) U) Initial program 100.0%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval30.5
Applied rewrites30.5%
Taylor expanded in K around 0
Applied rewrites20.9%
Applied rewrites20.9%
Applied rewrites64.9%
(FPCore (J l K U)
:precision binary64
(if (<= (+ (* (* J (- (exp l) (exp (- l)))) (cos (/ K 2.0))) U) 1e+306)
(fma (* 1.0 J) (* (fma 0.3333333333333333 (* l l) 2.0) l) U)
(fma
(* l 2.0)
(fma (* (fma (* K K) 0.0026041666666666665 -0.125) l) (* K K) l)
U)))
double code(double J, double l, double K, double U) {
double tmp;
if ((((J * (exp(l) - exp(-l))) * cos((K / 2.0))) + U) <= 1e+306) {
tmp = fma((1.0 * J), (fma(0.3333333333333333, (l * l), 2.0) * l), U);
} else {
tmp = fma((l * 2.0), fma((fma((K * K), 0.0026041666666666665, -0.125) * l), (K * K), l), U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (Float64(Float64(Float64(J * Float64(exp(l) - exp(Float64(-l)))) * cos(Float64(K / 2.0))) + U) <= 1e+306) tmp = fma(Float64(1.0 * J), Float64(fma(0.3333333333333333, Float64(l * l), 2.0) * l), U); else tmp = fma(Float64(l * 2.0), fma(Float64(fma(Float64(K * K), 0.0026041666666666665, -0.125) * l), Float64(K * K), l), U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[(N[(N[(J * N[(N[Exp[l], $MachinePrecision] - N[Exp[(-l)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], 1e+306], N[(N[(1.0 * J), $MachinePrecision] * N[(N[(0.3333333333333333 * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision], N[(N[(l * 2.0), $MachinePrecision] * N[(N[(N[(N[(K * K), $MachinePrecision] * 0.0026041666666666665 + -0.125), $MachinePrecision] * l), $MachinePrecision] * N[(K * K), $MachinePrecision] + l), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left(J \cdot \left(e^{\ell} - e^{-\ell}\right)\right) \cdot \cos \left(\frac{K}{2}\right) + U \leq 10^{+306}:\\
\;\;\;\;\mathsf{fma}\left(1 \cdot J, \mathsf{fma}\left(0.3333333333333333, \ell \cdot \ell, 2\right) \cdot \ell, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\ell \cdot 2, \mathsf{fma}\left(\mathsf{fma}\left(K \cdot K, 0.0026041666666666665, -0.125\right) \cdot \ell, K \cdot K, \ell\right), U\right)\\
\end{array}
\end{array}
if (+.f64 (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) U) < 1.00000000000000002e306Initial program 80.8%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6490.6
Applied rewrites90.6%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites90.6%
Taylor expanded in K around 0
Applied rewrites74.6%
if 1.00000000000000002e306 < (+.f64 (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) U) Initial program 100.0%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6429.4
Applied rewrites29.4%
Taylor expanded in K around 0
Applied rewrites36.0%
Applied rewrites82.2%
(FPCore (J l K U) :precision binary64 (if (<= (+ (* (* J (- (exp l) (exp (- l)))) (cos (/ K 2.0))) U) 1e+291) (fma (+ l l) J U) (fma (* l l) 4.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if ((((J * (exp(l) - exp(-l))) * cos((K / 2.0))) + U) <= 1e+291) {
tmp = fma((l + l), J, U);
} else {
tmp = fma((l * l), 4.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (Float64(Float64(Float64(J * Float64(exp(l) - exp(Float64(-l)))) * cos(Float64(K / 2.0))) + U) <= 1e+291) tmp = fma(Float64(l + l), J, U); else tmp = fma(Float64(l * l), 4.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[(N[(N[(J * N[(N[Exp[l], $MachinePrecision] - N[Exp[(-l)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], 1e+291], N[(N[(l + l), $MachinePrecision] * J + U), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * 4.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left(J \cdot \left(e^{\ell} - e^{-\ell}\right)\right) \cdot \cos \left(\frac{K}{2}\right) + U \leq 10^{+291}:\\
\;\;\;\;\mathsf{fma}\left(\ell + \ell, J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\ell \cdot \ell, 4, U\right)\\
\end{array}
\end{array}
if (+.f64 (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) U) < 9.9999999999999996e290Initial program 80.7%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval79.2
Applied rewrites79.2%
Taylor expanded in K around 0
Applied rewrites66.9%
Applied rewrites66.9%
if 9.9999999999999996e290 < (+.f64 (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) U) Initial program 100.0%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval30.5
Applied rewrites30.5%
Taylor expanded in K around 0
Applied rewrites20.9%
Applied rewrites20.9%
Applied rewrites64.9%
(FPCore (J l K U) :precision binary64 (if (<= (+ (* (* J (- (exp l) (exp (- l)))) (cos (/ K 2.0))) U) 1e+306) (fma (+ l l) J U) (fma (* J J) 4.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if ((((J * (exp(l) - exp(-l))) * cos((K / 2.0))) + U) <= 1e+306) {
tmp = fma((l + l), J, U);
} else {
tmp = fma((J * J), 4.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (Float64(Float64(Float64(J * Float64(exp(l) - exp(Float64(-l)))) * cos(Float64(K / 2.0))) + U) <= 1e+306) tmp = fma(Float64(l + l), J, U); else tmp = fma(Float64(J * J), 4.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[(N[(N[(J * N[(N[Exp[l], $MachinePrecision] - N[Exp[(-l)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], 1e+306], N[(N[(l + l), $MachinePrecision] * J + U), $MachinePrecision], N[(N[(J * J), $MachinePrecision] * 4.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left(J \cdot \left(e^{\ell} - e^{-\ell}\right)\right) \cdot \cos \left(\frac{K}{2}\right) + U \leq 10^{+306}:\\
\;\;\;\;\mathsf{fma}\left(\ell + \ell, J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(J \cdot J, 4, U\right)\\
\end{array}
\end{array}
if (+.f64 (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) U) < 1.00000000000000002e306Initial program 80.8%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval79.3
Applied rewrites79.3%
Taylor expanded in K around 0
Applied rewrites67.1%
Applied rewrites67.1%
if 1.00000000000000002e306 < (+.f64 (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) U) Initial program 100.0%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval29.4
Applied rewrites29.4%
Taylor expanded in K around 0
Applied rewrites19.6%
Applied rewrites19.6%
Applied rewrites27.0%
(FPCore (J l K U) :precision binary64 (if (<= (+ (* (* J (- (exp l) (exp (- l)))) (cos (/ K 2.0))) U) 1e+306) (fma (+ l l) J U) (fma (+ J J) J U)))
double code(double J, double l, double K, double U) {
double tmp;
if ((((J * (exp(l) - exp(-l))) * cos((K / 2.0))) + U) <= 1e+306) {
tmp = fma((l + l), J, U);
} else {
tmp = fma((J + J), J, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (Float64(Float64(Float64(J * Float64(exp(l) - exp(Float64(-l)))) * cos(Float64(K / 2.0))) + U) <= 1e+306) tmp = fma(Float64(l + l), J, U); else tmp = fma(Float64(J + J), J, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[(N[(N[(J * N[(N[Exp[l], $MachinePrecision] - N[Exp[(-l)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], 1e+306], N[(N[(l + l), $MachinePrecision] * J + U), $MachinePrecision], N[(N[(J + J), $MachinePrecision] * J + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left(J \cdot \left(e^{\ell} - e^{-\ell}\right)\right) \cdot \cos \left(\frac{K}{2}\right) + U \leq 10^{+306}:\\
\;\;\;\;\mathsf{fma}\left(\ell + \ell, J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(J + J, J, U\right)\\
\end{array}
\end{array}
if (+.f64 (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) U) < 1.00000000000000002e306Initial program 80.8%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval79.3
Applied rewrites79.3%
Taylor expanded in K around 0
Applied rewrites67.1%
Applied rewrites67.1%
if 1.00000000000000002e306 < (+.f64 (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) U) Initial program 100.0%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval29.4
Applied rewrites29.4%
Taylor expanded in K around 0
Applied rewrites19.6%
Applied rewrites19.6%
Applied rewrites27.0%
(FPCore (J l K U) :precision binary64 (if (<= (* (* J (- (exp l) (exp (- l)))) (cos (/ K 2.0))) 0.0) (fma l 2.0 U) (fma (+ J J) J U)))
double code(double J, double l, double K, double U) {
double tmp;
if (((J * (exp(l) - exp(-l))) * cos((K / 2.0))) <= 0.0) {
tmp = fma(l, 2.0, U);
} else {
tmp = fma((J + J), J, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (Float64(Float64(J * Float64(exp(l) - exp(Float64(-l)))) * cos(Float64(K / 2.0))) <= 0.0) tmp = fma(l, 2.0, U); else tmp = fma(Float64(J + J), J, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[(N[(J * N[(N[Exp[l], $MachinePrecision] - N[Exp[(-l)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 0.0], N[(l * 2.0 + U), $MachinePrecision], N[(N[(J + J), $MachinePrecision] * J + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\left(J \cdot \left(e^{\ell} - e^{-\ell}\right)\right) \cdot \cos \left(\frac{K}{2}\right) \leq 0:\\
\;\;\;\;\mathsf{fma}\left(\ell, 2, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(J + J, J, U\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) < 0.0Initial program 80.8%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval80.0
Applied rewrites80.0%
Taylor expanded in K around 0
Applied rewrites67.8%
Applied rewrites67.8%
Applied rewrites45.3%
if 0.0 < (*.f64 (*.f64 J (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l)))) (cos.f64 (/.f64 K #s(literal 2 binary64)))) Initial program 98.9%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval30.4
Applied rewrites30.4%
Taylor expanded in K around 0
Applied rewrites20.5%
Applied rewrites20.5%
Applied rewrites25.9%
(FPCore (J l K U)
:precision binary64
(let* ((t_0 (cos (/ K 2.0))))
(if (<= t_0 -0.91)
(fma (* (* (fma (* K K) -0.125 1.0) l) 2.0) J U)
(if (<= t_0 -0.02)
(fma (* 1.0 J) (* (fma (fabs (* 0.3333333333333333 l)) l 2.0) l) U)
(fma
(* 1.0 J)
(*
(fma
(fma 0.016666666666666666 (* l l) 0.3333333333333333)
(* l l)
2.0)
l)
U)))))
double code(double J, double l, double K, double U) {
double t_0 = cos((K / 2.0));
double tmp;
if (t_0 <= -0.91) {
tmp = fma(((fma((K * K), -0.125, 1.0) * l) * 2.0), J, U);
} else if (t_0 <= -0.02) {
tmp = fma((1.0 * J), (fma(fabs((0.3333333333333333 * l)), l, 2.0) * l), U);
} else {
tmp = fma((1.0 * J), (fma(fma(0.016666666666666666, (l * l), 0.3333333333333333), (l * l), 2.0) * l), U);
}
return tmp;
}
function code(J, l, K, U) t_0 = cos(Float64(K / 2.0)) tmp = 0.0 if (t_0 <= -0.91) tmp = fma(Float64(Float64(fma(Float64(K * K), -0.125, 1.0) * l) * 2.0), J, U); elseif (t_0 <= -0.02) tmp = fma(Float64(1.0 * J), Float64(fma(abs(Float64(0.3333333333333333 * l)), l, 2.0) * l), U); else tmp = fma(Float64(1.0 * J), Float64(fma(fma(0.016666666666666666, Float64(l * l), 0.3333333333333333), Float64(l * l), 2.0) * l), U); end return tmp end
code[J_, l_, K_, U_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t$95$0, -0.91], N[(N[(N[(N[(N[(K * K), $MachinePrecision] * -0.125 + 1.0), $MachinePrecision] * l), $MachinePrecision] * 2.0), $MachinePrecision] * J + U), $MachinePrecision], If[LessEqual[t$95$0, -0.02], N[(N[(1.0 * J), $MachinePrecision] * N[(N[(N[Abs[N[(0.3333333333333333 * l), $MachinePrecision]], $MachinePrecision] * l + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision], N[(N[(1.0 * J), $MachinePrecision] * N[(N[(N[(0.016666666666666666 * N[(l * l), $MachinePrecision] + 0.3333333333333333), $MachinePrecision] * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
\mathbf{if}\;t\_0 \leq -0.91:\\
\;\;\;\;\mathsf{fma}\left(\left(\mathsf{fma}\left(K \cdot K, -0.125, 1\right) \cdot \ell\right) \cdot 2, J, U\right)\\
\mathbf{elif}\;t\_0 \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(1 \cdot J, \mathsf{fma}\left(\left|0.3333333333333333 \cdot \ell\right|, \ell, 2\right) \cdot \ell, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 \cdot J, \mathsf{fma}\left(\mathsf{fma}\left(0.016666666666666666, \ell \cdot \ell, 0.3333333333333333\right), \ell \cdot \ell, 2\right) \cdot \ell, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.910000000000000031Initial program 85.4%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6456.5
Applied rewrites56.5%
Taylor expanded in K around 0
Applied rewrites56.8%
if -0.910000000000000031 < (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.0200000000000000004Initial program 89.7%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6490.0
Applied rewrites90.0%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites90.0%
Taylor expanded in K around 0
Applied rewrites49.4%
Applied rewrites68.7%
if -0.0200000000000000004 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.6%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6487.7
Applied rewrites87.7%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites87.7%
Taylor expanded in K around 0
Applied rewrites82.1%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6485.7
Applied rewrites85.7%
(FPCore (J l K U)
:precision binary64
(let* ((t_0 (cos (/ K 2.0))))
(if (<= t_0 -0.91)
(fma (* (* (fma (* K K) -0.125 1.0) l) 2.0) J U)
(if (<= t_0 -0.02)
(fma (* 1.0 J) (* (fma (fabs (* 0.3333333333333333 l)) l 2.0) l) U)
(fma (* 1.0 J) (* (fma 0.3333333333333333 (* l l) 2.0) l) U)))))
double code(double J, double l, double K, double U) {
double t_0 = cos((K / 2.0));
double tmp;
if (t_0 <= -0.91) {
tmp = fma(((fma((K * K), -0.125, 1.0) * l) * 2.0), J, U);
} else if (t_0 <= -0.02) {
tmp = fma((1.0 * J), (fma(fabs((0.3333333333333333 * l)), l, 2.0) * l), U);
} else {
tmp = fma((1.0 * J), (fma(0.3333333333333333, (l * l), 2.0) * l), U);
}
return tmp;
}
function code(J, l, K, U) t_0 = cos(Float64(K / 2.0)) tmp = 0.0 if (t_0 <= -0.91) tmp = fma(Float64(Float64(fma(Float64(K * K), -0.125, 1.0) * l) * 2.0), J, U); elseif (t_0 <= -0.02) tmp = fma(Float64(1.0 * J), Float64(fma(abs(Float64(0.3333333333333333 * l)), l, 2.0) * l), U); else tmp = fma(Float64(1.0 * J), Float64(fma(0.3333333333333333, Float64(l * l), 2.0) * l), U); end return tmp end
code[J_, l_, K_, U_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t$95$0, -0.91], N[(N[(N[(N[(N[(K * K), $MachinePrecision] * -0.125 + 1.0), $MachinePrecision] * l), $MachinePrecision] * 2.0), $MachinePrecision] * J + U), $MachinePrecision], If[LessEqual[t$95$0, -0.02], N[(N[(1.0 * J), $MachinePrecision] * N[(N[(N[Abs[N[(0.3333333333333333 * l), $MachinePrecision]], $MachinePrecision] * l + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision], N[(N[(1.0 * J), $MachinePrecision] * N[(N[(0.3333333333333333 * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
\mathbf{if}\;t\_0 \leq -0.91:\\
\;\;\;\;\mathsf{fma}\left(\left(\mathsf{fma}\left(K \cdot K, -0.125, 1\right) \cdot \ell\right) \cdot 2, J, U\right)\\
\mathbf{elif}\;t\_0 \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(1 \cdot J, \mathsf{fma}\left(\left|0.3333333333333333 \cdot \ell\right|, \ell, 2\right) \cdot \ell, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 \cdot J, \mathsf{fma}\left(0.3333333333333333, \ell \cdot \ell, 2\right) \cdot \ell, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.910000000000000031Initial program 85.4%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6456.5
Applied rewrites56.5%
Taylor expanded in K around 0
Applied rewrites56.8%
if -0.910000000000000031 < (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.0200000000000000004Initial program 89.7%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6490.0
Applied rewrites90.0%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites90.0%
Taylor expanded in K around 0
Applied rewrites49.4%
Applied rewrites68.7%
if -0.0200000000000000004 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.6%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6487.7
Applied rewrites87.7%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites87.7%
Taylor expanded in K around 0
Applied rewrites82.1%
(FPCore (J l K U) :precision binary64 (if (<= (cos (/ K 2.0)) -0.02) (fma (* (cos (* -0.5 K)) (* J (fma (* l l) 0.3333333333333333 2.0))) l U) (fma (* (* 1.0 (sinh l)) J) 2.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.02) {
tmp = fma((cos((-0.5 * K)) * (J * fma((l * l), 0.3333333333333333, 2.0))), l, U);
} else {
tmp = fma(((1.0 * sinh(l)) * J), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.02) tmp = fma(Float64(cos(Float64(-0.5 * K)) * Float64(J * fma(Float64(l * l), 0.3333333333333333, 2.0))), l, U); else tmp = fma(Float64(Float64(1.0 * sinh(l)) * J), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.02], N[(N[(N[Cos[N[(-0.5 * K), $MachinePrecision]], $MachinePrecision] * N[(J * N[(N[(l * l), $MachinePrecision] * 0.3333333333333333 + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * l + U), $MachinePrecision], N[(N[(N[(1.0 * N[Sinh[l], $MachinePrecision]), $MachinePrecision] * J), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(\cos \left(-0.5 \cdot K\right) \cdot \left(J \cdot \mathsf{fma}\left(\ell \cdot \ell, 0.3333333333333333, 2\right)\right), \ell, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(1 \cdot \sinh \ell\right) \cdot J, 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.0200000000000000004Initial program 88.3%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites84.9%
if -0.0200000000000000004 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.6%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in K around 0
Applied rewrites93.9%
lift-fma.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6493.9
lift-*.f64N/A
*-commutativeN/A
lower-*.f6493.9
Applied rewrites93.9%
(FPCore (J l K U) :precision binary64 (if (<= (cos (/ K 2.0)) -0.02) (fma (* (+ J J) l) (cos (* -0.5 K)) U) (fma (* (* 1.0 (sinh l)) J) 2.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.02) {
tmp = fma(((J + J) * l), cos((-0.5 * K)), U);
} else {
tmp = fma(((1.0 * sinh(l)) * J), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.02) tmp = fma(Float64(Float64(J + J) * l), cos(Float64(-0.5 * K)), U); else tmp = fma(Float64(Float64(1.0 * sinh(l)) * J), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.02], N[(N[(N[(J + J), $MachinePrecision] * l), $MachinePrecision] * N[Cos[N[(-0.5 * K), $MachinePrecision]], $MachinePrecision] + U), $MachinePrecision], N[(N[(N[(1.0 * N[Sinh[l], $MachinePrecision]), $MachinePrecision] * J), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(\left(J + J\right) \cdot \ell, \cos \left(-0.5 \cdot K\right), U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(1 \cdot \sinh \ell\right) \cdot J, 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.0200000000000000004Initial program 88.3%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval68.3
Applied rewrites68.3%
Applied rewrites68.3%
if -0.0200000000000000004 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.6%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in K around 0
Applied rewrites93.9%
lift-fma.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6493.9
lift-*.f64N/A
*-commutativeN/A
lower-*.f6493.9
Applied rewrites93.9%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) -0.02)
(fma
(* (fma (* K K) -0.125 1.0) J)
(* (fma 0.3333333333333333 (* l l) 2.0) l)
U)
(fma (* (* 1.0 (sinh l)) J) 2.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.02) {
tmp = fma((fma((K * K), -0.125, 1.0) * J), (fma(0.3333333333333333, (l * l), 2.0) * l), U);
} else {
tmp = fma(((1.0 * sinh(l)) * J), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.02) tmp = fma(Float64(fma(Float64(K * K), -0.125, 1.0) * J), Float64(fma(0.3333333333333333, Float64(l * l), 2.0) * l), U); else tmp = fma(Float64(Float64(1.0 * sinh(l)) * J), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.02], N[(N[(N[(N[(K * K), $MachinePrecision] * -0.125 + 1.0), $MachinePrecision] * J), $MachinePrecision] * N[(N[(0.3333333333333333 * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision], N[(N[(N[(1.0 * N[Sinh[l], $MachinePrecision]), $MachinePrecision] * J), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(K \cdot K, -0.125, 1\right) \cdot J, \mathsf{fma}\left(0.3333333333333333, \ell \cdot \ell, 2\right) \cdot \ell, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(1 \cdot \sinh \ell\right) \cdot J, 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.0200000000000000004Initial program 88.3%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6486.3
Applied rewrites86.3%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites86.3%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6462.2
Applied rewrites62.2%
if -0.0200000000000000004 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.6%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in K around 0
Applied rewrites93.9%
lift-fma.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6493.9
lift-*.f64N/A
*-commutativeN/A
lower-*.f6493.9
Applied rewrites93.9%
(FPCore (J l K U) :precision binary64 (fma (* (* (sinh l) (cos (* 0.5 K))) 2.0) J U))
double code(double J, double l, double K, double U) {
return fma(((sinh(l) * cos((0.5 * K))) * 2.0), J, U);
}
function code(J, l, K, U) return fma(Float64(Float64(sinh(l) * cos(Float64(0.5 * K))) * 2.0), J, U) end
code[J_, l_, K_, U_] := N[(N[(N[(N[Sinh[l], $MachinePrecision] * N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * 2.0), $MachinePrecision] * J + U), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\left(\sinh \ell \cdot \cos \left(0.5 \cdot K\right)\right) \cdot 2, J, U\right)
\end{array}
Initial program 85.6%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.9
Applied rewrites99.9%
(FPCore (J l K U)
:precision binary64
(if (<= K 5e-30)
(fma (* (* 1.0 (sinh l)) J) 2.0 U)
(+
(*
(*
J
(*
(fma
(fma
(fma 0.0003968253968253968 (* l l) 0.016666666666666666)
(* l l)
0.3333333333333333)
(* l l)
2.0)
l))
(cos (/ K 2.0)))
U)))
double code(double J, double l, double K, double U) {
double tmp;
if (K <= 5e-30) {
tmp = fma(((1.0 * sinh(l)) * J), 2.0, U);
} else {
tmp = ((J * (fma(fma(fma(0.0003968253968253968, (l * l), 0.016666666666666666), (l * l), 0.3333333333333333), (l * l), 2.0) * l)) * cos((K / 2.0))) + U;
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (K <= 5e-30) tmp = fma(Float64(Float64(1.0 * sinh(l)) * J), 2.0, U); else tmp = Float64(Float64(Float64(J * Float64(fma(fma(fma(0.0003968253968253968, Float64(l * l), 0.016666666666666666), Float64(l * l), 0.3333333333333333), Float64(l * l), 2.0) * l)) * cos(Float64(K / 2.0))) + U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[K, 5e-30], N[(N[(N[(1.0 * N[Sinh[l], $MachinePrecision]), $MachinePrecision] * J), $MachinePrecision] * 2.0 + U), $MachinePrecision], N[(N[(N[(J * N[(N[(N[(N[(0.0003968253968253968 * N[(l * l), $MachinePrecision] + 0.016666666666666666), $MachinePrecision] * N[(l * l), $MachinePrecision] + 0.3333333333333333), $MachinePrecision] * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;K \leq 5 \cdot 10^{-30}:\\
\;\;\;\;\mathsf{fma}\left(\left(1 \cdot \sinh \ell\right) \cdot J, 2, U\right)\\
\mathbf{else}:\\
\;\;\;\;\left(J \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.0003968253968253968, \ell \cdot \ell, 0.016666666666666666\right), \ell \cdot \ell, 0.3333333333333333\right), \ell \cdot \ell, 2\right) \cdot \ell\right)\right) \cdot \cos \left(\frac{K}{2}\right) + U\\
\end{array}
\end{array}
if K < 4.99999999999999972e-30Initial program 85.1%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in K around 0
Applied rewrites84.5%
lift-fma.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6484.5
lift-*.f64N/A
*-commutativeN/A
lower-*.f6484.5
Applied rewrites84.5%
if 4.99999999999999972e-30 < K Initial program 87.0%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6490.8
Applied rewrites90.8%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) -0.02)
(fma
(* (fma (* K K) -0.125 1.0) J)
(* (fma 0.3333333333333333 (* l l) 2.0) l)
U)
(fma
(* 1.0 J)
(*
(fma
(fma
(fma 0.0003968253968253968 (* l l) 0.016666666666666666)
(* l l)
0.3333333333333333)
(* l l)
2.0)
l)
U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.02) {
tmp = fma((fma((K * K), -0.125, 1.0) * J), (fma(0.3333333333333333, (l * l), 2.0) * l), U);
} else {
tmp = fma((1.0 * J), (fma(fma(fma(0.0003968253968253968, (l * l), 0.016666666666666666), (l * l), 0.3333333333333333), (l * l), 2.0) * l), U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.02) tmp = fma(Float64(fma(Float64(K * K), -0.125, 1.0) * J), Float64(fma(0.3333333333333333, Float64(l * l), 2.0) * l), U); else tmp = fma(Float64(1.0 * J), Float64(fma(fma(fma(0.0003968253968253968, Float64(l * l), 0.016666666666666666), Float64(l * l), 0.3333333333333333), Float64(l * l), 2.0) * l), U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.02], N[(N[(N[(N[(K * K), $MachinePrecision] * -0.125 + 1.0), $MachinePrecision] * J), $MachinePrecision] * N[(N[(0.3333333333333333 * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision], N[(N[(1.0 * J), $MachinePrecision] * N[(N[(N[(N[(0.0003968253968253968 * N[(l * l), $MachinePrecision] + 0.016666666666666666), $MachinePrecision] * N[(l * l), $MachinePrecision] + 0.3333333333333333), $MachinePrecision] * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(K \cdot K, -0.125, 1\right) \cdot J, \mathsf{fma}\left(0.3333333333333333, \ell \cdot \ell, 2\right) \cdot \ell, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 \cdot J, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.0003968253968253968, \ell \cdot \ell, 0.016666666666666666\right), \ell \cdot \ell, 0.3333333333333333\right), \ell \cdot \ell, 2\right) \cdot \ell, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.0200000000000000004Initial program 88.3%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6486.3
Applied rewrites86.3%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites86.3%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6462.2
Applied rewrites62.2%
if -0.0200000000000000004 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.6%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6487.7
Applied rewrites87.7%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites87.7%
Taylor expanded in K around 0
Applied rewrites82.1%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6487.7
Applied rewrites87.7%
(FPCore (J l K U)
:precision binary64
(if (<= K 5.2e-6)
(fma (* (* 1.0 (sinh l)) J) 2.0 U)
(+
(*
(*
J
(*
(fma (fma 0.016666666666666666 (* l l) 0.3333333333333333) (* l l) 2.0)
l))
(cos (/ K 2.0)))
U)))
double code(double J, double l, double K, double U) {
double tmp;
if (K <= 5.2e-6) {
tmp = fma(((1.0 * sinh(l)) * J), 2.0, U);
} else {
tmp = ((J * (fma(fma(0.016666666666666666, (l * l), 0.3333333333333333), (l * l), 2.0) * l)) * cos((K / 2.0))) + U;
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (K <= 5.2e-6) tmp = fma(Float64(Float64(1.0 * sinh(l)) * J), 2.0, U); else tmp = Float64(Float64(Float64(J * Float64(fma(fma(0.016666666666666666, Float64(l * l), 0.3333333333333333), Float64(l * l), 2.0) * l)) * cos(Float64(K / 2.0))) + U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[K, 5.2e-6], N[(N[(N[(1.0 * N[Sinh[l], $MachinePrecision]), $MachinePrecision] * J), $MachinePrecision] * 2.0 + U), $MachinePrecision], N[(N[(N[(J * N[(N[(N[(0.016666666666666666 * N[(l * l), $MachinePrecision] + 0.3333333333333333), $MachinePrecision] * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;K \leq 5.2 \cdot 10^{-6}:\\
\;\;\;\;\mathsf{fma}\left(\left(1 \cdot \sinh \ell\right) \cdot J, 2, U\right)\\
\mathbf{else}:\\
\;\;\;\;\left(J \cdot \left(\mathsf{fma}\left(\mathsf{fma}\left(0.016666666666666666, \ell \cdot \ell, 0.3333333333333333\right), \ell \cdot \ell, 2\right) \cdot \ell\right)\right) \cdot \cos \left(\frac{K}{2}\right) + U\\
\end{array}
\end{array}
if K < 5.20000000000000019e-6Initial program 85.4%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in K around 0
Applied rewrites84.8%
lift-fma.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f6484.8
lift-*.f64N/A
*-commutativeN/A
lower-*.f6484.8
Applied rewrites84.8%
if 5.20000000000000019e-6 < K Initial program 86.1%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6490.2
Applied rewrites90.2%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) -0.02)
(fma
(* (fma (* K K) -0.125 1.0) J)
(* (fma 0.3333333333333333 (* l l) 2.0) l)
U)
(fma
(* 1.0 J)
(*
(fma (fma 0.016666666666666666 (* l l) 0.3333333333333333) (* l l) 2.0)
l)
U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.02) {
tmp = fma((fma((K * K), -0.125, 1.0) * J), (fma(0.3333333333333333, (l * l), 2.0) * l), U);
} else {
tmp = fma((1.0 * J), (fma(fma(0.016666666666666666, (l * l), 0.3333333333333333), (l * l), 2.0) * l), U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.02) tmp = fma(Float64(fma(Float64(K * K), -0.125, 1.0) * J), Float64(fma(0.3333333333333333, Float64(l * l), 2.0) * l), U); else tmp = fma(Float64(1.0 * J), Float64(fma(fma(0.016666666666666666, Float64(l * l), 0.3333333333333333), Float64(l * l), 2.0) * l), U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.02], N[(N[(N[(N[(K * K), $MachinePrecision] * -0.125 + 1.0), $MachinePrecision] * J), $MachinePrecision] * N[(N[(0.3333333333333333 * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision], N[(N[(1.0 * J), $MachinePrecision] * N[(N[(N[(0.016666666666666666 * N[(l * l), $MachinePrecision] + 0.3333333333333333), $MachinePrecision] * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(K \cdot K, -0.125, 1\right) \cdot J, \mathsf{fma}\left(0.3333333333333333, \ell \cdot \ell, 2\right) \cdot \ell, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 \cdot J, \mathsf{fma}\left(\mathsf{fma}\left(0.016666666666666666, \ell \cdot \ell, 0.3333333333333333\right), \ell \cdot \ell, 2\right) \cdot \ell, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.0200000000000000004Initial program 88.3%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6486.3
Applied rewrites86.3%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites86.3%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6462.2
Applied rewrites62.2%
if -0.0200000000000000004 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.6%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6487.7
Applied rewrites87.7%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites87.7%
Taylor expanded in K around 0
Applied rewrites82.1%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6485.7
Applied rewrites85.7%
(FPCore (J l K U) :precision binary64 (if (<= K 20.0) (fma (* (* (sinh l) (fma (* K K) -0.125 1.0)) 2.0) J U) (+ (* (* J (* (fma (* l l) 0.3333333333333333 2.0) l)) (cos (/ K 2.0))) U)))
double code(double J, double l, double K, double U) {
double tmp;
if (K <= 20.0) {
tmp = fma(((sinh(l) * fma((K * K), -0.125, 1.0)) * 2.0), J, U);
} else {
tmp = ((J * (fma((l * l), 0.3333333333333333, 2.0) * l)) * cos((K / 2.0))) + U;
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (K <= 20.0) tmp = fma(Float64(Float64(sinh(l) * fma(Float64(K * K), -0.125, 1.0)) * 2.0), J, U); else tmp = Float64(Float64(Float64(J * Float64(fma(Float64(l * l), 0.3333333333333333, 2.0) * l)) * cos(Float64(K / 2.0))) + U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[K, 20.0], N[(N[(N[(N[Sinh[l], $MachinePrecision] * N[(N[(K * K), $MachinePrecision] * -0.125 + 1.0), $MachinePrecision]), $MachinePrecision] * 2.0), $MachinePrecision] * J + U), $MachinePrecision], N[(N[(N[(J * N[(N[(N[(l * l), $MachinePrecision] * 0.3333333333333333 + 2.0), $MachinePrecision] * l), $MachinePrecision]), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;K \leq 20:\\
\;\;\;\;\mathsf{fma}\left(\left(\sinh \ell \cdot \mathsf{fma}\left(K \cdot K, -0.125, 1\right)\right) \cdot 2, J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\left(J \cdot \left(\mathsf{fma}\left(\ell \cdot \ell, 0.3333333333333333, 2\right) \cdot \ell\right)\right) \cdot \cos \left(\frac{K}{2}\right) + U\\
\end{array}
\end{array}
if K < 20Initial program 85.2%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6474.4
Applied rewrites74.4%
if 20 < K Initial program 87.0%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6482.8
Applied rewrites82.8%
(FPCore (J l K U) :precision binary64 (if (<= (cos (/ K 2.0)) -0.02) (fma (* (* (fma (* K K) -0.125 1.0) l) 2.0) J U) (fma (* 1.0 J) (* (fma 0.3333333333333333 (* l l) 2.0) l) U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.02) {
tmp = fma(((fma((K * K), -0.125, 1.0) * l) * 2.0), J, U);
} else {
tmp = fma((1.0 * J), (fma(0.3333333333333333, (l * l), 2.0) * l), U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.02) tmp = fma(Float64(Float64(fma(Float64(K * K), -0.125, 1.0) * l) * 2.0), J, U); else tmp = fma(Float64(1.0 * J), Float64(fma(0.3333333333333333, Float64(l * l), 2.0) * l), U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.02], N[(N[(N[(N[(N[(K * K), $MachinePrecision] * -0.125 + 1.0), $MachinePrecision] * l), $MachinePrecision] * 2.0), $MachinePrecision] * J + U), $MachinePrecision], N[(N[(1.0 * J), $MachinePrecision] * N[(N[(0.3333333333333333 * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.02:\\
\;\;\;\;\mathsf{fma}\left(\left(\mathsf{fma}\left(K \cdot K, -0.125, 1\right) \cdot \ell\right) \cdot 2, J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 \cdot J, \mathsf{fma}\left(0.3333333333333333, \ell \cdot \ell, 2\right) \cdot \ell, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.0200000000000000004Initial program 88.3%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6468.2
Applied rewrites68.2%
Taylor expanded in K around 0
Applied rewrites58.1%
if -0.0200000000000000004 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.6%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6487.7
Applied rewrites87.7%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites87.7%
Taylor expanded in K around 0
Applied rewrites82.1%
(FPCore (J l K U) :precision binary64 (if (<= (cos (/ K 2.0)) -5e-310) (fma (* (* l l) 4.0) J U) (fma (* 1.0 J) (* (fma 0.3333333333333333 (* l l) 2.0) l) U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -5e-310) {
tmp = fma(((l * l) * 4.0), J, U);
} else {
tmp = fma((1.0 * J), (fma(0.3333333333333333, (l * l), 2.0) * l), U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -5e-310) tmp = fma(Float64(Float64(l * l) * 4.0), J, U); else tmp = fma(Float64(1.0 * J), Float64(fma(0.3333333333333333, Float64(l * l), 2.0) * l), U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -5e-310], N[(N[(N[(l * l), $MachinePrecision] * 4.0), $MachinePrecision] * J + U), $MachinePrecision], N[(N[(1.0 * J), $MachinePrecision] * N[(N[(0.3333333333333333 * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -5 \cdot 10^{-310}:\\
\;\;\;\;\mathsf{fma}\left(\left(\ell \cdot \ell\right) \cdot 4, J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1 \cdot J, \mathsf{fma}\left(0.3333333333333333, \ell \cdot \ell, 2\right) \cdot \ell, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -4.999999999999985e-310Initial program 88.3%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval68.3
Applied rewrites68.3%
Taylor expanded in K around 0
Applied rewrites39.8%
Applied rewrites39.8%
Applied rewrites56.4%
if -4.999999999999985e-310 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.6%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6487.7
Applied rewrites87.7%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites87.7%
Taylor expanded in K around 0
Applied rewrites82.1%
(FPCore (J l K U) :precision binary64 (if (<= K 20.0) (fma (* (* (sinh l) (fma (* K K) -0.125 1.0)) 2.0) J U) (fma (* (cos (* 0.5 K)) J) (* (fma 0.3333333333333333 (* l l) 2.0) l) U)))
double code(double J, double l, double K, double U) {
double tmp;
if (K <= 20.0) {
tmp = fma(((sinh(l) * fma((K * K), -0.125, 1.0)) * 2.0), J, U);
} else {
tmp = fma((cos((0.5 * K)) * J), (fma(0.3333333333333333, (l * l), 2.0) * l), U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (K <= 20.0) tmp = fma(Float64(Float64(sinh(l) * fma(Float64(K * K), -0.125, 1.0)) * 2.0), J, U); else tmp = fma(Float64(cos(Float64(0.5 * K)) * J), Float64(fma(0.3333333333333333, Float64(l * l), 2.0) * l), U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[K, 20.0], N[(N[(N[(N[Sinh[l], $MachinePrecision] * N[(N[(K * K), $MachinePrecision] * -0.125 + 1.0), $MachinePrecision]), $MachinePrecision] * 2.0), $MachinePrecision] * J + U), $MachinePrecision], N[(N[(N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision] * J), $MachinePrecision] * N[(N[(0.3333333333333333 * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision] * l), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;K \leq 20:\\
\;\;\;\;\mathsf{fma}\left(\left(\sinh \ell \cdot \mathsf{fma}\left(K \cdot K, -0.125, 1\right)\right) \cdot 2, J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\cos \left(0.5 \cdot K\right) \cdot J, \mathsf{fma}\left(0.3333333333333333, \ell \cdot \ell, 2\right) \cdot \ell, U\right)\\
\end{array}
\end{array}
if K < 20Initial program 85.2%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6474.4
Applied rewrites74.4%
if 20 < K Initial program 87.0%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6482.8
Applied rewrites82.8%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites82.8%
Taylor expanded in K around 0
lower-*.f6482.8
Applied rewrites82.8%
(FPCore (J l K U) :precision binary64 (if (<= K 20.0) (fma (* (* (sinh l) (fma (* K K) -0.125 1.0)) 2.0) J U) (fma (* (* (fma (* l l) 0.3333333333333333 2.0) l) (cos (* 0.5 K))) J U)))
double code(double J, double l, double K, double U) {
double tmp;
if (K <= 20.0) {
tmp = fma(((sinh(l) * fma((K * K), -0.125, 1.0)) * 2.0), J, U);
} else {
tmp = fma(((fma((l * l), 0.3333333333333333, 2.0) * l) * cos((0.5 * K))), J, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (K <= 20.0) tmp = fma(Float64(Float64(sinh(l) * fma(Float64(K * K), -0.125, 1.0)) * 2.0), J, U); else tmp = fma(Float64(Float64(fma(Float64(l * l), 0.3333333333333333, 2.0) * l) * cos(Float64(0.5 * K))), J, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[K, 20.0], N[(N[(N[(N[Sinh[l], $MachinePrecision] * N[(N[(K * K), $MachinePrecision] * -0.125 + 1.0), $MachinePrecision]), $MachinePrecision] * 2.0), $MachinePrecision] * J + U), $MachinePrecision], N[(N[(N[(N[(N[(l * l), $MachinePrecision] * 0.3333333333333333 + 2.0), $MachinePrecision] * l), $MachinePrecision] * N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * J + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;K \leq 20:\\
\;\;\;\;\mathsf{fma}\left(\left(\sinh \ell \cdot \mathsf{fma}\left(K \cdot K, -0.125, 1\right)\right) \cdot 2, J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(\mathsf{fma}\left(\ell \cdot \ell, 0.3333333333333333, 2\right) \cdot \ell\right) \cdot \cos \left(0.5 \cdot K\right), J, U\right)\\
\end{array}
\end{array}
if K < 20Initial program 85.2%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6474.4
Applied rewrites74.4%
if 20 < K Initial program 87.0%
Taylor expanded in l around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6482.8
Applied rewrites82.8%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-cos.f64N/A
lift-/.f64N/A
metadata-evalN/A
distribute-neg-frac2N/A
lift-/.f64N/A
cos-neg-revN/A
lift-cos.f64N/A
lift-*.f64N/A
Applied rewrites82.8%
Taylor expanded in K around 0
lower-*.f6482.8
Applied rewrites82.8%
lift-fma.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites82.8%
(FPCore (J l K U) :precision binary64 (fma l 2.0 U))
double code(double J, double l, double K, double U) {
return fma(l, 2.0, U);
}
function code(J, l, K, U) return fma(l, 2.0, U) end
code[J_, l_, K_, U_] := N[(l * 2.0 + U), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\ell, 2, U\right)
\end{array}
Initial program 85.6%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval66.9
Applied rewrites66.9%
Taylor expanded in K around 0
Applied rewrites55.2%
Applied rewrites55.2%
Applied rewrites34.1%
(FPCore (J l K U) :precision binary64 (fma J 2.0 U))
double code(double J, double l, double K, double U) {
return fma(J, 2.0, U);
}
function code(J, l, K, U) return fma(J, 2.0, U) end
code[J_, l_, K_, U_] := N[(J * 2.0 + U), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(J, 2, U\right)
\end{array}
Initial program 85.6%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval66.9
Applied rewrites66.9%
Taylor expanded in K around 0
Applied rewrites55.2%
Applied rewrites55.2%
Applied rewrites24.6%
herbie shell --seed 2024337
(FPCore (J l K U)
:name "Maksimov and Kolovsky, Equation (4)"
:precision binary64
(+ (* (* J (- (exp l) (exp (- l)))) (cos (/ K 2.0))) U))