
(FPCore (J K U) :precision binary64 (let* ((t_0 (cos (/ K 2.0)))) (* (* (* -2.0 J) t_0) (sqrt (+ 1.0 (pow (/ U (* (* 2.0 J) t_0)) 2.0))))))
double code(double J, double K, double U) {
double t_0 = cos((K / 2.0));
return ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U / ((2.0 * J) * t_0)), 2.0)));
}
real(8) function code(j, k, u)
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8), intent (in) :: u
real(8) :: t_0
t_0 = cos((k / 2.0d0))
code = (((-2.0d0) * j) * t_0) * sqrt((1.0d0 + ((u / ((2.0d0 * j) * t_0)) ** 2.0d0)))
end function
public static double code(double J, double K, double U) {
double t_0 = Math.cos((K / 2.0));
return ((-2.0 * J) * t_0) * Math.sqrt((1.0 + Math.pow((U / ((2.0 * J) * t_0)), 2.0)));
}
def code(J, K, U): t_0 = math.cos((K / 2.0)) return ((-2.0 * J) * t_0) * math.sqrt((1.0 + math.pow((U / ((2.0 * J) * t_0)), 2.0)))
function code(J, K, U) t_0 = cos(Float64(K / 2.0)) return Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) end
function tmp = code(J, K, U) t_0 = cos((K / 2.0)); tmp = ((-2.0 * J) * t_0) * sqrt((1.0 + ((U / ((2.0 * J) * t_0)) ^ 2.0))); end
code[J_, K_, U_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
\left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 13 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (J K U) :precision binary64 (let* ((t_0 (cos (/ K 2.0)))) (* (* (* -2.0 J) t_0) (sqrt (+ 1.0 (pow (/ U (* (* 2.0 J) t_0)) 2.0))))))
double code(double J, double K, double U) {
double t_0 = cos((K / 2.0));
return ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U / ((2.0 * J) * t_0)), 2.0)));
}
real(8) function code(j, k, u)
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8), intent (in) :: u
real(8) :: t_0
t_0 = cos((k / 2.0d0))
code = (((-2.0d0) * j) * t_0) * sqrt((1.0d0 + ((u / ((2.0d0 * j) * t_0)) ** 2.0d0)))
end function
public static double code(double J, double K, double U) {
double t_0 = Math.cos((K / 2.0));
return ((-2.0 * J) * t_0) * Math.sqrt((1.0 + Math.pow((U / ((2.0 * J) * t_0)), 2.0)));
}
def code(J, K, U): t_0 = math.cos((K / 2.0)) return ((-2.0 * J) * t_0) * math.sqrt((1.0 + math.pow((U / ((2.0 * J) * t_0)), 2.0)))
function code(J, K, U) t_0 = cos(Float64(K / 2.0)) return Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) end
function tmp = code(J, K, U) t_0 = cos((K / 2.0)); tmp = ((-2.0 * J) * t_0) * sqrt((1.0 + ((U / ((2.0 * J) * t_0)) ^ 2.0))); end
code[J_, K_, U_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
\left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}
\end{array}
\end{array}
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))))
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 5e+292)
(*
(* (* -2.0 J) (cos (* 0.5 K)))
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) (cos (* -0.5 K)))) 2.0))))
(* -1.0 (- U_m))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= 5e+292) {
tmp = ((-2.0 * J) * cos((0.5 * K))) * sqrt((1.0 + pow((U_m / ((2.0 * J) * cos((-0.5 * K)))), 2.0)));
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = Math.abs(U);
public static double code(double J, double K, double U_m) {
double t_0 = Math.cos((K / 2.0));
double t_1 = ((-2.0 * J) * t_0) * Math.sqrt((1.0 + Math.pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_1 <= -Double.POSITIVE_INFINITY) {
tmp = -U_m;
} else if (t_1 <= 5e+292) {
tmp = ((-2.0 * J) * Math.cos((0.5 * K))) * Math.sqrt((1.0 + Math.pow((U_m / ((2.0 * J) * Math.cos((-0.5 * K)))), 2.0)));
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = math.fabs(U) def code(J, K, U_m): t_0 = math.cos((K / 2.0)) t_1 = ((-2.0 * J) * t_0) * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J) * t_0)), 2.0))) tmp = 0 if t_1 <= -math.inf: tmp = -U_m elif t_1 <= 5e+292: tmp = ((-2.0 * J) * math.cos((0.5 * K))) * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J) * math.cos((-0.5 * K)))), 2.0))) else: tmp = -1.0 * -U_m return tmp
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= 5e+292) tmp = Float64(Float64(Float64(-2.0 * J) * cos(Float64(0.5 * K))) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * cos(Float64(-0.5 * K)))) ^ 2.0)))); else tmp = Float64(-1.0 * Float64(-U_m)); end return tmp end
U_m = abs(U); function tmp_2 = code(J, K, U_m) t_0 = cos((K / 2.0)); t_1 = ((-2.0 * J) * t_0) * sqrt((1.0 + ((U_m / ((2.0 * J) * t_0)) ^ 2.0))); tmp = 0.0; if (t_1 <= -Inf) tmp = -U_m; elseif (t_1 <= 5e+292) tmp = ((-2.0 * J) * cos((0.5 * K))) * sqrt((1.0 + ((U_m / ((2.0 * J) * cos((-0.5 * K)))) ^ 2.0))); else tmp = -1.0 * -U_m; end tmp_2 = tmp; end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, 5e+292], N[(N[(N[(-2.0 * J), $MachinePrecision] * N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * N[Cos[N[(-0.5 * K), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(-1.0 * (-U$95$m)), $MachinePrecision]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{+292}:\\
\;\;\;\;\left(\left(-2 \cdot J\right) \cdot \cos \left(0.5 \cdot K\right)\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot \cos \left(-0.5 \cdot K\right)}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -inf.0Initial program 6.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6449.0
Applied rewrites49.0%
if -inf.0 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < 4.9999999999999996e292Initial program 99.8%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.8
Applied rewrites99.8%
Taylor expanded in K around inf
metadata-evalN/A
distribute-lft-neg-inN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.8
Applied rewrites99.8%
if 4.9999999999999996e292 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 5.5%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f645.5
Applied rewrites5.5%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f644.1
Applied rewrites4.1%
Taylor expanded in U around -inf
Applied rewrites54.1%
Taylor expanded in J around 0
Applied rewrites54.1%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1 (* (* -2.0 J) t_0))
(t_2 (* t_1 (sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))))
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 -2e-118)
(* t_1 (sqrt (fma (* U_m (/ U_m (* J J))) 0.25 1.0)))
(if (<= t_2 -1e-218)
(fma
(* (* J J) -2.0)
(/ (- 0.5 (* -0.5 (cos (* 2.0 (* 0.5 K))))) U_m)
(- U_m))
(if (<= t_2 5e+292)
(* (* (* -2.0 J) (cos (* 0.5 K))) 1.0)
(* -1.0 (- U_m))))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = (-2.0 * J) * t_0;
double t_2 = t_1 * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= -2e-118) {
tmp = t_1 * sqrt(fma((U_m * (U_m / (J * J))), 0.25, 1.0));
} else if (t_2 <= -1e-218) {
tmp = fma(((J * J) * -2.0), ((0.5 - (-0.5 * cos((2.0 * (0.5 * K))))) / U_m), -U_m);
} else if (t_2 <= 5e+292) {
tmp = ((-2.0 * J) * cos((0.5 * K))) * 1.0;
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(-2.0 * J) * t_0) t_2 = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= -2e-118) tmp = Float64(t_1 * sqrt(fma(Float64(U_m * Float64(U_m / Float64(J * J))), 0.25, 1.0))); elseif (t_2 <= -1e-218) tmp = fma(Float64(Float64(J * J) * -2.0), Float64(Float64(0.5 - Float64(-0.5 * cos(Float64(2.0 * Float64(0.5 * K))))) / U_m), Float64(-U_m)); elseif (t_2 <= 5e+292) tmp = Float64(Float64(Float64(-2.0 * J) * cos(Float64(0.5 * K))) * 1.0); else tmp = Float64(-1.0 * Float64(-U_m)); end return tmp end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, -2e-118], N[(t$95$1 * N[Sqrt[N[(N[(U$95$m * N[(U$95$m / N[(J * J), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, -1e-218], N[(N[(N[(J * J), $MachinePrecision] * -2.0), $MachinePrecision] * N[(N[(0.5 - N[(-0.5 * N[Cos[N[(2.0 * N[(0.5 * K), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / U$95$m), $MachinePrecision] + (-U$95$m)), $MachinePrecision], If[LessEqual[t$95$2, 5e+292], N[(N[(N[(-2.0 * J), $MachinePrecision] * N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * 1.0), $MachinePrecision], N[(-1.0 * (-U$95$m)), $MachinePrecision]]]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(-2 \cdot J\right) \cdot t\_0\\
t_2 := t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq -2 \cdot 10^{-118}:\\
\;\;\;\;t\_1 \cdot \sqrt{\mathsf{fma}\left(U\_m \cdot \frac{U\_m}{J \cdot J}, 0.25, 1\right)}\\
\mathbf{elif}\;t\_2 \leq -1 \cdot 10^{-218}:\\
\;\;\;\;\mathsf{fma}\left(\left(J \cdot J\right) \cdot -2, \frac{0.5 - -0.5 \cdot \cos \left(2 \cdot \left(0.5 \cdot K\right)\right)}{U\_m}, -U\_m\right)\\
\mathbf{elif}\;t\_2 \leq 5 \cdot 10^{+292}:\\
\;\;\;\;\left(\left(-2 \cdot J\right) \cdot \cos \left(0.5 \cdot K\right)\right) \cdot 1\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -inf.0Initial program 6.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6449.0
Applied rewrites49.0%
if -inf.0 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1.99999999999999997e-118Initial program 99.7%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.7
Applied rewrites99.7%
Taylor expanded in K around 0
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6481.7
Applied rewrites81.7%
if -1.99999999999999997e-118 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1e-218Initial program 99.9%
Taylor expanded in J around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6421.1
Applied rewrites21.1%
Applied rewrites21.1%
if -1e-218 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < 4.9999999999999996e292Initial program 99.8%
Taylor expanded in J around inf
Applied rewrites63.0%
Taylor expanded in K around inf
metadata-evalN/A
distribute-lft-neg-inN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6463.0
Applied rewrites63.0%
if 4.9999999999999996e292 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 5.5%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f645.5
Applied rewrites5.5%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f644.1
Applied rewrites4.1%
Taylor expanded in U around -inf
Applied rewrites54.1%
Taylor expanded in J around 0
Applied rewrites54.1%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))))
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -2e-118)
(* (sqrt (fma (* U_m (/ (/ U_m J) J)) 0.25 1.0)) (* -2.0 J))
(if (<= t_1 -1e-218)
(fma
(* (* J J) -2.0)
(/ (- 0.5 (* -0.5 (cos (* 2.0 (* 0.5 K))))) U_m)
(- U_m))
(if (<= t_1 5e+292)
(* (* (* -2.0 J) (cos (* 0.5 K))) 1.0)
(* -1.0 (- U_m))))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -2e-118) {
tmp = sqrt(fma((U_m * ((U_m / J) / J)), 0.25, 1.0)) * (-2.0 * J);
} else if (t_1 <= -1e-218) {
tmp = fma(((J * J) * -2.0), ((0.5 - (-0.5 * cos((2.0 * (0.5 * K))))) / U_m), -U_m);
} else if (t_1 <= 5e+292) {
tmp = ((-2.0 * J) * cos((0.5 * K))) * 1.0;
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -2e-118) tmp = Float64(sqrt(fma(Float64(U_m * Float64(Float64(U_m / J) / J)), 0.25, 1.0)) * Float64(-2.0 * J)); elseif (t_1 <= -1e-218) tmp = fma(Float64(Float64(J * J) * -2.0), Float64(Float64(0.5 - Float64(-0.5 * cos(Float64(2.0 * Float64(0.5 * K))))) / U_m), Float64(-U_m)); elseif (t_1 <= 5e+292) tmp = Float64(Float64(Float64(-2.0 * J) * cos(Float64(0.5 * K))) * 1.0); else tmp = Float64(-1.0 * Float64(-U_m)); end return tmp end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -2e-118], N[(N[Sqrt[N[(N[(U$95$m * N[(N[(U$95$m / J), $MachinePrecision] / J), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision] * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -1e-218], N[(N[(N[(J * J), $MachinePrecision] * -2.0), $MachinePrecision] * N[(N[(0.5 - N[(-0.5 * N[Cos[N[(2.0 * N[(0.5 * K), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / U$95$m), $MachinePrecision] + (-U$95$m)), $MachinePrecision], If[LessEqual[t$95$1, 5e+292], N[(N[(N[(-2.0 * J), $MachinePrecision] * N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * 1.0), $MachinePrecision], N[(-1.0 * (-U$95$m)), $MachinePrecision]]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-118}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(U\_m \cdot \frac{\frac{U\_m}{J}}{J}, 0.25, 1\right)} \cdot \left(-2 \cdot J\right)\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-218}:\\
\;\;\;\;\mathsf{fma}\left(\left(J \cdot J\right) \cdot -2, \frac{0.5 - -0.5 \cdot \cos \left(2 \cdot \left(0.5 \cdot K\right)\right)}{U\_m}, -U\_m\right)\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{+292}:\\
\;\;\;\;\left(\left(-2 \cdot J\right) \cdot \cos \left(0.5 \cdot K\right)\right) \cdot 1\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -inf.0Initial program 6.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6449.0
Applied rewrites49.0%
if -inf.0 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1.99999999999999997e-118Initial program 99.7%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.7
Applied rewrites99.7%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6453.4
Applied rewrites53.4%
Applied rewrites59.6%
if -1.99999999999999997e-118 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1e-218Initial program 99.9%
Taylor expanded in J around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6421.1
Applied rewrites21.1%
Applied rewrites21.1%
if -1e-218 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < 4.9999999999999996e292Initial program 99.8%
Taylor expanded in J around inf
Applied rewrites63.0%
Taylor expanded in K around inf
metadata-evalN/A
distribute-lft-neg-inN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6463.0
Applied rewrites63.0%
if 4.9999999999999996e292 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 5.5%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f645.5
Applied rewrites5.5%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f644.1
Applied rewrites4.1%
Taylor expanded in U around -inf
Applied rewrites54.1%
Taylor expanded in J around 0
Applied rewrites54.1%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))))
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -2e-118)
(* (sqrt (fma (* U_m (/ (/ U_m J) J)) 0.25 1.0)) (* -2.0 J))
(if (<= t_1 -1e-218)
(- (* (/ (* J J) U_m) -2.0) U_m)
(if (<= t_1 5e+292)
(* (* (* -2.0 J) (cos (* 0.5 K))) 1.0)
(* -1.0 (- U_m))))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -2e-118) {
tmp = sqrt(fma((U_m * ((U_m / J) / J)), 0.25, 1.0)) * (-2.0 * J);
} else if (t_1 <= -1e-218) {
tmp = (((J * J) / U_m) * -2.0) - U_m;
} else if (t_1 <= 5e+292) {
tmp = ((-2.0 * J) * cos((0.5 * K))) * 1.0;
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -2e-118) tmp = Float64(sqrt(fma(Float64(U_m * Float64(Float64(U_m / J) / J)), 0.25, 1.0)) * Float64(-2.0 * J)); elseif (t_1 <= -1e-218) tmp = Float64(Float64(Float64(Float64(J * J) / U_m) * -2.0) - U_m); elseif (t_1 <= 5e+292) tmp = Float64(Float64(Float64(-2.0 * J) * cos(Float64(0.5 * K))) * 1.0); else tmp = Float64(-1.0 * Float64(-U_m)); end return tmp end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -2e-118], N[(N[Sqrt[N[(N[(U$95$m * N[(N[(U$95$m / J), $MachinePrecision] / J), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision] * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -1e-218], N[(N[(N[(N[(J * J), $MachinePrecision] / U$95$m), $MachinePrecision] * -2.0), $MachinePrecision] - U$95$m), $MachinePrecision], If[LessEqual[t$95$1, 5e+292], N[(N[(N[(-2.0 * J), $MachinePrecision] * N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * 1.0), $MachinePrecision], N[(-1.0 * (-U$95$m)), $MachinePrecision]]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-118}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(U\_m \cdot \frac{\frac{U\_m}{J}}{J}, 0.25, 1\right)} \cdot \left(-2 \cdot J\right)\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-218}:\\
\;\;\;\;\frac{J \cdot J}{U\_m} \cdot -2 - U\_m\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{+292}:\\
\;\;\;\;\left(\left(-2 \cdot J\right) \cdot \cos \left(0.5 \cdot K\right)\right) \cdot 1\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -inf.0Initial program 6.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6449.0
Applied rewrites49.0%
if -inf.0 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1.99999999999999997e-118Initial program 99.7%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.7
Applied rewrites99.7%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6453.4
Applied rewrites53.4%
Applied rewrites59.6%
if -1.99999999999999997e-118 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1e-218Initial program 99.9%
Taylor expanded in J around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6421.1
Applied rewrites21.1%
Taylor expanded in K around 0
Applied rewrites21.1%
if -1e-218 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < 4.9999999999999996e292Initial program 99.8%
Taylor expanded in J around inf
Applied rewrites63.0%
Taylor expanded in K around inf
metadata-evalN/A
distribute-lft-neg-inN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6463.0
Applied rewrites63.0%
if 4.9999999999999996e292 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 5.5%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f645.5
Applied rewrites5.5%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f644.1
Applied rewrites4.1%
Taylor expanded in U around -inf
Applied rewrites54.1%
Taylor expanded in J around 0
Applied rewrites54.1%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))))
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -2e-118)
(* (sqrt (fma (* U_m (/ (/ U_m J) J)) 0.25 1.0)) (* -2.0 J))
(if (<= t_1 -5e-289)
(- (* (/ (* J J) U_m) -2.0) U_m)
(* (- (* (/ (* (/ J U_m) J) U_m) -2.0) 1.0) (- U_m)))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -2e-118) {
tmp = sqrt(fma((U_m * ((U_m / J) / J)), 0.25, 1.0)) * (-2.0 * J);
} else if (t_1 <= -5e-289) {
tmp = (((J * J) / U_m) * -2.0) - U_m;
} else {
tmp = (((((J / U_m) * J) / U_m) * -2.0) - 1.0) * -U_m;
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -2e-118) tmp = Float64(sqrt(fma(Float64(U_m * Float64(Float64(U_m / J) / J)), 0.25, 1.0)) * Float64(-2.0 * J)); elseif (t_1 <= -5e-289) tmp = Float64(Float64(Float64(Float64(J * J) / U_m) * -2.0) - U_m); else tmp = Float64(Float64(Float64(Float64(Float64(Float64(J / U_m) * J) / U_m) * -2.0) - 1.0) * Float64(-U_m)); end return tmp end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -2e-118], N[(N[Sqrt[N[(N[(U$95$m * N[(N[(U$95$m / J), $MachinePrecision] / J), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision] * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -5e-289], N[(N[(N[(N[(J * J), $MachinePrecision] / U$95$m), $MachinePrecision] * -2.0), $MachinePrecision] - U$95$m), $MachinePrecision], N[(N[(N[(N[(N[(N[(J / U$95$m), $MachinePrecision] * J), $MachinePrecision] / U$95$m), $MachinePrecision] * -2.0), $MachinePrecision] - 1.0), $MachinePrecision] * (-U$95$m)), $MachinePrecision]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-118}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(U\_m \cdot \frac{\frac{U\_m}{J}}{J}, 0.25, 1\right)} \cdot \left(-2 \cdot J\right)\\
\mathbf{elif}\;t\_1 \leq -5 \cdot 10^{-289}:\\
\;\;\;\;\frac{J \cdot J}{U\_m} \cdot -2 - U\_m\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\frac{J}{U\_m} \cdot J}{U\_m} \cdot -2 - 1\right) \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -inf.0Initial program 6.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6449.0
Applied rewrites49.0%
if -inf.0 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1.99999999999999997e-118Initial program 99.7%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.7
Applied rewrites99.7%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6453.4
Applied rewrites53.4%
Applied rewrites59.6%
if -1.99999999999999997e-118 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -5.00000000000000029e-289Initial program 99.9%
Taylor expanded in J around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6418.2
Applied rewrites18.2%
Taylor expanded in K around 0
Applied rewrites18.1%
if -5.00000000000000029e-289 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 68.4%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6468.4
Applied rewrites68.4%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6431.6
Applied rewrites31.6%
Taylor expanded in U around -inf
Applied rewrites31.7%
Applied rewrites31.7%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))))
(if (<= t_1 -5e+285)
(- U_m)
(if (<= t_1 -2e-118)
(* (sqrt (fma (* U_m (/ U_m (* J J))) 0.25 1.0)) (* -2.0 J))
(if (<= t_1 -5e-289)
(- (* (/ (* J J) U_m) -2.0) U_m)
(* (- (* (/ (* (/ J U_m) J) U_m) -2.0) 1.0) (- U_m)))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_1 <= -5e+285) {
tmp = -U_m;
} else if (t_1 <= -2e-118) {
tmp = sqrt(fma((U_m * (U_m / (J * J))), 0.25, 1.0)) * (-2.0 * J);
} else if (t_1 <= -5e-289) {
tmp = (((J * J) / U_m) * -2.0) - U_m;
} else {
tmp = (((((J / U_m) * J) / U_m) * -2.0) - 1.0) * -U_m;
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= -5e+285) tmp = Float64(-U_m); elseif (t_1 <= -2e-118) tmp = Float64(sqrt(fma(Float64(U_m * Float64(U_m / Float64(J * J))), 0.25, 1.0)) * Float64(-2.0 * J)); elseif (t_1 <= -5e-289) tmp = Float64(Float64(Float64(Float64(J * J) / U_m) * -2.0) - U_m); else tmp = Float64(Float64(Float64(Float64(Float64(Float64(J / U_m) * J) / U_m) * -2.0) - 1.0) * Float64(-U_m)); end return tmp end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -5e+285], (-U$95$m), If[LessEqual[t$95$1, -2e-118], N[(N[Sqrt[N[(N[(U$95$m * N[(U$95$m / N[(J * J), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision] * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -5e-289], N[(N[(N[(N[(J * J), $MachinePrecision] / U$95$m), $MachinePrecision] * -2.0), $MachinePrecision] - U$95$m), $MachinePrecision], N[(N[(N[(N[(N[(N[(J / U$95$m), $MachinePrecision] * J), $MachinePrecision] / U$95$m), $MachinePrecision] * -2.0), $MachinePrecision] - 1.0), $MachinePrecision] * (-U$95$m)), $MachinePrecision]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+285}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-118}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(U\_m \cdot \frac{U\_m}{J \cdot J}, 0.25, 1\right)} \cdot \left(-2 \cdot J\right)\\
\mathbf{elif}\;t\_1 \leq -5 \cdot 10^{-289}:\\
\;\;\;\;\frac{J \cdot J}{U\_m} \cdot -2 - U\_m\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\frac{J}{U\_m} \cdot J}{U\_m} \cdot -2 - 1\right) \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -5.00000000000000016e285Initial program 15.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6446.9
Applied rewrites46.9%
if -5.00000000000000016e285 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1.99999999999999997e-118Initial program 99.7%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.7
Applied rewrites99.7%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6455.5
Applied rewrites55.5%
if -1.99999999999999997e-118 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -5.00000000000000029e-289Initial program 99.9%
Taylor expanded in J around 0
*-commutativeN/A
associate-/l*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-pow.f64N/A
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f6418.2
Applied rewrites18.2%
Taylor expanded in K around 0
Applied rewrites18.1%
if -5.00000000000000029e-289 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 68.4%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6468.4
Applied rewrites68.4%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6431.6
Applied rewrites31.6%
Taylor expanded in U around -inf
Applied rewrites31.7%
Applied rewrites31.7%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))))
(if (<= t_1 -5e+285)
(- U_m)
(if (<= t_1 -2e-102)
(fma (* U_m (/ U_m J)) -0.25 (* -2.0 J))
(if (<= t_1 -5e-289)
(- U_m)
(* (- (* (/ (* (/ J U_m) J) U_m) -2.0) 1.0) (- U_m)))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_1 <= -5e+285) {
tmp = -U_m;
} else if (t_1 <= -2e-102) {
tmp = fma((U_m * (U_m / J)), -0.25, (-2.0 * J));
} else if (t_1 <= -5e-289) {
tmp = -U_m;
} else {
tmp = (((((J / U_m) * J) / U_m) * -2.0) - 1.0) * -U_m;
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= -5e+285) tmp = Float64(-U_m); elseif (t_1 <= -2e-102) tmp = fma(Float64(U_m * Float64(U_m / J)), -0.25, Float64(-2.0 * J)); elseif (t_1 <= -5e-289) tmp = Float64(-U_m); else tmp = Float64(Float64(Float64(Float64(Float64(Float64(J / U_m) * J) / U_m) * -2.0) - 1.0) * Float64(-U_m)); end return tmp end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -5e+285], (-U$95$m), If[LessEqual[t$95$1, -2e-102], N[(N[(U$95$m * N[(U$95$m / J), $MachinePrecision]), $MachinePrecision] * -0.25 + N[(-2.0 * J), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -5e-289], (-U$95$m), N[(N[(N[(N[(N[(N[(J / U$95$m), $MachinePrecision] * J), $MachinePrecision] / U$95$m), $MachinePrecision] * -2.0), $MachinePrecision] - 1.0), $MachinePrecision] * (-U$95$m)), $MachinePrecision]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+285}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-102}:\\
\;\;\;\;\mathsf{fma}\left(U\_m \cdot \frac{U\_m}{J}, -0.25, -2 \cdot J\right)\\
\mathbf{elif}\;t\_1 \leq -5 \cdot 10^{-289}:\\
\;\;\;\;-U\_m\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\frac{J}{U\_m} \cdot J}{U\_m} \cdot -2 - 1\right) \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -5.00000000000000016e285 or -1.99999999999999987e-102 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -5.00000000000000029e-289Initial program 41.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6437.3
Applied rewrites37.3%
if -5.00000000000000016e285 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1.99999999999999987e-102Initial program 99.7%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.7
Applied rewrites99.7%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6456.2
Applied rewrites56.2%
Taylor expanded in U around 0
Applied rewrites41.3%
if -5.00000000000000029e-289 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 68.4%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6468.4
Applied rewrites68.4%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6431.6
Applied rewrites31.6%
Taylor expanded in U around -inf
Applied rewrites31.7%
Applied rewrites31.7%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))))
(if (<= t_1 -5e+285)
(- U_m)
(if (<= t_1 -2e-102)
(fma (* U_m (/ U_m J)) -0.25 (* -2.0 J))
(if (<= t_1 -5e-289) (- U_m) (* -1.0 (- U_m)))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_1 <= -5e+285) {
tmp = -U_m;
} else if (t_1 <= -2e-102) {
tmp = fma((U_m * (U_m / J)), -0.25, (-2.0 * J));
} else if (t_1 <= -5e-289) {
tmp = -U_m;
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= -5e+285) tmp = Float64(-U_m); elseif (t_1 <= -2e-102) tmp = fma(Float64(U_m * Float64(U_m / J)), -0.25, Float64(-2.0 * J)); elseif (t_1 <= -5e-289) tmp = Float64(-U_m); else tmp = Float64(-1.0 * Float64(-U_m)); end return tmp end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -5e+285], (-U$95$m), If[LessEqual[t$95$1, -2e-102], N[(N[(U$95$m * N[(U$95$m / J), $MachinePrecision]), $MachinePrecision] * -0.25 + N[(-2.0 * J), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -5e-289], (-U$95$m), N[(-1.0 * (-U$95$m)), $MachinePrecision]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+285}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-102}:\\
\;\;\;\;\mathsf{fma}\left(U\_m \cdot \frac{U\_m}{J}, -0.25, -2 \cdot J\right)\\
\mathbf{elif}\;t\_1 \leq -5 \cdot 10^{-289}:\\
\;\;\;\;-U\_m\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -5.00000000000000016e285 or -1.99999999999999987e-102 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -5.00000000000000029e-289Initial program 41.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6437.3
Applied rewrites37.3%
if -5.00000000000000016e285 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1.99999999999999987e-102Initial program 99.7%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.7
Applied rewrites99.7%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6456.2
Applied rewrites56.2%
Taylor expanded in U around 0
Applied rewrites41.3%
if -5.00000000000000029e-289 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 68.4%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6468.4
Applied rewrites68.4%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6431.6
Applied rewrites31.6%
Taylor expanded in U around -inf
Applied rewrites31.7%
Taylor expanded in J around 0
Applied rewrites31.9%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))))
(if (<= t_1 -5e+285)
(- U_m)
(if (<= t_1 -2e-102)
(* 1.0 (* -2.0 J))
(if (<= t_1 -5e-289) (- U_m) (* -1.0 (- U_m)))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_1 <= -5e+285) {
tmp = -U_m;
} else if (t_1 <= -2e-102) {
tmp = 1.0 * (-2.0 * J);
} else if (t_1 <= -5e-289) {
tmp = -U_m;
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = abs(u)
real(8) function code(j, k, u_m)
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8), intent (in) :: u_m
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = cos((k / 2.0d0))
t_1 = (((-2.0d0) * j) * t_0) * sqrt((1.0d0 + ((u_m / ((2.0d0 * j) * t_0)) ** 2.0d0)))
if (t_1 <= (-5d+285)) then
tmp = -u_m
else if (t_1 <= (-2d-102)) then
tmp = 1.0d0 * ((-2.0d0) * j)
else if (t_1 <= (-5d-289)) then
tmp = -u_m
else
tmp = (-1.0d0) * -u_m
end if
code = tmp
end function
U_m = Math.abs(U);
public static double code(double J, double K, double U_m) {
double t_0 = Math.cos((K / 2.0));
double t_1 = ((-2.0 * J) * t_0) * Math.sqrt((1.0 + Math.pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_1 <= -5e+285) {
tmp = -U_m;
} else if (t_1 <= -2e-102) {
tmp = 1.0 * (-2.0 * J);
} else if (t_1 <= -5e-289) {
tmp = -U_m;
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = math.fabs(U) def code(J, K, U_m): t_0 = math.cos((K / 2.0)) t_1 = ((-2.0 * J) * t_0) * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J) * t_0)), 2.0))) tmp = 0 if t_1 <= -5e+285: tmp = -U_m elif t_1 <= -2e-102: tmp = 1.0 * (-2.0 * J) elif t_1 <= -5e-289: tmp = -U_m else: tmp = -1.0 * -U_m return tmp
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= -5e+285) tmp = Float64(-U_m); elseif (t_1 <= -2e-102) tmp = Float64(1.0 * Float64(-2.0 * J)); elseif (t_1 <= -5e-289) tmp = Float64(-U_m); else tmp = Float64(-1.0 * Float64(-U_m)); end return tmp end
U_m = abs(U); function tmp_2 = code(J, K, U_m) t_0 = cos((K / 2.0)); t_1 = ((-2.0 * J) * t_0) * sqrt((1.0 + ((U_m / ((2.0 * J) * t_0)) ^ 2.0))); tmp = 0.0; if (t_1 <= -5e+285) tmp = -U_m; elseif (t_1 <= -2e-102) tmp = 1.0 * (-2.0 * J); elseif (t_1 <= -5e-289) tmp = -U_m; else tmp = -1.0 * -U_m; end tmp_2 = tmp; end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -5e+285], (-U$95$m), If[LessEqual[t$95$1, -2e-102], N[(1.0 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -5e-289], (-U$95$m), N[(-1.0 * (-U$95$m)), $MachinePrecision]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+285}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-102}:\\
\;\;\;\;1 \cdot \left(-2 \cdot J\right)\\
\mathbf{elif}\;t\_1 \leq -5 \cdot 10^{-289}:\\
\;\;\;\;-U\_m\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -5.00000000000000016e285 or -1.99999999999999987e-102 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -5.00000000000000029e-289Initial program 41.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6437.3
Applied rewrites37.3%
if -5.00000000000000016e285 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1.99999999999999987e-102Initial program 99.7%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.7
Applied rewrites99.7%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6456.2
Applied rewrites56.2%
Taylor expanded in J around inf
Applied rewrites41.2%
if -5.00000000000000029e-289 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 68.4%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6468.4
Applied rewrites68.4%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6431.6
Applied rewrites31.6%
Taylor expanded in U around -inf
Applied rewrites31.7%
Taylor expanded in J around 0
Applied rewrites31.9%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1 (* (* -2.0 J) t_0))
(t_2 (* t_1 (sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))))
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 5e+292)
(*
t_1
(sqrt
(+ 1.0 (pow (/ U_m (* (* 2.0 J) (fma -0.125 (* K K) 1.0))) 2.0))))
(* -1.0 (- U_m))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = (-2.0 * J) * t_0;
double t_2 = t_1 * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= 5e+292) {
tmp = t_1 * sqrt((1.0 + pow((U_m / ((2.0 * J) * fma(-0.125, (K * K), 1.0))), 2.0)));
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(-2.0 * J) * t_0) t_2 = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= 5e+292) tmp = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * fma(-0.125, Float64(K * K), 1.0))) ^ 2.0)))); else tmp = Float64(-1.0 * Float64(-U_m)); end return tmp end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, 5e+292], N[(t$95$1 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * N[(-0.125 * N[(K * K), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(-1.0 * (-U$95$m)), $MachinePrecision]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(-2 \cdot J\right) \cdot t\_0\\
t_2 := t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 5 \cdot 10^{+292}:\\
\;\;\;\;t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot \mathsf{fma}\left(-0.125, K \cdot K, 1\right)}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -inf.0Initial program 6.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6449.0
Applied rewrites49.0%
if -inf.0 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < 4.9999999999999996e292Initial program 99.8%
Taylor expanded in K around 0
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f6480.5
Applied rewrites80.5%
if 4.9999999999999996e292 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 5.5%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f645.5
Applied rewrites5.5%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f644.1
Applied rewrites4.1%
Taylor expanded in U around -inf
Applied rewrites54.1%
Taylor expanded in J around 0
Applied rewrites54.1%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1 (* (* -2.0 J) t_0))
(t_2 (* t_1 (sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))))
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 5e+292)
(* t_1 (sqrt (fma (* U_m (/ (/ U_m J) J)) 0.25 1.0)))
(* -1.0 (- U_m))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = (-2.0 * J) * t_0;
double t_2 = t_1 * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= 5e+292) {
tmp = t_1 * sqrt(fma((U_m * ((U_m / J) / J)), 0.25, 1.0));
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(-2.0 * J) * t_0) t_2 = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= 5e+292) tmp = Float64(t_1 * sqrt(fma(Float64(U_m * Float64(Float64(U_m / J) / J)), 0.25, 1.0))); else tmp = Float64(-1.0 * Float64(-U_m)); end return tmp end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, 5e+292], N[(t$95$1 * N[Sqrt[N[(N[(U$95$m * N[(N[(U$95$m / J), $MachinePrecision] / J), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(-1.0 * (-U$95$m)), $MachinePrecision]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(-2 \cdot J\right) \cdot t\_0\\
t_2 := t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 5 \cdot 10^{+292}:\\
\;\;\;\;t\_1 \cdot \sqrt{\mathsf{fma}\left(U\_m \cdot \frac{\frac{U\_m}{J}}{J}, 0.25, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -inf.0Initial program 6.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6449.0
Applied rewrites49.0%
if -inf.0 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < 4.9999999999999996e292Initial program 99.8%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6499.8
Applied rewrites99.8%
Taylor expanded in K around 0
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6471.5
Applied rewrites71.5%
Applied rewrites78.6%
if 4.9999999999999996e292 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 5.5%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f645.5
Applied rewrites5.5%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f644.1
Applied rewrites4.1%
Taylor expanded in U around -inf
Applied rewrites54.1%
Taylor expanded in J around 0
Applied rewrites54.1%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0))))
(if (<=
(*
(* (* -2.0 J) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J) t_0)) 2.0))))
-5e-289)
(- U_m)
(* -1.0 (- U_m)))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = cos((K / 2.0));
double tmp;
if ((((-2.0 * J) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J) * t_0)), 2.0)))) <= -5e-289) {
tmp = -U_m;
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = abs(u)
real(8) function code(j, k, u_m)
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8), intent (in) :: u_m
real(8) :: t_0
real(8) :: tmp
t_0 = cos((k / 2.0d0))
if (((((-2.0d0) * j) * t_0) * sqrt((1.0d0 + ((u_m / ((2.0d0 * j) * t_0)) ** 2.0d0)))) <= (-5d-289)) then
tmp = -u_m
else
tmp = (-1.0d0) * -u_m
end if
code = tmp
end function
U_m = Math.abs(U);
public static double code(double J, double K, double U_m) {
double t_0 = Math.cos((K / 2.0));
double tmp;
if ((((-2.0 * J) * t_0) * Math.sqrt((1.0 + Math.pow((U_m / ((2.0 * J) * t_0)), 2.0)))) <= -5e-289) {
tmp = -U_m;
} else {
tmp = -1.0 * -U_m;
}
return tmp;
}
U_m = math.fabs(U) def code(J, K, U_m): t_0 = math.cos((K / 2.0)) tmp = 0 if (((-2.0 * J) * t_0) * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J) * t_0)), 2.0)))) <= -5e-289: tmp = -U_m else: tmp = -1.0 * -U_m return tmp
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) tmp = 0.0 if (Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) <= -5e-289) tmp = Float64(-U_m); else tmp = Float64(-1.0 * Float64(-U_m)); end return tmp end
U_m = abs(U); function tmp_2 = code(J, K, U_m) t_0 = cos((K / 2.0)); tmp = 0.0; if ((((-2.0 * J) * t_0) * sqrt((1.0 + ((U_m / ((2.0 * J) * t_0)) ^ 2.0)))) <= -5e-289) tmp = -U_m; else tmp = -1.0 * -U_m; end tmp_2 = tmp; end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], -5e-289], (-U$95$m), N[(-1.0 * (-U$95$m)), $MachinePrecision]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
\mathbf{if}\;\left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}} \leq -5 \cdot 10^{-289}:\\
\;\;\;\;-U\_m\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot \left(-U\_m\right)\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -5.00000000000000029e-289Initial program 74.8%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6423.6
Applied rewrites23.6%
if -5.00000000000000029e-289 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 68.4%
Taylor expanded in K around inf
cos-neg-revN/A
lower-cos.f64N/A
distribute-lft-neg-inN/A
metadata-evalN/A
lower-*.f6468.4
Applied rewrites68.4%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6431.6
Applied rewrites31.6%
Taylor expanded in U around -inf
Applied rewrites31.7%
Taylor expanded in J around 0
Applied rewrites31.9%
U_m = (fabs.f64 U) (FPCore (J K U_m) :precision binary64 (- U_m))
U_m = fabs(U);
double code(double J, double K, double U_m) {
return -U_m;
}
U_m = abs(u)
real(8) function code(j, k, u_m)
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8), intent (in) :: u_m
code = -u_m
end function
U_m = Math.abs(U);
public static double code(double J, double K, double U_m) {
return -U_m;
}
U_m = math.fabs(U) def code(J, K, U_m): return -U_m
U_m = abs(U) function code(J, K, U_m) return Float64(-U_m) end
U_m = abs(U); function tmp = code(J, K, U_m) tmp = -U_m; end
U_m = N[Abs[U], $MachinePrecision] code[J_, K_, U$95$m_] := (-U$95$m)
\begin{array}{l}
U_m = \left|U\right|
\\
-U\_m
\end{array}
Initial program 71.9%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6426.9
Applied rewrites26.9%
herbie shell --seed 2024340
(FPCore (J K U)
:name "Maksimov and Kolovsky, Equation (3)"
:precision binary64
(* (* (* -2.0 J) (cos (/ K 2.0))) (sqrt (+ 1.0 (pow (/ U (* (* 2.0 J) (cos (/ K 2.0)))) 2.0)))))