
(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 14 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 (sqrt (+ 1.0 (pow (/ U_m (* t_0 (* J 2.0))) 2.0))))
(t_2 (* (* t_0 (* -2.0 J)) t_1)))
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 1e+304)
(* t_1 (* J (* -2.0 (cos (* K 0.5)))))
(fma 2.0 (/ (* J J) U_m) 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 = sqrt((1.0 + pow((U_m / (t_0 * (J * 2.0))), 2.0)));
double t_2 = (t_0 * (-2.0 * J)) * t_1;
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= 1e+304) {
tmp = t_1 * (J * (-2.0 * cos((K * 0.5))));
} else {
tmp = fma(2.0, ((J * J) / U_m), U_m);
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = sqrt(Float64(1.0 + (Float64(U_m / Float64(t_0 * Float64(J * 2.0))) ^ 2.0))) t_2 = Float64(Float64(t_0 * Float64(-2.0 * J)) * t_1) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= 1e+304) tmp = Float64(t_1 * Float64(J * Float64(-2.0 * cos(Float64(K * 0.5))))); else tmp = fma(2.0, Float64(Float64(J * J) / U_m), 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[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$0 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(t$95$0 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision]}, If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, 1e+304], N[(t$95$1 * N[(J * N[(-2.0 * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(J * J), $MachinePrecision] / U$95$m), $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 := \sqrt{1 + {\left(\frac{U\_m}{t\_0 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
t_2 := \left(t\_0 \cdot \left(-2 \cdot J\right)\right) \cdot t\_1\\
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 10^{+304}:\\
\;\;\;\;t\_1 \cdot \left(J \cdot \left(-2 \cdot \cos \left(K \cdot 0.5\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(2, \frac{J \cdot J}{U\_m}, 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.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6444.4
Simplified44.4%
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))))) < 9.9999999999999994e303Initial program 99.8%
lift-/.f64N/A
lift-cos.f64N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6499.8
lift-/.f64N/A
div-invN/A
metadata-evalN/A
lower-*.f6499.8
Applied egg-rr99.8%
if 9.9999999999999994e303 < (*.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.4%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Taylor expanded in K around 0
Simplified71.2%
Taylor expanded in J around 0
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Final simplification89.0%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0
(*
(* J (* -2.0 (cos (* K 0.5))))
(sqrt
(fma (* U_m U_m) (/ 0.25 (* (* J J) (fma 0.5 (cos K) 0.5))) 1.0))))
(t_1 (cos (/ K 2.0)))
(t_2
(*
(* t_1 (* -2.0 J))
(sqrt (+ 1.0 (pow (/ U_m (* t_1 (* J 2.0))) 2.0)))))
(t_3 (/ U_m (* -2.0 J)))
(t_4 (* (* -2.0 J) (sqrt (fma t_3 t_3 1.0)))))
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 -5e+202)
t_4
(if (<= t_2 -2e-21)
t_0
(if (<= t_2 -1e-306)
t_4
(if (<= t_2 2e-94)
(fma
(* (/ J U_m) (* U_m (- (- -1.0) (* -2.0 (* 0.5 (cos K))))))
(/ J U_m)
U_m)
(if (<= t_2 2e+299)
t_0
(fma U_m (/ (* (* J 2.0) (/ J U_m)) U_m) U_m)))))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = (J * (-2.0 * cos((K * 0.5)))) * sqrt(fma((U_m * U_m), (0.25 / ((J * J) * fma(0.5, cos(K), 0.5))), 1.0));
double t_1 = cos((K / 2.0));
double t_2 = (t_1 * (-2.0 * J)) * sqrt((1.0 + pow((U_m / (t_1 * (J * 2.0))), 2.0)));
double t_3 = U_m / (-2.0 * J);
double t_4 = (-2.0 * J) * sqrt(fma(t_3, t_3, 1.0));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= -5e+202) {
tmp = t_4;
} else if (t_2 <= -2e-21) {
tmp = t_0;
} else if (t_2 <= -1e-306) {
tmp = t_4;
} else if (t_2 <= 2e-94) {
tmp = fma(((J / U_m) * (U_m * (-(-1.0) - (-2.0 * (0.5 * cos(K)))))), (J / U_m), U_m);
} else if (t_2 <= 2e+299) {
tmp = t_0;
} else {
tmp = fma(U_m, (((J * 2.0) * (J / U_m)) / U_m), U_m);
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = Float64(Float64(J * Float64(-2.0 * cos(Float64(K * 0.5)))) * sqrt(fma(Float64(U_m * U_m), Float64(0.25 / Float64(Float64(J * J) * fma(0.5, cos(K), 0.5))), 1.0))) t_1 = cos(Float64(K / 2.0)) t_2 = Float64(Float64(t_1 * Float64(-2.0 * J)) * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_1 * Float64(J * 2.0))) ^ 2.0)))) t_3 = Float64(U_m / Float64(-2.0 * J)) t_4 = Float64(Float64(-2.0 * J) * sqrt(fma(t_3, t_3, 1.0))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= -5e+202) tmp = t_4; elseif (t_2 <= -2e-21) tmp = t_0; elseif (t_2 <= -1e-306) tmp = t_4; elseif (t_2 <= 2e-94) tmp = fma(Float64(Float64(J / U_m) * Float64(U_m * Float64(Float64(-(-1.0)) - Float64(-2.0 * Float64(0.5 * cos(K)))))), Float64(J / U_m), U_m); elseif (t_2 <= 2e+299) tmp = t_0; else tmp = fma(U_m, Float64(Float64(Float64(J * 2.0) * Float64(J / U_m)) / U_m), U_m); end return tmp end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[(N[(J * N[(-2.0 * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[(U$95$m * U$95$m), $MachinePrecision] * N[(0.25 / N[(N[(J * J), $MachinePrecision] * N[(0.5 * N[Cos[K], $MachinePrecision] + 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(t$95$1 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$1 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(U$95$m / N[(-2.0 * J), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(N[(-2.0 * J), $MachinePrecision] * N[Sqrt[N[(t$95$3 * t$95$3 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, -5e+202], t$95$4, If[LessEqual[t$95$2, -2e-21], t$95$0, If[LessEqual[t$95$2, -1e-306], t$95$4, If[LessEqual[t$95$2, 2e-94], N[(N[(N[(J / U$95$m), $MachinePrecision] * N[(U$95$m * N[((--1.0) - N[(-2.0 * N[(0.5 * N[Cos[K], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(J / U$95$m), $MachinePrecision] + U$95$m), $MachinePrecision], If[LessEqual[t$95$2, 2e+299], t$95$0, N[(U$95$m * N[(N[(N[(J * 2.0), $MachinePrecision] * N[(J / U$95$m), $MachinePrecision]), $MachinePrecision] / U$95$m), $MachinePrecision] + U$95$m), $MachinePrecision]]]]]]]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \left(J \cdot \left(-2 \cdot \cos \left(K \cdot 0.5\right)\right)\right) \cdot \sqrt{\mathsf{fma}\left(U\_m \cdot U\_m, \frac{0.25}{\left(J \cdot J\right) \cdot \mathsf{fma}\left(0.5, \cos K, 0.5\right)}, 1\right)}\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := \left(t\_1 \cdot \left(-2 \cdot J\right)\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_1 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
t_3 := \frac{U\_m}{-2 \cdot J}\\
t_4 := \left(-2 \cdot J\right) \cdot \sqrt{\mathsf{fma}\left(t\_3, t\_3, 1\right)}\\
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq -5 \cdot 10^{+202}:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;t\_2 \leq -2 \cdot 10^{-21}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_2 \leq -1 \cdot 10^{-306}:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{-94}:\\
\;\;\;\;\mathsf{fma}\left(\frac{J}{U\_m} \cdot \left(U\_m \cdot \left(\left(--1\right) - -2 \cdot \left(0.5 \cdot \cos K\right)\right)\right), \frac{J}{U\_m}, U\_m\right)\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+299}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(U\_m, \frac{\left(J \cdot 2\right) \cdot \frac{J}{U\_m}}{U\_m}, 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.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6444.4
Simplified44.4%
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.9999999999999999e202 or -1.99999999999999982e-21 < (*.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.00000000000000003e-306Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f6454.1
Simplified54.1%
Taylor expanded in K around 0
lower-*.f6461.5
Simplified61.5%
lift-*.f64N/A
lift-/.f64N/A
lift-pow.f64N/A
+-commutativeN/A
lift-pow.f64N/A
unpow2N/A
lift-/.f64N/A
lift-/.f64N/A
frac-2negN/A
distribute-frac-negN/A
frac-2negN/A
distribute-frac-negN/A
sqr-negN/A
lower-fma.f64N/A
Applied egg-rr61.5%
if -4.9999999999999999e202 < (*.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.99999999999999982e-21 or 1.9999999999999999e-94 < (*.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))))) < 2.0000000000000001e299Initial program 99.8%
lift-/.f64N/A
lift-cos.f64N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6499.8
lift-/.f64N/A
div-invN/A
metadata-evalN/A
lower-*.f6499.8
Applied egg-rr99.8%
Applied egg-rr94.3%
Taylor expanded in K around inf
lower-sqrt.f64N/A
+-commutativeN/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
lower-/.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f6486.2
Simplified86.2%
if -1.00000000000000003e-306 < (*.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.9999999999999999e-94Initial program 99.9%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f646.2
Simplified6.2%
Applied egg-rr13.8%
if 2.0000000000000001e299 < (*.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 11.1%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6467.2
Simplified67.2%
Taylor expanded in K around 0
Simplified67.2%
Taylor expanded in U around inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6467.4
Simplified67.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f6467.6
Applied egg-rr67.6%
Final simplification67.9%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0
(*
(* J (* -2.0 (cos (* K 0.5))))
(sqrt
(fma (/ (* U_m U_m) (* (* J J) (fma 0.5 (cos K) 0.5))) 0.25 1.0))))
(t_1 (cos (/ K 2.0)))
(t_2
(*
(* t_1 (* -2.0 J))
(sqrt (+ 1.0 (pow (/ U_m (* t_1 (* J 2.0))) 2.0)))))
(t_3 (/ U_m (* -2.0 J)))
(t_4 (* (* -2.0 J) (sqrt (fma t_3 t_3 1.0)))))
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 -5e+202)
t_4
(if (<= t_2 -5e-39)
t_0
(if (<= t_2 2e-185)
t_4
(if (<= t_2 2e+299)
t_0
(fma U_m (/ (* (* J 2.0) (/ J U_m)) U_m) U_m))))))))U_m = fabs(U);
double code(double J, double K, double U_m) {
double t_0 = (J * (-2.0 * cos((K * 0.5)))) * sqrt(fma(((U_m * U_m) / ((J * J) * fma(0.5, cos(K), 0.5))), 0.25, 1.0));
double t_1 = cos((K / 2.0));
double t_2 = (t_1 * (-2.0 * J)) * sqrt((1.0 + pow((U_m / (t_1 * (J * 2.0))), 2.0)));
double t_3 = U_m / (-2.0 * J);
double t_4 = (-2.0 * J) * sqrt(fma(t_3, t_3, 1.0));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= -5e+202) {
tmp = t_4;
} else if (t_2 <= -5e-39) {
tmp = t_0;
} else if (t_2 <= 2e-185) {
tmp = t_4;
} else if (t_2 <= 2e+299) {
tmp = t_0;
} else {
tmp = fma(U_m, (((J * 2.0) * (J / U_m)) / U_m), U_m);
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = Float64(Float64(J * Float64(-2.0 * cos(Float64(K * 0.5)))) * sqrt(fma(Float64(Float64(U_m * U_m) / Float64(Float64(J * J) * fma(0.5, cos(K), 0.5))), 0.25, 1.0))) t_1 = cos(Float64(K / 2.0)) t_2 = Float64(Float64(t_1 * Float64(-2.0 * J)) * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_1 * Float64(J * 2.0))) ^ 2.0)))) t_3 = Float64(U_m / Float64(-2.0 * J)) t_4 = Float64(Float64(-2.0 * J) * sqrt(fma(t_3, t_3, 1.0))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= -5e+202) tmp = t_4; elseif (t_2 <= -5e-39) tmp = t_0; elseif (t_2 <= 2e-185) tmp = t_4; elseif (t_2 <= 2e+299) tmp = t_0; else tmp = fma(U_m, Float64(Float64(Float64(J * 2.0) * Float64(J / U_m)) / U_m), U_m); end return tmp end
U_m = N[Abs[U], $MachinePrecision]
code[J_, K_, U$95$m_] := Block[{t$95$0 = N[(N[(J * N[(-2.0 * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[(N[(U$95$m * U$95$m), $MachinePrecision] / N[(N[(J * J), $MachinePrecision] * N[(0.5 * N[Cos[K], $MachinePrecision] + 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(t$95$1 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$1 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(U$95$m / N[(-2.0 * J), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(N[(-2.0 * J), $MachinePrecision] * N[Sqrt[N[(t$95$3 * t$95$3 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, -5e+202], t$95$4, If[LessEqual[t$95$2, -5e-39], t$95$0, If[LessEqual[t$95$2, 2e-185], t$95$4, If[LessEqual[t$95$2, 2e+299], t$95$0, N[(U$95$m * N[(N[(N[(J * 2.0), $MachinePrecision] * N[(J / U$95$m), $MachinePrecision]), $MachinePrecision] / U$95$m), $MachinePrecision] + U$95$m), $MachinePrecision]]]]]]]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
t_0 := \left(J \cdot \left(-2 \cdot \cos \left(K \cdot 0.5\right)\right)\right) \cdot \sqrt{\mathsf{fma}\left(\frac{U\_m \cdot U\_m}{\left(J \cdot J\right) \cdot \mathsf{fma}\left(0.5, \cos K, 0.5\right)}, 0.25, 1\right)}\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := \left(t\_1 \cdot \left(-2 \cdot J\right)\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_1 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
t_3 := \frac{U\_m}{-2 \cdot J}\\
t_4 := \left(-2 \cdot J\right) \cdot \sqrt{\mathsf{fma}\left(t\_3, t\_3, 1\right)}\\
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq -5 \cdot 10^{+202}:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;t\_2 \leq -5 \cdot 10^{-39}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{-185}:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+299}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(U\_m, \frac{\left(J \cdot 2\right) \cdot \frac{J}{U\_m}}{U\_m}, 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.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6444.4
Simplified44.4%
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.9999999999999999e202 or -4.9999999999999998e-39 < (*.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))))) < 2e-185Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f6455.1
Simplified55.1%
Taylor expanded in K around 0
lower-*.f6462.4
Simplified62.4%
lift-*.f64N/A
lift-/.f64N/A
lift-pow.f64N/A
+-commutativeN/A
lift-pow.f64N/A
unpow2N/A
lift-/.f64N/A
lift-/.f64N/A
frac-2negN/A
distribute-frac-negN/A
frac-2negN/A
distribute-frac-negN/A
sqr-negN/A
lower-fma.f64N/A
Applied egg-rr62.4%
if -4.9999999999999999e202 < (*.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.9999999999999998e-39 or 2e-185 < (*.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))))) < 2.0000000000000001e299Initial program 99.8%
lift-/.f64N/A
lift-cos.f64N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6499.8
lift-/.f64N/A
div-invN/A
metadata-evalN/A
lower-*.f6499.8
Applied egg-rr99.8%
Applied egg-rr95.1%
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
times-fracN/A
lift-*.f64N/A
associate-/r*N/A
associate-*r/N/A
lower-/.f64N/A
Applied egg-rr95.1%
Taylor expanded in K around inf
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f6483.7
Simplified83.7%
if 2.0000000000000001e299 < (*.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 11.1%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6467.2
Simplified67.2%
Taylor expanded in K around 0
Simplified67.2%
Taylor expanded in U around inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6467.4
Simplified67.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f6467.6
Applied egg-rr67.6%
Final simplification72.5%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1 (* t_0 (* -2.0 J)))
(t_2 (* t_1 (sqrt (+ 1.0 (pow (/ U_m (* t_0 (* J 2.0))) 2.0)))))
(t_3 (/ U_m (* -2.0 J))))
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 -1e-306)
(* (* -2.0 J) (sqrt (fma t_3 t_3 1.0)))
(if (<= t_2 2e-94)
(fma
(* (/ J U_m) (* U_m (- (- -1.0) (* -2.0 (* 0.5 (cos K))))))
(/ J U_m)
U_m)
(if (<= t_2 2e+299)
(* t_1 (sqrt (fma 0.25 (/ (* U_m U_m) (* J J)) 1.0)))
(fma U_m (/ (* (* J 2.0) (/ J U_m)) U_m) 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 = t_0 * (-2.0 * J);
double t_2 = t_1 * sqrt((1.0 + pow((U_m / (t_0 * (J * 2.0))), 2.0)));
double t_3 = U_m / (-2.0 * J);
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= -1e-306) {
tmp = (-2.0 * J) * sqrt(fma(t_3, t_3, 1.0));
} else if (t_2 <= 2e-94) {
tmp = fma(((J / U_m) * (U_m * (-(-1.0) - (-2.0 * (0.5 * cos(K)))))), (J / U_m), U_m);
} else if (t_2 <= 2e+299) {
tmp = t_1 * sqrt(fma(0.25, ((U_m * U_m) / (J * J)), 1.0));
} else {
tmp = fma(U_m, (((J * 2.0) * (J / U_m)) / U_m), U_m);
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(t_0 * Float64(-2.0 * J)) t_2 = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_0 * Float64(J * 2.0))) ^ 2.0)))) t_3 = Float64(U_m / Float64(-2.0 * J)) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= -1e-306) tmp = Float64(Float64(-2.0 * J) * sqrt(fma(t_3, t_3, 1.0))); elseif (t_2 <= 2e-94) tmp = fma(Float64(Float64(J / U_m) * Float64(U_m * Float64(Float64(-(-1.0)) - Float64(-2.0 * Float64(0.5 * cos(K)))))), Float64(J / U_m), U_m); elseif (t_2 <= 2e+299) tmp = Float64(t_1 * sqrt(fma(0.25, Float64(Float64(U_m * U_m) / Float64(J * J)), 1.0))); else tmp = fma(U_m, Float64(Float64(Float64(J * 2.0) * Float64(J / U_m)) / U_m), 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[(t$95$0 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$0 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(U$95$m / N[(-2.0 * J), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, -1e-306], N[(N[(-2.0 * J), $MachinePrecision] * N[Sqrt[N[(t$95$3 * t$95$3 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 2e-94], N[(N[(N[(J / U$95$m), $MachinePrecision] * N[(U$95$m * N[((--1.0) - N[(-2.0 * N[(0.5 * N[Cos[K], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(J / U$95$m), $MachinePrecision] + U$95$m), $MachinePrecision], If[LessEqual[t$95$2, 2e+299], N[(t$95$1 * N[Sqrt[N[(0.25 * N[(N[(U$95$m * U$95$m), $MachinePrecision] / N[(J * J), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(U$95$m * N[(N[(N[(J * 2.0), $MachinePrecision] * N[(J / U$95$m), $MachinePrecision]), $MachinePrecision] / U$95$m), $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 := t\_0 \cdot \left(-2 \cdot J\right)\\
t_2 := t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_0 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
t_3 := \frac{U\_m}{-2 \cdot J}\\
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq -1 \cdot 10^{-306}:\\
\;\;\;\;\left(-2 \cdot J\right) \cdot \sqrt{\mathsf{fma}\left(t\_3, t\_3, 1\right)}\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{-94}:\\
\;\;\;\;\mathsf{fma}\left(\frac{J}{U\_m} \cdot \left(U\_m \cdot \left(\left(--1\right) - -2 \cdot \left(0.5 \cdot \cos K\right)\right)\right), \frac{J}{U\_m}, U\_m\right)\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+299}:\\
\;\;\;\;t\_1 \cdot \sqrt{\mathsf{fma}\left(0.25, \frac{U\_m \cdot U\_m}{J \cdot J}, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(U\_m, \frac{\left(J \cdot 2\right) \cdot \frac{J}{U\_m}}{U\_m}, 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.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6444.4
Simplified44.4%
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.00000000000000003e-306Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f6446.9
Simplified46.9%
Taylor expanded in K around 0
lower-*.f6452.3
Simplified52.3%
lift-*.f64N/A
lift-/.f64N/A
lift-pow.f64N/A
+-commutativeN/A
lift-pow.f64N/A
unpow2N/A
lift-/.f64N/A
lift-/.f64N/A
frac-2negN/A
distribute-frac-negN/A
frac-2negN/A
distribute-frac-negN/A
sqr-negN/A
lower-fma.f64N/A
Applied egg-rr52.3%
if -1.00000000000000003e-306 < (*.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.9999999999999999e-94Initial program 99.9%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f646.2
Simplified6.2%
Applied egg-rr13.8%
if 1.9999999999999999e-94 < (*.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))))) < 2.0000000000000001e299Initial program 99.8%
Taylor expanded in K around 0
lower-sqrt.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6473.7
Simplified73.7%
if 2.0000000000000001e299 < (*.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 11.1%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6467.2
Simplified67.2%
Taylor expanded in K around 0
Simplified67.2%
Taylor expanded in U around inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
lower-*.f64N/A
unpow2N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6467.4
Simplified67.4%
lift-*.f64N/A
lift-/.f64N/A
associate-*r*N/A
lift-*.f64N/A
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower-/.f6467.6
Applied egg-rr67.6%
Final simplification56.6%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* t_0 (* -2.0 J))
(sqrt (+ 1.0 (pow (/ U_m (* t_0 (* J 2.0))) 2.0)))))
(t_2 (* -2.0 (cos (* K 0.5)))))
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -1e+135)
(*
J
(*
t_2
(sqrt
(fma
U_m
(/
U_m
(*
(+ 0.5 (* 0.5 (cos (* 2.0 (* K 0.5)))))
(* (* J 2.0) (* J 2.0))))
1.0))))
(if (<= t_1 1e+304)
(*
(* J t_2)
(sqrt
(+
1.0
(/
(* (/ U_m (* (* J 2.0) (fma 0.5 (cos K) 0.5))) (* U_m 0.5))
J))))
(fma 2.0 (/ (* J J) U_m) 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 = (t_0 * (-2.0 * J)) * sqrt((1.0 + pow((U_m / (t_0 * (J * 2.0))), 2.0)));
double t_2 = -2.0 * cos((K * 0.5));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -1e+135) {
tmp = J * (t_2 * sqrt(fma(U_m, (U_m / ((0.5 + (0.5 * cos((2.0 * (K * 0.5))))) * ((J * 2.0) * (J * 2.0)))), 1.0)));
} else if (t_1 <= 1e+304) {
tmp = (J * t_2) * sqrt((1.0 + (((U_m / ((J * 2.0) * fma(0.5, cos(K), 0.5))) * (U_m * 0.5)) / J)));
} else {
tmp = fma(2.0, ((J * J) / U_m), 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(t_0 * Float64(-2.0 * J)) * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_0 * Float64(J * 2.0))) ^ 2.0)))) t_2 = Float64(-2.0 * cos(Float64(K * 0.5))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -1e+135) tmp = Float64(J * Float64(t_2 * sqrt(fma(U_m, Float64(U_m / Float64(Float64(0.5 + Float64(0.5 * cos(Float64(2.0 * Float64(K * 0.5))))) * Float64(Float64(J * 2.0) * Float64(J * 2.0)))), 1.0)))); elseif (t_1 <= 1e+304) tmp = Float64(Float64(J * t_2) * sqrt(Float64(1.0 + Float64(Float64(Float64(U_m / Float64(Float64(J * 2.0) * fma(0.5, cos(K), 0.5))) * Float64(U_m * 0.5)) / J)))); else tmp = fma(2.0, Float64(Float64(J * J) / U_m), 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[(t$95$0 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$0 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(-2.0 * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -1e+135], N[(J * N[(t$95$2 * N[Sqrt[N[(U$95$m * N[(U$95$m / N[(N[(0.5 + N[(0.5 * N[Cos[N[(2.0 * N[(K * 0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(J * 2.0), $MachinePrecision] * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 1e+304], N[(N[(J * t$95$2), $MachinePrecision] * N[Sqrt[N[(1.0 + N[(N[(N[(U$95$m / N[(N[(J * 2.0), $MachinePrecision] * N[(0.5 * N[Cos[K], $MachinePrecision] + 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(U$95$m * 0.5), $MachinePrecision]), $MachinePrecision] / J), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(J * J), $MachinePrecision] / U$95$m), $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(t\_0 \cdot \left(-2 \cdot J\right)\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_0 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
t_2 := -2 \cdot \cos \left(K \cdot 0.5\right)\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{+135}:\\
\;\;\;\;J \cdot \left(t\_2 \cdot \sqrt{\mathsf{fma}\left(U\_m, \frac{U\_m}{\left(0.5 + 0.5 \cdot \cos \left(2 \cdot \left(K \cdot 0.5\right)\right)\right) \cdot \left(\left(J \cdot 2\right) \cdot \left(J \cdot 2\right)\right)}, 1\right)}\right)\\
\mathbf{elif}\;t\_1 \leq 10^{+304}:\\
\;\;\;\;\left(J \cdot t\_2\right) \cdot \sqrt{1 + \frac{\frac{U\_m}{\left(J \cdot 2\right) \cdot \mathsf{fma}\left(0.5, \cos K, 0.5\right)} \cdot \left(U\_m \cdot 0.5\right)}{J}}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(2, \frac{J \cdot J}{U\_m}, 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.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6444.4
Simplified44.4%
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))))) < -9.99999999999999962e134Initial program 99.7%
Applied egg-rr96.9%
if -9.99999999999999962e134 < (*.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))))) < 9.9999999999999994e303Initial program 99.8%
lift-/.f64N/A
lift-cos.f64N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6499.8
lift-/.f64N/A
div-invN/A
metadata-evalN/A
lower-*.f6499.8
Applied egg-rr99.8%
Applied egg-rr95.5%
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
times-fracN/A
lift-*.f64N/A
associate-/r*N/A
associate-*r/N/A
lower-/.f64N/A
Applied egg-rr95.5%
if 9.9999999999999994e303 < (*.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.4%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Taylor expanded in K around 0
Simplified71.2%
Taylor expanded in J around 0
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Final simplification86.0%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* t_0 (* -2.0 J))
(sqrt (+ 1.0 (pow (/ U_m (* t_0 (* J 2.0))) 2.0)))))
(t_2 (cos (* K 0.5))))
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -2e+205)
(*
-2.0
(*
(sqrt
(fma
U_m
(/
U_m
(*
(+ 0.5 (* 0.5 (cos (* 2.0 (* K 0.5)))))
(* (* J 2.0) (* J 2.0))))
1.0))
(* J t_2)))
(if (<= t_1 1e+304)
(*
(* J (* -2.0 t_2))
(sqrt
(+
1.0
(/
(* (/ U_m (* (* J 2.0) (fma 0.5 (cos K) 0.5))) (* U_m 0.5))
J))))
(fma 2.0 (/ (* J J) U_m) 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 = (t_0 * (-2.0 * J)) * sqrt((1.0 + pow((U_m / (t_0 * (J * 2.0))), 2.0)));
double t_2 = cos((K * 0.5));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -2e+205) {
tmp = -2.0 * (sqrt(fma(U_m, (U_m / ((0.5 + (0.5 * cos((2.0 * (K * 0.5))))) * ((J * 2.0) * (J * 2.0)))), 1.0)) * (J * t_2));
} else if (t_1 <= 1e+304) {
tmp = (J * (-2.0 * t_2)) * sqrt((1.0 + (((U_m / ((J * 2.0) * fma(0.5, cos(K), 0.5))) * (U_m * 0.5)) / J)));
} else {
tmp = fma(2.0, ((J * J) / U_m), 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(t_0 * Float64(-2.0 * J)) * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_0 * Float64(J * 2.0))) ^ 2.0)))) t_2 = cos(Float64(K * 0.5)) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -2e+205) tmp = Float64(-2.0 * Float64(sqrt(fma(U_m, Float64(U_m / Float64(Float64(0.5 + Float64(0.5 * cos(Float64(2.0 * Float64(K * 0.5))))) * Float64(Float64(J * 2.0) * Float64(J * 2.0)))), 1.0)) * Float64(J * t_2))); elseif (t_1 <= 1e+304) tmp = Float64(Float64(J * Float64(-2.0 * t_2)) * sqrt(Float64(1.0 + Float64(Float64(Float64(U_m / Float64(Float64(J * 2.0) * fma(0.5, cos(K), 0.5))) * Float64(U_m * 0.5)) / J)))); else tmp = fma(2.0, Float64(Float64(J * J) / U_m), 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[(t$95$0 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$0 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -2e+205], N[(-2.0 * N[(N[Sqrt[N[(U$95$m * N[(U$95$m / N[(N[(0.5 + N[(0.5 * N[Cos[N[(2.0 * N[(K * 0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(J * 2.0), $MachinePrecision] * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision] * N[(J * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 1e+304], N[(N[(J * N[(-2.0 * t$95$2), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[(N[(N[(U$95$m / N[(N[(J * 2.0), $MachinePrecision] * N[(0.5 * N[Cos[K], $MachinePrecision] + 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(U$95$m * 0.5), $MachinePrecision]), $MachinePrecision] / J), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(J * J), $MachinePrecision] / U$95$m), $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(t\_0 \cdot \left(-2 \cdot J\right)\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_0 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
t_2 := \cos \left(K \cdot 0.5\right)\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{+205}:\\
\;\;\;\;-2 \cdot \left(\sqrt{\mathsf{fma}\left(U\_m, \frac{U\_m}{\left(0.5 + 0.5 \cdot \cos \left(2 \cdot \left(K \cdot 0.5\right)\right)\right) \cdot \left(\left(J \cdot 2\right) \cdot \left(J \cdot 2\right)\right)}, 1\right)} \cdot \left(J \cdot t\_2\right)\right)\\
\mathbf{elif}\;t\_1 \leq 10^{+304}:\\
\;\;\;\;\left(J \cdot \left(-2 \cdot t\_2\right)\right) \cdot \sqrt{1 + \frac{\frac{U\_m}{\left(J \cdot 2\right) \cdot \mathsf{fma}\left(0.5, \cos K, 0.5\right)} \cdot \left(U\_m \cdot 0.5\right)}{J}}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(2, \frac{J \cdot J}{U\_m}, 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.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6444.4
Simplified44.4%
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))))) < -2.00000000000000003e205Initial program 99.7%
Applied egg-rr95.4%
if -2.00000000000000003e205 < (*.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))))) < 9.9999999999999994e303Initial program 99.8%
lift-/.f64N/A
lift-cos.f64N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6499.8
lift-/.f64N/A
div-invN/A
metadata-evalN/A
lower-*.f6499.8
Applied egg-rr99.8%
Applied egg-rr95.8%
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
times-fracN/A
lift-*.f64N/A
associate-/r*N/A
associate-*r/N/A
lower-/.f64N/A
Applied egg-rr95.8%
if 9.9999999999999994e303 < (*.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.4%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Taylor expanded in K around 0
Simplified71.2%
Taylor expanded in J around 0
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Final simplification86.0%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* t_0 (* -2.0 J))
(sqrt (+ 1.0 (pow (/ U_m (* t_0 (* J 2.0))) 2.0)))))
(t_2 (/ U_m (* -2.0 J))))
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -5e-245)
(* (* -2.0 J) (sqrt (fma t_2 t_2 1.0)))
(if (<= t_1 1e+304)
(* (cos (* K 0.5)) (* -2.0 J))
(fma 2.0 (/ (* J J) U_m) 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 = (t_0 * (-2.0 * J)) * sqrt((1.0 + pow((U_m / (t_0 * (J * 2.0))), 2.0)));
double t_2 = U_m / (-2.0 * J);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -5e-245) {
tmp = (-2.0 * J) * sqrt(fma(t_2, t_2, 1.0));
} else if (t_1 <= 1e+304) {
tmp = cos((K * 0.5)) * (-2.0 * J);
} else {
tmp = fma(2.0, ((J * J) / U_m), 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(t_0 * Float64(-2.0 * J)) * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_0 * Float64(J * 2.0))) ^ 2.0)))) t_2 = Float64(U_m / Float64(-2.0 * J)) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -5e-245) tmp = Float64(Float64(-2.0 * J) * sqrt(fma(t_2, t_2, 1.0))); elseif (t_1 <= 1e+304) tmp = Float64(cos(Float64(K * 0.5)) * Float64(-2.0 * J)); else tmp = fma(2.0, Float64(Float64(J * J) / U_m), 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[(t$95$0 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$0 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(U$95$m / N[(-2.0 * J), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -5e-245], N[(N[(-2.0 * J), $MachinePrecision] * N[Sqrt[N[(t$95$2 * t$95$2 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 1e+304], N[(N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision] * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(J * J), $MachinePrecision] / U$95$m), $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(t\_0 \cdot \left(-2 \cdot J\right)\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_0 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
t_2 := \frac{U\_m}{-2 \cdot J}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -5 \cdot 10^{-245}:\\
\;\;\;\;\left(-2 \cdot J\right) \cdot \sqrt{\mathsf{fma}\left(t\_2, t\_2, 1\right)}\\
\mathbf{elif}\;t\_1 \leq 10^{+304}:\\
\;\;\;\;\cos \left(K \cdot 0.5\right) \cdot \left(-2 \cdot J\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(2, \frac{J \cdot J}{U\_m}, 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.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6444.4
Simplified44.4%
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.9999999999999997e-245Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f6447.5
Simplified47.5%
Taylor expanded in K around 0
lower-*.f6453.0
Simplified53.0%
lift-*.f64N/A
lift-/.f64N/A
lift-pow.f64N/A
+-commutativeN/A
lift-pow.f64N/A
unpow2N/A
lift-/.f64N/A
lift-/.f64N/A
frac-2negN/A
distribute-frac-negN/A
frac-2negN/A
distribute-frac-negN/A
sqr-negN/A
lower-fma.f64N/A
Applied egg-rr53.0%
if -4.9999999999999997e-245 < (*.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))))) < 9.9999999999999994e303Initial program 99.8%
Taylor expanded in J around inf
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6466.2
Simplified66.2%
if 9.9999999999999994e303 < (*.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.4%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Taylor expanded in K around 0
Simplified71.2%
Taylor expanded in J around 0
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Final simplification59.4%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* t_0 (* -2.0 J))
(sqrt (+ 1.0 (pow (/ U_m (* t_0 (* J 2.0))) 2.0))))))
(if (<= t_1 -2e+276)
(- U_m)
(if (<= t_1 -4e-154)
(* (* -2.0 J) (sqrt (fma 0.25 (/ (* U_m U_m) (* J J)) 1.0)))
(if (<= t_1 -1e-306) (- U_m) (fma 2.0 (/ (* J J) U_m) 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 = (t_0 * (-2.0 * J)) * sqrt((1.0 + pow((U_m / (t_0 * (J * 2.0))), 2.0)));
double tmp;
if (t_1 <= -2e+276) {
tmp = -U_m;
} else if (t_1 <= -4e-154) {
tmp = (-2.0 * J) * sqrt(fma(0.25, ((U_m * U_m) / (J * J)), 1.0));
} else if (t_1 <= -1e-306) {
tmp = -U_m;
} else {
tmp = fma(2.0, ((J * J) / U_m), 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(t_0 * Float64(-2.0 * J)) * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_0 * Float64(J * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_1 <= -2e+276) tmp = Float64(-U_m); elseif (t_1 <= -4e-154) tmp = Float64(Float64(-2.0 * J) * sqrt(fma(0.25, Float64(Float64(U_m * U_m) / Float64(J * J)), 1.0))); elseif (t_1 <= -1e-306) tmp = Float64(-U_m); else tmp = fma(2.0, Float64(Float64(J * J) / U_m), 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[(t$95$0 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$0 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2e+276], (-U$95$m), If[LessEqual[t$95$1, -4e-154], N[(N[(-2.0 * J), $MachinePrecision] * N[Sqrt[N[(0.25 * N[(N[(U$95$m * U$95$m), $MachinePrecision] / N[(J * J), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -1e-306], (-U$95$m), N[(2.0 * N[(N[(J * J), $MachinePrecision] / U$95$m), $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(t\_0 \cdot \left(-2 \cdot J\right)\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_0 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{+276}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -4 \cdot 10^{-154}:\\
\;\;\;\;\left(-2 \cdot J\right) \cdot \sqrt{\mathsf{fma}\left(0.25, \frac{U\_m \cdot U\_m}{J \cdot J}, 1\right)}\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-306}:\\
\;\;\;\;-U\_m\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(2, \frac{J \cdot J}{U\_m}, 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))))) < -2.0000000000000001e276 or -3.9999999999999999e-154 < (*.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.00000000000000003e-306Initial program 25.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6437.4
Simplified37.4%
if -2.0000000000000001e276 < (*.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))))) < -3.9999999999999999e-154Initial program 99.8%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6440.9
Simplified40.9%
if -1.00000000000000003e-306 < (*.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 77.6%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6429.6
Simplified29.6%
Taylor expanded in K around 0
Simplified29.6%
Taylor expanded in J around 0
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6431.1
Simplified31.1%
Final simplification35.3%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* t_0 (* -2.0 J))
(sqrt (+ 1.0 (pow (/ U_m (* t_0 (* J 2.0))) 2.0))))))
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -2e-104)
(* -2.0 J)
(if (<= t_1 -1e-306) (- U_m) (fma 2.0 (/ (* J J) U_m) 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 = (t_0 * (-2.0 * J)) * sqrt((1.0 + pow((U_m / (t_0 * (J * 2.0))), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -2e-104) {
tmp = -2.0 * J;
} else if (t_1 <= -1e-306) {
tmp = -U_m;
} else {
tmp = fma(2.0, ((J * J) / U_m), 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(t_0 * Float64(-2.0 * J)) * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_0 * Float64(J * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -2e-104) tmp = Float64(-2.0 * J); elseif (t_1 <= -1e-306) tmp = Float64(-U_m); else tmp = fma(2.0, Float64(Float64(J * J) / U_m), 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[(t$95$0 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$0 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -2e-104], N[(-2.0 * J), $MachinePrecision], If[LessEqual[t$95$1, -1e-306], (-U$95$m), N[(2.0 * N[(N[(J * J), $MachinePrecision] / U$95$m), $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(t\_0 \cdot \left(-2 \cdot J\right)\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_0 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-104}:\\
\;\;\;\;-2 \cdot J\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-306}:\\
\;\;\;\;-U\_m\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(2, \frac{J \cdot J}{U\_m}, 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.0 or -1.99999999999999985e-104 < (*.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.00000000000000003e-306Initial program 27.5%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6436.8
Simplified36.8%
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.99999999999999985e-104Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f6447.7
Simplified47.7%
Taylor expanded in J around inf
*-commutativeN/A
lower-*.f6434.4
Simplified34.4%
if -1.00000000000000003e-306 < (*.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 77.6%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6429.6
Simplified29.6%
Taylor expanded in K around 0
Simplified29.6%
Taylor expanded in J around 0
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6431.1
Simplified31.1%
Final simplification33.1%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1 (* t_0 (* -2.0 J)))
(t_2 (* t_1 (sqrt (+ 1.0 (pow (/ U_m (* t_0 (* J 2.0))) 2.0))))))
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 1e+304)
(*
t_1
(sqrt
(+
1.0
(*
(/ U_m (* J 2.0))
(/ U_m (* (* J 2.0) (+ 0.5 (* 0.5 (cos (* 2.0 (* K 0.5)))))))))))
(fma 2.0 (/ (* J J) U_m) 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 = t_0 * (-2.0 * J);
double t_2 = t_1 * sqrt((1.0 + pow((U_m / (t_0 * (J * 2.0))), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= 1e+304) {
tmp = t_1 * sqrt((1.0 + ((U_m / (J * 2.0)) * (U_m / ((J * 2.0) * (0.5 + (0.5 * cos((2.0 * (K * 0.5))))))))));
} else {
tmp = fma(2.0, ((J * J) / U_m), U_m);
}
return tmp;
}
U_m = abs(U) function code(J, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(t_0 * Float64(-2.0 * J)) t_2 = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_0 * Float64(J * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= 1e+304) tmp = Float64(t_1 * sqrt(Float64(1.0 + Float64(Float64(U_m / Float64(J * 2.0)) * Float64(U_m / Float64(Float64(J * 2.0) * Float64(0.5 + Float64(0.5 * cos(Float64(2.0 * Float64(K * 0.5))))))))))); else tmp = fma(2.0, Float64(Float64(J * J) / U_m), 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[(t$95$0 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$0 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, 1e+304], N[(t$95$1 * N[Sqrt[N[(1.0 + N[(N[(U$95$m / N[(J * 2.0), $MachinePrecision]), $MachinePrecision] * N[(U$95$m / N[(N[(J * 2.0), $MachinePrecision] * N[(0.5 + N[(0.5 * N[Cos[N[(2.0 * N[(K * 0.5), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(J * J), $MachinePrecision] / U$95$m), $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 := t\_0 \cdot \left(-2 \cdot J\right)\\
t_2 := t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_0 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 10^{+304}:\\
\;\;\;\;t\_1 \cdot \sqrt{1 + \frac{U\_m}{J \cdot 2} \cdot \frac{U\_m}{\left(J \cdot 2\right) \cdot \left(0.5 + 0.5 \cdot \cos \left(2 \cdot \left(K \cdot 0.5\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(2, \frac{J \cdot J}{U\_m}, 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.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6444.4
Simplified44.4%
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))))) < 9.9999999999999994e303Initial program 99.8%
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lift-/.f64N/A
unpow2N/A
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
associate-*l/N/A
lift-/.f64N/A
associate-*r/N/A
Applied egg-rr99.6%
if 9.9999999999999994e303 < (*.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.4%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Taylor expanded in K around 0
Simplified71.2%
Taylor expanded in J around 0
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Final simplification88.9%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* t_0 (* -2.0 J))
(sqrt (+ 1.0 (pow (/ U_m (* t_0 (* J 2.0))) 2.0))))))
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 1e+304)
(*
(* J (* -2.0 (cos (* K 0.5))))
(sqrt
(+
1.0
(/ (* (/ U_m (* (* J 2.0) (fma 0.5 (cos K) 0.5))) (* U_m 0.5)) J))))
(fma 2.0 (/ (* J J) U_m) 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 = (t_0 * (-2.0 * J)) * sqrt((1.0 + pow((U_m / (t_0 * (J * 2.0))), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= 1e+304) {
tmp = (J * (-2.0 * cos((K * 0.5)))) * sqrt((1.0 + (((U_m / ((J * 2.0) * fma(0.5, cos(K), 0.5))) * (U_m * 0.5)) / J)));
} else {
tmp = fma(2.0, ((J * J) / U_m), 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(t_0 * Float64(-2.0 * J)) * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_0 * Float64(J * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= 1e+304) tmp = Float64(Float64(J * Float64(-2.0 * cos(Float64(K * 0.5)))) * sqrt(Float64(1.0 + Float64(Float64(Float64(U_m / Float64(Float64(J * 2.0) * fma(0.5, cos(K), 0.5))) * Float64(U_m * 0.5)) / J)))); else tmp = fma(2.0, Float64(Float64(J * J) / U_m), 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[(t$95$0 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$0 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, 1e+304], N[(N[(J * N[(-2.0 * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[(N[(N[(U$95$m / N[(N[(J * 2.0), $MachinePrecision] * N[(0.5 * N[Cos[K], $MachinePrecision] + 0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(U$95$m * 0.5), $MachinePrecision]), $MachinePrecision] / J), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(J * J), $MachinePrecision] / U$95$m), $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(t\_0 \cdot \left(-2 \cdot J\right)\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_0 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq 10^{+304}:\\
\;\;\;\;\left(J \cdot \left(-2 \cdot \cos \left(K \cdot 0.5\right)\right)\right) \cdot \sqrt{1 + \frac{\frac{U\_m}{\left(J \cdot 2\right) \cdot \mathsf{fma}\left(0.5, \cos K, 0.5\right)} \cdot \left(U\_m \cdot 0.5\right)}{J}}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(2, \frac{J \cdot J}{U\_m}, 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.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6444.4
Simplified44.4%
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))))) < 9.9999999999999994e303Initial program 99.8%
lift-/.f64N/A
lift-cos.f64N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6499.8
lift-/.f64N/A
div-invN/A
metadata-evalN/A
lower-*.f6499.8
Applied egg-rr99.8%
Applied egg-rr94.3%
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lift-+.f64N/A
lift-*.f64N/A
times-fracN/A
lift-*.f64N/A
associate-/r*N/A
associate-*r/N/A
lower-/.f64N/A
Applied egg-rr94.3%
if 9.9999999999999994e303 < (*.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.4%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Taylor expanded in K around 0
Simplified71.2%
Taylor expanded in J around 0
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6471.2
Simplified71.2%
Final simplification84.9%
U_m = (fabs.f64 U)
(FPCore (J K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* t_0 (* -2.0 J))
(sqrt (+ 1.0 (pow (/ U_m (* t_0 (* J 2.0))) 2.0)))))
(t_2 (/ U_m (* -2.0 J))))
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -1e-306)
(* (* -2.0 J) (sqrt (fma t_2 t_2 1.0)))
(fma 2.0 (/ (* J J) U_m) 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 = (t_0 * (-2.0 * J)) * sqrt((1.0 + pow((U_m / (t_0 * (J * 2.0))), 2.0)));
double t_2 = U_m / (-2.0 * J);
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -1e-306) {
tmp = (-2.0 * J) * sqrt(fma(t_2, t_2, 1.0));
} else {
tmp = fma(2.0, ((J * J) / U_m), 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(t_0 * Float64(-2.0 * J)) * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_0 * Float64(J * 2.0))) ^ 2.0)))) t_2 = Float64(U_m / Float64(-2.0 * J)) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -1e-306) tmp = Float64(Float64(-2.0 * J) * sqrt(fma(t_2, t_2, 1.0))); else tmp = fma(2.0, Float64(Float64(J * J) / U_m), 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[(t$95$0 * N[(-2.0 * J), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$0 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(U$95$m / N[(-2.0 * J), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -1e-306], N[(N[(-2.0 * J), $MachinePrecision] * N[Sqrt[N[(t$95$2 * t$95$2 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(J * J), $MachinePrecision] / U$95$m), $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(t\_0 \cdot \left(-2 \cdot J\right)\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_0 \cdot \left(J \cdot 2\right)}\right)}^{2}}\\
t_2 := \frac{U\_m}{-2 \cdot J}\\
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-306}:\\
\;\;\;\;\left(-2 \cdot J\right) \cdot \sqrt{\mathsf{fma}\left(t\_2, t\_2, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(2, \frac{J \cdot J}{U\_m}, 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.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6444.4
Simplified44.4%
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.00000000000000003e-306Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f6446.9
Simplified46.9%
Taylor expanded in K around 0
lower-*.f6452.3
Simplified52.3%
lift-*.f64N/A
lift-/.f64N/A
lift-pow.f64N/A
+-commutativeN/A
lift-pow.f64N/A
unpow2N/A
lift-/.f64N/A
lift-/.f64N/A
frac-2negN/A
distribute-frac-negN/A
frac-2negN/A
distribute-frac-negN/A
sqr-negN/A
lower-fma.f64N/A
Applied egg-rr52.3%
if -1.00000000000000003e-306 < (*.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 77.6%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6429.6
Simplified29.6%
Taylor expanded in K around 0
Simplified29.6%
Taylor expanded in J around 0
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f6431.1
Simplified31.1%
Final simplification40.4%
U_m = (fabs.f64 U) (FPCore (J K U_m) :precision binary64 (if (<= U_m 270000000000.0) (* -2.0 J) (- U_m)))
U_m = fabs(U);
double code(double J, double K, double U_m) {
double tmp;
if (U_m <= 270000000000.0) {
tmp = -2.0 * J;
} else {
tmp = -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) :: tmp
if (u_m <= 270000000000.0d0) then
tmp = (-2.0d0) * j
else
tmp = -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 tmp;
if (U_m <= 270000000000.0) {
tmp = -2.0 * J;
} else {
tmp = -U_m;
}
return tmp;
}
U_m = math.fabs(U) def code(J, K, U_m): tmp = 0 if U_m <= 270000000000.0: tmp = -2.0 * J else: tmp = -U_m return tmp
U_m = abs(U) function code(J, K, U_m) tmp = 0.0 if (U_m <= 270000000000.0) tmp = Float64(-2.0 * J); else tmp = Float64(-U_m); end return tmp end
U_m = abs(U); function tmp_2 = code(J, K, U_m) tmp = 0.0; if (U_m <= 270000000000.0) tmp = -2.0 * J; else tmp = -U_m; end tmp_2 = tmp; end
U_m = N[Abs[U], $MachinePrecision] code[J_, K_, U$95$m_] := If[LessEqual[U$95$m, 270000000000.0], N[(-2.0 * J), $MachinePrecision], (-U$95$m)]
\begin{array}{l}
U_m = \left|U\right|
\\
\begin{array}{l}
\mathbf{if}\;U\_m \leq 270000000000:\\
\;\;\;\;-2 \cdot J\\
\mathbf{else}:\\
\;\;\;\;-U\_m\\
\end{array}
\end{array}
if U < 2.7e11Initial program 84.7%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f6444.6
Simplified44.6%
Taylor expanded in J around inf
*-commutativeN/A
lower-*.f6431.4
Simplified31.4%
if 2.7e11 < U Initial program 49.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6433.0
Simplified33.0%
Final simplification31.8%
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 75.9%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6420.9
Simplified20.9%
herbie shell --seed 2024219
(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)))))