\[\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}
\]
↓
\[\begin{array}{l}
t_1 := \sqrt{n \cdot 2}\\
\mathbf{if}\;n \leq -1 \cdot 10^{-310}:\\
\;\;\;\;\sqrt{2 \cdot \left(n \cdot \left(U \cdot t\right) + \frac{\ell \cdot -2 + \left(n \cdot \ell\right) \cdot \frac{U* - U}{Om}}{\frac{Om}{n \cdot \left(U \cdot \ell\right)}}\right)}\\
\mathbf{elif}\;n \leq 1.7 \cdot 10^{-267}:\\
\;\;\;\;\sqrt{U \cdot \left(\frac{\ell}{Om} \cdot \mathsf{fma}\left(\frac{\ell}{Om}, n \cdot \left(U* - U\right), \ell \cdot -2\right)\right) + U \cdot t} \cdot t_1\\
\mathbf{elif}\;n \leq 3.3 \cdot 10^{-217}:\\
\;\;\;\;\sqrt{2 \cdot \left(U \cdot \left(n \cdot t\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;t_1 \cdot \sqrt{U \cdot \mathsf{fma}\left(\frac{\ell}{Om}, \mathsf{fma}\left(\ell, -2, n \cdot \left(\left(U* - U\right) \cdot \frac{\ell}{Om}\right)\right), t\right)}\\
\end{array}
\]
double code(double n, double U, double t, double l, double Om, double U_42_) {
return sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l * l) / Om))) - ((n * pow((l / Om), 2.0)) * (U - U_42_)))));
}
↓
double code(double n, double U, double t, double l, double Om, double U_42_) {
double t_1 = sqrt((n * 2.0));
double tmp;
if (n <= -1e-310) {
tmp = sqrt((2.0 * ((n * (U * t)) + (((l * -2.0) + ((n * l) * ((U_42_ - U) / Om))) / (Om / (n * (U * l)))))));
} else if (n <= 1.7e-267) {
tmp = sqrt(((U * ((l / Om) * fma((l / Om), (n * (U_42_ - U)), (l * -2.0)))) + (U * t))) * t_1;
} else if (n <= 3.3e-217) {
tmp = sqrt((2.0 * (U * (n * t))));
} else {
tmp = t_1 * sqrt((U * fma((l / Om), fma(l, -2.0, (n * ((U_42_ - U) * (l / Om)))), t)));
}
return tmp;
}
function code(n, U, t, l, Om, U_42_)
return sqrt(Float64(Float64(Float64(2.0 * n) * U) * Float64(Float64(t - Float64(2.0 * Float64(Float64(l * l) / Om))) - Float64(Float64(n * (Float64(l / Om) ^ 2.0)) * Float64(U - U_42_)))))
end
↓
function code(n, U, t, l, Om, U_42_)
t_1 = sqrt(Float64(n * 2.0))
tmp = 0.0
if (n <= -1e-310)
tmp = sqrt(Float64(2.0 * Float64(Float64(n * Float64(U * t)) + Float64(Float64(Float64(l * -2.0) + Float64(Float64(n * l) * Float64(Float64(U_42_ - U) / Om))) / Float64(Om / Float64(n * Float64(U * l)))))));
elseif (n <= 1.7e-267)
tmp = Float64(sqrt(Float64(Float64(U * Float64(Float64(l / Om) * fma(Float64(l / Om), Float64(n * Float64(U_42_ - U)), Float64(l * -2.0)))) + Float64(U * t))) * t_1);
elseif (n <= 3.3e-217)
tmp = sqrt(Float64(2.0 * Float64(U * Float64(n * t))));
else
tmp = Float64(t_1 * sqrt(Float64(U * fma(Float64(l / Om), fma(l, -2.0, Float64(n * Float64(Float64(U_42_ - U) * Float64(l / Om)))), t))));
end
return tmp
end
code[n_, U_, t_, l_, Om_, U$42$_] := N[Sqrt[N[(N[(N[(2.0 * n), $MachinePrecision] * U), $MachinePrecision] * N[(N[(t - N[(2.0 * N[(N[(l * l), $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(n * N[Power[N[(l / Om), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(U - U$42$), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
↓
code[n_, U_, t_, l_, Om_, U$42$_] := Block[{t$95$1 = N[Sqrt[N[(n * 2.0), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[n, -1e-310], N[Sqrt[N[(2.0 * N[(N[(n * N[(U * t), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(l * -2.0), $MachinePrecision] + N[(N[(n * l), $MachinePrecision] * N[(N[(U$42$ - U), $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(Om / N[(n * N[(U * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[n, 1.7e-267], N[(N[Sqrt[N[(N[(U * N[(N[(l / Om), $MachinePrecision] * N[(N[(l / Om), $MachinePrecision] * N[(n * N[(U$42$ - U), $MachinePrecision]), $MachinePrecision] + N[(l * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(U * t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * t$95$1), $MachinePrecision], If[LessEqual[n, 3.3e-217], N[Sqrt[N[(2.0 * N[(U * N[(n * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[(t$95$1 * N[Sqrt[N[(U * N[(N[(l / Om), $MachinePrecision] * N[(l * -2.0 + N[(n * N[(N[(U$42$ - U), $MachinePrecision] * N[(l / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]]
\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}
↓
\begin{array}{l}
t_1 := \sqrt{n \cdot 2}\\
\mathbf{if}\;n \leq -1 \cdot 10^{-310}:\\
\;\;\;\;\sqrt{2 \cdot \left(n \cdot \left(U \cdot t\right) + \frac{\ell \cdot -2 + \left(n \cdot \ell\right) \cdot \frac{U* - U}{Om}}{\frac{Om}{n \cdot \left(U \cdot \ell\right)}}\right)}\\
\mathbf{elif}\;n \leq 1.7 \cdot 10^{-267}:\\
\;\;\;\;\sqrt{U \cdot \left(\frac{\ell}{Om} \cdot \mathsf{fma}\left(\frac{\ell}{Om}, n \cdot \left(U* - U\right), \ell \cdot -2\right)\right) + U \cdot t} \cdot t_1\\
\mathbf{elif}\;n \leq 3.3 \cdot 10^{-217}:\\
\;\;\;\;\sqrt{2 \cdot \left(U \cdot \left(n \cdot t\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;t_1 \cdot \sqrt{U \cdot \mathsf{fma}\left(\frac{\ell}{Om}, \mathsf{fma}\left(\ell, -2, n \cdot \left(\left(U* - U\right) \cdot \frac{\ell}{Om}\right)\right), t\right)}\\
\end{array}