
(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)));
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(j, k, u)
use fmin_fmax_functions
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}
Herbie found 18 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)));
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(j, k, u)
use fmin_fmax_functions
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)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (* 0.5 K)))
(t_1 (cos (/ K 2.0)))
(t_2
(*
(* (* -2.0 J_m) t_1)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_1)) 2.0))))))
(*
J_s
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 2e+298)
(*
(* (* -2.0 J_m) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (+ J_m J_m) t_0)) 2.0))))
U_m)))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((0.5 * K));
double t_1 = cos((K / 2.0));
double t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_1)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= 2e+298) {
tmp = ((-2.0 * J_m) * t_0) * sqrt((1.0 + pow((U_m / ((J_m + J_m) * t_0)), 2.0)));
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = Math.abs(U);
J\_m = Math.abs(J);
J\_s = Math.copySign(1.0, J);
public static double code(double J_s, double J_m, double K, double U_m) {
double t_0 = Math.cos((0.5 * K));
double t_1 = Math.cos((K / 2.0));
double t_2 = ((-2.0 * J_m) * t_1) * Math.sqrt((1.0 + Math.pow((U_m / ((2.0 * J_m) * t_1)), 2.0)));
double tmp;
if (t_2 <= -Double.POSITIVE_INFINITY) {
tmp = -U_m;
} else if (t_2 <= 2e+298) {
tmp = ((-2.0 * J_m) * t_0) * Math.sqrt((1.0 + Math.pow((U_m / ((J_m + J_m) * t_0)), 2.0)));
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = math.fabs(U) J\_m = math.fabs(J) J\_s = math.copysign(1.0, J) def code(J_s, J_m, K, U_m): t_0 = math.cos((0.5 * K)) t_1 = math.cos((K / 2.0)) t_2 = ((-2.0 * J_m) * t_1) * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J_m) * t_1)), 2.0))) tmp = 0 if t_2 <= -math.inf: tmp = -U_m elif t_2 <= 2e+298: tmp = ((-2.0 * J_m) * t_0) * math.sqrt((1.0 + math.pow((U_m / ((J_m + J_m) * t_0)), 2.0))) else: tmp = U_m return J_s * tmp
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(0.5 * K)) t_1 = cos(Float64(K / 2.0)) t_2 = Float64(Float64(Float64(-2.0 * J_m) * t_1) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_1)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= 2e+298) tmp = Float64(Float64(Float64(-2.0 * J_m) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(J_m + J_m) * t_0)) ^ 2.0)))); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = abs(U); J\_m = abs(J); J\_s = sign(J) * abs(1.0); function tmp_2 = code(J_s, J_m, K, U_m) t_0 = cos((0.5 * K)); t_1 = cos((K / 2.0)); t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + ((U_m / ((2.0 * J_m) * t_1)) ^ 2.0))); tmp = 0.0; if (t_2 <= -Inf) tmp = -U_m; elseif (t_2 <= 2e+298) tmp = ((-2.0 * J_m) * t_0) * sqrt((1.0 + ((U_m / ((J_m + J_m) * t_0)) ^ 2.0))); else tmp = U_m; end tmp_2 = J_s * tmp; end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, 2e+298], N[(N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(J$95$m + J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], U$95$m]]), $MachinePrecision]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(0.5 \cdot K\right)\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := \left(\left(-2 \cdot J\_m\right) \cdot t\_1\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_1}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+298}:\\
\;\;\;\;\left(\left(-2 \cdot J\_m\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(J\_m + J\_m\right) \cdot t\_0}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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.9999999999999999e298Initial program 99.8%
Taylor expanded in K around 0
lower-*.f6499.8
Applied rewrites99.8%
Taylor expanded in K around 0
lower-*.f6499.8
Applied rewrites99.8%
lift-*.f64N/A
count-2-revN/A
lift-+.f6499.8
Applied rewrites99.8%
if 1.9999999999999999e298 < (*.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 12.1%
Taylor expanded in U around -inf
Applied rewrites94.5%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (* 0.5 K)))
(t_1 (cos (/ K 2.0)))
(t_2
(*
(* (* -2.0 J_m) t_1)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_1)) 2.0))))))
(*
J_s
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 2e+298)
(*
(* J_m -2.0)
(* t_0 (sqrt (+ (pow (/ U_m (* (+ J_m J_m) t_0)) 2.0) 1.0))))
U_m)))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((0.5 * K));
double t_1 = cos((K / 2.0));
double t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_1)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= 2e+298) {
tmp = (J_m * -2.0) * (t_0 * sqrt((pow((U_m / ((J_m + J_m) * t_0)), 2.0) + 1.0)));
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = Math.abs(U);
J\_m = Math.abs(J);
J\_s = Math.copySign(1.0, J);
public static double code(double J_s, double J_m, double K, double U_m) {
double t_0 = Math.cos((0.5 * K));
double t_1 = Math.cos((K / 2.0));
double t_2 = ((-2.0 * J_m) * t_1) * Math.sqrt((1.0 + Math.pow((U_m / ((2.0 * J_m) * t_1)), 2.0)));
double tmp;
if (t_2 <= -Double.POSITIVE_INFINITY) {
tmp = -U_m;
} else if (t_2 <= 2e+298) {
tmp = (J_m * -2.0) * (t_0 * Math.sqrt((Math.pow((U_m / ((J_m + J_m) * t_0)), 2.0) + 1.0)));
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = math.fabs(U) J\_m = math.fabs(J) J\_s = math.copysign(1.0, J) def code(J_s, J_m, K, U_m): t_0 = math.cos((0.5 * K)) t_1 = math.cos((K / 2.0)) t_2 = ((-2.0 * J_m) * t_1) * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J_m) * t_1)), 2.0))) tmp = 0 if t_2 <= -math.inf: tmp = -U_m elif t_2 <= 2e+298: tmp = (J_m * -2.0) * (t_0 * math.sqrt((math.pow((U_m / ((J_m + J_m) * t_0)), 2.0) + 1.0))) else: tmp = U_m return J_s * tmp
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(0.5 * K)) t_1 = cos(Float64(K / 2.0)) t_2 = Float64(Float64(Float64(-2.0 * J_m) * t_1) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_1)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= 2e+298) tmp = Float64(Float64(J_m * -2.0) * Float64(t_0 * sqrt(Float64((Float64(U_m / Float64(Float64(J_m + J_m) * t_0)) ^ 2.0) + 1.0)))); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = abs(U); J\_m = abs(J); J\_s = sign(J) * abs(1.0); function tmp_2 = code(J_s, J_m, K, U_m) t_0 = cos((0.5 * K)); t_1 = cos((K / 2.0)); t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + ((U_m / ((2.0 * J_m) * t_1)) ^ 2.0))); tmp = 0.0; if (t_2 <= -Inf) tmp = -U_m; elseif (t_2 <= 2e+298) tmp = (J_m * -2.0) * (t_0 * sqrt((((U_m / ((J_m + J_m) * t_0)) ^ 2.0) + 1.0))); else tmp = U_m; end tmp_2 = J_s * tmp; end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, 2e+298], N[(N[(J$95$m * -2.0), $MachinePrecision] * N[(t$95$0 * N[Sqrt[N[(N[Power[N[(U$95$m / N[(N[(J$95$m + J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], U$95$m]]), $MachinePrecision]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(0.5 \cdot K\right)\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := \left(\left(-2 \cdot J\_m\right) \cdot t\_1\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_1}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+298}:\\
\;\;\;\;\left(J\_m \cdot -2\right) \cdot \left(t\_0 \cdot \sqrt{{\left(\frac{U\_m}{\left(J\_m + J\_m\right) \cdot t\_0}\right)}^{2} + 1}\right)\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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.9999999999999999e298Initial program 99.8%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
lift-sqrt.f64N/A
lift-+.f64N/A
lift-pow.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
Applied rewrites97.6%
Taylor expanded in K around 0
lift-*.f6497.6
Applied rewrites97.6%
Taylor expanded in K around 0
lift-*.f6497.6
Applied rewrites97.6%
lift-cosh.f64N/A
lift-asinh.f64N/A
cosh-asinhN/A
lower-sqrt.f64N/A
unpow2N/A
lift-pow.f64N/A
lower-+.f6499.7
Applied rewrites99.7%
if 1.9999999999999999e298 < (*.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 12.1%
Taylor expanded in U around -inf
Applied rewrites94.5%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (* 0.5 K)))
(t_1 (cos (/ K 2.0)))
(t_2
(*
(* (* -2.0 J_m) t_1)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_1)) 2.0))))))
(*
J_s
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 2e+298)
(* (* (* -2.0 J_m) t_0) (cosh (asinh (/ U_m (* (+ J_m J_m) t_0)))))
U_m)))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((0.5 * K));
double t_1 = cos((K / 2.0));
double t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_1)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= 2e+298) {
tmp = ((-2.0 * J_m) * t_0) * cosh(asinh((U_m / ((J_m + J_m) * t_0))));
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = math.fabs(U) J\_m = math.fabs(J) J\_s = math.copysign(1.0, J) def code(J_s, J_m, K, U_m): t_0 = math.cos((0.5 * K)) t_1 = math.cos((K / 2.0)) t_2 = ((-2.0 * J_m) * t_1) * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J_m) * t_1)), 2.0))) tmp = 0 if t_2 <= -math.inf: tmp = -U_m elif t_2 <= 2e+298: tmp = ((-2.0 * J_m) * t_0) * math.cosh(math.asinh((U_m / ((J_m + J_m) * t_0)))) else: tmp = U_m return J_s * tmp
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(0.5 * K)) t_1 = cos(Float64(K / 2.0)) t_2 = Float64(Float64(Float64(-2.0 * J_m) * t_1) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_1)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= 2e+298) tmp = Float64(Float64(Float64(-2.0 * J_m) * t_0) * cosh(asinh(Float64(U_m / Float64(Float64(J_m + J_m) * t_0))))); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = abs(U); J\_m = abs(J); J\_s = sign(J) * abs(1.0); function tmp_2 = code(J_s, J_m, K, U_m) t_0 = cos((0.5 * K)); t_1 = cos((K / 2.0)); t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + ((U_m / ((2.0 * J_m) * t_1)) ^ 2.0))); tmp = 0.0; if (t_2 <= -Inf) tmp = -U_m; elseif (t_2 <= 2e+298) tmp = ((-2.0 * J_m) * t_0) * cosh(asinh((U_m / ((J_m + J_m) * t_0)))); else tmp = U_m; end tmp_2 = J_s * tmp; end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, 2e+298], N[(N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Cosh[N[ArcSinh[N[(U$95$m / N[(N[(J$95$m + J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], U$95$m]]), $MachinePrecision]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(0.5 \cdot K\right)\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := \left(\left(-2 \cdot J\_m\right) \cdot t\_1\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_1}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+298}:\\
\;\;\;\;\left(\left(-2 \cdot J\_m\right) \cdot t\_0\right) \cdot \cosh \sinh^{-1} \left(\frac{U\_m}{\left(J\_m + J\_m\right) \cdot t\_0}\right)\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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.9999999999999999e298Initial program 99.8%
lift-sqrt.f64N/A
lift-+.f64N/A
lift-pow.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
+-commutativeN/A
unpow2N/A
cosh-asinh-revN/A
lower-cosh.f64N/A
lower-asinh.f64N/A
Applied rewrites97.6%
Taylor expanded in K around inf
lift-cos.f64N/A
lift-*.f6497.6
Applied rewrites97.6%
Taylor expanded in K around inf
lift-cos.f64N/A
lift-*.f6497.6
Applied rewrites97.6%
if 1.9999999999999999e298 < (*.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 12.1%
Taylor expanded in U around -inf
Applied rewrites94.5%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (* 0.5 K)))
(t_1 (cos (/ K 2.0)))
(t_2
(*
(* (* -2.0 J_m) t_1)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_1)) 2.0))))))
(*
J_s
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 2e+298)
(* (* J_m -2.0) (* t_0 (cosh (asinh (/ U_m (* (+ J_m J_m) t_0))))))
U_m)))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((0.5 * K));
double t_1 = cos((K / 2.0));
double t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_1)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= 2e+298) {
tmp = (J_m * -2.0) * (t_0 * cosh(asinh((U_m / ((J_m + J_m) * t_0)))));
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = math.fabs(U) J\_m = math.fabs(J) J\_s = math.copysign(1.0, J) def code(J_s, J_m, K, U_m): t_0 = math.cos((0.5 * K)) t_1 = math.cos((K / 2.0)) t_2 = ((-2.0 * J_m) * t_1) * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J_m) * t_1)), 2.0))) tmp = 0 if t_2 <= -math.inf: tmp = -U_m elif t_2 <= 2e+298: tmp = (J_m * -2.0) * (t_0 * math.cosh(math.asinh((U_m / ((J_m + J_m) * t_0))))) else: tmp = U_m return J_s * tmp
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(0.5 * K)) t_1 = cos(Float64(K / 2.0)) t_2 = Float64(Float64(Float64(-2.0 * J_m) * t_1) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_1)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= 2e+298) tmp = Float64(Float64(J_m * -2.0) * Float64(t_0 * cosh(asinh(Float64(U_m / Float64(Float64(J_m + J_m) * t_0)))))); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = abs(U); J\_m = abs(J); J\_s = sign(J) * abs(1.0); function tmp_2 = code(J_s, J_m, K, U_m) t_0 = cos((0.5 * K)); t_1 = cos((K / 2.0)); t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + ((U_m / ((2.0 * J_m) * t_1)) ^ 2.0))); tmp = 0.0; if (t_2 <= -Inf) tmp = -U_m; elseif (t_2 <= 2e+298) tmp = (J_m * -2.0) * (t_0 * cosh(asinh((U_m / ((J_m + J_m) * t_0))))); else tmp = U_m; end tmp_2 = J_s * tmp; end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, 2e+298], N[(N[(J$95$m * -2.0), $MachinePrecision] * N[(t$95$0 * N[Cosh[N[ArcSinh[N[(U$95$m / N[(N[(J$95$m + J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], U$95$m]]), $MachinePrecision]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(0.5 \cdot K\right)\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := \left(\left(-2 \cdot J\_m\right) \cdot t\_1\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_1}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+298}:\\
\;\;\;\;\left(J\_m \cdot -2\right) \cdot \left(t\_0 \cdot \cosh \sinh^{-1} \left(\frac{U\_m}{\left(J\_m + J\_m\right) \cdot t\_0}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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.9999999999999999e298Initial program 99.8%
lift-*.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
lift-sqrt.f64N/A
lift-+.f64N/A
lift-pow.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
Applied rewrites97.6%
Taylor expanded in K around 0
lift-*.f6497.6
Applied rewrites97.6%
Taylor expanded in K around 0
lift-*.f6497.6
Applied rewrites97.6%
if 1.9999999999999999e298 < (*.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 12.1%
Taylor expanded in U around -inf
Applied rewrites94.5%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1 (* (* -2.0 J_m) t_0))
(t_2 (* t_1 (sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_0)) 2.0)))))
(t_3
(*
(* (* (cos (* 0.5 K)) J_m) -2.0)
(sqrt
(fma
(/
(* U_m U_m)
(* (+ 0.5 (* 0.5 (cos (* 2.0 (* 0.5 K))))) (* J_m J_m)))
0.25
1.0)))))
(*
J_s
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 -4e+159)
(* t_1 (sqrt (+ 1.0 (pow (/ U_m (+ J_m J_m)) 2.0))))
(if (<= t_2 -5e-91)
t_3
(if (<= t_2 5e-185)
(* (* (sqrt (fma (* (/ U_m J_m) (/ U_m J_m)) 0.25 1.0)) J_m) -2.0)
(if (<= t_2 2e+298) t_3 U_m))))))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = (-2.0 * J_m) * t_0;
double t_2 = t_1 * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double t_3 = ((cos((0.5 * K)) * J_m) * -2.0) * sqrt(fma(((U_m * U_m) / ((0.5 + (0.5 * cos((2.0 * (0.5 * K))))) * (J_m * J_m))), 0.25, 1.0));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= -4e+159) {
tmp = t_1 * sqrt((1.0 + pow((U_m / (J_m + J_m)), 2.0)));
} else if (t_2 <= -5e-91) {
tmp = t_3;
} else if (t_2 <= 5e-185) {
tmp = (sqrt(fma(((U_m / J_m) * (U_m / J_m)), 0.25, 1.0)) * J_m) * -2.0;
} else if (t_2 <= 2e+298) {
tmp = t_3;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(-2.0 * J_m) * t_0) t_2 = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_0)) ^ 2.0)))) t_3 = Float64(Float64(Float64(cos(Float64(0.5 * K)) * J_m) * -2.0) * sqrt(fma(Float64(Float64(U_m * U_m) / Float64(Float64(0.5 + Float64(0.5 * cos(Float64(2.0 * Float64(0.5 * K))))) * Float64(J_m * J_m))), 0.25, 1.0))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= -4e+159) tmp = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(J_m + J_m)) ^ 2.0)))); elseif (t_2 <= -5e-91) tmp = t_3; elseif (t_2 <= 5e-185) tmp = Float64(Float64(sqrt(fma(Float64(Float64(U_m / J_m) * Float64(U_m / J_m)), 0.25, 1.0)) * J_m) * -2.0); elseif (t_2 <= 2e+298) tmp = t_3; else tmp = U_m; end return Float64(J_s * tmp) end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, 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$95$m), $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$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision] * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision] * N[Sqrt[N[(N[(N[(U$95$m * U$95$m), $MachinePrecision] / N[(N[(0.5 + N[(0.5 * N[Cos[N[(2.0 * N[(0.5 * K), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(J$95$m * J$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, -4e+159], N[(t$95$1 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(J$95$m + J$95$m), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, -5e-91], t$95$3, If[LessEqual[t$95$2, 5e-185], N[(N[(N[Sqrt[N[(N[(N[(U$95$m / J$95$m), $MachinePrecision] * N[(U$95$m / J$95$m), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision] * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision], If[LessEqual[t$95$2, 2e+298], t$95$3, U$95$m]]]]]), $MachinePrecision]]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(-2 \cdot J\_m\right) \cdot t\_0\\
t_2 := t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_0}\right)}^{2}}\\
t_3 := \left(\left(\cos \left(0.5 \cdot K\right) \cdot J\_m\right) \cdot -2\right) \cdot \sqrt{\mathsf{fma}\left(\frac{U\_m \cdot U\_m}{\left(0.5 + 0.5 \cdot \cos \left(2 \cdot \left(0.5 \cdot K\right)\right)\right) \cdot \left(J\_m \cdot J\_m\right)}, 0.25, 1\right)}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq -4 \cdot 10^{+159}:\\
\;\;\;\;t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{J\_m + J\_m}\right)}^{2}}\\
\mathbf{elif}\;t\_2 \leq -5 \cdot 10^{-91}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t\_2 \leq 5 \cdot 10^{-185}:\\
\;\;\;\;\left(\sqrt{\mathsf{fma}\left(\frac{U\_m}{J\_m} \cdot \frac{U\_m}{J\_m}, 0.25, 1\right)} \cdot J\_m\right) \cdot -2\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+298}:\\
\;\;\;\;t\_3\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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))))) < -3.9999999999999997e159Initial program 99.8%
Taylor expanded in K around 0
count-2-revN/A
lower-+.f6494.3
Applied rewrites94.3%
if -3.9999999999999997e159 < (*.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.99999999999999997e-91 or 5.0000000000000003e-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))))) < 1.9999999999999999e298Initial program 99.8%
Taylor expanded in K around inf
associate-*r*N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
Applied rewrites89.1%
if -4.99999999999999997e-91 < (*.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.0000000000000003e-185Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6421.4
Applied rewrites21.4%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6477.4
Applied rewrites77.4%
if 1.9999999999999999e298 < (*.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 99.8%
Taylor expanded in U around -inf
Applied rewrites14.7%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1 (* (* -2.0 J_m) t_0))
(t_2 (* t_1 (sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_0)) 2.0))))))
(*
J_s
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 2e+298)
(* t_1 (sqrt (+ 1.0 (pow (/ U_m (+ J_m J_m)) 2.0))))
U_m)))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = (-2.0 * J_m) * t_0;
double t_2 = t_1 * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= 2e+298) {
tmp = t_1 * sqrt((1.0 + pow((U_m / (J_m + J_m)), 2.0)));
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = Math.abs(U);
J\_m = Math.abs(J);
J\_s = Math.copySign(1.0, J);
public static double code(double J_s, double J_m, double K, double U_m) {
double t_0 = Math.cos((K / 2.0));
double t_1 = (-2.0 * J_m) * t_0;
double t_2 = t_1 * Math.sqrt((1.0 + Math.pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double tmp;
if (t_2 <= -Double.POSITIVE_INFINITY) {
tmp = -U_m;
} else if (t_2 <= 2e+298) {
tmp = t_1 * Math.sqrt((1.0 + Math.pow((U_m / (J_m + J_m)), 2.0)));
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = math.fabs(U) J\_m = math.fabs(J) J\_s = math.copysign(1.0, J) def code(J_s, J_m, K, U_m): t_0 = math.cos((K / 2.0)) t_1 = (-2.0 * J_m) * t_0 t_2 = t_1 * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J_m) * t_0)), 2.0))) tmp = 0 if t_2 <= -math.inf: tmp = -U_m elif t_2 <= 2e+298: tmp = t_1 * math.sqrt((1.0 + math.pow((U_m / (J_m + J_m)), 2.0))) else: tmp = U_m return J_s * tmp
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(-2.0 * J_m) * t_0) t_2 = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= 2e+298) tmp = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(J_m + J_m)) ^ 2.0)))); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = abs(U); J\_m = abs(J); J\_s = sign(J) * abs(1.0); function tmp_2 = code(J_s, J_m, K, U_m) t_0 = cos((K / 2.0)); t_1 = (-2.0 * J_m) * t_0; t_2 = t_1 * sqrt((1.0 + ((U_m / ((2.0 * J_m) * t_0)) ^ 2.0))); tmp = 0.0; if (t_2 <= -Inf) tmp = -U_m; elseif (t_2 <= 2e+298) tmp = t_1 * sqrt((1.0 + ((U_m / (J_m + J_m)) ^ 2.0))); else tmp = U_m; end tmp_2 = J_s * tmp; end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, 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$95$m), $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$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, 2e+298], N[(t$95$1 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(J$95$m + J$95$m), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], U$95$m]]), $MachinePrecision]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(-2 \cdot J\_m\right) \cdot t\_0\\
t_2 := t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_0}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+298}:\\
\;\;\;\;t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{J\_m + J\_m}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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.9999999999999999e298Initial program 99.8%
Taylor expanded in K around 0
count-2-revN/A
lower-+.f6487.2
Applied rewrites87.2%
if 1.9999999999999999e298 < (*.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 12.1%
Taylor expanded in U around -inf
Applied rewrites94.5%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1 (* (* -2.0 J_m) t_0))
(t_2 (* t_1 (sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_0)) 2.0))))))
(*
J_s
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 2e+298) (* t_1 (cosh (asinh (/ U_m (+ J_m J_m))))) U_m)))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = (-2.0 * J_m) * t_0;
double t_2 = t_1 * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= 2e+298) {
tmp = t_1 * cosh(asinh((U_m / (J_m + J_m))));
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = math.fabs(U) J\_m = math.fabs(J) J\_s = math.copysign(1.0, J) def code(J_s, J_m, K, U_m): t_0 = math.cos((K / 2.0)) t_1 = (-2.0 * J_m) * t_0 t_2 = t_1 * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J_m) * t_0)), 2.0))) tmp = 0 if t_2 <= -math.inf: tmp = -U_m elif t_2 <= 2e+298: tmp = t_1 * math.cosh(math.asinh((U_m / (J_m + J_m)))) else: tmp = U_m return J_s * tmp
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(-2.0 * J_m) * t_0) t_2 = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= 2e+298) tmp = Float64(t_1 * cosh(asinh(Float64(U_m / Float64(J_m + J_m))))); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = abs(U); J\_m = abs(J); J\_s = sign(J) * abs(1.0); function tmp_2 = code(J_s, J_m, K, U_m) t_0 = cos((K / 2.0)); t_1 = (-2.0 * J_m) * t_0; t_2 = t_1 * sqrt((1.0 + ((U_m / ((2.0 * J_m) * t_0)) ^ 2.0))); tmp = 0.0; if (t_2 <= -Inf) tmp = -U_m; elseif (t_2 <= 2e+298) tmp = t_1 * cosh(asinh((U_m / (J_m + J_m)))); else tmp = U_m; end tmp_2 = J_s * tmp; end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, 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$95$m), $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$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, 2e+298], N[(t$95$1 * N[Cosh[N[ArcSinh[N[(U$95$m / N[(J$95$m + J$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], U$95$m]]), $MachinePrecision]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(-2 \cdot J\_m\right) \cdot t\_0\\
t_2 := t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_0}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+298}:\\
\;\;\;\;t\_1 \cdot \cosh \sinh^{-1} \left(\frac{U\_m}{J\_m + J\_m}\right)\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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.9999999999999999e298Initial program 99.8%
lift-sqrt.f64N/A
lift-+.f64N/A
lift-pow.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-cos.f64N/A
+-commutativeN/A
unpow2N/A
cosh-asinh-revN/A
lower-cosh.f64N/A
lower-asinh.f64N/A
Applied rewrites97.6%
Taylor expanded in K around 0
count-2-revN/A
count-2-revN/A
lift-+.f6486.2
Applied rewrites86.2%
if 1.9999999999999999e298 < (*.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 12.1%
Taylor expanded in U around -inf
Applied rewrites94.5%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1 (* (* -2.0 J_m) t_0))
(t_2 (* t_1 (sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_0)) 2.0))))))
(*
J_s
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 -2e+227)
(fma (* U_m (/ U_m J_m)) -0.25 (* (* (cos (* 0.5 K)) J_m) -2.0))
(if (<= t_2 5e-185)
(* (* (sqrt (fma (* (/ U_m J_m) (/ U_m J_m)) 0.25 1.0)) J_m) -2.0)
(if (<= t_2 2e+298)
(* t_1 (sqrt (fma (/ (* U_m U_m) (* J_m J_m)) 0.25 1.0)))
U_m)))))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = (-2.0 * J_m) * t_0;
double t_2 = t_1 * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= -2e+227) {
tmp = fma((U_m * (U_m / J_m)), -0.25, ((cos((0.5 * K)) * J_m) * -2.0));
} else if (t_2 <= 5e-185) {
tmp = (sqrt(fma(((U_m / J_m) * (U_m / J_m)), 0.25, 1.0)) * J_m) * -2.0;
} else if (t_2 <= 2e+298) {
tmp = t_1 * sqrt(fma(((U_m * U_m) / (J_m * J_m)), 0.25, 1.0));
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(-2.0 * J_m) * t_0) t_2 = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= -2e+227) tmp = fma(Float64(U_m * Float64(U_m / J_m)), -0.25, Float64(Float64(cos(Float64(0.5 * K)) * J_m) * -2.0)); elseif (t_2 <= 5e-185) tmp = Float64(Float64(sqrt(fma(Float64(Float64(U_m / J_m) * Float64(U_m / J_m)), 0.25, 1.0)) * J_m) * -2.0); elseif (t_2 <= 2e+298) tmp = Float64(t_1 * sqrt(fma(Float64(Float64(U_m * U_m) / Float64(J_m * J_m)), 0.25, 1.0))); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, 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$95$m), $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$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, -2e+227], N[(N[(U$95$m * N[(U$95$m / J$95$m), $MachinePrecision]), $MachinePrecision] * -0.25 + N[(N[(N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision] * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 5e-185], N[(N[(N[Sqrt[N[(N[(N[(U$95$m / J$95$m), $MachinePrecision] * N[(U$95$m / J$95$m), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision] * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision], If[LessEqual[t$95$2, 2e+298], N[(t$95$1 * N[Sqrt[N[(N[(N[(U$95$m * U$95$m), $MachinePrecision] / N[(J$95$m * J$95$m), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], U$95$m]]]]), $MachinePrecision]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(-2 \cdot J\_m\right) \cdot t\_0\\
t_2 := t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_0}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq -2 \cdot 10^{+227}:\\
\;\;\;\;\mathsf{fma}\left(U\_m \cdot \frac{U\_m}{J\_m}, -0.25, \left(\cos \left(0.5 \cdot K\right) \cdot J\_m\right) \cdot -2\right)\\
\mathbf{elif}\;t\_2 \leq 5 \cdot 10^{-185}:\\
\;\;\;\;\left(\sqrt{\mathsf{fma}\left(\frac{U\_m}{J\_m} \cdot \frac{U\_m}{J\_m}, 0.25, 1\right)} \cdot J\_m\right) \cdot -2\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+298}:\\
\;\;\;\;t\_1 \cdot \sqrt{\mathsf{fma}\left(\frac{U\_m \cdot U\_m}{J\_m \cdot J\_m}, 0.25, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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))))) < -2.0000000000000002e227Initial program 99.8%
Taylor expanded in U around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites65.6%
lift-*.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lower-/.f6479.4
Applied rewrites79.4%
Taylor expanded in K around 0
Applied rewrites79.3%
if -2.0000000000000002e227 < (*.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.0000000000000003e-185Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6459.7
Applied rewrites59.7%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6479.4
Applied rewrites79.4%
if 5.0000000000000003e-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))))) < 1.9999999999999999e298Initial program 99.7%
Taylor expanded in K around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6467.6
Applied rewrites67.6%
if 1.9999999999999999e298 < (*.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 12.1%
Taylor expanded in U around -inf
Applied rewrites94.5%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (* (* (cos (* 0.5 K)) J_m) -2.0))
(t_1 (cos (/ K 2.0)))
(t_2
(*
(* (* -2.0 J_m) t_1)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_1)) 2.0))))))
(*
J_s
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 -2e+227)
(fma (* U_m (/ U_m J_m)) -0.25 t_0)
(if (<= t_2 -1e-268)
(* (* (sqrt (fma (* (/ U_m J_m) (/ U_m J_m)) 0.25 1.0)) J_m) -2.0)
(if (<= t_2 2e+298) t_0 U_m)))))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = (cos((0.5 * K)) * J_m) * -2.0;
double t_1 = cos((K / 2.0));
double t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_1)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= -2e+227) {
tmp = fma((U_m * (U_m / J_m)), -0.25, t_0);
} else if (t_2 <= -1e-268) {
tmp = (sqrt(fma(((U_m / J_m) * (U_m / J_m)), 0.25, 1.0)) * J_m) * -2.0;
} else if (t_2 <= 2e+298) {
tmp = t_0;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = Float64(Float64(cos(Float64(0.5 * K)) * J_m) * -2.0) t_1 = cos(Float64(K / 2.0)) t_2 = Float64(Float64(Float64(-2.0 * J_m) * t_1) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_1)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= -2e+227) tmp = fma(Float64(U_m * Float64(U_m / J_m)), -0.25, t_0); elseif (t_2 <= -1e-268) tmp = Float64(Float64(sqrt(fma(Float64(Float64(U_m / J_m) * Float64(U_m / J_m)), 0.25, 1.0)) * J_m) * -2.0); elseif (t_2 <= 2e+298) tmp = t_0; else tmp = U_m; end return Float64(J_s * tmp) end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, K_, U$95$m_] := Block[{t$95$0 = N[(N[(N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision] * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, -2e+227], N[(N[(U$95$m * N[(U$95$m / J$95$m), $MachinePrecision]), $MachinePrecision] * -0.25 + t$95$0), $MachinePrecision], If[LessEqual[t$95$2, -1e-268], N[(N[(N[Sqrt[N[(N[(N[(U$95$m / J$95$m), $MachinePrecision] * N[(U$95$m / J$95$m), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision] * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision], If[LessEqual[t$95$2, 2e+298], t$95$0, U$95$m]]]]), $MachinePrecision]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \left(\cos \left(0.5 \cdot K\right) \cdot J\_m\right) \cdot -2\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := \left(\left(-2 \cdot J\_m\right) \cdot t\_1\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_1}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq -2 \cdot 10^{+227}:\\
\;\;\;\;\mathsf{fma}\left(U\_m \cdot \frac{U\_m}{J\_m}, -0.25, t\_0\right)\\
\mathbf{elif}\;t\_2 \leq -1 \cdot 10^{-268}:\\
\;\;\;\;\left(\sqrt{\mathsf{fma}\left(\frac{U\_m}{J\_m} \cdot \frac{U\_m}{J\_m}, 0.25, 1\right)} \cdot J\_m\right) \cdot -2\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+298}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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))))) < -2.0000000000000002e227Initial program 99.8%
Taylor expanded in U around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites65.6%
lift-*.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lower-/.f6479.4
Applied rewrites79.4%
Taylor expanded in K around 0
Applied rewrites79.3%
if -2.0000000000000002e227 < (*.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.99999999999999958e-269Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6461.6
Applied rewrites61.6%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6481.2
Applied rewrites81.2%
if -9.99999999999999958e-269 < (*.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.9999999999999999e298Initial program 99.7%
Taylor expanded in J around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f6469.5
Applied rewrites69.5%
if 1.9999999999999999e298 < (*.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 12.1%
Taylor expanded in U around -inf
Applied rewrites94.5%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J_m) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_0)) 2.0)))))
(t_2 (* (* (cos (* 0.5 K)) J_m) -2.0)))
(*
J_s
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -2e+227)
t_2
(if (<= t_1 -1e-268)
(* (* (sqrt (fma (* (/ U_m J_m) (/ U_m J_m)) 0.25 1.0)) J_m) -2.0)
(if (<= t_1 2e+298) t_2 U_m)))))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J_m) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double t_2 = (cos((0.5 * K)) * J_m) * -2.0;
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -2e+227) {
tmp = t_2;
} else if (t_1 <= -1e-268) {
tmp = (sqrt(fma(((U_m / J_m) * (U_m / J_m)), 0.25, 1.0)) * J_m) * -2.0;
} else if (t_1 <= 2e+298) {
tmp = t_2;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J_m) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_0)) ^ 2.0)))) t_2 = Float64(Float64(cos(Float64(0.5 * K)) * J_m) * -2.0) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -2e+227) tmp = t_2; elseif (t_1 <= -1e-268) tmp = Float64(Float64(sqrt(fma(Float64(Float64(U_m / J_m) * Float64(U_m / J_m)), 0.25, 1.0)) * J_m) * -2.0); elseif (t_1 <= 2e+298) tmp = t_2; else tmp = U_m; end return Float64(J_s * tmp) end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, 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$95$m), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[Cos[N[(0.5 * K), $MachinePrecision]], $MachinePrecision] * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -2e+227], t$95$2, If[LessEqual[t$95$1, -1e-268], N[(N[(N[Sqrt[N[(N[(N[(U$95$m / J$95$m), $MachinePrecision] * N[(U$95$m / J$95$m), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision] * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision], If[LessEqual[t$95$1, 2e+298], t$95$2, U$95$m]]]]), $MachinePrecision]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\_m\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_0}\right)}^{2}}\\
t_2 := \left(\cos \left(0.5 \cdot K\right) \cdot J\_m\right) \cdot -2\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{+227}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-268}:\\
\;\;\;\;\left(\sqrt{\mathsf{fma}\left(\frac{U\_m}{J\_m} \cdot \frac{U\_m}{J\_m}, 0.25, 1\right)} \cdot J\_m\right) \cdot -2\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{+298}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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))))) < -2.0000000000000002e227 or -9.99999999999999958e-269 < (*.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.9999999999999999e298Initial program 99.7%
Taylor expanded in J around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f6473.1
Applied rewrites73.1%
if -2.0000000000000002e227 < (*.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.99999999999999958e-269Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6461.6
Applied rewrites61.6%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6481.2
Applied rewrites81.2%
if 1.9999999999999999e298 < (*.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 99.7%
Taylor expanded in U around -inf
Applied rewrites20.7%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J_m) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_0)) 2.0))))))
(*
J_s
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -1e-268)
(* (* (sqrt (fma (* (/ U_m J_m) (/ U_m J_m)) 0.25 1.0)) J_m) -2.0)
U_m)))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J_m) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -1e-268) {
tmp = (sqrt(fma(((U_m / J_m) * (U_m / J_m)), 0.25, 1.0)) * J_m) * -2.0;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J_m) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -1e-268) tmp = Float64(Float64(sqrt(fma(Float64(Float64(U_m / J_m) * Float64(U_m / J_m)), 0.25, 1.0)) * J_m) * -2.0); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, 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$95$m), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -1e-268], N[(N[(N[Sqrt[N[(N[(N[(U$95$m / J$95$m), $MachinePrecision] * N[(U$95$m / J$95$m), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision] * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision], U$95$m]]), $MachinePrecision]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\_m\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_0}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-268}:\\
\;\;\;\;\left(\sqrt{\mathsf{fma}\left(\frac{U\_m}{J\_m} \cdot \frac{U\_m}{J\_m}, 0.25, 1\right)} \cdot J\_m\right) \cdot -2\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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))))) < -9.99999999999999958e-269Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6458.2
Applied rewrites58.2%
lift-*.f64N/A
lift-*.f64N/A
lift-/.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lower-/.f6480.6
Applied rewrites80.6%
if -9.99999999999999958e-269 < (*.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 74.9%
Taylor expanded in U around -inf
Applied rewrites50.9%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J_m) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_0)) 2.0))))))
(*
J_s
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -1e+225)
(fma (* (/ U_m 1.0) (/ U_m J_m)) -0.25 (* (* 1.0 J_m) -2.0))
(if (<= t_1 -1e+218)
(* (* (* 0.5 (/ U_m J_m)) J_m) -2.0)
(if (<= t_1 -2e-103)
(* (* (sqrt (fma (/ (* U_m U_m) (* J_m J_m)) 0.25 1.0)) J_m) -2.0)
(if (<= t_1 -1e-268)
(* (fma 0.5 U_m (/ (* J_m J_m) U_m)) -2.0)
U_m))))))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J_m) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -1e+225) {
tmp = fma(((U_m / 1.0) * (U_m / J_m)), -0.25, ((1.0 * J_m) * -2.0));
} else if (t_1 <= -1e+218) {
tmp = ((0.5 * (U_m / J_m)) * J_m) * -2.0;
} else if (t_1 <= -2e-103) {
tmp = (sqrt(fma(((U_m * U_m) / (J_m * J_m)), 0.25, 1.0)) * J_m) * -2.0;
} else if (t_1 <= -1e-268) {
tmp = fma(0.5, U_m, ((J_m * J_m) / U_m)) * -2.0;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J_m) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -1e+225) tmp = fma(Float64(Float64(U_m / 1.0) * Float64(U_m / J_m)), -0.25, Float64(Float64(1.0 * J_m) * -2.0)); elseif (t_1 <= -1e+218) tmp = Float64(Float64(Float64(0.5 * Float64(U_m / J_m)) * J_m) * -2.0); elseif (t_1 <= -2e-103) tmp = Float64(Float64(sqrt(fma(Float64(Float64(U_m * U_m) / Float64(J_m * J_m)), 0.25, 1.0)) * J_m) * -2.0); elseif (t_1 <= -1e-268) tmp = Float64(fma(0.5, U_m, Float64(Float64(J_m * J_m) / U_m)) * -2.0); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, 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$95$m), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -1e+225], N[(N[(N[(U$95$m / 1.0), $MachinePrecision] * N[(U$95$m / J$95$m), $MachinePrecision]), $MachinePrecision] * -0.25 + N[(N[(1.0 * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -1e+218], N[(N[(N[(0.5 * N[(U$95$m / J$95$m), $MachinePrecision]), $MachinePrecision] * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision], If[LessEqual[t$95$1, -2e-103], N[(N[(N[Sqrt[N[(N[(N[(U$95$m * U$95$m), $MachinePrecision] / N[(J$95$m * J$95$m), $MachinePrecision]), $MachinePrecision] * 0.25 + 1.0), $MachinePrecision]], $MachinePrecision] * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision], If[LessEqual[t$95$1, -1e-268], N[(N[(0.5 * U$95$m + N[(N[(J$95$m * J$95$m), $MachinePrecision] / U$95$m), $MachinePrecision]), $MachinePrecision] * -2.0), $MachinePrecision], U$95$m]]]]]), $MachinePrecision]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\_m\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_0}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{+225}:\\
\;\;\;\;\mathsf{fma}\left(\frac{U\_m}{1} \cdot \frac{U\_m}{J\_m}, -0.25, \left(1 \cdot J\_m\right) \cdot -2\right)\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{+218}:\\
\;\;\;\;\left(\left(0.5 \cdot \frac{U\_m}{J\_m}\right) \cdot J\_m\right) \cdot -2\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-103}:\\
\;\;\;\;\left(\sqrt{\mathsf{fma}\left(\frac{U\_m \cdot U\_m}{J\_m \cdot J\_m}, 0.25, 1\right)} \cdot J\_m\right) \cdot -2\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-268}:\\
\;\;\;\;\mathsf{fma}\left(0.5, U\_m, \frac{J\_m \cdot J\_m}{U\_m}\right) \cdot -2\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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))))) < -9.99999999999999928e224Initial program 99.8%
Taylor expanded in U around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites65.6%
lift-*.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lower-/.f6479.6
Applied rewrites79.6%
Taylor expanded in K around 0
Applied rewrites79.4%
Taylor expanded in K around 0
Applied rewrites58.8%
if -9.99999999999999928e224 < (*.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.00000000000000008e218Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6450.4
Applied rewrites50.4%
Taylor expanded in J around 0
lower-*.f64N/A
lower-/.f6422.2
Applied rewrites22.2%
if -1.00000000000000008e218 < (*.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.99999999999999992e-103Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.4
Applied rewrites71.4%
if -1.99999999999999992e-103 < (*.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.99999999999999958e-269Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6420.8
Applied rewrites20.8%
Taylor expanded in J around 0
lower-fma.f64N/A
lower-/.f64N/A
pow2N/A
lift-*.f6453.2
Applied rewrites53.2%
if -9.99999999999999958e-269 < (*.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 74.9%
Taylor expanded in U around -inf
Applied rewrites50.9%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (* (fma 0.5 U_m (/ (* J_m J_m) U_m)) -2.0))
(t_1 (cos (/ K 2.0)))
(t_2
(*
(* (* -2.0 J_m) t_1)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_1)) 2.0))))))
(*
J_s
(if (<= t_2 (- INFINITY))
(- U_m)
(if (<= t_2 -4e+117)
(* J_m -2.0)
(if (<= t_2 -5e+71)
t_0
(if (<= t_2 -2e-103)
(fma (/ (* U_m U_m) J_m) -0.25 (* J_m -2.0))
(if (<= t_2 -1e-268) t_0 U_m))))))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = fma(0.5, U_m, ((J_m * J_m) / U_m)) * -2.0;
double t_1 = cos((K / 2.0));
double t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_1)), 2.0)));
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_2 <= -4e+117) {
tmp = J_m * -2.0;
} else if (t_2 <= -5e+71) {
tmp = t_0;
} else if (t_2 <= -2e-103) {
tmp = fma(((U_m * U_m) / J_m), -0.25, (J_m * -2.0));
} else if (t_2 <= -1e-268) {
tmp = t_0;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = Float64(fma(0.5, U_m, Float64(Float64(J_m * J_m) / U_m)) * -2.0) t_1 = cos(Float64(K / 2.0)) t_2 = Float64(Float64(Float64(-2.0 * J_m) * t_1) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_1)) ^ 2.0)))) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_2 <= -4e+117) tmp = Float64(J_m * -2.0); elseif (t_2 <= -5e+71) tmp = t_0; elseif (t_2 <= -2e-103) tmp = fma(Float64(Float64(U_m * U_m) / J_m), -0.25, Float64(J_m * -2.0)); elseif (t_2 <= -1e-268) tmp = t_0; else tmp = U_m; end return Float64(J_s * tmp) end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, K_, U$95$m_] := Block[{t$95$0 = N[(N[(0.5 * U$95$m + N[(N[(J$95$m * J$95$m), $MachinePrecision] / U$95$m), $MachinePrecision]), $MachinePrecision] * -2.0), $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, (-Infinity)], (-U$95$m), If[LessEqual[t$95$2, -4e+117], N[(J$95$m * -2.0), $MachinePrecision], If[LessEqual[t$95$2, -5e+71], t$95$0, If[LessEqual[t$95$2, -2e-103], N[(N[(N[(U$95$m * U$95$m), $MachinePrecision] / J$95$m), $MachinePrecision] * -0.25 + N[(J$95$m * -2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, -1e-268], t$95$0, U$95$m]]]]]), $MachinePrecision]]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(0.5, U\_m, \frac{J\_m \cdot J\_m}{U\_m}\right) \cdot -2\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := \left(\left(-2 \cdot J\_m\right) \cdot t\_1\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_1}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq -4 \cdot 10^{+117}:\\
\;\;\;\;J\_m \cdot -2\\
\mathbf{elif}\;t\_2 \leq -5 \cdot 10^{+71}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;t\_2 \leq -2 \cdot 10^{-103}:\\
\;\;\;\;\mathsf{fma}\left(\frac{U\_m \cdot U\_m}{J\_m}, -0.25, J\_m \cdot -2\right)\\
\mathbf{elif}\;t\_2 \leq -1 \cdot 10^{-268}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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))))) < -4.0000000000000002e117Initial program 99.8%
Taylor expanded in J around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f6479.1
Applied rewrites79.1%
Taylor expanded in K around 0
Applied rewrites58.4%
if -4.0000000000000002e117 < (*.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.99999999999999972e71 or -1.99999999999999992e-103 < (*.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.99999999999999958e-269Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6443.9
Applied rewrites43.9%
Taylor expanded in J around 0
lower-fma.f64N/A
lower-/.f64N/A
pow2N/A
lift-*.f6444.3
Applied rewrites44.3%
if -4.99999999999999972e71 < (*.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.99999999999999992e-103Initial program 99.8%
Taylor expanded in U around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites68.8%
Taylor expanded in K around 0
Applied rewrites68.7%
Taylor expanded in K around 0
Applied rewrites51.0%
if -9.99999999999999958e-269 < (*.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 99.8%
Taylor expanded in U around -inf
Applied rewrites2.9%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J_m) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_0)) 2.0))))))
(*
J_s
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -2e-103)
(fma (* (/ U_m 1.0) (/ U_m J_m)) -0.25 (* (* 1.0 J_m) -2.0))
(if (<= t_1 -1e-268)
(* (fma 0.5 U_m (/ (* J_m J_m) U_m)) -2.0)
U_m))))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J_m) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -2e-103) {
tmp = fma(((U_m / 1.0) * (U_m / J_m)), -0.25, ((1.0 * J_m) * -2.0));
} else if (t_1 <= -1e-268) {
tmp = fma(0.5, U_m, ((J_m * J_m) / U_m)) * -2.0;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J_m) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -2e-103) tmp = fma(Float64(Float64(U_m / 1.0) * Float64(U_m / J_m)), -0.25, Float64(Float64(1.0 * J_m) * -2.0)); elseif (t_1 <= -1e-268) tmp = Float64(fma(0.5, U_m, Float64(Float64(J_m * J_m) / U_m)) * -2.0); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, 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$95$m), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -2e-103], N[(N[(N[(U$95$m / 1.0), $MachinePrecision] * N[(U$95$m / J$95$m), $MachinePrecision]), $MachinePrecision] * -0.25 + N[(N[(1.0 * J$95$m), $MachinePrecision] * -2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, -1e-268], N[(N[(0.5 * U$95$m + N[(N[(J$95$m * J$95$m), $MachinePrecision] / U$95$m), $MachinePrecision]), $MachinePrecision] * -2.0), $MachinePrecision], U$95$m]]]), $MachinePrecision]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\_m\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_0}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-103}:\\
\;\;\;\;\mathsf{fma}\left(\frac{U\_m}{1} \cdot \frac{U\_m}{J\_m}, -0.25, \left(1 \cdot J\_m\right) \cdot -2\right)\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-268}:\\
\;\;\;\;\mathsf{fma}\left(0.5, U\_m, \frac{J\_m \cdot J\_m}{U\_m}\right) \cdot -2\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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.99999999999999992e-103Initial program 99.8%
Taylor expanded in U around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
Applied rewrites70.7%
lift-*.f64N/A
lift-/.f64N/A
lift-*.f64N/A
lift-*.f64N/A
lift-cos.f64N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
lower-/.f6475.5
Applied rewrites75.5%
Taylor expanded in K around 0
Applied rewrites75.4%
Taylor expanded in K around 0
Applied rewrites55.7%
if -1.99999999999999992e-103 < (*.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.99999999999999958e-269Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6420.8
Applied rewrites20.8%
Taylor expanded in J around 0
lower-fma.f64N/A
lower-/.f64N/A
pow2N/A
lift-*.f6453.2
Applied rewrites53.2%
if -9.99999999999999958e-269 < (*.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 74.9%
Taylor expanded in U around -inf
Applied rewrites50.9%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J_m) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_0)) 2.0))))))
(*
J_s
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -2e-103)
(* J_m -2.0)
(if (<= t_1 -1e-268)
(* (fma 0.5 U_m (/ (* J_m J_m) U_m)) -2.0)
U_m))))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J_m) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -2e-103) {
tmp = J_m * -2.0;
} else if (t_1 <= -1e-268) {
tmp = fma(0.5, U_m, ((J_m * J_m) / U_m)) * -2.0;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J_m) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -2e-103) tmp = Float64(J_m * -2.0); elseif (t_1 <= -1e-268) tmp = Float64(fma(0.5, U_m, Float64(Float64(J_m * J_m) / U_m)) * -2.0); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, 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$95$m), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -2e-103], N[(J$95$m * -2.0), $MachinePrecision], If[LessEqual[t$95$1, -1e-268], N[(N[(0.5 * U$95$m + N[(N[(J$95$m * J$95$m), $MachinePrecision] / U$95$m), $MachinePrecision]), $MachinePrecision] * -2.0), $MachinePrecision], U$95$m]]]), $MachinePrecision]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\_m\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_0}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-103}:\\
\;\;\;\;J\_m \cdot -2\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-268}:\\
\;\;\;\;\mathsf{fma}\left(0.5, U\_m, \frac{J\_m \cdot J\_m}{U\_m}\right) \cdot -2\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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 5.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.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.99999999999999992e-103Initial program 99.8%
Taylor expanded in J around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f6475.1
Applied rewrites75.1%
Taylor expanded in K around 0
Applied rewrites55.5%
if -1.99999999999999992e-103 < (*.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.99999999999999958e-269Initial program 99.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6420.8
Applied rewrites20.8%
Taylor expanded in J around 0
lower-fma.f64N/A
lower-/.f64N/A
pow2N/A
lift-*.f6453.2
Applied rewrites53.2%
if -9.99999999999999958e-269 < (*.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 74.9%
Taylor expanded in U around -inf
Applied rewrites50.9%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1
(*
(* (* -2.0 J_m) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_0)) 2.0))))))
(*
J_s
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 -2e-103) (* J_m -2.0) (if (<= t_1 -1e-268) (- U_m) U_m))))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((K / 2.0));
double t_1 = ((-2.0 * J_m) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= -2e-103) {
tmp = J_m * -2.0;
} else if (t_1 <= -1e-268) {
tmp = -U_m;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = Math.abs(U);
J\_m = Math.abs(J);
J\_s = Math.copySign(1.0, J);
public static double code(double J_s, double J_m, double K, double U_m) {
double t_0 = Math.cos((K / 2.0));
double t_1 = ((-2.0 * J_m) * t_0) * Math.sqrt((1.0 + Math.pow((U_m / ((2.0 * J_m) * t_0)), 2.0)));
double tmp;
if (t_1 <= -Double.POSITIVE_INFINITY) {
tmp = -U_m;
} else if (t_1 <= -2e-103) {
tmp = J_m * -2.0;
} else if (t_1 <= -1e-268) {
tmp = -U_m;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = math.fabs(U) J\_m = math.fabs(J) J\_s = math.copysign(1.0, J) def code(J_s, J_m, K, U_m): t_0 = math.cos((K / 2.0)) t_1 = ((-2.0 * J_m) * t_0) * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J_m) * t_0)), 2.0))) tmp = 0 if t_1 <= -math.inf: tmp = -U_m elif t_1 <= -2e-103: tmp = J_m * -2.0 elif t_1 <= -1e-268: tmp = -U_m else: tmp = U_m return J_s * tmp
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(K / 2.0)) t_1 = Float64(Float64(Float64(-2.0 * J_m) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_0)) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= -2e-103) tmp = Float64(J_m * -2.0); elseif (t_1 <= -1e-268) tmp = Float64(-U_m); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = abs(U); J\_m = abs(J); J\_s = sign(J) * abs(1.0); function tmp_2 = code(J_s, J_m, K, U_m) t_0 = cos((K / 2.0)); t_1 = ((-2.0 * J_m) * t_0) * sqrt((1.0 + ((U_m / ((2.0 * J_m) * t_0)) ^ 2.0))); tmp = 0.0; if (t_1 <= -Inf) tmp = -U_m; elseif (t_1 <= -2e-103) tmp = J_m * -2.0; elseif (t_1 <= -1e-268) tmp = -U_m; else tmp = U_m; end tmp_2 = J_s * tmp; end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, 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$95$m), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$1, (-Infinity)], (-U$95$m), If[LessEqual[t$95$1, -2e-103], N[(J$95$m * -2.0), $MachinePrecision], If[LessEqual[t$95$1, -1e-268], (-U$95$m), U$95$m]]]), $MachinePrecision]]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
t_1 := \left(\left(-2 \cdot J\_m\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_0}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-103}:\\
\;\;\;\;J\_m \cdot -2\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-268}:\\
\;\;\;\;-U\_m\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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.99999999999999992e-103 < (*.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.99999999999999958e-269Initial program 29.8%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6487.9
Applied rewrites87.9%
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.99999999999999992e-103Initial program 99.8%
Taylor expanded in J around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f6475.1
Applied rewrites75.1%
Taylor expanded in K around 0
Applied rewrites55.5%
if -9.99999999999999958e-269 < (*.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 29.8%
Taylor expanded in U around -inf
Applied rewrites1.5%
U_m = (fabs.f64 U)
J\_m = (fabs.f64 J)
J\_s = (copysign.f64 #s(literal 1 binary64) J)
(FPCore (J_s J_m K U_m)
:precision binary64
(let* ((t_0 (cos (/ K 2.0))))
(*
J_s
(if (<=
(*
(* (* -2.0 J_m) t_0)
(sqrt (+ 1.0 (pow (/ U_m (* (* 2.0 J_m) t_0)) 2.0))))
-1e-268)
(- U_m)
U_m))))U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
double t_0 = cos((K / 2.0));
double tmp;
if ((((-2.0 * J_m) * t_0) * sqrt((1.0 + pow((U_m / ((2.0 * J_m) * t_0)), 2.0)))) <= -1e-268) {
tmp = -U_m;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = private
J\_m = private
J\_s = private
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(j_s, j_m, k, u_m)
use fmin_fmax_functions
real(8), intent (in) :: j_s
real(8), intent (in) :: j_m
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_m) * t_0) * sqrt((1.0d0 + ((u_m / ((2.0d0 * j_m) * t_0)) ** 2.0d0)))) <= (-1d-268)) then
tmp = -u_m
else
tmp = u_m
end if
code = j_s * tmp
end function
U_m = Math.abs(U);
J\_m = Math.abs(J);
J\_s = Math.copySign(1.0, J);
public static double code(double J_s, double J_m, double K, double U_m) {
double t_0 = Math.cos((K / 2.0));
double tmp;
if ((((-2.0 * J_m) * t_0) * Math.sqrt((1.0 + Math.pow((U_m / ((2.0 * J_m) * t_0)), 2.0)))) <= -1e-268) {
tmp = -U_m;
} else {
tmp = U_m;
}
return J_s * tmp;
}
U_m = math.fabs(U) J\_m = math.fabs(J) J\_s = math.copysign(1.0, J) def code(J_s, J_m, K, U_m): t_0 = math.cos((K / 2.0)) tmp = 0 if (((-2.0 * J_m) * t_0) * math.sqrt((1.0 + math.pow((U_m / ((2.0 * J_m) * t_0)), 2.0)))) <= -1e-268: tmp = -U_m else: tmp = U_m return J_s * tmp
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) t_0 = cos(Float64(K / 2.0)) tmp = 0.0 if (Float64(Float64(Float64(-2.0 * J_m) * t_0) * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(2.0 * J_m) * t_0)) ^ 2.0)))) <= -1e-268) tmp = Float64(-U_m); else tmp = U_m; end return Float64(J_s * tmp) end
U_m = abs(U); J\_m = abs(J); J\_s = sign(J) * abs(1.0); function tmp_2 = code(J_s, J_m, K, U_m) t_0 = cos((K / 2.0)); tmp = 0.0; if ((((-2.0 * J_m) * t_0) * sqrt((1.0 + ((U_m / ((2.0 * J_m) * t_0)) ^ 2.0)))) <= -1e-268) tmp = -U_m; else tmp = U_m; end tmp_2 = J_s * tmp; end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, K_, U$95$m_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, N[(J$95$s * If[LessEqual[N[(N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(2.0 * J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], -1e-268], (-U$95$m), U$95$m]), $MachinePrecision]]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;\left(\left(-2 \cdot J\_m\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(2 \cdot J\_m\right) \cdot t\_0}\right)}^{2}} \leq -1 \cdot 10^{-268}:\\
\;\;\;\;-U\_m\\
\mathbf{else}:\\
\;\;\;\;U\_m\\
\end{array}
\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))))) < -9.99999999999999958e-269Initial program 72.5%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6452.0
Applied rewrites52.0%
if -9.99999999999999958e-269 < (*.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 74.9%
Taylor expanded in U around -inf
Applied rewrites50.9%
U_m = (fabs.f64 U) J\_m = (fabs.f64 J) J\_s = (copysign.f64 #s(literal 1 binary64) J) (FPCore (J_s J_m K U_m) :precision binary64 (* J_s U_m))
U_m = fabs(U);
J\_m = fabs(J);
J\_s = copysign(1.0, J);
double code(double J_s, double J_m, double K, double U_m) {
return J_s * U_m;
}
U_m = private
J\_m = private
J\_s = private
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(j_s, j_m, k, u_m)
use fmin_fmax_functions
real(8), intent (in) :: j_s
real(8), intent (in) :: j_m
real(8), intent (in) :: k
real(8), intent (in) :: u_m
code = j_s * u_m
end function
U_m = Math.abs(U);
J\_m = Math.abs(J);
J\_s = Math.copySign(1.0, J);
public static double code(double J_s, double J_m, double K, double U_m) {
return J_s * U_m;
}
U_m = math.fabs(U) J\_m = math.fabs(J) J\_s = math.copysign(1.0, J) def code(J_s, J_m, K, U_m): return J_s * U_m
U_m = abs(U) J\_m = abs(J) J\_s = copysign(1.0, J) function code(J_s, J_m, K, U_m) return Float64(J_s * U_m) end
U_m = abs(U); J\_m = abs(J); J\_s = sign(J) * abs(1.0); function tmp = code(J_s, J_m, K, U_m) tmp = J_s * U_m; end
U_m = N[Abs[U], $MachinePrecision]
J\_m = N[Abs[J], $MachinePrecision]
J\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[J]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[J$95$s_, J$95$m_, K_, U$95$m_] := N[(J$95$s * U$95$m), $MachinePrecision]
\begin{array}{l}
U_m = \left|U\right|
\\
J\_m = \left|J\right|
\\
J\_s = \mathsf{copysign}\left(1, J\right)
\\
J\_s \cdot U\_m
\end{array}
Initial program 73.1%
Taylor expanded in U around -inf
Applied rewrites14.0%
herbie shell --seed 2025120
(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)))))