\[\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(\frac{K}{2}\right)\\
t_1 := J \cdot \left(\mathsf{hypot}\left(1, \frac{U}{t_0 \cdot \left(J \cdot 2\right)}\right) \cdot \left(t_0 \cdot -2\right)\right)\\
\mathbf{if}\;U \leq -2.1046685566989484 \cdot 10^{+186}:\\
\;\;\;\;U\\
\mathbf{elif}\;U \leq 3.4376729941037334 \cdot 10^{+247}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;U \leq 1.822574863810491 \cdot 10^{+286}:\\
\;\;\;\;U\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\]
(FPCore (J K U)
:precision binary64
(*
(* (* -2.0 J) (cos (/ K 2.0)))
(sqrt (+ 1.0 (pow (/ U (* (* 2.0 J) (cos (/ K 2.0)))) 2.0)))))
↓
(FPCore (J K U)
:precision binary64
(let* ((t_0 (cos (/ K 2.0)))
(t_1 (* J (* (hypot 1.0 (/ U (* t_0 (* J 2.0)))) (* t_0 -2.0)))))
(if (<= U -2.1046685566989484e+186)
U
(if (<= U 3.4376729941037334e+247)
t_1
(if (<= U 1.822574863810491e+286) U t_1)))))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 / 2.0));
double t_1 = J * (hypot(1.0, (U / (t_0 * (J * 2.0)))) * (t_0 * -2.0));
double tmp;
if (U <= -2.1046685566989484e+186) {
tmp = U;
} else if (U <= 3.4376729941037334e+247) {
tmp = t_1;
} else if (U <= 1.822574863810491e+286) {
tmp = U;
} else {
tmp = t_1;
}
return tmp;
}
public static double code(double J, double K, double U) {
return ((-2.0 * J) * Math.cos((K / 2.0))) * Math.sqrt((1.0 + Math.pow((U / ((2.0 * J) * Math.cos((K / 2.0)))), 2.0)));
}
↓
public static double code(double J, double K, double U) {
double t_0 = Math.cos((K / 2.0));
double t_1 = J * (Math.hypot(1.0, (U / (t_0 * (J * 2.0)))) * (t_0 * -2.0));
double tmp;
if (U <= -2.1046685566989484e+186) {
tmp = U;
} else if (U <= 3.4376729941037334e+247) {
tmp = t_1;
} else if (U <= 1.822574863810491e+286) {
tmp = U;
} else {
tmp = t_1;
}
return tmp;
}
def code(J, K, U):
return ((-2.0 * J) * math.cos((K / 2.0))) * math.sqrt((1.0 + math.pow((U / ((2.0 * J) * math.cos((K / 2.0)))), 2.0)))
↓
def code(J, K, U):
t_0 = math.cos((K / 2.0))
t_1 = J * (math.hypot(1.0, (U / (t_0 * (J * 2.0)))) * (t_0 * -2.0))
tmp = 0
if U <= -2.1046685566989484e+186:
tmp = U
elif U <= 3.4376729941037334e+247:
tmp = t_1
elif U <= 1.822574863810491e+286:
tmp = U
else:
tmp = t_1
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 / 2.0))
t_1 = Float64(J * Float64(hypot(1.0, Float64(U / Float64(t_0 * Float64(J * 2.0)))) * Float64(t_0 * -2.0)))
tmp = 0.0
if (U <= -2.1046685566989484e+186)
tmp = U;
elseif (U <= 3.4376729941037334e+247)
tmp = t_1;
elseif (U <= 1.822574863810491e+286)
tmp = U;
else
tmp = t_1;
end
return tmp
end
function tmp = code(J, K, U)
tmp = ((-2.0 * J) * cos((K / 2.0))) * sqrt((1.0 + ((U / ((2.0 * J) * cos((K / 2.0)))) ^ 2.0)));
end
↓
function tmp_2 = code(J, K, U)
t_0 = cos((K / 2.0));
t_1 = J * (hypot(1.0, (U / (t_0 * (J * 2.0)))) * (t_0 * -2.0));
tmp = 0.0;
if (U <= -2.1046685566989484e+186)
tmp = U;
elseif (U <= 3.4376729941037334e+247)
tmp = t_1;
elseif (U <= 1.822574863810491e+286)
tmp = U;
else
tmp = t_1;
end
tmp_2 = 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 / 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[(J * N[(N[Sqrt[1.0 ^ 2 + N[(U / N[(t$95$0 * N[(J * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] ^ 2], $MachinePrecision] * N[(t$95$0 * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[U, -2.1046685566989484e+186], U, If[LessEqual[U, 3.4376729941037334e+247], t$95$1, If[LessEqual[U, 1.822574863810491e+286], U, t$95$1]]]]]
\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(\frac{K}{2}\right)\\
t_1 := J \cdot \left(\mathsf{hypot}\left(1, \frac{U}{t_0 \cdot \left(J \cdot 2\right)}\right) \cdot \left(t_0 \cdot -2\right)\right)\\
\mathbf{if}\;U \leq -2.1046685566989484 \cdot 10^{+186}:\\
\;\;\;\;U\\
\mathbf{elif}\;U \leq 3.4376729941037334 \cdot 10^{+247}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;U \leq 1.822574863810491 \cdot 10^{+286}:\\
\;\;\;\;U\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
