
(FPCore (J K U) :precision binary64 (let* ((t_0 (cos (/ K 2.0)))) (* (* (* -2.0 J) t_0) (sqrt (+ 1.0 (pow (/ U (* (* 2.0 J) t_0)) 2.0))))))
double code(double J, double K, double U) {
double t_0 = cos((K / 2.0));
return ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U / ((2.0 * J) * t_0)), 2.0)));
}
real(8) function code(j, k, u)
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8), intent (in) :: u
real(8) :: t_0
t_0 = cos((k / 2.0d0))
code = (((-2.0d0) * j) * t_0) * sqrt((1.0d0 + ((u / ((2.0d0 * j) * t_0)) ** 2.0d0)))
end function
public static double code(double J, double K, double U) {
double t_0 = Math.cos((K / 2.0));
return ((-2.0 * J) * t_0) * Math.sqrt((1.0 + Math.pow((U / ((2.0 * J) * t_0)), 2.0)));
}
def code(J, K, U): t_0 = math.cos((K / 2.0)) return ((-2.0 * J) * t_0) * math.sqrt((1.0 + math.pow((U / ((2.0 * J) * t_0)), 2.0)))
function code(J, K, U) t_0 = cos(Float64(K / 2.0)) return Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) end
function tmp = code(J, K, U) t_0 = cos((K / 2.0)); tmp = ((-2.0 * J) * t_0) * sqrt((1.0 + ((U / ((2.0 * J) * t_0)) ^ 2.0))); end
code[J_, K_, U_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
\left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 16 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (J K U) :precision binary64 (let* ((t_0 (cos (/ K 2.0)))) (* (* (* -2.0 J) t_0) (sqrt (+ 1.0 (pow (/ U (* (* 2.0 J) t_0)) 2.0))))))
double code(double J, double K, double U) {
double t_0 = cos((K / 2.0));
return ((-2.0 * J) * t_0) * sqrt((1.0 + pow((U / ((2.0 * J) * t_0)), 2.0)));
}
real(8) function code(j, k, u)
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8), intent (in) :: u
real(8) :: t_0
t_0 = cos((k / 2.0d0))
code = (((-2.0d0) * j) * t_0) * sqrt((1.0d0 + ((u / ((2.0d0 * j) * t_0)) ** 2.0d0)))
end function
public static double code(double J, double K, double U) {
double t_0 = Math.cos((K / 2.0));
return ((-2.0 * J) * t_0) * Math.sqrt((1.0 + Math.pow((U / ((2.0 * J) * t_0)), 2.0)));
}
def code(J, K, U): t_0 = math.cos((K / 2.0)) return ((-2.0 * J) * t_0) * math.sqrt((1.0 + math.pow((U / ((2.0 * J) * t_0)), 2.0)))
function code(J, K, U) t_0 = cos(Float64(K / 2.0)) return Float64(Float64(Float64(-2.0 * J) * t_0) * sqrt(Float64(1.0 + (Float64(U / Float64(Float64(2.0 * J) * t_0)) ^ 2.0)))) end
function tmp = code(J, K, U) t_0 = cos((K / 2.0)); tmp = ((-2.0 * J) * t_0) * sqrt((1.0 + ((U / ((2.0 * J) * t_0)) ^ 2.0))); end
code[J_, K_, U_] := Block[{t$95$0 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, N[(N[(N[(-2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U / N[(N[(2.0 * J), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos \left(\frac{K}{2}\right)\\
\left(\left(-2 \cdot J\right) \cdot t\_0\right) \cdot \sqrt{1 + {\left(\frac{U}{\left(2 \cdot J\right) \cdot t\_0}\right)}^{2}}
\end{array}
\end{array}
U_m = (fabs.f64 U)
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 (* t_0 (* J_m 2.0))) 2.0))))))
(*
J_s
(if (<= t_1 (- INFINITY))
(- U_m)
(if (<= t_1 2e+286)
t_1
(*
(- U_m)
(fma
-2.0
(* (pow (cos (* K 0.5)) 2.0) (/ (* J_m J_m) (* U_m U_m)))
-1.0)))))))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 / (t_0 * (J_m * 2.0))), 2.0)));
double tmp;
if (t_1 <= -((double) INFINITY)) {
tmp = -U_m;
} else if (t_1 <= 2e+286) {
tmp = t_1;
} else {
tmp = -U_m * fma(-2.0, (pow(cos((K * 0.5)), 2.0) * ((J_m * J_m) / (U_m * U_m))), -1.0);
}
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(t_0 * Float64(J_m * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_1 <= Float64(-Inf)) tmp = Float64(-U_m); elseif (t_1 <= 2e+286) tmp = t_1; else tmp = Float64(Float64(-U_m) * fma(-2.0, Float64((cos(Float64(K * 0.5)) ^ 2.0) * Float64(Float64(J_m * J_m) / Float64(U_m * U_m))), -1.0)); 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[(t$95$0 * N[(J$95$m * 2.0), $MachinePrecision]), $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+286], t$95$1, N[((-U$95$m) * N[(-2.0 * N[(N[Power[N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision], 2.0], $MachinePrecision] * N[(N[(J$95$m * J$95$m), $MachinePrecision] / N[(U$95$m * U$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]]]), $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}{t\_0 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -\infty:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{+286}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\left(-U\_m\right) \cdot \mathsf{fma}\left(-2, {\cos \left(K \cdot 0.5\right)}^{2} \cdot \frac{J\_m \cdot J\_m}{U\_m \cdot U\_m}, -1\right)\\
\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.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
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.00000000000000007e286Initial program 99.9%
if 2.00000000000000007e286 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 11.5%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6448.7
Applied rewrites48.7%
Final simplification85.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 0.5)))
(t_1 (cos (/ K 2.0)))
(t_2 (* (* -2.0 J_m) t_1))
(t_3 (* t_2 (sqrt (+ 1.0 (pow (/ U_m (* t_1 (* J_m 2.0))) 2.0))))))
(*
J_s
(if (<= t_3 -2e+295)
(- U_m)
(if (<= t_3 2e-63)
(* t_2 (sqrt (+ 1.0 (pow (/ U_m (* (* J_m 2.0) 1.0)) 2.0))))
(if (<= t_3 5e+265)
(*
J_m
(*
(* -2.0 t_0)
(sqrt
(fma
U_m
(/ U_m (* (fma 0.5 (cos K) 0.5) (* (* J_m J_m) 4.0)))
1.0))))
(/ (* (- U_m) t_0) (sqrt (fma (cos K) 0.5 0.5)))))))))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 * 0.5));
double t_1 = cos((K / 2.0));
double t_2 = (-2.0 * J_m) * t_1;
double t_3 = t_2 * sqrt((1.0 + pow((U_m / (t_1 * (J_m * 2.0))), 2.0)));
double tmp;
if (t_3 <= -2e+295) {
tmp = -U_m;
} else if (t_3 <= 2e-63) {
tmp = t_2 * sqrt((1.0 + pow((U_m / ((J_m * 2.0) * 1.0)), 2.0)));
} else if (t_3 <= 5e+265) {
tmp = J_m * ((-2.0 * t_0) * sqrt(fma(U_m, (U_m / (fma(0.5, cos(K), 0.5) * ((J_m * J_m) * 4.0))), 1.0)));
} else {
tmp = (-U_m * t_0) / sqrt(fma(cos(K), 0.5, 0.5));
}
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 * 0.5)) t_1 = cos(Float64(K / 2.0)) t_2 = Float64(Float64(-2.0 * J_m) * t_1) t_3 = Float64(t_2 * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_1 * Float64(J_m * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_3 <= -2e+295) tmp = Float64(-U_m); elseif (t_3 <= 2e-63) tmp = Float64(t_2 * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(J_m * 2.0) * 1.0)) ^ 2.0)))); elseif (t_3 <= 5e+265) tmp = Float64(J_m * Float64(Float64(-2.0 * t_0) * sqrt(fma(U_m, Float64(U_m / Float64(fma(0.5, cos(K), 0.5) * Float64(Float64(J_m * J_m) * 4.0))), 1.0)))); else tmp = Float64(Float64(Float64(-U_m) * t_0) / sqrt(fma(cos(K), 0.5, 0.5))); 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 * 0.5), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(t$95$2 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$1 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$3, -2e+295], (-U$95$m), If[LessEqual[t$95$3, 2e-63], N[(t$95$2 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(J$95$m * 2.0), $MachinePrecision] * 1.0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 5e+265], N[(J$95$m * N[(N[(-2.0 * t$95$0), $MachinePrecision] * N[Sqrt[N[(U$95$m * N[(U$95$m / N[(N[(0.5 * N[Cos[K], $MachinePrecision] + 0.5), $MachinePrecision] * N[(N[(J$95$m * J$95$m), $MachinePrecision] * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[((-U$95$m) * t$95$0), $MachinePrecision] / N[Sqrt[N[(N[Cos[K], $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]), $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(K \cdot 0.5\right)\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := \left(-2 \cdot J\_m\right) \cdot t\_1\\
t_3 := t\_2 \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_1 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_3 \leq -2 \cdot 10^{+295}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_3 \leq 2 \cdot 10^{-63}:\\
\;\;\;\;t\_2 \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(J\_m \cdot 2\right) \cdot 1}\right)}^{2}}\\
\mathbf{elif}\;t\_3 \leq 5 \cdot 10^{+265}:\\
\;\;\;\;J\_m \cdot \left(\left(-2 \cdot t\_0\right) \cdot \sqrt{\mathsf{fma}\left(U\_m, \frac{U\_m}{\mathsf{fma}\left(0.5, \cos K, 0.5\right) \cdot \left(\left(J\_m \cdot J\_m\right) \cdot 4\right)}, 1\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-U\_m\right) \cdot t\_0}{\sqrt{\mathsf{fma}\left(\cos K, 0.5, 0.5\right)}}\\
\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))))) < -2e295Initial program 5.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
if -2e295 < (*.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.00000000000000013e-63Initial program 99.9%
Taylor expanded in K around 0
Applied rewrites84.4%
if 2.00000000000000013e-63 < (*.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.0000000000000002e265Initial program 99.8%
Applied rewrites94.5%
Taylor expanded in K around inf
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6494.5
Applied rewrites94.5%
if 5.0000000000000002e265 < (*.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 18.9%
lift-pow.f64N/A
unpow2N/A
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
associate-*l/N/A
lift-/.f64N/A
associate-*r/N/A
associate-/r*N/A
lower-/.f64N/A
Applied rewrites13.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f6450.2
Applied rewrites50.2%
Applied rewrites50.2%
Final simplification77.6%
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 0.5)))
(t_1 (cos (/ K 2.0)))
(t_2 (* (* -2.0 J_m) t_1))
(t_3 (* t_2 (sqrt (+ 1.0 (pow (/ U_m (* t_1 (* J_m 2.0))) 2.0))))))
(*
J_s
(if (<= t_3 -2e+295)
(- U_m)
(if (<= t_3 2e-63)
(* t_2 (sqrt (+ 1.0 (pow (/ U_m (* (* J_m 2.0) 1.0)) 2.0))))
(if (<= t_3 4e+231)
(*
J_m
(*
(* -2.0 t_0)
(sqrt
(fma 0.5 (/ (* U_m U_m) (* J_m (fma J_m (cos K) J_m))) 1.0))))
(/ (* (- U_m) t_0) (sqrt (fma (cos K) 0.5 0.5)))))))))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 * 0.5));
double t_1 = cos((K / 2.0));
double t_2 = (-2.0 * J_m) * t_1;
double t_3 = t_2 * sqrt((1.0 + pow((U_m / (t_1 * (J_m * 2.0))), 2.0)));
double tmp;
if (t_3 <= -2e+295) {
tmp = -U_m;
} else if (t_3 <= 2e-63) {
tmp = t_2 * sqrt((1.0 + pow((U_m / ((J_m * 2.0) * 1.0)), 2.0)));
} else if (t_3 <= 4e+231) {
tmp = J_m * ((-2.0 * t_0) * sqrt(fma(0.5, ((U_m * U_m) / (J_m * fma(J_m, cos(K), J_m))), 1.0)));
} else {
tmp = (-U_m * t_0) / sqrt(fma(cos(K), 0.5, 0.5));
}
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 * 0.5)) t_1 = cos(Float64(K / 2.0)) t_2 = Float64(Float64(-2.0 * J_m) * t_1) t_3 = Float64(t_2 * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_1 * Float64(J_m * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_3 <= -2e+295) tmp = Float64(-U_m); elseif (t_3 <= 2e-63) tmp = Float64(t_2 * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(J_m * 2.0) * 1.0)) ^ 2.0)))); elseif (t_3 <= 4e+231) tmp = Float64(J_m * Float64(Float64(-2.0 * t_0) * sqrt(fma(0.5, Float64(Float64(U_m * U_m) / Float64(J_m * fma(J_m, cos(K), J_m))), 1.0)))); else tmp = Float64(Float64(Float64(-U_m) * t_0) / sqrt(fma(cos(K), 0.5, 0.5))); 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 * 0.5), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(t$95$2 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$1 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$3, -2e+295], (-U$95$m), If[LessEqual[t$95$3, 2e-63], N[(t$95$2 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(J$95$m * 2.0), $MachinePrecision] * 1.0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 4e+231], N[(J$95$m * N[(N[(-2.0 * t$95$0), $MachinePrecision] * N[Sqrt[N[(0.5 * N[(N[(U$95$m * U$95$m), $MachinePrecision] / N[(J$95$m * N[(J$95$m * N[Cos[K], $MachinePrecision] + J$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[((-U$95$m) * t$95$0), $MachinePrecision] / N[Sqrt[N[(N[Cos[K], $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]), $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(K \cdot 0.5\right)\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := \left(-2 \cdot J\_m\right) \cdot t\_1\\
t_3 := t\_2 \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_1 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_3 \leq -2 \cdot 10^{+295}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_3 \leq 2 \cdot 10^{-63}:\\
\;\;\;\;t\_2 \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(J\_m \cdot 2\right) \cdot 1}\right)}^{2}}\\
\mathbf{elif}\;t\_3 \leq 4 \cdot 10^{+231}:\\
\;\;\;\;J\_m \cdot \left(\left(-2 \cdot t\_0\right) \cdot \sqrt{\mathsf{fma}\left(0.5, \frac{U\_m \cdot U\_m}{J\_m \cdot \mathsf{fma}\left(J\_m, \cos K, J\_m\right)}, 1\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-U\_m\right) \cdot t\_0}{\sqrt{\mathsf{fma}\left(\cos K, 0.5, 0.5\right)}}\\
\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))))) < -2e295Initial program 5.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
if -2e295 < (*.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.00000000000000013e-63Initial program 99.9%
Taylor expanded in K around 0
Applied rewrites84.4%
if 2.00000000000000013e-63 < (*.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.0000000000000002e231Initial program 99.8%
Applied rewrites98.1%
lift-fma.f64N/A
Applied rewrites99.8%
Taylor expanded in U around 0
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f6491.8
Applied rewrites91.8%
if 4.0000000000000002e231 < (*.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 37.8%
lift-pow.f64N/A
unpow2N/A
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
associate-*l/N/A
lift-/.f64N/A
associate-*r/N/A
associate-/r*N/A
lower-/.f64N/A
Applied rewrites31.8%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f6443.2
Applied rewrites43.2%
Applied rewrites43.2%
Final simplification73.8%
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 (* t_0 (* J_m 2.0))) 2.0))))))
(*
J_s
(if (<= t_2 -2e+295)
(- U_m)
(if (<= t_2 -5e+59)
(* t_1 (fma U_m (* U_m (/ 0.125 (* J_m J_m))) 1.0))
(if (<= t_2 -2e-168)
(* (sqrt (fma 0.25 (/ (/ (* U_m U_m) J_m) J_m) 1.0)) (* -2.0 J_m))
(* (* -2.0 J_m) (cos (* K 0.5)))))))))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 / (t_0 * (J_m * 2.0))), 2.0)));
double tmp;
if (t_2 <= -2e+295) {
tmp = -U_m;
} else if (t_2 <= -5e+59) {
tmp = t_1 * fma(U_m, (U_m * (0.125 / (J_m * J_m))), 1.0);
} else if (t_2 <= -2e-168) {
tmp = sqrt(fma(0.25, (((U_m * U_m) / J_m) / J_m), 1.0)) * (-2.0 * J_m);
} else {
tmp = (-2.0 * J_m) * cos((K * 0.5));
}
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(t_0 * Float64(J_m * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_2 <= -2e+295) tmp = Float64(-U_m); elseif (t_2 <= -5e+59) tmp = Float64(t_1 * fma(U_m, Float64(U_m * Float64(0.125 / Float64(J_m * J_m))), 1.0)); elseif (t_2 <= -2e-168) tmp = Float64(sqrt(fma(0.25, Float64(Float64(Float64(U_m * U_m) / J_m) / J_m), 1.0)) * Float64(-2.0 * J_m)); else tmp = Float64(Float64(-2.0 * J_m) * cos(Float64(K * 0.5))); 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[(t$95$0 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, -2e+295], (-U$95$m), If[LessEqual[t$95$2, -5e+59], N[(t$95$1 * N[(U$95$m * N[(U$95$m * N[(0.125 / N[(J$95$m * J$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, -2e-168], N[(N[Sqrt[N[(0.25 * N[(N[(N[(U$95$m * U$95$m), $MachinePrecision] / J$95$m), $MachinePrecision] / J$95$m), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision] * N[(-2.0 * J$95$m), $MachinePrecision]), $MachinePrecision], N[(N[(-2.0 * J$95$m), $MachinePrecision] * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]), $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}{t\_0 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -2 \cdot 10^{+295}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq -5 \cdot 10^{+59}:\\
\;\;\;\;t\_1 \cdot \mathsf{fma}\left(U\_m, U\_m \cdot \frac{0.125}{J\_m \cdot J\_m}, 1\right)\\
\mathbf{elif}\;t\_2 \leq -2 \cdot 10^{-168}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.25, \frac{\frac{U\_m \cdot U\_m}{J\_m}}{J\_m}, 1\right)} \cdot \left(-2 \cdot J\_m\right)\\
\mathbf{else}:\\
\;\;\;\;\left(-2 \cdot J\_m\right) \cdot \cos \left(K \cdot 0.5\right)\\
\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))))) < -2e295Initial program 5.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
if -2e295 < (*.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.9999999999999997e59Initial program 99.9%
Taylor expanded in K around 0
lower-sqrt.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6474.0
Applied rewrites74.0%
Taylor expanded in U around 0
Applied rewrites76.5%
if -4.9999999999999997e59 < (*.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.0000000000000001e-168Initial program 100.0%
Taylor expanded in K around 0
lower-sqrt.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6473.2
Applied rewrites73.2%
Applied rewrites78.3%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f6460.9
Applied rewrites60.9%
if -2.0000000000000001e-168 < (*.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 78.1%
Taylor expanded in J around inf
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-*.f6451.3
Applied rewrites51.3%
Final simplification57.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 (* t_0 (* J_m 2.0))) 2.0)))))
(t_2 (* (* -2.0 J_m) (cos (* K 0.5)))))
(*
J_s
(if (<= t_1 -2e+295)
(- U_m)
(if (<= t_1 -5e+59)
t_2
(if (<= t_1 -2e-168)
(* (sqrt (fma 0.25 (/ (/ (* U_m U_m) J_m) J_m) 1.0)) (* -2.0 J_m))
t_2))))))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 / (t_0 * (J_m * 2.0))), 2.0)));
double t_2 = (-2.0 * J_m) * cos((K * 0.5));
double tmp;
if (t_1 <= -2e+295) {
tmp = -U_m;
} else if (t_1 <= -5e+59) {
tmp = t_2;
} else if (t_1 <= -2e-168) {
tmp = sqrt(fma(0.25, (((U_m * U_m) / J_m) / J_m), 1.0)) * (-2.0 * J_m);
} else {
tmp = t_2;
}
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(t_0 * Float64(J_m * 2.0))) ^ 2.0)))) t_2 = Float64(Float64(-2.0 * J_m) * cos(Float64(K * 0.5))) tmp = 0.0 if (t_1 <= -2e+295) tmp = Float64(-U_m); elseif (t_1 <= -5e+59) tmp = t_2; elseif (t_1 <= -2e-168) tmp = Float64(sqrt(fma(0.25, Float64(Float64(Float64(U_m * U_m) / J_m) / J_m), 1.0)) * Float64(-2.0 * J_m)); else tmp = t_2; 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[(t$95$0 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(-2.0 * J$95$m), $MachinePrecision] * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$1, -2e+295], (-U$95$m), If[LessEqual[t$95$1, -5e+59], t$95$2, If[LessEqual[t$95$1, -2e-168], N[(N[Sqrt[N[(0.25 * N[(N[(N[(U$95$m * U$95$m), $MachinePrecision] / J$95$m), $MachinePrecision] / J$95$m), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision] * N[(-2.0 * J$95$m), $MachinePrecision]), $MachinePrecision], t$95$2]]]), $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}{t\_0 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
t_2 := \left(-2 \cdot J\_m\right) \cdot \cos \left(K \cdot 0.5\right)\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{+295}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -5 \cdot 10^{+59}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-168}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.25, \frac{\frac{U\_m \cdot U\_m}{J\_m}}{J\_m}, 1\right)} \cdot \left(-2 \cdot J\_m\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\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))))) < -2e295Initial program 5.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
if -2e295 < (*.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.9999999999999997e59 or -2.0000000000000001e-168 < (*.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 83.8%
Taylor expanded in J around inf
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-*.f6457.9
Applied rewrites57.9%
if -4.9999999999999997e59 < (*.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.0000000000000001e-168Initial program 100.0%
Taylor expanded in K around 0
lower-sqrt.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6473.2
Applied rewrites73.2%
Applied rewrites78.3%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f6460.9
Applied rewrites60.9%
Final simplification57.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 0.5)))
(t_1 (cos (/ K 2.0)))
(t_2
(*
(* (* -2.0 J_m) t_1)
(sqrt (+ 1.0 (pow (/ U_m (* t_1 (* J_m 2.0))) 2.0))))))
(*
J_s
(if (<= t_2 -2e+295)
(- U_m)
(if (<= t_2 2e+286)
(*
J_m
(*
(sqrt (fma (/ U_m (fma (cos K) J_m J_m)) (/ U_m (* J_m 2.0)) 1.0))
(* -2.0 t_0)))
(*
(- U_m)
(fma -2.0 (* (pow t_0 2.0) (/ (* J_m J_m) (* U_m U_m))) -1.0)))))))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 * 0.5));
double t_1 = cos((K / 2.0));
double t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + pow((U_m / (t_1 * (J_m * 2.0))), 2.0)));
double tmp;
if (t_2 <= -2e+295) {
tmp = -U_m;
} else if (t_2 <= 2e+286) {
tmp = J_m * (sqrt(fma((U_m / fma(cos(K), J_m, J_m)), (U_m / (J_m * 2.0)), 1.0)) * (-2.0 * t_0));
} else {
tmp = -U_m * fma(-2.0, (pow(t_0, 2.0) * ((J_m * J_m) / (U_m * U_m))), -1.0);
}
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 * 0.5)) 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(t_1 * Float64(J_m * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_2 <= -2e+295) tmp = Float64(-U_m); elseif (t_2 <= 2e+286) tmp = Float64(J_m * Float64(sqrt(fma(Float64(U_m / fma(cos(K), J_m, J_m)), Float64(U_m / Float64(J_m * 2.0)), 1.0)) * Float64(-2.0 * t_0))); else tmp = Float64(Float64(-U_m) * fma(-2.0, Float64((t_0 ^ 2.0) * Float64(Float64(J_m * J_m) / Float64(U_m * U_m))), -1.0)); 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 * 0.5), $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[(t$95$1 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, -2e+295], (-U$95$m), If[LessEqual[t$95$2, 2e+286], N[(J$95$m * N[(N[Sqrt[N[(N[(U$95$m / N[(N[Cos[K], $MachinePrecision] * J$95$m + J$95$m), $MachinePrecision]), $MachinePrecision] * N[(U$95$m / N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision] * N[(-2.0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[((-U$95$m) * N[(-2.0 * N[(N[Power[t$95$0, 2.0], $MachinePrecision] * N[(N[(J$95$m * J$95$m), $MachinePrecision] / N[(U$95$m * U$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]]]), $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(K \cdot 0.5\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}{t\_1 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -2 \cdot 10^{+295}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 2 \cdot 10^{+286}:\\
\;\;\;\;J\_m \cdot \left(\sqrt{\mathsf{fma}\left(\frac{U\_m}{\mathsf{fma}\left(\cos K, J\_m, J\_m\right)}, \frac{U\_m}{J\_m \cdot 2}, 1\right)} \cdot \left(-2 \cdot t\_0\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(-U\_m\right) \cdot \mathsf{fma}\left(-2, {t\_0}^{2} \cdot \frac{J\_m \cdot J\_m}{U\_m \cdot U\_m}, -1\right)\\
\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))))) < -2e295Initial program 5.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
if -2e295 < (*.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.00000000000000007e286Initial program 99.9%
Applied rewrites80.6%
lift-fma.f64N/A
Applied rewrites99.7%
lift-*.f64N/A
*-rgt-identity99.7
lift-*.f64N/A
*-rgt-identity99.7
lift-*.f64N/A
*-rgt-identity99.7
Applied rewrites99.7%
if 2.00000000000000007e286 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 11.5%
Taylor expanded in U around -inf
associate-*r*N/A
lower-*.f64N/A
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-pow.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6448.7
Applied rewrites48.7%
Final simplification85.8%
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 (* t_0 (* J_m 2.0))) 2.0))))))
(*
J_s
(if (<= t_2 -2e+295)
(- U_m)
(if (<= t_2 5e+265)
(* t_1 (sqrt (+ 1.0 (pow (/ U_m (* (* J_m 2.0) 1.0)) 2.0))))
(/ (* (- U_m) (cos (* K 0.5))) (sqrt (fma (cos K) 0.5 0.5))))))))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 / (t_0 * (J_m * 2.0))), 2.0)));
double tmp;
if (t_2 <= -2e+295) {
tmp = -U_m;
} else if (t_2 <= 5e+265) {
tmp = t_1 * sqrt((1.0 + pow((U_m / ((J_m * 2.0) * 1.0)), 2.0)));
} else {
tmp = (-U_m * cos((K * 0.5))) / sqrt(fma(cos(K), 0.5, 0.5));
}
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(t_0 * Float64(J_m * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_2 <= -2e+295) tmp = Float64(-U_m); elseif (t_2 <= 5e+265) tmp = Float64(t_1 * sqrt(Float64(1.0 + (Float64(U_m / Float64(Float64(J_m * 2.0) * 1.0)) ^ 2.0)))); else tmp = Float64(Float64(Float64(-U_m) * cos(Float64(K * 0.5))) / sqrt(fma(cos(K), 0.5, 0.5))); 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[(t$95$0 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, -2e+295], (-U$95$m), If[LessEqual[t$95$2, 5e+265], N[(t$95$1 * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(N[(J$95$m * 2.0), $MachinePrecision] * 1.0), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(N[((-U$95$m) * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sqrt[N[(N[Cos[K], $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]), $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}{t\_0 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -2 \cdot 10^{+295}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq 5 \cdot 10^{+265}:\\
\;\;\;\;t\_1 \cdot \sqrt{1 + {\left(\frac{U\_m}{\left(J\_m \cdot 2\right) \cdot 1}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-U\_m\right) \cdot \cos \left(K \cdot 0.5\right)}{\sqrt{\mathsf{fma}\left(\cos K, 0.5, 0.5\right)}}\\
\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))))) < -2e295Initial program 5.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
if -2e295 < (*.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.0000000000000002e265Initial program 99.9%
Taylor expanded in K around 0
Applied rewrites84.1%
if 5.0000000000000002e265 < (*.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 18.9%
lift-pow.f64N/A
unpow2N/A
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
associate-*l/N/A
lift-/.f64N/A
associate-*r/N/A
associate-/r*N/A
lower-/.f64N/A
Applied rewrites13.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f6450.2
Applied rewrites50.2%
Applied rewrites50.2%
Final simplification74.6%
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 0.5)))
(t_1 (/ U_m (* J_m 2.0)))
(t_2 (cos (/ K 2.0)))
(t_3
(*
(* (* -2.0 J_m) t_2)
(sqrt (+ 1.0 (pow (/ U_m (* t_2 (* J_m 2.0))) 2.0))))))
(*
J_s
(if (<= t_3 -2e+295)
(- U_m)
(if (<= t_3 5e+265)
(* J_m (* (* -2.0 t_0) (sqrt (fma t_1 t_1 1.0))))
(/ (* (- U_m) t_0) (sqrt (fma (cos K) 0.5 0.5))))))))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 * 0.5));
double t_1 = U_m / (J_m * 2.0);
double t_2 = cos((K / 2.0));
double t_3 = ((-2.0 * J_m) * t_2) * sqrt((1.0 + pow((U_m / (t_2 * (J_m * 2.0))), 2.0)));
double tmp;
if (t_3 <= -2e+295) {
tmp = -U_m;
} else if (t_3 <= 5e+265) {
tmp = J_m * ((-2.0 * t_0) * sqrt(fma(t_1, t_1, 1.0)));
} else {
tmp = (-U_m * t_0) / sqrt(fma(cos(K), 0.5, 0.5));
}
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 * 0.5)) t_1 = Float64(U_m / Float64(J_m * 2.0)) t_2 = cos(Float64(K / 2.0)) t_3 = Float64(Float64(Float64(-2.0 * J_m) * t_2) * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_2 * Float64(J_m * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_3 <= -2e+295) tmp = Float64(-U_m); elseif (t_3 <= 5e+265) tmp = Float64(J_m * Float64(Float64(-2.0 * t_0) * sqrt(fma(t_1, t_1, 1.0)))); else tmp = Float64(Float64(Float64(-U_m) * t_0) / sqrt(fma(cos(K), 0.5, 0.5))); 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 * 0.5), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(U$95$m / N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$2), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$2 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$3, -2e+295], (-U$95$m), If[LessEqual[t$95$3, 5e+265], N[(J$95$m * N[(N[(-2.0 * t$95$0), $MachinePrecision] * N[Sqrt[N[(t$95$1 * t$95$1 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[((-U$95$m) * t$95$0), $MachinePrecision] / N[Sqrt[N[(N[Cos[K], $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]), $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(K \cdot 0.5\right)\\
t_1 := \frac{U\_m}{J\_m \cdot 2}\\
t_2 := \cos \left(\frac{K}{2}\right)\\
t_3 := \left(\left(-2 \cdot J\_m\right) \cdot t\_2\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_2 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_3 \leq -2 \cdot 10^{+295}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_3 \leq 5 \cdot 10^{+265}:\\
\;\;\;\;J\_m \cdot \left(\left(-2 \cdot t\_0\right) \cdot \sqrt{\mathsf{fma}\left(t\_1, t\_1, 1\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(-U\_m\right) \cdot t\_0}{\sqrt{\mathsf{fma}\left(\cos K, 0.5, 0.5\right)}}\\
\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))))) < -2e295Initial program 5.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
if -2e295 < (*.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.0000000000000002e265Initial program 99.9%
Applied rewrites81.9%
lift-fma.f64N/A
Applied rewrites99.7%
Taylor expanded in K around 0
lower-*.f6484.1
Applied rewrites84.1%
if 5.0000000000000002e265 < (*.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 18.9%
lift-pow.f64N/A
unpow2N/A
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
associate-*l/N/A
lift-/.f64N/A
associate-*r/N/A
associate-/r*N/A
lower-/.f64N/A
Applied rewrites13.7%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f6450.2
Applied rewrites50.2%
Applied rewrites50.2%
Final simplification74.6%
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 0.5)))
(t_1 (cos (/ K 2.0)))
(t_2
(*
(* (* -2.0 J_m) t_1)
(sqrt (+ 1.0 (pow (/ U_m (* t_1 (* J_m 2.0))) 2.0))))))
(*
J_s
(if (<= t_2 -2e+295)
(- U_m)
(if (<= t_2 -1e-157)
(*
J_m
(* (* -2.0 t_0) (sqrt (fma U_m (/ U_m (* (* J_m J_m) 4.0)) 1.0))))
(* (* -2.0 J_m) t_0))))))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 * 0.5));
double t_1 = cos((K / 2.0));
double t_2 = ((-2.0 * J_m) * t_1) * sqrt((1.0 + pow((U_m / (t_1 * (J_m * 2.0))), 2.0)));
double tmp;
if (t_2 <= -2e+295) {
tmp = -U_m;
} else if (t_2 <= -1e-157) {
tmp = J_m * ((-2.0 * t_0) * sqrt(fma(U_m, (U_m / ((J_m * J_m) * 4.0)), 1.0)));
} else {
tmp = (-2.0 * J_m) * t_0;
}
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 * 0.5)) 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(t_1 * Float64(J_m * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_2 <= -2e+295) tmp = Float64(-U_m); elseif (t_2 <= -1e-157) tmp = Float64(J_m * Float64(Float64(-2.0 * t_0) * sqrt(fma(U_m, Float64(U_m / Float64(Float64(J_m * J_m) * 4.0)), 1.0)))); else tmp = Float64(Float64(-2.0 * J_m) * t_0); 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 * 0.5), $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[(t$95$1 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$2, -2e+295], (-U$95$m), If[LessEqual[t$95$2, -1e-157], N[(J$95$m * N[(N[(-2.0 * t$95$0), $MachinePrecision] * N[Sqrt[N[(U$95$m * N[(U$95$m / N[(N[(J$95$m * J$95$m), $MachinePrecision] * 4.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(-2.0 * J$95$m), $MachinePrecision] * t$95$0), $MachinePrecision]]]), $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(K \cdot 0.5\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}{t\_1 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_2 \leq -2 \cdot 10^{+295}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_2 \leq -1 \cdot 10^{-157}:\\
\;\;\;\;J\_m \cdot \left(\left(-2 \cdot t\_0\right) \cdot \sqrt{\mathsf{fma}\left(U\_m, \frac{U\_m}{\left(J\_m \cdot J\_m\right) \cdot 4}, 1\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(-2 \cdot J\_m\right) \cdot t\_0\\
\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))))) < -2e295Initial program 5.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
if -2e295 < (*.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.99999999999999943e-158Initial program 99.9%
Applied rewrites88.8%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6480.7
Applied rewrites80.7%
if -9.99999999999999943e-158 < (*.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 78.4%
Taylor expanded in J around inf
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-*.f6450.7
Applied rewrites50.7%
Final simplification60.6%
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 (* t_0 (* J_m 2.0))) 2.0))))))
(*
J_s
(if (<= t_1 -2e+270)
(- U_m)
(if (<= t_1 -5e-102)
(* (sqrt (fma 0.25 (/ (/ (* U_m U_m) J_m) J_m) 1.0)) (* -2.0 J_m))
(*
(fma K (* K (* J_m 0.25)) (* -2.0 J_m))
(sqrt (fma 0.25 (/ (* U_m (/ U_m J_m)) J_m) 1.0))))))))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 / (t_0 * (J_m * 2.0))), 2.0)));
double tmp;
if (t_1 <= -2e+270) {
tmp = -U_m;
} else if (t_1 <= -5e-102) {
tmp = sqrt(fma(0.25, (((U_m * U_m) / J_m) / J_m), 1.0)) * (-2.0 * J_m);
} else {
tmp = fma(K, (K * (J_m * 0.25)), (-2.0 * J_m)) * sqrt(fma(0.25, ((U_m * (U_m / J_m)) / J_m), 1.0));
}
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(t_0 * Float64(J_m * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_1 <= -2e+270) tmp = Float64(-U_m); elseif (t_1 <= -5e-102) tmp = Float64(sqrt(fma(0.25, Float64(Float64(Float64(U_m * U_m) / J_m) / J_m), 1.0)) * Float64(-2.0 * J_m)); else tmp = Float64(fma(K, Float64(K * Float64(J_m * 0.25)), Float64(-2.0 * J_m)) * sqrt(fma(0.25, Float64(Float64(U_m * Float64(U_m / J_m)) / J_m), 1.0))); 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[(t$95$0 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$1, -2e+270], (-U$95$m), If[LessEqual[t$95$1, -5e-102], N[(N[Sqrt[N[(0.25 * N[(N[(N[(U$95$m * U$95$m), $MachinePrecision] / J$95$m), $MachinePrecision] / J$95$m), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision] * N[(-2.0 * J$95$m), $MachinePrecision]), $MachinePrecision], N[(N[(K * N[(K * N[(J$95$m * 0.25), $MachinePrecision]), $MachinePrecision] + N[(-2.0 * J$95$m), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(0.25 * N[(N[(U$95$m * N[(U$95$m / J$95$m), $MachinePrecision]), $MachinePrecision] / J$95$m), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]), $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}{t\_0 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{+270}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -5 \cdot 10^{-102}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.25, \frac{\frac{U\_m \cdot U\_m}{J\_m}}{J\_m}, 1\right)} \cdot \left(-2 \cdot J\_m\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(K, K \cdot \left(J\_m \cdot 0.25\right), -2 \cdot J\_m\right) \cdot \sqrt{\mathsf{fma}\left(0.25, \frac{U\_m \cdot \frac{U\_m}{J\_m}}{J\_m}, 1\right)}\\
\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))))) < -2.0000000000000001e270Initial program 15.1%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6453.7
Applied rewrites53.7%
if -2.0000000000000001e270 < (*.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.00000000000000026e-102Initial program 99.9%
Taylor expanded in K around 0
lower-sqrt.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6477.5
Applied rewrites77.5%
Applied rewrites78.7%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f6453.6
Applied rewrites53.6%
if -5.00000000000000026e-102 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) Initial program 78.9%
Taylor expanded in K around 0
lower-sqrt.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6447.3
Applied rewrites47.3%
Taylor expanded in K around 0
associate-*r*N/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6426.9
Applied rewrites26.9%
Applied rewrites35.5%
Final simplification43.8%
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 (* t_0 (* J_m 2.0))) 2.0))))))
(*
J_s
(if (<= t_1 -2e+270)
(- U_m)
(if (<= t_1 -2e-240)
(* (sqrt (fma 0.25 (/ (/ (* U_m U_m) J_m) J_m) 1.0)) (* -2.0 J_m))
(* (fma -0.125 (* U_m (* K K)) U_m) (- 1.0)))))))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 / (t_0 * (J_m * 2.0))), 2.0)));
double tmp;
if (t_1 <= -2e+270) {
tmp = -U_m;
} else if (t_1 <= -2e-240) {
tmp = sqrt(fma(0.25, (((U_m * U_m) / J_m) / J_m), 1.0)) * (-2.0 * J_m);
} else {
tmp = fma(-0.125, (U_m * (K * K)), U_m) * -1.0;
}
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(t_0 * Float64(J_m * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_1 <= -2e+270) tmp = Float64(-U_m); elseif (t_1 <= -2e-240) tmp = Float64(sqrt(fma(0.25, Float64(Float64(Float64(U_m * U_m) / J_m) / J_m), 1.0)) * Float64(-2.0 * J_m)); else tmp = Float64(fma(-0.125, Float64(U_m * Float64(K * K)), U_m) * Float64(-1.0)); 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[(t$95$0 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$1, -2e+270], (-U$95$m), If[LessEqual[t$95$1, -2e-240], N[(N[Sqrt[N[(0.25 * N[(N[(N[(U$95$m * U$95$m), $MachinePrecision] / J$95$m), $MachinePrecision] / J$95$m), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision] * N[(-2.0 * J$95$m), $MachinePrecision]), $MachinePrecision], N[(N[(-0.125 * N[(U$95$m * N[(K * K), $MachinePrecision]), $MachinePrecision] + U$95$m), $MachinePrecision] * (-1.0)), $MachinePrecision]]]), $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}{t\_0 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{+270}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -2 \cdot 10^{-240}:\\
\;\;\;\;\sqrt{\mathsf{fma}\left(0.25, \frac{\frac{U\_m \cdot U\_m}{J\_m}}{J\_m}, 1\right)} \cdot \left(-2 \cdot J\_m\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-0.125, U\_m \cdot \left(K \cdot K\right), U\_m\right) \cdot \left(-1\right)\\
\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))))) < -2.0000000000000001e270Initial program 15.1%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6453.7
Applied rewrites53.7%
if -2.0000000000000001e270 < (*.f64 (*.f64 (*.f64 #s(literal -2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64)))) (sqrt.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 U (*.f64 (*.f64 #s(literal 2 binary64) J) (cos.f64 (/.f64 K #s(literal 2 binary64))))) #s(literal 2 binary64))))) < -1.9999999999999999e-240Initial program 99.9%
Taylor expanded in K around 0
lower-sqrt.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6468.3
Applied rewrites68.3%
Applied rewrites76.2%
Taylor expanded in K around 0
*-commutativeN/A
lower-*.f6451.6
Applied rewrites51.6%
if -1.9999999999999999e-240 < (*.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 76.9%
lift-pow.f64N/A
unpow2N/A
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
associate-*l/N/A
lift-/.f64N/A
associate-*r/N/A
associate-/r*N/A
lower-/.f64N/A
Applied rewrites73.9%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f6426.6
Applied rewrites26.6%
Taylor expanded in K around 0
Applied rewrites18.2%
Taylor expanded in K around 0
Applied rewrites16.4%
Final simplification34.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 (* t_0 (* J_m 2.0))) 2.0))))))
(*
J_s
(if (<= t_1 -5e+259)
(- U_m)
(if (<= t_1 -1e-128)
(* (* -2.0 J_m) (sqrt (fma 0.25 (/ (* U_m U_m) (* J_m J_m)) 1.0)))
(* (fma -0.125 (* U_m (* K K)) U_m) (- 1.0)))))))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 / (t_0 * (J_m * 2.0))), 2.0)));
double tmp;
if (t_1 <= -5e+259) {
tmp = -U_m;
} else if (t_1 <= -1e-128) {
tmp = (-2.0 * J_m) * sqrt(fma(0.25, ((U_m * U_m) / (J_m * J_m)), 1.0));
} else {
tmp = fma(-0.125, (U_m * (K * K)), U_m) * -1.0;
}
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(t_0 * Float64(J_m * 2.0))) ^ 2.0)))) tmp = 0.0 if (t_1 <= -5e+259) tmp = Float64(-U_m); elseif (t_1 <= -1e-128) tmp = Float64(Float64(-2.0 * J_m) * sqrt(fma(0.25, Float64(Float64(U_m * U_m) / Float64(J_m * J_m)), 1.0))); else tmp = Float64(fma(-0.125, Float64(U_m * Float64(K * K)), U_m) * Float64(-1.0)); 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[(t$95$0 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, N[(J$95$s * If[LessEqual[t$95$1, -5e+259], (-U$95$m), If[LessEqual[t$95$1, -1e-128], N[(N[(-2.0 * J$95$m), $MachinePrecision] * N[Sqrt[N[(0.25 * N[(N[(U$95$m * U$95$m), $MachinePrecision] / N[(J$95$m * J$95$m), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(N[(-0.125 * N[(U$95$m * N[(K * K), $MachinePrecision]), $MachinePrecision] + U$95$m), $MachinePrecision] * (-1.0)), $MachinePrecision]]]), $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}{t\_0 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}}\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+259}:\\
\;\;\;\;-U\_m\\
\mathbf{elif}\;t\_1 \leq -1 \cdot 10^{-128}:\\
\;\;\;\;\left(-2 \cdot J\_m\right) \cdot \sqrt{\mathsf{fma}\left(0.25, \frac{U\_m \cdot U\_m}{J\_m \cdot J\_m}, 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-0.125, U\_m \cdot \left(K \cdot K\right), U\_m\right) \cdot \left(-1\right)\\
\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))))) < -5.00000000000000033e259Initial program 23.2%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6450.6
Applied rewrites50.6%
if -5.00000000000000033e259 < (*.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.00000000000000005e-128Initial program 99.9%
Taylor expanded in K around 0
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f64N/A
lower-*.f6453.5
Applied rewrites53.5%
if -1.00000000000000005e-128 < (*.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 78.7%
lift-pow.f64N/A
unpow2N/A
lift-/.f64N/A
lift-*.f64N/A
associate-/r*N/A
associate-*l/N/A
lift-/.f64N/A
associate-*r/N/A
associate-/r*N/A
lower-/.f64N/A
Applied rewrites75.9%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-sqrt.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-*.f64N/A
lower-cos.f64N/A
lower-*.f6426.6
Applied rewrites26.6%
Taylor expanded in K around 0
Applied rewrites17.7%
Taylor expanded in K around 0
Applied rewrites15.6%
Final simplification32.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 (/ U_m (* J_m 2.0))) (t_1 (cos (/ K 2.0))))
(*
J_s
(if (<=
(*
(* (* -2.0 J_m) t_1)
(sqrt (+ 1.0 (pow (/ U_m (* t_1 (* J_m 2.0))) 2.0))))
-2e+295)
(- U_m)
(* J_m (* (* -2.0 (cos (* K 0.5))) (sqrt (fma t_0 t_0 1.0))))))))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 = U_m / (J_m * 2.0);
double t_1 = cos((K / 2.0));
double tmp;
if ((((-2.0 * J_m) * t_1) * sqrt((1.0 + pow((U_m / (t_1 * (J_m * 2.0))), 2.0)))) <= -2e+295) {
tmp = -U_m;
} else {
tmp = J_m * ((-2.0 * cos((K * 0.5))) * sqrt(fma(t_0, t_0, 1.0)));
}
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(U_m / Float64(J_m * 2.0)) t_1 = cos(Float64(K / 2.0)) tmp = 0.0 if (Float64(Float64(Float64(-2.0 * J_m) * t_1) * sqrt(Float64(1.0 + (Float64(U_m / Float64(t_1 * Float64(J_m * 2.0))) ^ 2.0)))) <= -2e+295) tmp = Float64(-U_m); else tmp = Float64(J_m * Float64(Float64(-2.0 * cos(Float64(K * 0.5))) * sqrt(fma(t_0, t_0, 1.0)))); 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[(U$95$m / N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = 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$1), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U$95$m / N[(t$95$1 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], -2e+295], (-U$95$m), N[(J$95$m * N[(N[(-2.0 * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$0 * t$95$0 + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $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 := \frac{U\_m}{J\_m \cdot 2}\\
t_1 := \cos \left(\frac{K}{2}\right)\\
J\_s \cdot \begin{array}{l}
\mathbf{if}\;\left(\left(-2 \cdot J\_m\right) \cdot t\_1\right) \cdot \sqrt{1 + {\left(\frac{U\_m}{t\_1 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}} \leq -2 \cdot 10^{+295}:\\
\;\;\;\;-U\_m\\
\mathbf{else}:\\
\;\;\;\;J\_m \cdot \left(\left(-2 \cdot \cos \left(K \cdot 0.5\right)\right) \cdot \sqrt{\mathsf{fma}\left(t\_0, t\_0, 1\right)}\right)\\
\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))))) < -2e295Initial program 5.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
if -2e295 < (*.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 86.2%
Applied rewrites69.6%
lift-fma.f64N/A
Applied rewrites86.0%
Taylor expanded in K around 0
lower-*.f6472.1
Applied rewrites72.1%
Final simplification69.3%
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 (* t_0 (* J_m 2.0))) 2.0))))
-2e+295)
(- U_m)
(*
(* -2.0 (* J_m (cos (* K 0.5))))
(sqrt (fma 0.25 (/ (/ (* U_m U_m) J_m) J_m) 1.0)))))))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 / (t_0 * (J_m * 2.0))), 2.0)))) <= -2e+295) {
tmp = -U_m;
} else {
tmp = (-2.0 * (J_m * cos((K * 0.5)))) * sqrt(fma(0.25, (((U_m * U_m) / J_m) / J_m), 1.0));
}
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(t_0 * Float64(J_m * 2.0))) ^ 2.0)))) <= -2e+295) tmp = Float64(-U_m); else tmp = Float64(Float64(-2.0 * Float64(J_m * cos(Float64(K * 0.5)))) * sqrt(fma(0.25, Float64(Float64(Float64(U_m * U_m) / J_m) / J_m), 1.0))); 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]}, 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[(t$95$0 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], -2e+295], (-U$95$m), N[(N[(-2.0 * N[(J$95$m * N[Cos[N[(K * 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(0.25 * N[(N[(N[(U$95$m * U$95$m), $MachinePrecision] / J$95$m), $MachinePrecision] / J$95$m), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]), $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}{t\_0 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}} \leq -2 \cdot 10^{+295}:\\
\;\;\;\;-U\_m\\
\mathbf{else}:\\
\;\;\;\;\left(-2 \cdot \left(J\_m \cdot \cos \left(K \cdot 0.5\right)\right)\right) \cdot \sqrt{\mathsf{fma}\left(0.25, \frac{\frac{U\_m \cdot U\_m}{J\_m}}{J\_m}, 1\right)}\\
\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))))) < -2e295Initial program 5.0%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6455.2
Applied rewrites55.2%
if -2e295 < (*.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 86.2%
Taylor expanded in K around 0
lower-sqrt.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6457.4
Applied rewrites57.4%
Applied rewrites63.7%
lift-*.f64N/A
lift-*.f64N/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6463.7
lift-/.f64N/A
div-invN/A
metadata-evalN/A
lift-*.f6463.7
Applied rewrites63.7%
Final simplification62.3%
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 (* t_0 (* J_m 2.0))) 2.0))))
-2e+270)
(- U_m)
(* (fma K (* K (* J_m 0.25)) (* -2.0 J_m)) 1.0)))))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 / (t_0 * (J_m * 2.0))), 2.0)))) <= -2e+270) {
tmp = -U_m;
} else {
tmp = fma(K, (K * (J_m * 0.25)), (-2.0 * J_m)) * 1.0;
}
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(t_0 * Float64(J_m * 2.0))) ^ 2.0)))) <= -2e+270) tmp = Float64(-U_m); else tmp = Float64(fma(K, Float64(K * Float64(J_m * 0.25)), Float64(-2.0 * J_m)) * 1.0); 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]}, 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[(t$95$0 * N[(J$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], -2e+270], (-U$95$m), N[(N[(K * N[(K * N[(J$95$m * 0.25), $MachinePrecision]), $MachinePrecision] + N[(-2.0 * J$95$m), $MachinePrecision]), $MachinePrecision] * 1.0), $MachinePrecision]]), $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}{t\_0 \cdot \left(J\_m \cdot 2\right)}\right)}^{2}} \leq -2 \cdot 10^{+270}:\\
\;\;\;\;-U\_m\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(K, K \cdot \left(J\_m \cdot 0.25\right), -2 \cdot J\_m\right) \cdot 1\\
\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))))) < -2.0000000000000001e270Initial program 15.1%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6453.7
Applied rewrites53.7%
if -2.0000000000000001e270 < (*.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 85.9%
Taylor expanded in K around 0
lower-sqrt.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
unpow2N/A
lower-*.f6457.4
Applied rewrites57.4%
Taylor expanded in K around 0
associate-*r*N/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6433.2
Applied rewrites33.2%
Taylor expanded in U around 0
Applied rewrites30.9%
Final simplification35.1%
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 = abs(u)
J\_m = abs(j)
J\_s = copysign(1.0d0, j)
real(8) function code(j_s, j_m, k, u_m)
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 * Float64(-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 \left(-U\_m\right)
\end{array}
Initial program 72.9%
Taylor expanded in J around 0
mul-1-negN/A
lower-neg.f6428.3
Applied rewrites28.3%
herbie shell --seed 2024223
(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)))))