
(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 21 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 (* (* J (cos (* K 0.5))) (sinh l)) 2.0 U))
double code(double J, double l, double K, double U) {
return fma(((J * cos((K * 0.5))) * sinh(l)), 2.0, U);
}
function code(J, l, K, U) return fma(Float64(Float64(J * cos(Float64(K * 0.5))) * sinh(l)), 2.0, U) end
code[J_, l_, K_, U_] := N[(N[(N[(J * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Sinh[l], $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\left(J \cdot \cos \left(K \cdot 0.5\right)\right) \cdot \sinh \ell, 2, U\right)
\end{array}
Initial program 86.3%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
(FPCore (J l K U)
:precision binary64
(let* ((t_0 (- (exp l) (exp (- l)))))
(if (<= t_0 -0.0005)
(fma (* J l) (fma -0.25 (* K K) 2.0) U)
(if (<= t_0 2e-13)
(fma J (* l 2.0) U)
(fma (* l (fma (* K K) -0.25 2.0)) J U)))))
double code(double J, double l, double K, double U) {
double t_0 = exp(l) - exp(-l);
double tmp;
if (t_0 <= -0.0005) {
tmp = fma((J * l), fma(-0.25, (K * K), 2.0), U);
} else if (t_0 <= 2e-13) {
tmp = fma(J, (l * 2.0), U);
} else {
tmp = fma((l * fma((K * K), -0.25, 2.0)), J, U);
}
return tmp;
}
function code(J, l, K, U) t_0 = Float64(exp(l) - exp(Float64(-l))) tmp = 0.0 if (t_0 <= -0.0005) tmp = fma(Float64(J * l), fma(-0.25, Float64(K * K), 2.0), U); elseif (t_0 <= 2e-13) tmp = fma(J, Float64(l * 2.0), U); else tmp = fma(Float64(l * fma(Float64(K * K), -0.25, 2.0)), J, U); end return tmp end
code[J_, l_, K_, U_] := Block[{t$95$0 = N[(N[Exp[l], $MachinePrecision] - N[Exp[(-l)], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -0.0005], N[(N[(J * l), $MachinePrecision] * N[(-0.25 * N[(K * K), $MachinePrecision] + 2.0), $MachinePrecision] + U), $MachinePrecision], If[LessEqual[t$95$0, 2e-13], N[(J * N[(l * 2.0), $MachinePrecision] + U), $MachinePrecision], N[(N[(l * N[(N[(K * K), $MachinePrecision] * -0.25 + 2.0), $MachinePrecision]), $MachinePrecision] * J + U), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{\ell} - e^{-\ell}\\
\mathbf{if}\;t\_0 \leq -0.0005:\\
\;\;\;\;\mathsf{fma}\left(J \cdot \ell, \mathsf{fma}\left(-0.25, K \cdot K, 2\right), U\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-13}:\\
\;\;\;\;\mathsf{fma}\left(J, \ell \cdot 2, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\ell \cdot \mathsf{fma}\left(K \cdot K, -0.25, 2\right), J, U\right)\\
\end{array}
\end{array}
if (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l))) < -5.0000000000000001e-4Initial program 99.8%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6418.0
Applied rewrites18.0%
Taylor expanded in K around 0
Applied rewrites35.4%
if -5.0000000000000001e-4 < (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l))) < 2.0000000000000001e-13Initial program 72.2%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6499.9
Applied rewrites99.9%
Taylor expanded in K around 0
Applied rewrites89.0%
if 2.0000000000000001e-13 < (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l))) Initial program 100.0%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6435.8
Applied rewrites35.8%
Applied rewrites35.8%
Taylor expanded in K around 0
Applied rewrites43.4%
(FPCore (J l K U)
:precision binary64
(let* ((t_0 (- (exp l) (exp (- l))))
(t_1 (fma (* J l) (fma -0.25 (* K K) 2.0) U)))
(if (<= t_0 -0.0005) t_1 (if (<= t_0 2e-13) (fma J (* l 2.0) U) t_1))))
double code(double J, double l, double K, double U) {
double t_0 = exp(l) - exp(-l);
double t_1 = fma((J * l), fma(-0.25, (K * K), 2.0), U);
double tmp;
if (t_0 <= -0.0005) {
tmp = t_1;
} else if (t_0 <= 2e-13) {
tmp = fma(J, (l * 2.0), U);
} else {
tmp = t_1;
}
return tmp;
}
function code(J, l, K, U) t_0 = Float64(exp(l) - exp(Float64(-l))) t_1 = fma(Float64(J * l), fma(-0.25, Float64(K * K), 2.0), U) tmp = 0.0 if (t_0 <= -0.0005) tmp = t_1; elseif (t_0 <= 2e-13) tmp = fma(J, Float64(l * 2.0), U); else tmp = t_1; end return tmp end
code[J_, l_, K_, U_] := Block[{t$95$0 = N[(N[Exp[l], $MachinePrecision] - N[Exp[(-l)], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(J * l), $MachinePrecision] * N[(-0.25 * N[(K * K), $MachinePrecision] + 2.0), $MachinePrecision] + U), $MachinePrecision]}, If[LessEqual[t$95$0, -0.0005], t$95$1, If[LessEqual[t$95$0, 2e-13], N[(J * N[(l * 2.0), $MachinePrecision] + U), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := e^{\ell} - e^{-\ell}\\
t_1 := \mathsf{fma}\left(J \cdot \ell, \mathsf{fma}\left(-0.25, K \cdot K, 2\right), U\right)\\
\mathbf{if}\;t\_0 \leq -0.0005:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{-13}:\\
\;\;\;\;\mathsf{fma}\left(J, \ell \cdot 2, U\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l))) < -5.0000000000000001e-4 or 2.0000000000000001e-13 < (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l))) Initial program 99.9%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6426.0
Applied rewrites26.0%
Taylor expanded in K around 0
Applied rewrites38.2%
if -5.0000000000000001e-4 < (-.f64 (exp.f64 l) (exp.f64 (neg.f64 l))) < 2.0000000000000001e-13Initial program 72.2%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6499.9
Applied rewrites99.9%
Taylor expanded in K around 0
Applied rewrites89.0%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) 0.999)
(fma
(*
(* J (cos (* K 0.5)))
(fma
(fma
(* l l)
(fma (* l l) 0.0001984126984126984 0.008333333333333333)
0.16666666666666666)
(* l (* l l))
l))
2.0
U)
(fma (* (sinh l) (* J 1.0)) 2.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= 0.999) {
tmp = fma(((J * cos((K * 0.5))) * fma(fma((l * l), fma((l * l), 0.0001984126984126984, 0.008333333333333333), 0.16666666666666666), (l * (l * l)), l)), 2.0, U);
} else {
tmp = fma((sinh(l) * (J * 1.0)), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= 0.999) tmp = fma(Float64(Float64(J * cos(Float64(K * 0.5))) * fma(fma(Float64(l * l), fma(Float64(l * l), 0.0001984126984126984, 0.008333333333333333), 0.16666666666666666), Float64(l * Float64(l * l)), l)), 2.0, U); else tmp = fma(Float64(sinh(l) * Float64(J * 1.0)), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], 0.999], N[(N[(N[(J * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(l * l), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] * 0.0001984126984126984 + 0.008333333333333333), $MachinePrecision] + 0.16666666666666666), $MachinePrecision] * N[(l * N[(l * l), $MachinePrecision]), $MachinePrecision] + l), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision], N[(N[(N[Sinh[l], $MachinePrecision] * N[(J * 1.0), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq 0.999:\\
\;\;\;\;\mathsf{fma}\left(\left(J \cdot \cos \left(K \cdot 0.5\right)\right) \cdot \mathsf{fma}\left(\mathsf{fma}\left(\ell \cdot \ell, \mathsf{fma}\left(\ell \cdot \ell, 0.0001984126984126984, 0.008333333333333333\right), 0.16666666666666666\right), \ell \cdot \left(\ell \cdot \ell\right), \ell\right), 2, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\sinh \ell \cdot \left(J \cdot 1\right), 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < 0.998999999999999999Initial program 87.3%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in l around 0
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow3N/A
*-lft-identityN/A
lower-fma.f64N/A
Applied rewrites93.6%
if 0.998999999999999999 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 85.3%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in K around 0
Applied rewrites99.7%
Final simplification96.9%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) 0.999)
(fma
(*
(* J (cos (* K 0.5)))
(fma
(* l l)
(* l (fma (* l l) 0.008333333333333333 0.16666666666666666))
l))
2.0
U)
(fma (* (sinh l) (* J 1.0)) 2.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= 0.999) {
tmp = fma(((J * cos((K * 0.5))) * fma((l * l), (l * fma((l * l), 0.008333333333333333, 0.16666666666666666)), l)), 2.0, U);
} else {
tmp = fma((sinh(l) * (J * 1.0)), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= 0.999) tmp = fma(Float64(Float64(J * cos(Float64(K * 0.5))) * fma(Float64(l * l), Float64(l * fma(Float64(l * l), 0.008333333333333333, 0.16666666666666666)), l)), 2.0, U); else tmp = fma(Float64(sinh(l) * Float64(J * 1.0)), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], 0.999], N[(N[(N[(J * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] * N[(l * N[(N[(l * l), $MachinePrecision] * 0.008333333333333333 + 0.16666666666666666), $MachinePrecision]), $MachinePrecision] + l), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision], N[(N[(N[Sinh[l], $MachinePrecision] * N[(J * 1.0), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq 0.999:\\
\;\;\;\;\mathsf{fma}\left(\left(J \cdot \cos \left(K \cdot 0.5\right)\right) \cdot \mathsf{fma}\left(\ell \cdot \ell, \ell \cdot \mathsf{fma}\left(\ell \cdot \ell, 0.008333333333333333, 0.16666666666666666\right), \ell\right), 2, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\sinh \ell \cdot \left(J \cdot 1\right), 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < 0.998999999999999999Initial program 87.3%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in l around 0
+-commutativeN/A
distribute-rgt-inN/A
*-lft-identityN/A
associate-*l*N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6491.7
Applied rewrites91.7%
if 0.998999999999999999 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 85.3%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in K around 0
Applied rewrites99.7%
Final simplification96.0%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) 0.999)
(fma
(cos (* K 0.5))
(*
l
(*
J
(fma (* l l) (fma l (* l 0.016666666666666666) 0.3333333333333333) 2.0)))
U)
(fma (* (sinh l) (* J 1.0)) 2.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= 0.999) {
tmp = fma(cos((K * 0.5)), (l * (J * fma((l * l), fma(l, (l * 0.016666666666666666), 0.3333333333333333), 2.0))), U);
} else {
tmp = fma((sinh(l) * (J * 1.0)), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= 0.999) tmp = fma(cos(Float64(K * 0.5)), Float64(l * Float64(J * fma(Float64(l * l), fma(l, Float64(l * 0.016666666666666666), 0.3333333333333333), 2.0))), U); else tmp = fma(Float64(sinh(l) * Float64(J * 1.0)), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], 0.999], N[(N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision] * N[(l * N[(J * N[(N[(l * l), $MachinePrecision] * N[(l * N[(l * 0.016666666666666666), $MachinePrecision] + 0.3333333333333333), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], N[(N[(N[Sinh[l], $MachinePrecision] * N[(J * 1.0), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq 0.999:\\
\;\;\;\;\mathsf{fma}\left(\cos \left(K \cdot 0.5\right), \ell \cdot \left(J \cdot \mathsf{fma}\left(\ell \cdot \ell, \mathsf{fma}\left(\ell, \ell \cdot 0.016666666666666666, 0.3333333333333333\right), 2\right)\right), U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\sinh \ell \cdot \left(J \cdot 1\right), 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < 0.998999999999999999Initial program 87.3%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites90.9%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6490.9
Applied rewrites90.9%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
distribute-lft-inN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
Applied rewrites91.7%
if 0.998999999999999999 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 85.3%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in K around 0
Applied rewrites99.7%
Final simplification96.0%
(FPCore (J l K U) :precision binary64 (if (<= (cos (/ K 2.0)) 0.4) (fma (* J (cos (* K 0.5))) (* l (fma l (* l 0.3333333333333333) 2.0)) U) (fma (* (sinh l) (* J 1.0)) 2.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= 0.4) {
tmp = fma((J * cos((K * 0.5))), (l * fma(l, (l * 0.3333333333333333), 2.0)), U);
} else {
tmp = fma((sinh(l) * (J * 1.0)), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= 0.4) tmp = fma(Float64(J * cos(Float64(K * 0.5))), Float64(l * fma(l, Float64(l * 0.3333333333333333), 2.0)), U); else tmp = fma(Float64(sinh(l) * Float64(J * 1.0)), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], 0.4], N[(N[(J * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(l * N[(l * N[(l * 0.3333333333333333), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], N[(N[(N[Sinh[l], $MachinePrecision] * N[(J * 1.0), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq 0.4:\\
\;\;\;\;\mathsf{fma}\left(J \cdot \cos \left(K \cdot 0.5\right), \ell \cdot \mathsf{fma}\left(\ell, \ell \cdot 0.3333333333333333, 2\right), U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\sinh \ell \cdot \left(J \cdot 1\right), 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < 0.40000000000000002Initial program 90.4%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6488.3
Applied rewrites88.3%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
associate-*r*N/A
lift-/.f64N/A
div-invN/A
metadata-evalN/A
lift-*.f64N/A
*-commutativeN/A
lift-*.f64N/A
lower-fma.f6488.3
Applied rewrites88.3%
if 0.40000000000000002 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.5%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in K around 0
Applied rewrites96.5%
Final simplification94.0%
(FPCore (J l K U) :precision binary64 (if (<= (cos (/ K 2.0)) -0.05) (fma l (* (cos (* K 0.5)) (* J (fma 0.3333333333333333 (* l l) 2.0))) U) (fma (* (sinh l) (* J 1.0)) 2.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.05) {
tmp = fma(l, (cos((K * 0.5)) * (J * fma(0.3333333333333333, (l * l), 2.0))), U);
} else {
tmp = fma((sinh(l) * (J * 1.0)), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.05) tmp = fma(l, Float64(cos(Float64(K * 0.5)) * Float64(J * fma(0.3333333333333333, Float64(l * l), 2.0))), U); else tmp = fma(Float64(sinh(l) * Float64(J * 1.0)), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.05], N[(l * N[(N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision] * N[(J * N[(0.3333333333333333 * N[(l * l), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], N[(N[(N[Sinh[l], $MachinePrecision] * N[(J * 1.0), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.05:\\
\;\;\;\;\mathsf{fma}\left(\ell, \cos \left(K \cdot 0.5\right) \cdot \left(J \cdot \mathsf{fma}\left(0.3333333333333333, \ell \cdot \ell, 2\right)\right), U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\sinh \ell \cdot \left(J \cdot 1\right), 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.050000000000000003Initial program 91.4%
Taylor expanded in l around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites87.2%
if -0.050000000000000003 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.4%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in K around 0
Applied rewrites95.7%
Final simplification93.4%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) -0.05)
(fma
(fma K (* K -0.125) 1.0)
(* l (* l (* J (* l (* (* l l) 0.016666666666666666)))))
U)
(fma (* (sinh l) (* J 1.0)) 2.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.05) {
tmp = fma(fma(K, (K * -0.125), 1.0), (l * (l * (J * (l * ((l * l) * 0.016666666666666666))))), U);
} else {
tmp = fma((sinh(l) * (J * 1.0)), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.05) tmp = fma(fma(K, Float64(K * -0.125), 1.0), Float64(l * Float64(l * Float64(J * Float64(l * Float64(Float64(l * l) * 0.016666666666666666))))), U); else tmp = fma(Float64(sinh(l) * Float64(J * 1.0)), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.05], N[(N[(K * N[(K * -0.125), $MachinePrecision] + 1.0), $MachinePrecision] * N[(l * N[(l * N[(J * N[(l * N[(N[(l * l), $MachinePrecision] * 0.016666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], N[(N[(N[Sinh[l], $MachinePrecision] * N[(J * 1.0), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.05:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(K, K \cdot -0.125, 1\right), \ell \cdot \left(\ell \cdot \left(J \cdot \left(\ell \cdot \left(\left(\ell \cdot \ell\right) \cdot 0.016666666666666666\right)\right)\right)\right), U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\sinh \ell \cdot \left(J \cdot 1\right), 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.050000000000000003Initial program 91.4%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites92.7%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6492.7
Applied rewrites92.7%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f6458.1
Applied rewrites58.1%
Taylor expanded in l around inf
Applied rewrites67.6%
if -0.050000000000000003 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.4%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in K around 0
Applied rewrites95.7%
Final simplification88.2%
(FPCore (J l K U) :precision binary64 (fma (* (cos (* K 0.5)) (* (sinh l) 2.0)) J U))
double code(double J, double l, double K, double U) {
return fma((cos((K * 0.5)) * (sinh(l) * 2.0)), J, U);
}
function code(J, l, K, U) return fma(Float64(cos(Float64(K * 0.5)) * Float64(sinh(l) * 2.0)), J, U) end
code[J_, l_, K_, U_] := N[(N[(N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision] * N[(N[Sinh[l], $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision] * J + U), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\cos \left(K \cdot 0.5\right) \cdot \left(\sinh \ell \cdot 2\right), J, U\right)
\end{array}
Initial program 86.3%
lift-+.f64N/A
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites99.9%
Final simplification99.9%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) -0.205)
(fma
(fma K (* K -0.125) 1.0)
(* l (* l (* J (* l (* (* l l) 0.016666666666666666)))))
U)
(fma
(*
(fma
(fma
(* l l)
(fma (* l l) 0.0001984126984126984 0.008333333333333333)
0.16666666666666666)
(* l (* l l))
l)
(* J 1.0))
2.0
U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.205) {
tmp = fma(fma(K, (K * -0.125), 1.0), (l * (l * (J * (l * ((l * l) * 0.016666666666666666))))), U);
} else {
tmp = fma((fma(fma((l * l), fma((l * l), 0.0001984126984126984, 0.008333333333333333), 0.16666666666666666), (l * (l * l)), l) * (J * 1.0)), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.205) tmp = fma(fma(K, Float64(K * -0.125), 1.0), Float64(l * Float64(l * Float64(J * Float64(l * Float64(Float64(l * l) * 0.016666666666666666))))), U); else tmp = fma(Float64(fma(fma(Float64(l * l), fma(Float64(l * l), 0.0001984126984126984, 0.008333333333333333), 0.16666666666666666), Float64(l * Float64(l * l)), l) * Float64(J * 1.0)), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.205], N[(N[(K * N[(K * -0.125), $MachinePrecision] + 1.0), $MachinePrecision] * N[(l * N[(l * N[(J * N[(l * N[(N[(l * l), $MachinePrecision] * 0.016666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], N[(N[(N[(N[(N[(l * l), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] * 0.0001984126984126984 + 0.008333333333333333), $MachinePrecision] + 0.16666666666666666), $MachinePrecision] * N[(l * N[(l * l), $MachinePrecision]), $MachinePrecision] + l), $MachinePrecision] * N[(J * 1.0), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.205:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(K, K \cdot -0.125, 1\right), \ell \cdot \left(\ell \cdot \left(J \cdot \left(\ell \cdot \left(\left(\ell \cdot \ell\right) \cdot 0.016666666666666666\right)\right)\right)\right), U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\ell \cdot \ell, \mathsf{fma}\left(\ell \cdot \ell, 0.0001984126984126984, 0.008333333333333333\right), 0.16666666666666666\right), \ell \cdot \left(\ell \cdot \ell\right), \ell\right) \cdot \left(J \cdot 1\right), 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.204999999999999988Initial program 90.9%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites93.6%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6493.6
Applied rewrites93.6%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f6459.8
Applied rewrites59.8%
Taylor expanded in l around inf
Applied rewrites68.5%
if -0.204999999999999988 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.7%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in l around 0
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*l*N/A
unpow2N/A
unpow3N/A
*-lft-identityN/A
lower-fma.f64N/A
Applied rewrites91.1%
Taylor expanded in K around 0
Applied rewrites86.5%
Final simplification81.9%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) -0.205)
(fma
(fma K (* K -0.125) 1.0)
(* l (* l (* J (* l (* (* l l) 0.016666666666666666)))))
U)
(fma
(*
(* J 1.0)
(fma
(* l l)
(* l (fma (* l l) 0.008333333333333333 0.16666666666666666))
l))
2.0
U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.205) {
tmp = fma(fma(K, (K * -0.125), 1.0), (l * (l * (J * (l * ((l * l) * 0.016666666666666666))))), U);
} else {
tmp = fma(((J * 1.0) * fma((l * l), (l * fma((l * l), 0.008333333333333333, 0.16666666666666666)), l)), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.205) tmp = fma(fma(K, Float64(K * -0.125), 1.0), Float64(l * Float64(l * Float64(J * Float64(l * Float64(Float64(l * l) * 0.016666666666666666))))), U); else tmp = fma(Float64(Float64(J * 1.0) * fma(Float64(l * l), Float64(l * fma(Float64(l * l), 0.008333333333333333, 0.16666666666666666)), l)), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.205], N[(N[(K * N[(K * -0.125), $MachinePrecision] + 1.0), $MachinePrecision] * N[(l * N[(l * N[(J * N[(l * N[(N[(l * l), $MachinePrecision] * 0.016666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], N[(N[(N[(J * 1.0), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] * N[(l * N[(N[(l * l), $MachinePrecision] * 0.008333333333333333 + 0.16666666666666666), $MachinePrecision]), $MachinePrecision] + l), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.205:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(K, K \cdot -0.125, 1\right), \ell \cdot \left(\ell \cdot \left(J \cdot \left(\ell \cdot \left(\left(\ell \cdot \ell\right) \cdot 0.016666666666666666\right)\right)\right)\right), U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(J \cdot 1\right) \cdot \mathsf{fma}\left(\ell \cdot \ell, \ell \cdot \mathsf{fma}\left(\ell \cdot \ell, 0.008333333333333333, 0.16666666666666666\right), \ell\right), 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.204999999999999988Initial program 90.9%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites93.6%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6493.6
Applied rewrites93.6%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f6459.8
Applied rewrites59.8%
Taylor expanded in l around inf
Applied rewrites68.5%
if -0.204999999999999988 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.7%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in K around 0
Applied rewrites94.8%
Taylor expanded in l around 0
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
*-lft-identityN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6485.9
Applied rewrites85.9%
Final simplification81.5%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) -0.15)
(fma
(fma K (* K -0.125) 1.0)
(* l (* J (fma l (* l 0.3333333333333333) 2.0)))
U)
(fma
(*
(* J 1.0)
(fma
(* l l)
(* l (fma (* l l) 0.008333333333333333 0.16666666666666666))
l))
2.0
U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.15) {
tmp = fma(fma(K, (K * -0.125), 1.0), (l * (J * fma(l, (l * 0.3333333333333333), 2.0))), U);
} else {
tmp = fma(((J * 1.0) * fma((l * l), (l * fma((l * l), 0.008333333333333333, 0.16666666666666666)), l)), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.15) tmp = fma(fma(K, Float64(K * -0.125), 1.0), Float64(l * Float64(J * fma(l, Float64(l * 0.3333333333333333), 2.0))), U); else tmp = fma(Float64(Float64(J * 1.0) * fma(Float64(l * l), Float64(l * fma(Float64(l * l), 0.008333333333333333, 0.16666666666666666)), l)), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.15], N[(N[(K * N[(K * -0.125), $MachinePrecision] + 1.0), $MachinePrecision] * N[(l * N[(J * N[(l * N[(l * 0.3333333333333333), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], N[(N[(N[(J * 1.0), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] * N[(l * N[(N[(l * l), $MachinePrecision] * 0.008333333333333333 + 0.16666666666666666), $MachinePrecision]), $MachinePrecision] + l), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.15:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(K, K \cdot -0.125, 1\right), \ell \cdot \left(J \cdot \mathsf{fma}\left(\ell, \ell \cdot 0.3333333333333333, 2\right)\right), U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(J \cdot 1\right) \cdot \mathsf{fma}\left(\ell \cdot \ell, \ell \cdot \mathsf{fma}\left(\ell \cdot \ell, 0.008333333333333333, 0.16666666666666666\right), \ell\right), 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.149999999999999994Initial program 91.2%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites93.8%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6493.8
Applied rewrites93.8%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f6459.6
Applied rewrites59.6%
Taylor expanded in l around 0
lower-*.f64N/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6459.6
Applied rewrites59.6%
if -0.149999999999999994 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.5%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in K around 0
Applied rewrites95.3%
Taylor expanded in l around 0
+-commutativeN/A
distribute-rgt-inN/A
associate-*l*N/A
*-lft-identityN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6486.2
Applied rewrites86.2%
Final simplification79.3%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) -0.15)
(fma
(fma K (* K -0.125) 1.0)
(* l (* J (fma l (* l 0.3333333333333333) 2.0)))
U)
(fma
1.0
(*
l
(*
J
(fma l (* l (fma (* l l) 0.016666666666666666 0.3333333333333333)) 2.0)))
U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.15) {
tmp = fma(fma(K, (K * -0.125), 1.0), (l * (J * fma(l, (l * 0.3333333333333333), 2.0))), U);
} else {
tmp = fma(1.0, (l * (J * fma(l, (l * fma((l * l), 0.016666666666666666, 0.3333333333333333)), 2.0))), U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.15) tmp = fma(fma(K, Float64(K * -0.125), 1.0), Float64(l * Float64(J * fma(l, Float64(l * 0.3333333333333333), 2.0))), U); else tmp = fma(1.0, Float64(l * Float64(J * fma(l, Float64(l * fma(Float64(l * l), 0.016666666666666666, 0.3333333333333333)), 2.0))), U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.15], N[(N[(K * N[(K * -0.125), $MachinePrecision] + 1.0), $MachinePrecision] * N[(l * N[(J * N[(l * N[(l * 0.3333333333333333), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], N[(1.0 * N[(l * N[(J * N[(l * N[(l * N[(N[(l * l), $MachinePrecision] * 0.016666666666666666 + 0.3333333333333333), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.15:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(K, K \cdot -0.125, 1\right), \ell \cdot \left(J \cdot \mathsf{fma}\left(\ell, \ell \cdot 0.3333333333333333, 2\right)\right), U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1, \ell \cdot \left(J \cdot \mathsf{fma}\left(\ell, \ell \cdot \mathsf{fma}\left(\ell \cdot \ell, 0.016666666666666666, 0.3333333333333333\right), 2\right)\right), U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.149999999999999994Initial program 91.2%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites93.8%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6493.8
Applied rewrites93.8%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f6459.6
Applied rewrites59.6%
Taylor expanded in l around 0
lower-*.f64N/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6459.6
Applied rewrites59.6%
if -0.149999999999999994 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.5%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites86.5%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6486.5
Applied rewrites86.5%
Taylor expanded in K around 0
Applied rewrites82.3%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-inN/A
lower-*.f64N/A
+-commutativeN/A
Applied rewrites84.3%
(FPCore (J l K U)
:precision binary64
(if (<= (cos (/ K 2.0)) -0.15)
(fma (* l (fma (* K K) -0.25 2.0)) J U)
(fma
1.0
(*
l
(*
J
(fma l (* l (fma (* l l) 0.016666666666666666 0.3333333333333333)) 2.0)))
U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.15) {
tmp = fma((l * fma((K * K), -0.25, 2.0)), J, U);
} else {
tmp = fma(1.0, (l * (J * fma(l, (l * fma((l * l), 0.016666666666666666, 0.3333333333333333)), 2.0))), U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.15) tmp = fma(Float64(l * fma(Float64(K * K), -0.25, 2.0)), J, U); else tmp = fma(1.0, Float64(l * Float64(J * fma(l, Float64(l * fma(Float64(l * l), 0.016666666666666666, 0.3333333333333333)), 2.0))), U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.15], N[(N[(l * N[(N[(K * K), $MachinePrecision] * -0.25 + 2.0), $MachinePrecision]), $MachinePrecision] * J + U), $MachinePrecision], N[(1.0 * N[(l * N[(J * N[(l * N[(l * N[(N[(l * l), $MachinePrecision] * 0.016666666666666666 + 0.3333333333333333), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.15:\\
\;\;\;\;\mathsf{fma}\left(\ell \cdot \mathsf{fma}\left(K \cdot K, -0.25, 2\right), J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1, \ell \cdot \left(J \cdot \mathsf{fma}\left(\ell, \ell \cdot \mathsf{fma}\left(\ell \cdot \ell, 0.016666666666666666, 0.3333333333333333\right), 2\right)\right), U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.149999999999999994Initial program 91.2%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6460.0
Applied rewrites60.0%
Applied rewrites60.1%
Taylor expanded in K around 0
Applied rewrites55.4%
if -0.149999999999999994 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.5%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites86.5%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6486.5
Applied rewrites86.5%
Taylor expanded in K around 0
Applied rewrites82.3%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-inN/A
lower-*.f64N/A
+-commutativeN/A
Applied rewrites84.3%
(FPCore (J l K U) :precision binary64 (if (<= (cos (/ K 2.0)) -0.15) (fma (* l (fma (* K K) -0.25 2.0)) J U) (fma (* (* J 1.0) (fma l (* (* l l) 0.16666666666666666) l)) 2.0 U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.15) {
tmp = fma((l * fma((K * K), -0.25, 2.0)), J, U);
} else {
tmp = fma(((J * 1.0) * fma(l, ((l * l) * 0.16666666666666666), l)), 2.0, U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.15) tmp = fma(Float64(l * fma(Float64(K * K), -0.25, 2.0)), J, U); else tmp = fma(Float64(Float64(J * 1.0) * fma(l, Float64(Float64(l * l) * 0.16666666666666666), l)), 2.0, U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.15], N[(N[(l * N[(N[(K * K), $MachinePrecision] * -0.25 + 2.0), $MachinePrecision]), $MachinePrecision] * J + U), $MachinePrecision], N[(N[(N[(J * 1.0), $MachinePrecision] * N[(l * N[(N[(l * l), $MachinePrecision] * 0.16666666666666666), $MachinePrecision] + l), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.15:\\
\;\;\;\;\mathsf{fma}\left(\ell \cdot \mathsf{fma}\left(K \cdot K, -0.25, 2\right), J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\left(J \cdot 1\right) \cdot \mathsf{fma}\left(\ell, \left(\ell \cdot \ell\right) \cdot 0.16666666666666666, \ell\right), 2, U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.149999999999999994Initial program 91.2%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6460.0
Applied rewrites60.0%
Applied rewrites60.1%
Taylor expanded in K around 0
Applied rewrites55.4%
if -0.149999999999999994 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.5%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in K around 0
Applied rewrites95.3%
Taylor expanded in l around 0
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6480.8
Applied rewrites80.8%
(FPCore (J l K U) :precision binary64 (if (<= (cos (/ K 2.0)) -0.15) (fma (* l (fma (* K K) -0.25 2.0)) J U) (fma 1.0 (* l (fma l (* J (* l 0.3333333333333333)) (* J 2.0))) U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.15) {
tmp = fma((l * fma((K * K), -0.25, 2.0)), J, U);
} else {
tmp = fma(1.0, (l * fma(l, (J * (l * 0.3333333333333333)), (J * 2.0))), U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.15) tmp = fma(Float64(l * fma(Float64(K * K), -0.25, 2.0)), J, U); else tmp = fma(1.0, Float64(l * fma(l, Float64(J * Float64(l * 0.3333333333333333)), Float64(J * 2.0))), U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.15], N[(N[(l * N[(N[(K * K), $MachinePrecision] * -0.25 + 2.0), $MachinePrecision]), $MachinePrecision] * J + U), $MachinePrecision], N[(1.0 * N[(l * N[(l * N[(J * N[(l * 0.3333333333333333), $MachinePrecision]), $MachinePrecision] + N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.15:\\
\;\;\;\;\mathsf{fma}\left(\ell \cdot \mathsf{fma}\left(K \cdot K, -0.25, 2\right), J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1, \ell \cdot \mathsf{fma}\left(\ell, J \cdot \left(\ell \cdot 0.3333333333333333\right), J \cdot 2\right), U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.149999999999999994Initial program 91.2%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6460.0
Applied rewrites60.0%
Applied rewrites60.1%
Taylor expanded in K around 0
Applied rewrites55.4%
if -0.149999999999999994 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.5%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites86.5%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6486.5
Applied rewrites86.5%
Taylor expanded in K around 0
Applied rewrites82.3%
Taylor expanded in l around 0
Applied rewrites77.9%
(FPCore (J l K U) :precision binary64 (if (<= (cos (/ K 2.0)) -0.15) (fma (* l (fma (* K K) -0.25 2.0)) J U) (fma 1.0 (* l (* J (fma l (* l 0.3333333333333333) 2.0))) U)))
double code(double J, double l, double K, double U) {
double tmp;
if (cos((K / 2.0)) <= -0.15) {
tmp = fma((l * fma((K * K), -0.25, 2.0)), J, U);
} else {
tmp = fma(1.0, (l * (J * fma(l, (l * 0.3333333333333333), 2.0))), U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (cos(Float64(K / 2.0)) <= -0.15) tmp = fma(Float64(l * fma(Float64(K * K), -0.25, 2.0)), J, U); else tmp = fma(1.0, Float64(l * Float64(J * fma(l, Float64(l * 0.3333333333333333), 2.0))), U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision], -0.15], N[(N[(l * N[(N[(K * K), $MachinePrecision] * -0.25 + 2.0), $MachinePrecision]), $MachinePrecision] * J + U), $MachinePrecision], N[(1.0 * N[(l * N[(J * N[(l * N[(l * 0.3333333333333333), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\cos \left(\frac{K}{2}\right) \leq -0.15:\\
\;\;\;\;\mathsf{fma}\left(\ell \cdot \mathsf{fma}\left(K \cdot K, -0.25, 2\right), J, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(1, \ell \cdot \left(J \cdot \mathsf{fma}\left(\ell, \ell \cdot 0.3333333333333333, 2\right)\right), U\right)\\
\end{array}
\end{array}
if (cos.f64 (/.f64 K #s(literal 2 binary64))) < -0.149999999999999994Initial program 91.2%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6460.0
Applied rewrites60.0%
Applied rewrites60.1%
Taylor expanded in K around 0
Applied rewrites55.4%
if -0.149999999999999994 < (cos.f64 (/.f64 K #s(literal 2 binary64))) Initial program 84.5%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites86.5%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6486.5
Applied rewrites86.5%
Taylor expanded in K around 0
Applied rewrites82.3%
Taylor expanded in l around 0
lower-*.f64N/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6477.9
Applied rewrites77.9%
(FPCore (J l K U)
:precision binary64
(if (<= l -3.2)
(fma (* (sinh l) (* J (fma -0.125 (* K K) 1.0))) 2.0 U)
(if (<= l 9000000000.0)
(fma (* (* J (cos (* K 0.5))) l) 2.0 U)
(if (<= l 9.5e+120)
(fma (* (sinh l) (* J 1.0)) 2.0 U)
(fma
(fma K (* K -0.125) 1.0)
(* l (* J (fma l (* l 0.3333333333333333) 2.0)))
U)))))
double code(double J, double l, double K, double U) {
double tmp;
if (l <= -3.2) {
tmp = fma((sinh(l) * (J * fma(-0.125, (K * K), 1.0))), 2.0, U);
} else if (l <= 9000000000.0) {
tmp = fma(((J * cos((K * 0.5))) * l), 2.0, U);
} else if (l <= 9.5e+120) {
tmp = fma((sinh(l) * (J * 1.0)), 2.0, U);
} else {
tmp = fma(fma(K, (K * -0.125), 1.0), (l * (J * fma(l, (l * 0.3333333333333333), 2.0))), U);
}
return tmp;
}
function code(J, l, K, U) tmp = 0.0 if (l <= -3.2) tmp = fma(Float64(sinh(l) * Float64(J * fma(-0.125, Float64(K * K), 1.0))), 2.0, U); elseif (l <= 9000000000.0) tmp = fma(Float64(Float64(J * cos(Float64(K * 0.5))) * l), 2.0, U); elseif (l <= 9.5e+120) tmp = fma(Float64(sinh(l) * Float64(J * 1.0)), 2.0, U); else tmp = fma(fma(K, Float64(K * -0.125), 1.0), Float64(l * Float64(J * fma(l, Float64(l * 0.3333333333333333), 2.0))), U); end return tmp end
code[J_, l_, K_, U_] := If[LessEqual[l, -3.2], N[(N[(N[Sinh[l], $MachinePrecision] * N[(J * N[(-0.125 * N[(K * K), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision], If[LessEqual[l, 9000000000.0], N[(N[(N[(J * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * l), $MachinePrecision] * 2.0 + U), $MachinePrecision], If[LessEqual[l, 9.5e+120], N[(N[(N[Sinh[l], $MachinePrecision] * N[(J * 1.0), $MachinePrecision]), $MachinePrecision] * 2.0 + U), $MachinePrecision], N[(N[(K * N[(K * -0.125), $MachinePrecision] + 1.0), $MachinePrecision] * N[(l * N[(J * N[(l * N[(l * 0.3333333333333333), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\ell \leq -3.2:\\
\;\;\;\;\mathsf{fma}\left(\sinh \ell \cdot \left(J \cdot \mathsf{fma}\left(-0.125, K \cdot K, 1\right)\right), 2, U\right)\\
\mathbf{elif}\;\ell \leq 9000000000:\\
\;\;\;\;\mathsf{fma}\left(\left(J \cdot \cos \left(K \cdot 0.5\right)\right) \cdot \ell, 2, U\right)\\
\mathbf{elif}\;\ell \leq 9.5 \cdot 10^{+120}:\\
\;\;\;\;\mathsf{fma}\left(\sinh \ell \cdot \left(J \cdot 1\right), 2, U\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(K, K \cdot -0.125, 1\right), \ell \cdot \left(J \cdot \mathsf{fma}\left(\ell, \ell \cdot 0.3333333333333333, 2\right)\right), U\right)\\
\end{array}
\end{array}
if l < -3.2000000000000002Initial program 100.0%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in K around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6480.0
Applied rewrites80.0%
if -3.2000000000000002 < l < 9e9Initial program 72.9%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6499.1
Applied rewrites99.1%
Applied rewrites99.1%
if 9e9 < l < 9.5e120Initial program 100.0%
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
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
Applied rewrites100.0%
Taylor expanded in K around 0
Applied rewrites90.5%
if 9.5e120 < l Initial program 100.0%
Taylor expanded in l around 0
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
Applied rewrites100.0%
lift-+.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64100.0
Applied rewrites100.0%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f6488.6
Applied rewrites88.6%
Taylor expanded in l around 0
lower-*.f64N/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6488.6
Applied rewrites88.6%
Final simplification91.7%
(FPCore (J l K U) :precision binary64 (fma J (* l 2.0) U))
double code(double J, double l, double K, double U) {
return fma(J, (l * 2.0), U);
}
function code(J, l, K, U) return fma(J, Float64(l * 2.0), U) end
code[J_, l_, K_, U_] := N[(J * N[(l * 2.0), $MachinePrecision] + U), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(J, \ell \cdot 2, U\right)
\end{array}
Initial program 86.3%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6462.4
Applied rewrites62.4%
Taylor expanded in K around 0
Applied rewrites54.3%
(FPCore (J l K U) :precision binary64 (* l (* J 2.0)))
double code(double J, double l, double K, double U) {
return l * (J * 2.0);
}
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 = l * (j * 2.0d0)
end function
public static double code(double J, double l, double K, double U) {
return l * (J * 2.0);
}
def code(J, l, K, U): return l * (J * 2.0)
function code(J, l, K, U) return Float64(l * Float64(J * 2.0)) end
function tmp = code(J, l, K, U) tmp = l * (J * 2.0); end
code[J_, l_, K_, U_] := N[(l * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\ell \cdot \left(J \cdot 2\right)
\end{array}
Initial program 86.3%
Taylor expanded in l around 0
+-commutativeN/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
associate-*l*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6462.4
Applied rewrites62.4%
Taylor expanded in K around 0
Applied rewrites54.3%
Taylor expanded in J around inf
Applied rewrites20.0%
herbie shell --seed 2024222
(FPCore (J l K U)
:name "Maksimov and Kolovsky, Equation (4)"
:precision binary64
(+ (* (* J (- (exp l) (exp (- l)))) (cos (/ K 2.0))) U))