\[\left(\left(-2 \cdot J\right) \cdot \cos \left(\frac{K}{2}\right)\right) \cdot \sqrt{1 + {\left(\frac{U}{\left(2 \cdot J\right) \cdot \cos \left(\frac{K}{2}\right)}\right)}^{2}}
\]
↓
\[\begin{array}{l}
t_0 := \cos \left(K \cdot 0.5\right)\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := J \cdot \left(\mathsf{hypot}\left(1, \frac{U}{t_1 \cdot \left(J \cdot 2\right)}\right) \cdot \left(-2 \cdot t_1\right)\right)\\
\mathbf{if}\;J \leq -4.4204213601496506 \cdot 10^{-120}:\\
\;\;\;\;\left(\left(J \cdot -2\right) \cdot t_1\right) \cdot \mathsf{hypot}\left(1, U \cdot \frac{0.5}{J \cdot t_0}\right)\\
\mathbf{elif}\;J \leq -1.5110225112256388 \cdot 10^{-202}:\\
\;\;\;\;\mathsf{fma}\left(2, \frac{J}{\frac{U}{J}} \cdot {t_0}^{2}, U\right)\\
\mathbf{elif}\;J \leq 1.8809044213424615 \cdot 10^{-271}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;J \leq 3.658807514656197 \cdot 10^{-134}:\\
\;\;\;\;U\\
\mathbf{else}:\\
\;\;\;\;t_2\\
\end{array}
\]
double code(double J, double K, double U) {
return ((-2.0 * J) * cos((K / 2.0))) * sqrt((1.0 + pow((U / ((2.0 * J) * cos((K / 2.0)))), 2.0)));
}
↓
double code(double J, double K, double U) {
double t_0 = cos((K * 0.5));
double t_1 = cos((K / 2.0));
double t_2 = J * (hypot(1.0, (U / (t_1 * (J * 2.0)))) * (-2.0 * t_1));
double tmp;
if (J <= -4.4204213601496506e-120) {
tmp = ((J * -2.0) * t_1) * hypot(1.0, (U * (0.5 / (J * t_0))));
} else if (J <= -1.5110225112256388e-202) {
tmp = fma(2.0, ((J / (U / J)) * pow(t_0, 2.0)), U);
} else if (J <= 1.8809044213424615e-271) {
tmp = t_2;
} else if (J <= 3.658807514656197e-134) {
tmp = U;
} else {
tmp = t_2;
}
return tmp;
}
function code(J, K, U)
return Float64(Float64(Float64(-2.0 * J) * cos(Float64(K / 2.0))) * sqrt(Float64(1.0 + (Float64(U / Float64(Float64(2.0 * J) * cos(Float64(K / 2.0)))) ^ 2.0))))
end
↓
function code(J, K, U)
t_0 = cos(Float64(K * 0.5))
t_1 = cos(Float64(K / 2.0))
t_2 = Float64(J * Float64(hypot(1.0, Float64(U / Float64(t_1 * Float64(J * 2.0)))) * Float64(-2.0 * t_1)))
tmp = 0.0
if (J <= -4.4204213601496506e-120)
tmp = Float64(Float64(Float64(J * -2.0) * t_1) * hypot(1.0, Float64(U * Float64(0.5 / Float64(J * t_0)))));
elseif (J <= -1.5110225112256388e-202)
tmp = fma(2.0, Float64(Float64(J / Float64(U / J)) * (t_0 ^ 2.0)), U);
elseif (J <= 1.8809044213424615e-271)
tmp = t_2;
elseif (J <= 3.658807514656197e-134)
tmp = U;
else
tmp = t_2;
end
return tmp
end
code[J_, K_, U_] := N[(N[(N[(-2.0 * J), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 + N[Power[N[(U / N[(N[(2.0 * J), $MachinePrecision] * N[Cos[N[(K / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
↓
code[J_, K_, U_] := 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[(J * N[(N[Sqrt[1.0 ^ 2 + N[(U / N[(t$95$1 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] ^ 2], $MachinePrecision] * N[(-2.0 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[J, -4.4204213601496506e-120], N[(N[(N[(J * -2.0), $MachinePrecision] * t$95$1), $MachinePrecision] * N[Sqrt[1.0 ^ 2 + N[(U * N[(0.5 / N[(J * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] ^ 2], $MachinePrecision]), $MachinePrecision], If[LessEqual[J, -1.5110225112256388e-202], N[(2.0 * N[(N[(J / N[(U / J), $MachinePrecision]), $MachinePrecision] * N[Power[t$95$0, 2.0], $MachinePrecision]), $MachinePrecision] + U), $MachinePrecision], If[LessEqual[J, 1.8809044213424615e-271], t$95$2, If[LessEqual[J, 3.658807514656197e-134], U, t$95$2]]]]]]]
\left(\left(-2 \cdot J\right) \cdot \cos \left(\frac{K}{2}\right)\right) \cdot \sqrt{1 + {\left(\frac{U}{\left(2 \cdot J\right) \cdot \cos \left(\frac{K}{2}\right)}\right)}^{2}}
↓
\begin{array}{l}
t_0 := \cos \left(K \cdot 0.5\right)\\
t_1 := \cos \left(\frac{K}{2}\right)\\
t_2 := J \cdot \left(\mathsf{hypot}\left(1, \frac{U}{t_1 \cdot \left(J \cdot 2\right)}\right) \cdot \left(-2 \cdot t_1\right)\right)\\
\mathbf{if}\;J \leq -4.4204213601496506 \cdot 10^{-120}:\\
\;\;\;\;\left(\left(J \cdot -2\right) \cdot t_1\right) \cdot \mathsf{hypot}\left(1, U \cdot \frac{0.5}{J \cdot t_0}\right)\\
\mathbf{elif}\;J \leq -1.5110225112256388 \cdot 10^{-202}:\\
\;\;\;\;\mathsf{fma}\left(2, \frac{J}{\frac{U}{J}} \cdot {t_0}^{2}, U\right)\\
\mathbf{elif}\;J \leq 1.8809044213424615 \cdot 10^{-271}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;J \leq 3.658807514656197 \cdot 10^{-134}:\\
\;\;\;\;U\\
\mathbf{else}:\\
\;\;\;\;t_2\\
\end{array}
