
(FPCore (n U t l Om U*) :precision binary64 (sqrt (* (* (* 2.0 n) U) (- (- t (* 2.0 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2.0)) (- U U*))))))
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_)))));
}
real(8) function code(n, u, t, l, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: om
real(8), intent (in) :: u_42
code = sqrt((((2.0d0 * n) * u) * ((t - (2.0d0 * ((l * l) / om))) - ((n * ((l / om) ** 2.0d0)) * (u - u_42)))))
end function
public static double code(double n, double U, double t, double l, double Om, double U_42_) {
return Math.sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l * l) / Om))) - ((n * Math.pow((l / Om), 2.0)) * (U - U_42_)))));
}
def code(n, U, t, l, Om, U_42_): return math.sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l * l) / Om))) - ((n * math.pow((l / Om), 2.0)) * (U - U_42_)))))
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 tmp = code(n, U, t, l, Om, U_42_) tmp = sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l * l) / Om))) - ((n * ((l / Om) ^ 2.0)) * (U - U_42_))))); 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]
\begin{array}{l}
\\
\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)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 23 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (n U t l Om U*) :precision binary64 (sqrt (* (* (* 2.0 n) U) (- (- t (* 2.0 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2.0)) (- U U*))))))
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_)))));
}
real(8) function code(n, u, t, l, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: om
real(8), intent (in) :: u_42
code = sqrt((((2.0d0 * n) * u) * ((t - (2.0d0 * ((l * l) / om))) - ((n * ((l / om) ** 2.0d0)) * (u - u_42)))))
end function
public static double code(double n, double U, double t, double l, double Om, double U_42_) {
return Math.sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l * l) / Om))) - ((n * Math.pow((l / Om), 2.0)) * (U - U_42_)))));
}
def code(n, U, t, l, Om, U_42_): return math.sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l * l) / Om))) - ((n * math.pow((l / Om), 2.0)) * (U - U_42_)))))
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 tmp = code(n, U, t, l, Om, U_42_) tmp = sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l * l) / Om))) - ((n * ((l / Om) ^ 2.0)) * (U - U_42_))))); 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]
\begin{array}{l}
\\
\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)}
\end{array}
l_m = (fabs.f64 l)
(FPCore (n U t l_m Om U*)
:precision binary64
(let* ((t_1 (pow (/ l_m Om) 2.0))
(t_2
(sqrt
(*
(* (* 2.0 n) U)
(+ (- t (* 2.0 (/ (* l_m l_m) Om))) (* (* n t_1) (- U* U)))))))
(if (<= t_2 0.0)
(*
(sqrt (* 2.0 n))
(pow
(* U (- t (fma 2.0 (* l_m (/ l_m Om)) (* n (* t_1 (- U U*))))))
0.5))
(if (<= t_2 INFINITY)
(sqrt
(*
(* 2.0 (* n U))
(+
t
(-
(* (* n (pow (/ Om l_m) -2.0)) (- U* U))
(* (/ l_m Om) (* 2.0 l_m))))))
(pow
(*
(exp
(*
0.25
(+
(log
(*
U
(* n (- (* 2.0 (/ -1.0 Om)) (/ (* n (- U U*)) (pow Om 2.0))))))
(* -2.0 (log (/ 1.0 l_m))))))
(sqrt (sqrt 2.0)))
2.0)))))l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double t_1 = pow((l_m / Om), 2.0);
double t_2 = sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l_m * l_m) / Om))) + ((n * t_1) * (U_42_ - U)))));
double tmp;
if (t_2 <= 0.0) {
tmp = sqrt((2.0 * n)) * pow((U * (t - fma(2.0, (l_m * (l_m / Om)), (n * (t_1 * (U - U_42_)))))), 0.5);
} else if (t_2 <= ((double) INFINITY)) {
tmp = sqrt(((2.0 * (n * U)) * (t + (((n * pow((Om / l_m), -2.0)) * (U_42_ - U)) - ((l_m / Om) * (2.0 * l_m))))));
} else {
tmp = pow((exp((0.25 * (log((U * (n * ((2.0 * (-1.0 / Om)) - ((n * (U - U_42_)) / pow(Om, 2.0)))))) + (-2.0 * log((1.0 / l_m)))))) * sqrt(sqrt(2.0))), 2.0);
}
return tmp;
}
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) t_1 = Float64(l_m / Om) ^ 2.0 t_2 = sqrt(Float64(Float64(Float64(2.0 * n) * U) * Float64(Float64(t - Float64(2.0 * Float64(Float64(l_m * l_m) / Om))) + Float64(Float64(n * t_1) * Float64(U_42_ - U))))) tmp = 0.0 if (t_2 <= 0.0) tmp = Float64(sqrt(Float64(2.0 * n)) * (Float64(U * Float64(t - fma(2.0, Float64(l_m * Float64(l_m / Om)), Float64(n * Float64(t_1 * Float64(U - U_42_)))))) ^ 0.5)); elseif (t_2 <= Inf) tmp = sqrt(Float64(Float64(2.0 * Float64(n * U)) * Float64(t + Float64(Float64(Float64(n * (Float64(Om / l_m) ^ -2.0)) * Float64(U_42_ - U)) - Float64(Float64(l_m / Om) * Float64(2.0 * l_m)))))); else tmp = Float64(exp(Float64(0.25 * Float64(log(Float64(U * Float64(n * Float64(Float64(2.0 * Float64(-1.0 / Om)) - Float64(Float64(n * Float64(U - U_42_)) / (Om ^ 2.0)))))) + Float64(-2.0 * log(Float64(1.0 / l_m)))))) * sqrt(sqrt(2.0))) ^ 2.0; end return tmp end
l_m = N[Abs[l], $MachinePrecision]
code[n_, U_, t_, l$95$m_, Om_, U$42$_] := Block[{t$95$1 = N[Power[N[(l$95$m / Om), $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[Sqrt[N[(N[(N[(2.0 * n), $MachinePrecision] * U), $MachinePrecision] * N[(N[(t - N[(2.0 * N[(N[(l$95$m * l$95$m), $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(n * t$95$1), $MachinePrecision] * N[(U$42$ - U), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t$95$2, 0.0], N[(N[Sqrt[N[(2.0 * n), $MachinePrecision]], $MachinePrecision] * N[Power[N[(U * N[(t - N[(2.0 * N[(l$95$m * N[(l$95$m / Om), $MachinePrecision]), $MachinePrecision] + N[(n * N[(t$95$1 * N[(U - U$42$), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, Infinity], N[Sqrt[N[(N[(2.0 * N[(n * U), $MachinePrecision]), $MachinePrecision] * N[(t + N[(N[(N[(n * N[Power[N[(Om / l$95$m), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(U$42$ - U), $MachinePrecision]), $MachinePrecision] - N[(N[(l$95$m / Om), $MachinePrecision] * N[(2.0 * l$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Power[N[(N[Exp[N[(0.25 * N[(N[Log[N[(U * N[(n * N[(N[(2.0 * N[(-1.0 / Om), $MachinePrecision]), $MachinePrecision] - N[(N[(n * N[(U - U$42$), $MachinePrecision]), $MachinePrecision] / N[Power[Om, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + N[(-2.0 * N[Log[N[(1.0 / l$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[Sqrt[N[Sqrt[2.0], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]]]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
t_1 := {\left(\frac{l\_m}{Om}\right)}^{2}\\
t_2 := \sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{l\_m \cdot l\_m}{Om}\right) + \left(n \cdot t\_1\right) \cdot \left(U* - U\right)\right)}\\
\mathbf{if}\;t\_2 \leq 0:\\
\;\;\;\;\sqrt{2 \cdot n} \cdot {\left(U \cdot \left(t - \mathsf{fma}\left(2, l\_m \cdot \frac{l\_m}{Om}, n \cdot \left(t\_1 \cdot \left(U - U*\right)\right)\right)\right)\right)}^{0.5}\\
\mathbf{elif}\;t\_2 \leq \infty:\\
\;\;\;\;\sqrt{\left(2 \cdot \left(n \cdot U\right)\right) \cdot \left(t + \left(\left(n \cdot {\left(\frac{Om}{l\_m}\right)}^{-2}\right) \cdot \left(U* - U\right) - \frac{l\_m}{Om} \cdot \left(2 \cdot l\_m\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;{\left(e^{0.25 \cdot \left(\log \left(U \cdot \left(n \cdot \left(2 \cdot \frac{-1}{Om} - \frac{n \cdot \left(U - U*\right)}{{Om}^{2}}\right)\right)\right) + -2 \cdot \log \left(\frac{1}{l\_m}\right)\right)} \cdot \sqrt{\sqrt{2}}\right)}^{2}\\
\end{array}
\end{array}
if (sqrt.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) n) U) (-.f64 (-.f64 t (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 l l) Om))) (*.f64 (*.f64 n (pow.f64 (/.f64 l Om) #s(literal 2 binary64))) (-.f64 U U*))))) < 0.0Initial program 16.2%
associate-*r/16.2%
*-commutative16.2%
add-sqr-sqrt14.1%
associate-*r*14.1%
Applied egg-rr14.1%
pow1/214.1%
associate-*l*21.6%
unpow-prod-down38.5%
pow1/238.5%
Applied egg-rr46.5%
if 0.0 < (sqrt.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) n) U) (-.f64 (-.f64 t (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 l l) Om))) (*.f64 (*.f64 n (pow.f64 (/.f64 l Om) #s(literal 2 binary64))) (-.f64 U U*))))) < +inf.0Initial program 66.1%
Simplified73.7%
associate-*r*73.7%
fma-define73.8%
associate-*r*72.7%
Applied egg-rr72.7%
clear-num72.7%
inv-pow72.7%
Applied egg-rr72.7%
unpow-172.7%
Simplified72.7%
fma-undefine72.7%
clear-num72.7%
associate-*r*73.7%
clear-num73.7%
inv-pow73.7%
pow-pow73.8%
metadata-eval73.8%
Applied egg-rr73.8%
if +inf.0 < (sqrt.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) n) U) (-.f64 (-.f64 t (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 l l) Om))) (*.f64 (*.f64 n (pow.f64 (/.f64 l Om) #s(literal 2 binary64))) (-.f64 U U*))))) Initial program 0.0%
Simplified7.4%
associate-*r*7.4%
fma-define22.4%
associate-*r*22.5%
Applied egg-rr22.5%
clear-num22.5%
inv-pow22.5%
Applied egg-rr22.5%
unpow-122.5%
Simplified22.5%
add-sqr-sqrt22.5%
pow222.5%
Applied egg-rr22.3%
Taylor expanded in l around inf 26.5%
Final simplification63.4%
l_m = (fabs.f64 l)
(FPCore (n U t l_m Om U*)
:precision binary64
(let* ((t_1 (pow (/ l_m Om) 2.0))
(t_2
(sqrt
(*
(* (* 2.0 n) U)
(+ (- t (* 2.0 (/ (* l_m l_m) Om))) (* (* n t_1) (- U* U)))))))
(if (<= t_2 0.0)
(*
(sqrt (* 2.0 n))
(pow
(* U (- t (fma 2.0 (* l_m (/ l_m Om)) (* n (* t_1 (- U U*))))))
0.5))
(if (<= t_2 INFINITY)
(sqrt
(*
(* 2.0 (* n U))
(+
t
(-
(* (* n (pow (/ Om l_m) -2.0)) (- U* U))
(* (/ l_m Om) (* 2.0 l_m))))))
(* (* (/ 1.0 Om) (* l_m (* n (sqrt 2.0)))) (sqrt (* U U*)))))))l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double t_1 = pow((l_m / Om), 2.0);
double t_2 = sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l_m * l_m) / Om))) + ((n * t_1) * (U_42_ - U)))));
double tmp;
if (t_2 <= 0.0) {
tmp = sqrt((2.0 * n)) * pow((U * (t - fma(2.0, (l_m * (l_m / Om)), (n * (t_1 * (U - U_42_)))))), 0.5);
} else if (t_2 <= ((double) INFINITY)) {
tmp = sqrt(((2.0 * (n * U)) * (t + (((n * pow((Om / l_m), -2.0)) * (U_42_ - U)) - ((l_m / Om) * (2.0 * l_m))))));
} else {
tmp = ((1.0 / Om) * (l_m * (n * sqrt(2.0)))) * sqrt((U * U_42_));
}
return tmp;
}
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) t_1 = Float64(l_m / Om) ^ 2.0 t_2 = sqrt(Float64(Float64(Float64(2.0 * n) * U) * Float64(Float64(t - Float64(2.0 * Float64(Float64(l_m * l_m) / Om))) + Float64(Float64(n * t_1) * Float64(U_42_ - U))))) tmp = 0.0 if (t_2 <= 0.0) tmp = Float64(sqrt(Float64(2.0 * n)) * (Float64(U * Float64(t - fma(2.0, Float64(l_m * Float64(l_m / Om)), Float64(n * Float64(t_1 * Float64(U - U_42_)))))) ^ 0.5)); elseif (t_2 <= Inf) tmp = sqrt(Float64(Float64(2.0 * Float64(n * U)) * Float64(t + Float64(Float64(Float64(n * (Float64(Om / l_m) ^ -2.0)) * Float64(U_42_ - U)) - Float64(Float64(l_m / Om) * Float64(2.0 * l_m)))))); else tmp = Float64(Float64(Float64(1.0 / Om) * Float64(l_m * Float64(n * sqrt(2.0)))) * sqrt(Float64(U * U_42_))); end return tmp end
l_m = N[Abs[l], $MachinePrecision]
code[n_, U_, t_, l$95$m_, Om_, U$42$_] := Block[{t$95$1 = N[Power[N[(l$95$m / Om), $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[Sqrt[N[(N[(N[(2.0 * n), $MachinePrecision] * U), $MachinePrecision] * N[(N[(t - N[(2.0 * N[(N[(l$95$m * l$95$m), $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(n * t$95$1), $MachinePrecision] * N[(U$42$ - U), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t$95$2, 0.0], N[(N[Sqrt[N[(2.0 * n), $MachinePrecision]], $MachinePrecision] * N[Power[N[(U * N[(t - N[(2.0 * N[(l$95$m * N[(l$95$m / Om), $MachinePrecision]), $MachinePrecision] + N[(n * N[(t$95$1 * N[(U - U$42$), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, Infinity], N[Sqrt[N[(N[(2.0 * N[(n * U), $MachinePrecision]), $MachinePrecision] * N[(t + N[(N[(N[(n * N[Power[N[(Om / l$95$m), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(U$42$ - U), $MachinePrecision]), $MachinePrecision] - N[(N[(l$95$m / Om), $MachinePrecision] * N[(2.0 * l$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[(N[(N[(1.0 / Om), $MachinePrecision] * N[(l$95$m * N[(n * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(U * U$42$), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
t_1 := {\left(\frac{l\_m}{Om}\right)}^{2}\\
t_2 := \sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{l\_m \cdot l\_m}{Om}\right) + \left(n \cdot t\_1\right) \cdot \left(U* - U\right)\right)}\\
\mathbf{if}\;t\_2 \leq 0:\\
\;\;\;\;\sqrt{2 \cdot n} \cdot {\left(U \cdot \left(t - \mathsf{fma}\left(2, l\_m \cdot \frac{l\_m}{Om}, n \cdot \left(t\_1 \cdot \left(U - U*\right)\right)\right)\right)\right)}^{0.5}\\
\mathbf{elif}\;t\_2 \leq \infty:\\
\;\;\;\;\sqrt{\left(2 \cdot \left(n \cdot U\right)\right) \cdot \left(t + \left(\left(n \cdot {\left(\frac{Om}{l\_m}\right)}^{-2}\right) \cdot \left(U* - U\right) - \frac{l\_m}{Om} \cdot \left(2 \cdot l\_m\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{1}{Om} \cdot \left(l\_m \cdot \left(n \cdot \sqrt{2}\right)\right)\right) \cdot \sqrt{U \cdot U*}\\
\end{array}
\end{array}
if (sqrt.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) n) U) (-.f64 (-.f64 t (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 l l) Om))) (*.f64 (*.f64 n (pow.f64 (/.f64 l Om) #s(literal 2 binary64))) (-.f64 U U*))))) < 0.0Initial program 16.2%
associate-*r/16.2%
*-commutative16.2%
add-sqr-sqrt14.1%
associate-*r*14.1%
Applied egg-rr14.1%
pow1/214.1%
associate-*l*21.6%
unpow-prod-down38.5%
pow1/238.5%
Applied egg-rr46.5%
if 0.0 < (sqrt.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) n) U) (-.f64 (-.f64 t (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 l l) Om))) (*.f64 (*.f64 n (pow.f64 (/.f64 l Om) #s(literal 2 binary64))) (-.f64 U U*))))) < +inf.0Initial program 66.1%
Simplified73.7%
associate-*r*73.7%
fma-define73.8%
associate-*r*72.7%
Applied egg-rr72.7%
clear-num72.7%
inv-pow72.7%
Applied egg-rr72.7%
unpow-172.7%
Simplified72.7%
fma-undefine72.7%
clear-num72.7%
associate-*r*73.7%
clear-num73.7%
inv-pow73.7%
pow-pow73.8%
metadata-eval73.8%
Applied egg-rr73.8%
if +inf.0 < (sqrt.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) n) U) (-.f64 (-.f64 t (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 l l) Om))) (*.f64 (*.f64 n (pow.f64 (/.f64 l Om) #s(literal 2 binary64))) (-.f64 U U*))))) Initial program 0.0%
Simplified8.0%
Taylor expanded in U* around inf 32.6%
div-inv32.6%
Applied egg-rr32.6%
Final simplification64.4%
l_m = (fabs.f64 l)
(FPCore (n U t l_m Om U*)
:precision binary64
(let* ((t_1
(sqrt
(*
(* (* 2.0 n) U)
(+
(- t (* 2.0 (/ (* l_m l_m) Om)))
(* (* n (pow (/ l_m Om) 2.0)) (- U* U)))))))
(if (<= t_1 0.0)
(* (sqrt (* 2.0 n)) (sqrt (* U t)))
(if (<= t_1 INFINITY)
(sqrt
(*
(* 2.0 (* n U))
(+
t
(-
(* (* n (pow (/ Om l_m) -2.0)) (- U* U))
(* (/ l_m Om) (* 2.0 l_m))))))
(* (* (/ 1.0 Om) (* l_m (* n (sqrt 2.0)))) (sqrt (* U U*)))))))l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double t_1 = sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l_m * l_m) / Om))) + ((n * pow((l_m / Om), 2.0)) * (U_42_ - U)))));
double tmp;
if (t_1 <= 0.0) {
tmp = sqrt((2.0 * n)) * sqrt((U * t));
} else if (t_1 <= ((double) INFINITY)) {
tmp = sqrt(((2.0 * (n * U)) * (t + (((n * pow((Om / l_m), -2.0)) * (U_42_ - U)) - ((l_m / Om) * (2.0 * l_m))))));
} else {
tmp = ((1.0 / Om) * (l_m * (n * sqrt(2.0)))) * sqrt((U * U_42_));
}
return tmp;
}
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double t_1 = Math.sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l_m * l_m) / Om))) + ((n * Math.pow((l_m / Om), 2.0)) * (U_42_ - U)))));
double tmp;
if (t_1 <= 0.0) {
tmp = Math.sqrt((2.0 * n)) * Math.sqrt((U * t));
} else if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = Math.sqrt(((2.0 * (n * U)) * (t + (((n * Math.pow((Om / l_m), -2.0)) * (U_42_ - U)) - ((l_m / Om) * (2.0 * l_m))))));
} else {
tmp = ((1.0 / Om) * (l_m * (n * Math.sqrt(2.0)))) * Math.sqrt((U * U_42_));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): t_1 = math.sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l_m * l_m) / Om))) + ((n * math.pow((l_m / Om), 2.0)) * (U_42_ - U))))) tmp = 0 if t_1 <= 0.0: tmp = math.sqrt((2.0 * n)) * math.sqrt((U * t)) elif t_1 <= math.inf: tmp = math.sqrt(((2.0 * (n * U)) * (t + (((n * math.pow((Om / l_m), -2.0)) * (U_42_ - U)) - ((l_m / Om) * (2.0 * l_m)))))) else: tmp = ((1.0 / Om) * (l_m * (n * math.sqrt(2.0)))) * math.sqrt((U * U_42_)) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) t_1 = sqrt(Float64(Float64(Float64(2.0 * n) * U) * Float64(Float64(t - Float64(2.0 * Float64(Float64(l_m * l_m) / Om))) + Float64(Float64(n * (Float64(l_m / Om) ^ 2.0)) * Float64(U_42_ - U))))) tmp = 0.0 if (t_1 <= 0.0) tmp = Float64(sqrt(Float64(2.0 * n)) * sqrt(Float64(U * t))); elseif (t_1 <= Inf) tmp = sqrt(Float64(Float64(2.0 * Float64(n * U)) * Float64(t + Float64(Float64(Float64(n * (Float64(Om / l_m) ^ -2.0)) * Float64(U_42_ - U)) - Float64(Float64(l_m / Om) * Float64(2.0 * l_m)))))); else tmp = Float64(Float64(Float64(1.0 / Om) * Float64(l_m * Float64(n * sqrt(2.0)))) * sqrt(Float64(U * U_42_))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) t_1 = sqrt((((2.0 * n) * U) * ((t - (2.0 * ((l_m * l_m) / Om))) + ((n * ((l_m / Om) ^ 2.0)) * (U_42_ - U))))); tmp = 0.0; if (t_1 <= 0.0) tmp = sqrt((2.0 * n)) * sqrt((U * t)); elseif (t_1 <= Inf) tmp = sqrt(((2.0 * (n * U)) * (t + (((n * ((Om / l_m) ^ -2.0)) * (U_42_ - U)) - ((l_m / Om) * (2.0 * l_m)))))); else tmp = ((1.0 / Om) * (l_m * (n * sqrt(2.0)))) * sqrt((U * U_42_)); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision]
code[n_, U_, t_, l$95$m_, Om_, U$42$_] := Block[{t$95$1 = N[Sqrt[N[(N[(N[(2.0 * n), $MachinePrecision] * U), $MachinePrecision] * N[(N[(t - N[(2.0 * N[(N[(l$95$m * l$95$m), $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(n * N[Power[N[(l$95$m / Om), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(U$42$ - U), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[t$95$1, 0.0], N[(N[Sqrt[N[(2.0 * n), $MachinePrecision]], $MachinePrecision] * N[Sqrt[N[(U * t), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, Infinity], N[Sqrt[N[(N[(2.0 * N[(n * U), $MachinePrecision]), $MachinePrecision] * N[(t + N[(N[(N[(n * N[Power[N[(Om / l$95$m), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(U$42$ - U), $MachinePrecision]), $MachinePrecision] - N[(N[(l$95$m / Om), $MachinePrecision] * N[(2.0 * l$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[(N[(N[(1.0 / Om), $MachinePrecision] * N[(l$95$m * N[(n * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(U * U$42$), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
t_1 := \sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{l\_m \cdot l\_m}{Om}\right) + \left(n \cdot {\left(\frac{l\_m}{Om}\right)}^{2}\right) \cdot \left(U* - U\right)\right)}\\
\mathbf{if}\;t\_1 \leq 0:\\
\;\;\;\;\sqrt{2 \cdot n} \cdot \sqrt{U \cdot t}\\
\mathbf{elif}\;t\_1 \leq \infty:\\
\;\;\;\;\sqrt{\left(2 \cdot \left(n \cdot U\right)\right) \cdot \left(t + \left(\left(n \cdot {\left(\frac{Om}{l\_m}\right)}^{-2}\right) \cdot \left(U* - U\right) - \frac{l\_m}{Om} \cdot \left(2 \cdot l\_m\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{1}{Om} \cdot \left(l\_m \cdot \left(n \cdot \sqrt{2}\right)\right)\right) \cdot \sqrt{U \cdot U*}\\
\end{array}
\end{array}
if (sqrt.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) n) U) (-.f64 (-.f64 t (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 l l) Om))) (*.f64 (*.f64 n (pow.f64 (/.f64 l Om) #s(literal 2 binary64))) (-.f64 U U*))))) < 0.0Initial program 16.2%
Simplified30.9%
Taylor expanded in t around inf 27.2%
sqrt-prod39.1%
Applied egg-rr39.1%
*-commutative39.1%
Simplified39.1%
if 0.0 < (sqrt.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) n) U) (-.f64 (-.f64 t (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 l l) Om))) (*.f64 (*.f64 n (pow.f64 (/.f64 l Om) #s(literal 2 binary64))) (-.f64 U U*))))) < +inf.0Initial program 66.1%
Simplified73.7%
associate-*r*73.7%
fma-define73.8%
associate-*r*72.7%
Applied egg-rr72.7%
clear-num72.7%
inv-pow72.7%
Applied egg-rr72.7%
unpow-172.7%
Simplified72.7%
fma-undefine72.7%
clear-num72.7%
associate-*r*73.7%
clear-num73.7%
inv-pow73.7%
pow-pow73.8%
metadata-eval73.8%
Applied egg-rr73.8%
if +inf.0 < (sqrt.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) n) U) (-.f64 (-.f64 t (*.f64 #s(literal 2 binary64) (/.f64 (*.f64 l l) Om))) (*.f64 (*.f64 n (pow.f64 (/.f64 l Om) #s(literal 2 binary64))) (-.f64 U U*))))) Initial program 0.0%
Simplified8.0%
Taylor expanded in U* around inf 32.6%
div-inv32.6%
Applied egg-rr32.6%
Final simplification63.6%
l_m = (fabs.f64 l)
(FPCore (n U t l_m Om U*)
:precision binary64
(let* ((t_1
(sqrt
(* (* (* 2.0 n) U) (+ t (* (* n (pow (/ l_m Om) 2.0)) (- U* U)))))))
(if (<= U* -4e+67)
t_1
(if (<= U* 2.3e-252)
(sqrt
(+ (* -4.0 (/ (* U (* n (pow l_m 2.0))) Om)) (* 2.0 (* U (* n t)))))
(if (<= U* 4e-114)
(sqrt (* (* 2.0 n) (* U (- t (* 2.0 (/ (pow l_m 2.0) Om))))))
t_1)))))l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double t_1 = sqrt((((2.0 * n) * U) * (t + ((n * pow((l_m / Om), 2.0)) * (U_42_ - U)))));
double tmp;
if (U_42_ <= -4e+67) {
tmp = t_1;
} else if (U_42_ <= 2.3e-252) {
tmp = sqrt(((-4.0 * ((U * (n * pow(l_m, 2.0))) / Om)) + (2.0 * (U * (n * t)))));
} else if (U_42_ <= 4e-114) {
tmp = sqrt(((2.0 * n) * (U * (t - (2.0 * (pow(l_m, 2.0) / Om))))));
} else {
tmp = t_1;
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: t_1
real(8) :: tmp
t_1 = sqrt((((2.0d0 * n) * u) * (t + ((n * ((l_m / om) ** 2.0d0)) * (u_42 - u)))))
if (u_42 <= (-4d+67)) then
tmp = t_1
else if (u_42 <= 2.3d-252) then
tmp = sqrt((((-4.0d0) * ((u * (n * (l_m ** 2.0d0))) / om)) + (2.0d0 * (u * (n * t)))))
else if (u_42 <= 4d-114) then
tmp = sqrt(((2.0d0 * n) * (u * (t - (2.0d0 * ((l_m ** 2.0d0) / om))))))
else
tmp = t_1
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double t_1 = Math.sqrt((((2.0 * n) * U) * (t + ((n * Math.pow((l_m / Om), 2.0)) * (U_42_ - U)))));
double tmp;
if (U_42_ <= -4e+67) {
tmp = t_1;
} else if (U_42_ <= 2.3e-252) {
tmp = Math.sqrt(((-4.0 * ((U * (n * Math.pow(l_m, 2.0))) / Om)) + (2.0 * (U * (n * t)))));
} else if (U_42_ <= 4e-114) {
tmp = Math.sqrt(((2.0 * n) * (U * (t - (2.0 * (Math.pow(l_m, 2.0) / Om))))));
} else {
tmp = t_1;
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): t_1 = math.sqrt((((2.0 * n) * U) * (t + ((n * math.pow((l_m / Om), 2.0)) * (U_42_ - U))))) tmp = 0 if U_42_ <= -4e+67: tmp = t_1 elif U_42_ <= 2.3e-252: tmp = math.sqrt(((-4.0 * ((U * (n * math.pow(l_m, 2.0))) / Om)) + (2.0 * (U * (n * t))))) elif U_42_ <= 4e-114: tmp = math.sqrt(((2.0 * n) * (U * (t - (2.0 * (math.pow(l_m, 2.0) / Om)))))) else: tmp = t_1 return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) t_1 = sqrt(Float64(Float64(Float64(2.0 * n) * U) * Float64(t + Float64(Float64(n * (Float64(l_m / Om) ^ 2.0)) * Float64(U_42_ - U))))) tmp = 0.0 if (U_42_ <= -4e+67) tmp = t_1; elseif (U_42_ <= 2.3e-252) tmp = sqrt(Float64(Float64(-4.0 * Float64(Float64(U * Float64(n * (l_m ^ 2.0))) / Om)) + Float64(2.0 * Float64(U * Float64(n * t))))); elseif (U_42_ <= 4e-114) tmp = sqrt(Float64(Float64(2.0 * n) * Float64(U * Float64(t - Float64(2.0 * Float64((l_m ^ 2.0) / Om)))))); else tmp = t_1; end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) t_1 = sqrt((((2.0 * n) * U) * (t + ((n * ((l_m / Om) ^ 2.0)) * (U_42_ - U))))); tmp = 0.0; if (U_42_ <= -4e+67) tmp = t_1; elseif (U_42_ <= 2.3e-252) tmp = sqrt(((-4.0 * ((U * (n * (l_m ^ 2.0))) / Om)) + (2.0 * (U * (n * t))))); elseif (U_42_ <= 4e-114) tmp = sqrt(((2.0 * n) * (U * (t - (2.0 * ((l_m ^ 2.0) / Om)))))); else tmp = t_1; end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision]
code[n_, U_, t_, l$95$m_, Om_, U$42$_] := Block[{t$95$1 = N[Sqrt[N[(N[(N[(2.0 * n), $MachinePrecision] * U), $MachinePrecision] * N[(t + N[(N[(n * N[Power[N[(l$95$m / Om), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(U$42$ - U), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[U$42$, -4e+67], t$95$1, If[LessEqual[U$42$, 2.3e-252], N[Sqrt[N[(N[(-4.0 * N[(N[(U * N[(n * N[Power[l$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision] + N[(2.0 * N[(U * N[(n * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[U$42$, 4e-114], N[Sqrt[N[(N[(2.0 * n), $MachinePrecision] * N[(U * N[(t - N[(2.0 * N[(N[Power[l$95$m, 2.0], $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], t$95$1]]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
t_1 := \sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(t + \left(n \cdot {\left(\frac{l\_m}{Om}\right)}^{2}\right) \cdot \left(U* - U\right)\right)}\\
\mathbf{if}\;U* \leq -4 \cdot 10^{+67}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;U* \leq 2.3 \cdot 10^{-252}:\\
\;\;\;\;\sqrt{-4 \cdot \frac{U \cdot \left(n \cdot {l\_m}^{2}\right)}{Om} + 2 \cdot \left(U \cdot \left(n \cdot t\right)\right)}\\
\mathbf{elif}\;U* \leq 4 \cdot 10^{-114}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(t - 2 \cdot \frac{{l\_m}^{2}}{Om}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if U* < -3.99999999999999993e67 or 4.0000000000000002e-114 < U* Initial program 51.5%
associate-*r/56.6%
*-commutative56.6%
add-sqr-sqrt32.7%
associate-*r*32.7%
Applied egg-rr32.7%
Taylor expanded in t around inf 56.2%
if -3.99999999999999993e67 < U* < 2.2999999999999998e-252Initial program 54.3%
Simplified60.7%
Taylor expanded in Om around inf 57.8%
if 2.2999999999999998e-252 < U* < 4.0000000000000002e-114Initial program 34.6%
Simplified51.8%
Taylor expanded in n around 0 57.6%
Final simplification56.9%
l_m = (fabs.f64 l)
(FPCore (n U t l_m Om U*)
:precision binary64
(let* ((t_1 (* n (pow (/ l_m Om) 2.0))))
(if (<= l_m 2.9e+270)
(sqrt
(*
(* 2.0 (* n U))
(+ t (- (* t_1 (- U* U)) (* 2.0 (* l_m (/ l_m Om)))))))
(sqrt (* (* (* 2.0 n) U) (* t_1 U*))))))l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double t_1 = n * pow((l_m / Om), 2.0);
double tmp;
if (l_m <= 2.9e+270) {
tmp = sqrt(((2.0 * (n * U)) * (t + ((t_1 * (U_42_ - U)) - (2.0 * (l_m * (l_m / Om)))))));
} else {
tmp = sqrt((((2.0 * n) * U) * (t_1 * U_42_)));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: t_1
real(8) :: tmp
t_1 = n * ((l_m / om) ** 2.0d0)
if (l_m <= 2.9d+270) then
tmp = sqrt(((2.0d0 * (n * u)) * (t + ((t_1 * (u_42 - u)) - (2.0d0 * (l_m * (l_m / om)))))))
else
tmp = sqrt((((2.0d0 * n) * u) * (t_1 * u_42)))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double t_1 = n * Math.pow((l_m / Om), 2.0);
double tmp;
if (l_m <= 2.9e+270) {
tmp = Math.sqrt(((2.0 * (n * U)) * (t + ((t_1 * (U_42_ - U)) - (2.0 * (l_m * (l_m / Om)))))));
} else {
tmp = Math.sqrt((((2.0 * n) * U) * (t_1 * U_42_)));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): t_1 = n * math.pow((l_m / Om), 2.0) tmp = 0 if l_m <= 2.9e+270: tmp = math.sqrt(((2.0 * (n * U)) * (t + ((t_1 * (U_42_ - U)) - (2.0 * (l_m * (l_m / Om))))))) else: tmp = math.sqrt((((2.0 * n) * U) * (t_1 * U_42_))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) t_1 = Float64(n * (Float64(l_m / Om) ^ 2.0)) tmp = 0.0 if (l_m <= 2.9e+270) tmp = sqrt(Float64(Float64(2.0 * Float64(n * U)) * Float64(t + Float64(Float64(t_1 * Float64(U_42_ - U)) - Float64(2.0 * Float64(l_m * Float64(l_m / Om))))))); else tmp = sqrt(Float64(Float64(Float64(2.0 * n) * U) * Float64(t_1 * U_42_))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) t_1 = n * ((l_m / Om) ^ 2.0); tmp = 0.0; if (l_m <= 2.9e+270) tmp = sqrt(((2.0 * (n * U)) * (t + ((t_1 * (U_42_ - U)) - (2.0 * (l_m * (l_m / Om))))))); else tmp = sqrt((((2.0 * n) * U) * (t_1 * U_42_))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision]
code[n_, U_, t_, l$95$m_, Om_, U$42$_] := Block[{t$95$1 = N[(n * N[Power[N[(l$95$m / Om), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[l$95$m, 2.9e+270], N[Sqrt[N[(N[(2.0 * N[(n * U), $MachinePrecision]), $MachinePrecision] * N[(t + N[(N[(t$95$1 * N[(U$42$ - U), $MachinePrecision]), $MachinePrecision] - N[(2.0 * N[(l$95$m * N[(l$95$m / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Sqrt[N[(N[(N[(2.0 * n), $MachinePrecision] * U), $MachinePrecision] * N[(t$95$1 * U$42$), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
t_1 := n \cdot {\left(\frac{l\_m}{Om}\right)}^{2}\\
\mathbf{if}\;l\_m \leq 2.9 \cdot 10^{+270}:\\
\;\;\;\;\sqrt{\left(2 \cdot \left(n \cdot U\right)\right) \cdot \left(t + \left(t\_1 \cdot \left(U* - U\right) - 2 \cdot \left(l\_m \cdot \frac{l\_m}{Om}\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(t\_1 \cdot U*\right)}\\
\end{array}
\end{array}
if l < 2.8999999999999999e270Initial program 51.9%
Simplified58.9%
if 2.8999999999999999e270 < l Initial program 12.7%
associate-*r/12.7%
*-commutative12.7%
add-sqr-sqrt12.7%
associate-*r*12.7%
Applied egg-rr12.7%
Taylor expanded in U* around inf 22.8%
associate-/l*22.8%
*-commutative22.8%
associate-/l*22.8%
unpow222.8%
unpow222.8%
times-frac34.6%
unpow234.6%
Simplified34.6%
Final simplification58.1%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (if (or (<= n -9e-113) (not (<= n 2.7e-96))) (sqrt (* (* (* 2.0 n) U) (+ t (* (* n (pow (/ l_m Om) 2.0)) (- U* U))))) (sqrt (* 2.0 (* U (* n (- t (* 2.0 (/ (pow l_m 2.0) Om)))))))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if ((n <= -9e-113) || !(n <= 2.7e-96)) {
tmp = sqrt((((2.0 * n) * U) * (t + ((n * pow((l_m / Om), 2.0)) * (U_42_ - U)))));
} else {
tmp = sqrt((2.0 * (U * (n * (t - (2.0 * (pow(l_m, 2.0) / Om)))))));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if ((n <= (-9d-113)) .or. (.not. (n <= 2.7d-96))) then
tmp = sqrt((((2.0d0 * n) * u) * (t + ((n * ((l_m / om) ** 2.0d0)) * (u_42 - u)))))
else
tmp = sqrt((2.0d0 * (u * (n * (t - (2.0d0 * ((l_m ** 2.0d0) / om)))))))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if ((n <= -9e-113) || !(n <= 2.7e-96)) {
tmp = Math.sqrt((((2.0 * n) * U) * (t + ((n * Math.pow((l_m / Om), 2.0)) * (U_42_ - U)))));
} else {
tmp = Math.sqrt((2.0 * (U * (n * (t - (2.0 * (Math.pow(l_m, 2.0) / Om)))))));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if (n <= -9e-113) or not (n <= 2.7e-96): tmp = math.sqrt((((2.0 * n) * U) * (t + ((n * math.pow((l_m / Om), 2.0)) * (U_42_ - U))))) else: tmp = math.sqrt((2.0 * (U * (n * (t - (2.0 * (math.pow(l_m, 2.0) / Om))))))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if ((n <= -9e-113) || !(n <= 2.7e-96)) tmp = sqrt(Float64(Float64(Float64(2.0 * n) * U) * Float64(t + Float64(Float64(n * (Float64(l_m / Om) ^ 2.0)) * Float64(U_42_ - U))))); else tmp = sqrt(Float64(2.0 * Float64(U * Float64(n * Float64(t - Float64(2.0 * Float64((l_m ^ 2.0) / Om))))))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if ((n <= -9e-113) || ~((n <= 2.7e-96))) tmp = sqrt((((2.0 * n) * U) * (t + ((n * ((l_m / Om) ^ 2.0)) * (U_42_ - U))))); else tmp = sqrt((2.0 * (U * (n * (t - (2.0 * ((l_m ^ 2.0) / Om))))))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[Or[LessEqual[n, -9e-113], N[Not[LessEqual[n, 2.7e-96]], $MachinePrecision]], N[Sqrt[N[(N[(N[(2.0 * n), $MachinePrecision] * U), $MachinePrecision] * N[(t + N[(N[(n * N[Power[N[(l$95$m / Om), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(U$42$ - U), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Sqrt[N[(2.0 * N[(U * N[(n * N[(t - N[(2.0 * N[(N[Power[l$95$m, 2.0], $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;n \leq -9 \cdot 10^{-113} \lor \neg \left(n \leq 2.7 \cdot 10^{-96}\right):\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(t + \left(n \cdot {\left(\frac{l\_m}{Om}\right)}^{2}\right) \cdot \left(U* - U\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2 \cdot \left(U \cdot \left(n \cdot \left(t - 2 \cdot \frac{{l\_m}^{2}}{Om}\right)\right)\right)}\\
\end{array}
\end{array}
if n < -9.0000000000000002e-113 or 2.7e-96 < n Initial program 53.8%
associate-*r/60.7%
*-commutative60.7%
add-sqr-sqrt36.3%
associate-*r*36.3%
Applied egg-rr36.3%
Taylor expanded in t around inf 58.7%
if -9.0000000000000002e-113 < n < 2.7e-96Initial program 44.0%
Simplified44.2%
Taylor expanded in n around 0 48.5%
Final simplification55.3%
l_m = (fabs.f64 l)
(FPCore (n U t l_m Om U*)
:precision binary64
(if (<= Om -4.2e-120)
(pow (* (* n t) (* 2.0 U)) 0.5)
(if (<= Om 4.9e-49)
(/ (* (* l_m (* n (sqrt 2.0))) (sqrt (* U U*))) Om)
(sqrt (fabs (* 2.0 (* n (* U t))))))))l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (Om <= -4.2e-120) {
tmp = pow(((n * t) * (2.0 * U)), 0.5);
} else if (Om <= 4.9e-49) {
tmp = ((l_m * (n * sqrt(2.0))) * sqrt((U * U_42_))) / Om;
} else {
tmp = sqrt(fabs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (om <= (-4.2d-120)) then
tmp = ((n * t) * (2.0d0 * u)) ** 0.5d0
else if (om <= 4.9d-49) then
tmp = ((l_m * (n * sqrt(2.0d0))) * sqrt((u * u_42))) / om
else
tmp = sqrt(abs((2.0d0 * (n * (u * t)))))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (Om <= -4.2e-120) {
tmp = Math.pow(((n * t) * (2.0 * U)), 0.5);
} else if (Om <= 4.9e-49) {
tmp = ((l_m * (n * Math.sqrt(2.0))) * Math.sqrt((U * U_42_))) / Om;
} else {
tmp = Math.sqrt(Math.abs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if Om <= -4.2e-120: tmp = math.pow(((n * t) * (2.0 * U)), 0.5) elif Om <= 4.9e-49: tmp = ((l_m * (n * math.sqrt(2.0))) * math.sqrt((U * U_42_))) / Om else: tmp = math.sqrt(math.fabs((2.0 * (n * (U * t))))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (Om <= -4.2e-120) tmp = Float64(Float64(n * t) * Float64(2.0 * U)) ^ 0.5; elseif (Om <= 4.9e-49) tmp = Float64(Float64(Float64(l_m * Float64(n * sqrt(2.0))) * sqrt(Float64(U * U_42_))) / Om); else tmp = sqrt(abs(Float64(2.0 * Float64(n * Float64(U * t))))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (Om <= -4.2e-120) tmp = ((n * t) * (2.0 * U)) ^ 0.5; elseif (Om <= 4.9e-49) tmp = ((l_m * (n * sqrt(2.0))) * sqrt((U * U_42_))) / Om; else tmp = sqrt(abs((2.0 * (n * (U * t))))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[Om, -4.2e-120], N[Power[N[(N[(n * t), $MachinePrecision] * N[(2.0 * U), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision], If[LessEqual[Om, 4.9e-49], N[(N[(N[(l$95$m * N[(n * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(U * U$42$), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / Om), $MachinePrecision], N[Sqrt[N[Abs[N[(2.0 * N[(n * N[(U * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;Om \leq -4.2 \cdot 10^{-120}:\\
\;\;\;\;{\left(\left(n \cdot t\right) \cdot \left(2 \cdot U\right)\right)}^{0.5}\\
\mathbf{elif}\;Om \leq 4.9 \cdot 10^{-49}:\\
\;\;\;\;\frac{\left(l\_m \cdot \left(n \cdot \sqrt{2}\right)\right) \cdot \sqrt{U \cdot U*}}{Om}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left|2 \cdot \left(n \cdot \left(U \cdot t\right)\right)\right|}\\
\end{array}
\end{array}
if Om < -4.2000000000000001e-120Initial program 54.4%
Simplified60.3%
Taylor expanded in t around inf 44.0%
pow1/246.0%
associate-*r*46.0%
Applied egg-rr46.0%
if -4.2000000000000001e-120 < Om < 4.9000000000000002e-49Initial program 43.1%
Simplified43.2%
Taylor expanded in U* around inf 36.4%
associate-*l/37.6%
Applied egg-rr37.6%
if 4.9000000000000002e-49 < Om Initial program 52.5%
Simplified60.9%
Taylor expanded in t around inf 44.5%
add-sqr-sqrt44.5%
pow1/244.5%
pow1/245.9%
pow-prod-down31.2%
pow231.2%
Applied egg-rr31.2%
unpow1/231.2%
associate-*l*31.2%
associate-*r*29.1%
*-commutative29.1%
associate-*r*30.3%
unpow230.3%
rem-sqrt-square42.9%
associate-*r*42.9%
*-commutative42.9%
associate-*r*46.4%
Simplified46.4%
Final simplification43.7%
l_m = (fabs.f64 l)
(FPCore (n U t l_m Om U*)
:precision binary64
(if (<= Om -3.5e-120)
(pow (* (* n t) (* 2.0 U)) 0.5)
(if (<= Om 5.4e-49)
(* (sqrt (* U U*)) (/ (* l_m (* n (sqrt 2.0))) Om))
(sqrt (fabs (* 2.0 (* n (* U t))))))))l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (Om <= -3.5e-120) {
tmp = pow(((n * t) * (2.0 * U)), 0.5);
} else if (Om <= 5.4e-49) {
tmp = sqrt((U * U_42_)) * ((l_m * (n * sqrt(2.0))) / Om);
} else {
tmp = sqrt(fabs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (om <= (-3.5d-120)) then
tmp = ((n * t) * (2.0d0 * u)) ** 0.5d0
else if (om <= 5.4d-49) then
tmp = sqrt((u * u_42)) * ((l_m * (n * sqrt(2.0d0))) / om)
else
tmp = sqrt(abs((2.0d0 * (n * (u * t)))))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (Om <= -3.5e-120) {
tmp = Math.pow(((n * t) * (2.0 * U)), 0.5);
} else if (Om <= 5.4e-49) {
tmp = Math.sqrt((U * U_42_)) * ((l_m * (n * Math.sqrt(2.0))) / Om);
} else {
tmp = Math.sqrt(Math.abs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if Om <= -3.5e-120: tmp = math.pow(((n * t) * (2.0 * U)), 0.5) elif Om <= 5.4e-49: tmp = math.sqrt((U * U_42_)) * ((l_m * (n * math.sqrt(2.0))) / Om) else: tmp = math.sqrt(math.fabs((2.0 * (n * (U * t))))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (Om <= -3.5e-120) tmp = Float64(Float64(n * t) * Float64(2.0 * U)) ^ 0.5; elseif (Om <= 5.4e-49) tmp = Float64(sqrt(Float64(U * U_42_)) * Float64(Float64(l_m * Float64(n * sqrt(2.0))) / Om)); else tmp = sqrt(abs(Float64(2.0 * Float64(n * Float64(U * t))))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (Om <= -3.5e-120) tmp = ((n * t) * (2.0 * U)) ^ 0.5; elseif (Om <= 5.4e-49) tmp = sqrt((U * U_42_)) * ((l_m * (n * sqrt(2.0))) / Om); else tmp = sqrt(abs((2.0 * (n * (U * t))))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[Om, -3.5e-120], N[Power[N[(N[(n * t), $MachinePrecision] * N[(2.0 * U), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision], If[LessEqual[Om, 5.4e-49], N[(N[Sqrt[N[(U * U$42$), $MachinePrecision]], $MachinePrecision] * N[(N[(l$95$m * N[(n * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision], N[Sqrt[N[Abs[N[(2.0 * N[(n * N[(U * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;Om \leq -3.5 \cdot 10^{-120}:\\
\;\;\;\;{\left(\left(n \cdot t\right) \cdot \left(2 \cdot U\right)\right)}^{0.5}\\
\mathbf{elif}\;Om \leq 5.4 \cdot 10^{-49}:\\
\;\;\;\;\sqrt{U \cdot U*} \cdot \frac{l\_m \cdot \left(n \cdot \sqrt{2}\right)}{Om}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left|2 \cdot \left(n \cdot \left(U \cdot t\right)\right)\right|}\\
\end{array}
\end{array}
if Om < -3.5e-120Initial program 54.4%
Simplified60.3%
Taylor expanded in t around inf 44.0%
pow1/246.0%
associate-*r*46.0%
Applied egg-rr46.0%
if -3.5e-120 < Om < 5.3999999999999999e-49Initial program 43.1%
Simplified43.2%
Taylor expanded in U* around inf 36.4%
if 5.3999999999999999e-49 < Om Initial program 52.5%
Simplified60.9%
Taylor expanded in t around inf 44.5%
add-sqr-sqrt44.5%
pow1/244.5%
pow1/245.9%
pow-prod-down31.2%
pow231.2%
Applied egg-rr31.2%
unpow1/231.2%
associate-*l*31.2%
associate-*r*29.1%
*-commutative29.1%
associate-*r*30.3%
unpow230.3%
rem-sqrt-square42.9%
associate-*r*42.9%
*-commutative42.9%
associate-*r*46.4%
Simplified46.4%
Final simplification43.3%
l_m = (fabs.f64 l)
(FPCore (n U t l_m Om U*)
:precision binary64
(if (<= Om -4.7e-122)
(pow (* (* n t) (* 2.0 U)) 0.5)
(if (<= Om 1.3e-48)
(* (sqrt (* U U*)) (* (/ (sqrt 2.0) Om) (* n l_m)))
(sqrt (fabs (* 2.0 (* n (* U t))))))))l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (Om <= -4.7e-122) {
tmp = pow(((n * t) * (2.0 * U)), 0.5);
} else if (Om <= 1.3e-48) {
tmp = sqrt((U * U_42_)) * ((sqrt(2.0) / Om) * (n * l_m));
} else {
tmp = sqrt(fabs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (om <= (-4.7d-122)) then
tmp = ((n * t) * (2.0d0 * u)) ** 0.5d0
else if (om <= 1.3d-48) then
tmp = sqrt((u * u_42)) * ((sqrt(2.0d0) / om) * (n * l_m))
else
tmp = sqrt(abs((2.0d0 * (n * (u * t)))))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (Om <= -4.7e-122) {
tmp = Math.pow(((n * t) * (2.0 * U)), 0.5);
} else if (Om <= 1.3e-48) {
tmp = Math.sqrt((U * U_42_)) * ((Math.sqrt(2.0) / Om) * (n * l_m));
} else {
tmp = Math.sqrt(Math.abs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if Om <= -4.7e-122: tmp = math.pow(((n * t) * (2.0 * U)), 0.5) elif Om <= 1.3e-48: tmp = math.sqrt((U * U_42_)) * ((math.sqrt(2.0) / Om) * (n * l_m)) else: tmp = math.sqrt(math.fabs((2.0 * (n * (U * t))))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (Om <= -4.7e-122) tmp = Float64(Float64(n * t) * Float64(2.0 * U)) ^ 0.5; elseif (Om <= 1.3e-48) tmp = Float64(sqrt(Float64(U * U_42_)) * Float64(Float64(sqrt(2.0) / Om) * Float64(n * l_m))); else tmp = sqrt(abs(Float64(2.0 * Float64(n * Float64(U * t))))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (Om <= -4.7e-122) tmp = ((n * t) * (2.0 * U)) ^ 0.5; elseif (Om <= 1.3e-48) tmp = sqrt((U * U_42_)) * ((sqrt(2.0) / Om) * (n * l_m)); else tmp = sqrt(abs((2.0 * (n * (U * t))))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[Om, -4.7e-122], N[Power[N[(N[(n * t), $MachinePrecision] * N[(2.0 * U), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision], If[LessEqual[Om, 1.3e-48], N[(N[Sqrt[N[(U * U$42$), $MachinePrecision]], $MachinePrecision] * N[(N[(N[Sqrt[2.0], $MachinePrecision] / Om), $MachinePrecision] * N[(n * l$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Sqrt[N[Abs[N[(2.0 * N[(n * N[(U * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;Om \leq -4.7 \cdot 10^{-122}:\\
\;\;\;\;{\left(\left(n \cdot t\right) \cdot \left(2 \cdot U\right)\right)}^{0.5}\\
\mathbf{elif}\;Om \leq 1.3 \cdot 10^{-48}:\\
\;\;\;\;\sqrt{U \cdot U*} \cdot \left(\frac{\sqrt{2}}{Om} \cdot \left(n \cdot l\_m\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left|2 \cdot \left(n \cdot \left(U \cdot t\right)\right)\right|}\\
\end{array}
\end{array}
if Om < -4.6999999999999999e-122Initial program 54.4%
Simplified60.3%
Taylor expanded in t around inf 44.0%
pow1/246.0%
associate-*r*46.0%
Applied egg-rr46.0%
if -4.6999999999999999e-122 < Om < 1.29999999999999994e-48Initial program 43.1%
Simplified43.2%
Taylor expanded in U* around inf 36.4%
associate-/l*34.5%
Simplified34.5%
Taylor expanded in l around 0 36.4%
associate-/l*34.5%
associate-*r/34.5%
associate-*r*36.3%
*-commutative36.3%
Simplified36.3%
if 1.29999999999999994e-48 < Om Initial program 52.5%
Simplified60.9%
Taylor expanded in t around inf 44.5%
add-sqr-sqrt44.5%
pow1/244.5%
pow1/245.9%
pow-prod-down31.2%
pow231.2%
Applied egg-rr31.2%
unpow1/231.2%
associate-*l*31.2%
associate-*r*29.1%
*-commutative29.1%
associate-*r*30.3%
unpow230.3%
rem-sqrt-square42.9%
associate-*r*42.9%
*-commutative42.9%
associate-*r*46.4%
Simplified46.4%
Final simplification43.3%
l_m = (fabs.f64 l)
(FPCore (n U t l_m Om U*)
:precision binary64
(if (<= Om -5.3e-120)
(pow (* (* n t) (* 2.0 U)) 0.5)
(if (<= Om 4.8e-49)
(* (* l_m (* n (sqrt 2.0))) (/ (sqrt (* U U*)) Om))
(sqrt (fabs (* 2.0 (* n (* U t))))))))l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (Om <= -5.3e-120) {
tmp = pow(((n * t) * (2.0 * U)), 0.5);
} else if (Om <= 4.8e-49) {
tmp = (l_m * (n * sqrt(2.0))) * (sqrt((U * U_42_)) / Om);
} else {
tmp = sqrt(fabs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (om <= (-5.3d-120)) then
tmp = ((n * t) * (2.0d0 * u)) ** 0.5d0
else if (om <= 4.8d-49) then
tmp = (l_m * (n * sqrt(2.0d0))) * (sqrt((u * u_42)) / om)
else
tmp = sqrt(abs((2.0d0 * (n * (u * t)))))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (Om <= -5.3e-120) {
tmp = Math.pow(((n * t) * (2.0 * U)), 0.5);
} else if (Om <= 4.8e-49) {
tmp = (l_m * (n * Math.sqrt(2.0))) * (Math.sqrt((U * U_42_)) / Om);
} else {
tmp = Math.sqrt(Math.abs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if Om <= -5.3e-120: tmp = math.pow(((n * t) * (2.0 * U)), 0.5) elif Om <= 4.8e-49: tmp = (l_m * (n * math.sqrt(2.0))) * (math.sqrt((U * U_42_)) / Om) else: tmp = math.sqrt(math.fabs((2.0 * (n * (U * t))))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (Om <= -5.3e-120) tmp = Float64(Float64(n * t) * Float64(2.0 * U)) ^ 0.5; elseif (Om <= 4.8e-49) tmp = Float64(Float64(l_m * Float64(n * sqrt(2.0))) * Float64(sqrt(Float64(U * U_42_)) / Om)); else tmp = sqrt(abs(Float64(2.0 * Float64(n * Float64(U * t))))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (Om <= -5.3e-120) tmp = ((n * t) * (2.0 * U)) ^ 0.5; elseif (Om <= 4.8e-49) tmp = (l_m * (n * sqrt(2.0))) * (sqrt((U * U_42_)) / Om); else tmp = sqrt(abs((2.0 * (n * (U * t))))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[Om, -5.3e-120], N[Power[N[(N[(n * t), $MachinePrecision] * N[(2.0 * U), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision], If[LessEqual[Om, 4.8e-49], N[(N[(l$95$m * N[(n * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sqrt[N[(U * U$42$), $MachinePrecision]], $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision], N[Sqrt[N[Abs[N[(2.0 * N[(n * N[(U * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;Om \leq -5.3 \cdot 10^{-120}:\\
\;\;\;\;{\left(\left(n \cdot t\right) \cdot \left(2 \cdot U\right)\right)}^{0.5}\\
\mathbf{elif}\;Om \leq 4.8 \cdot 10^{-49}:\\
\;\;\;\;\left(l\_m \cdot \left(n \cdot \sqrt{2}\right)\right) \cdot \frac{\sqrt{U \cdot U*}}{Om}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left|2 \cdot \left(n \cdot \left(U \cdot t\right)\right)\right|}\\
\end{array}
\end{array}
if Om < -5.29999999999999997e-120Initial program 54.4%
Simplified60.3%
Taylor expanded in t around inf 44.0%
pow1/246.0%
associate-*r*46.0%
Applied egg-rr46.0%
if -5.29999999999999997e-120 < Om < 4.79999999999999985e-49Initial program 43.1%
Simplified43.2%
Taylor expanded in U* around inf 32.5%
Taylor expanded in n around 0 36.4%
associate-*l/37.6%
associate-/l*36.3%
Simplified36.3%
if 4.79999999999999985e-49 < Om Initial program 52.5%
Simplified60.9%
Taylor expanded in t around inf 44.5%
add-sqr-sqrt44.5%
pow1/244.5%
pow1/245.9%
pow-prod-down31.2%
pow231.2%
Applied egg-rr31.2%
unpow1/231.2%
associate-*l*31.2%
associate-*r*29.1%
*-commutative29.1%
associate-*r*30.3%
unpow230.3%
rem-sqrt-square42.9%
associate-*r*42.9%
*-commutative42.9%
associate-*r*46.4%
Simplified46.4%
Final simplification43.3%
l_m = (fabs.f64 l)
(FPCore (n U t l_m Om U*)
:precision binary64
(if (<= l_m 2.25e+84)
(sqrt (fabs (* 2.0 (* n (* U t)))))
(if (<= l_m 2.6e+132)
(sqrt (* -4.0 (/ (* U (* n (pow l_m 2.0))) Om)))
(* l_m (* (sqrt (* U U*)) (* n (/ (sqrt 2.0) Om)))))))l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 2.25e+84) {
tmp = sqrt(fabs((2.0 * (n * (U * t)))));
} else if (l_m <= 2.6e+132) {
tmp = sqrt((-4.0 * ((U * (n * pow(l_m, 2.0))) / Om)));
} else {
tmp = l_m * (sqrt((U * U_42_)) * (n * (sqrt(2.0) / Om)));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (l_m <= 2.25d+84) then
tmp = sqrt(abs((2.0d0 * (n * (u * t)))))
else if (l_m <= 2.6d+132) then
tmp = sqrt(((-4.0d0) * ((u * (n * (l_m ** 2.0d0))) / om)))
else
tmp = l_m * (sqrt((u * u_42)) * (n * (sqrt(2.0d0) / om)))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 2.25e+84) {
tmp = Math.sqrt(Math.abs((2.0 * (n * (U * t)))));
} else if (l_m <= 2.6e+132) {
tmp = Math.sqrt((-4.0 * ((U * (n * Math.pow(l_m, 2.0))) / Om)));
} else {
tmp = l_m * (Math.sqrt((U * U_42_)) * (n * (Math.sqrt(2.0) / Om)));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if l_m <= 2.25e+84: tmp = math.sqrt(math.fabs((2.0 * (n * (U * t))))) elif l_m <= 2.6e+132: tmp = math.sqrt((-4.0 * ((U * (n * math.pow(l_m, 2.0))) / Om))) else: tmp = l_m * (math.sqrt((U * U_42_)) * (n * (math.sqrt(2.0) / Om))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (l_m <= 2.25e+84) tmp = sqrt(abs(Float64(2.0 * Float64(n * Float64(U * t))))); elseif (l_m <= 2.6e+132) tmp = sqrt(Float64(-4.0 * Float64(Float64(U * Float64(n * (l_m ^ 2.0))) / Om))); else tmp = Float64(l_m * Float64(sqrt(Float64(U * U_42_)) * Float64(n * Float64(sqrt(2.0) / Om)))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (l_m <= 2.25e+84) tmp = sqrt(abs((2.0 * (n * (U * t))))); elseif (l_m <= 2.6e+132) tmp = sqrt((-4.0 * ((U * (n * (l_m ^ 2.0))) / Om))); else tmp = l_m * (sqrt((U * U_42_)) * (n * (sqrt(2.0) / Om))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[l$95$m, 2.25e+84], N[Sqrt[N[Abs[N[(2.0 * N[(n * N[(U * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision], If[LessEqual[l$95$m, 2.6e+132], N[Sqrt[N[(-4.0 * N[(N[(U * N[(n * N[Power[l$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[(l$95$m * N[(N[Sqrt[N[(U * U$42$), $MachinePrecision]], $MachinePrecision] * N[(n * N[(N[Sqrt[2.0], $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;l\_m \leq 2.25 \cdot 10^{+84}:\\
\;\;\;\;\sqrt{\left|2 \cdot \left(n \cdot \left(U \cdot t\right)\right)\right|}\\
\mathbf{elif}\;l\_m \leq 2.6 \cdot 10^{+132}:\\
\;\;\;\;\sqrt{-4 \cdot \frac{U \cdot \left(n \cdot {l\_m}^{2}\right)}{Om}}\\
\mathbf{else}:\\
\;\;\;\;l\_m \cdot \left(\sqrt{U \cdot U*} \cdot \left(n \cdot \frac{\sqrt{2}}{Om}\right)\right)\\
\end{array}
\end{array}
if l < 2.2499999999999999e84Initial program 56.6%
Simplified58.6%
Taylor expanded in t around inf 39.6%
add-sqr-sqrt39.6%
pow1/239.6%
pow1/240.2%
pow-prod-down31.2%
pow231.2%
Applied egg-rr31.2%
unpow1/231.2%
associate-*l*31.2%
associate-*r*31.2%
*-commutative31.2%
associate-*r*31.2%
unpow231.2%
rem-sqrt-square41.7%
associate-*r*41.3%
*-commutative41.3%
associate-*r*41.0%
Simplified41.0%
if 2.2499999999999999e84 < l < 2.6e132Initial program 30.7%
Simplified30.7%
Taylor expanded in Om around inf 30.7%
Taylor expanded in l around inf 30.7%
if 2.6e132 < l Initial program 23.6%
Simplified44.6%
associate-*r*44.6%
fma-define47.0%
associate-*r*47.0%
Applied egg-rr47.0%
clear-num46.9%
inv-pow46.9%
Applied egg-rr46.9%
unpow-146.9%
Simplified46.9%
Taylor expanded in U* around inf 25.5%
associate-/l*25.6%
associate-*r/25.6%
associate-*r*29.5%
*-commutative29.5%
Simplified29.5%
Final simplification38.9%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (if (<= l_m 1.32e+154) (sqrt (* (* 2.0 n) (* U (- t (* 2.0 (/ (pow l_m 2.0) Om)))))) (sqrt (* (pow (/ l_m Om) 2.0) (* (* (* 2.0 n) U) (* n U*))))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 1.32e+154) {
tmp = sqrt(((2.0 * n) * (U * (t - (2.0 * (pow(l_m, 2.0) / Om))))));
} else {
tmp = sqrt((pow((l_m / Om), 2.0) * (((2.0 * n) * U) * (n * U_42_))));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (l_m <= 1.32d+154) then
tmp = sqrt(((2.0d0 * n) * (u * (t - (2.0d0 * ((l_m ** 2.0d0) / om))))))
else
tmp = sqrt((((l_m / om) ** 2.0d0) * (((2.0d0 * n) * u) * (n * u_42))))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 1.32e+154) {
tmp = Math.sqrt(((2.0 * n) * (U * (t - (2.0 * (Math.pow(l_m, 2.0) / Om))))));
} else {
tmp = Math.sqrt((Math.pow((l_m / Om), 2.0) * (((2.0 * n) * U) * (n * U_42_))));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if l_m <= 1.32e+154: tmp = math.sqrt(((2.0 * n) * (U * (t - (2.0 * (math.pow(l_m, 2.0) / Om)))))) else: tmp = math.sqrt((math.pow((l_m / Om), 2.0) * (((2.0 * n) * U) * (n * U_42_)))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (l_m <= 1.32e+154) tmp = sqrt(Float64(Float64(2.0 * n) * Float64(U * Float64(t - Float64(2.0 * Float64((l_m ^ 2.0) / Om)))))); else tmp = sqrt(Float64((Float64(l_m / Om) ^ 2.0) * Float64(Float64(Float64(2.0 * n) * U) * Float64(n * U_42_)))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (l_m <= 1.32e+154) tmp = sqrt(((2.0 * n) * (U * (t - (2.0 * ((l_m ^ 2.0) / Om)))))); else tmp = sqrt((((l_m / Om) ^ 2.0) * (((2.0 * n) * U) * (n * U_42_)))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[l$95$m, 1.32e+154], N[Sqrt[N[(N[(2.0 * n), $MachinePrecision] * N[(U * N[(t - N[(2.0 * N[(N[Power[l$95$m, 2.0], $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Sqrt[N[(N[Power[N[(l$95$m / Om), $MachinePrecision], 2.0], $MachinePrecision] * N[(N[(N[(2.0 * n), $MachinePrecision] * U), $MachinePrecision] * N[(n * U$42$), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;l\_m \leq 1.32 \cdot 10^{+154}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(t - 2 \cdot \frac{{l\_m}^{2}}{Om}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{{\left(\frac{l\_m}{Om}\right)}^{2} \cdot \left(\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(n \cdot U*\right)\right)}\\
\end{array}
\end{array}
if l < 1.31999999999999998e154Initial program 55.8%
Simplified58.1%
Taylor expanded in n around 0 49.2%
if 1.31999999999999998e154 < l Initial program 21.8%
associate-*r/46.0%
*-commutative46.0%
add-sqr-sqrt46.0%
associate-*r*46.0%
Applied egg-rr46.0%
Taylor expanded in U* around inf 23.6%
associate-/l*23.6%
*-commutative23.6%
associate-/l*23.6%
unpow223.6%
unpow223.6%
times-frac32.3%
unpow232.3%
Simplified32.3%
*-un-lft-identity32.3%
*-commutative32.3%
associate-*r*32.3%
Applied egg-rr32.3%
*-lft-identity32.3%
associate-*r*29.5%
*-commutative29.5%
*-commutative29.5%
Simplified29.5%
Final simplification46.1%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (if (<= l_m 3.35e+143) (sqrt (* (* 2.0 n) (* U (- t (* 2.0 (/ (pow l_m 2.0) Om)))))) (sqrt (* (* (* 2.0 n) U) (* (pow (/ l_m Om) 2.0) (* n U*))))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 3.35e+143) {
tmp = sqrt(((2.0 * n) * (U * (t - (2.0 * (pow(l_m, 2.0) / Om))))));
} else {
tmp = sqrt((((2.0 * n) * U) * (pow((l_m / Om), 2.0) * (n * U_42_))));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (l_m <= 3.35d+143) then
tmp = sqrt(((2.0d0 * n) * (u * (t - (2.0d0 * ((l_m ** 2.0d0) / om))))))
else
tmp = sqrt((((2.0d0 * n) * u) * (((l_m / om) ** 2.0d0) * (n * u_42))))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 3.35e+143) {
tmp = Math.sqrt(((2.0 * n) * (U * (t - (2.0 * (Math.pow(l_m, 2.0) / Om))))));
} else {
tmp = Math.sqrt((((2.0 * n) * U) * (Math.pow((l_m / Om), 2.0) * (n * U_42_))));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if l_m <= 3.35e+143: tmp = math.sqrt(((2.0 * n) * (U * (t - (2.0 * (math.pow(l_m, 2.0) / Om)))))) else: tmp = math.sqrt((((2.0 * n) * U) * (math.pow((l_m / Om), 2.0) * (n * U_42_)))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (l_m <= 3.35e+143) tmp = sqrt(Float64(Float64(2.0 * n) * Float64(U * Float64(t - Float64(2.0 * Float64((l_m ^ 2.0) / Om)))))); else tmp = sqrt(Float64(Float64(Float64(2.0 * n) * U) * Float64((Float64(l_m / Om) ^ 2.0) * Float64(n * U_42_)))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (l_m <= 3.35e+143) tmp = sqrt(((2.0 * n) * (U * (t - (2.0 * ((l_m ^ 2.0) / Om)))))); else tmp = sqrt((((2.0 * n) * U) * (((l_m / Om) ^ 2.0) * (n * U_42_)))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[l$95$m, 3.35e+143], N[Sqrt[N[(N[(2.0 * n), $MachinePrecision] * N[(U * N[(t - N[(2.0 * N[(N[Power[l$95$m, 2.0], $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Sqrt[N[(N[(N[(2.0 * n), $MachinePrecision] * U), $MachinePrecision] * N[(N[Power[N[(l$95$m / Om), $MachinePrecision], 2.0], $MachinePrecision] * N[(n * U$42$), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;l\_m \leq 3.35 \cdot 10^{+143}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(t - 2 \cdot \frac{{l\_m}^{2}}{Om}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left({\left(\frac{l\_m}{Om}\right)}^{2} \cdot \left(n \cdot U*\right)\right)}\\
\end{array}
\end{array}
if l < 3.3500000000000001e143Initial program 56.3%
Simplified58.6%
Taylor expanded in n around 0 49.6%
if 3.3500000000000001e143 < l Initial program 21.0%
associate-*r/44.0%
*-commutative44.0%
add-sqr-sqrt44.0%
associate-*r*44.0%
Applied egg-rr44.0%
Taylor expanded in U* around inf 22.7%
associate-/l*22.7%
*-commutative22.7%
associate-/l*22.7%
unpow222.7%
unpow222.7%
times-frac31.0%
unpow231.0%
Simplified31.0%
pow131.0%
associate-*r*31.0%
Applied egg-rr31.0%
unpow131.0%
*-commutative31.0%
*-commutative31.0%
Simplified31.0%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (if (<= l_m 3.35e+143) (sqrt (* (* 2.0 n) (* U (- t (* 2.0 (/ (pow l_m 2.0) Om)))))) (sqrt (* (* (* 2.0 n) U) (* (* n (pow (/ l_m Om) 2.0)) U*)))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 3.35e+143) {
tmp = sqrt(((2.0 * n) * (U * (t - (2.0 * (pow(l_m, 2.0) / Om))))));
} else {
tmp = sqrt((((2.0 * n) * U) * ((n * pow((l_m / Om), 2.0)) * U_42_)));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (l_m <= 3.35d+143) then
tmp = sqrt(((2.0d0 * n) * (u * (t - (2.0d0 * ((l_m ** 2.0d0) / om))))))
else
tmp = sqrt((((2.0d0 * n) * u) * ((n * ((l_m / om) ** 2.0d0)) * u_42)))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 3.35e+143) {
tmp = Math.sqrt(((2.0 * n) * (U * (t - (2.0 * (Math.pow(l_m, 2.0) / Om))))));
} else {
tmp = Math.sqrt((((2.0 * n) * U) * ((n * Math.pow((l_m / Om), 2.0)) * U_42_)));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if l_m <= 3.35e+143: tmp = math.sqrt(((2.0 * n) * (U * (t - (2.0 * (math.pow(l_m, 2.0) / Om)))))) else: tmp = math.sqrt((((2.0 * n) * U) * ((n * math.pow((l_m / Om), 2.0)) * U_42_))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (l_m <= 3.35e+143) tmp = sqrt(Float64(Float64(2.0 * n) * Float64(U * Float64(t - Float64(2.0 * Float64((l_m ^ 2.0) / Om)))))); else tmp = sqrt(Float64(Float64(Float64(2.0 * n) * U) * Float64(Float64(n * (Float64(l_m / Om) ^ 2.0)) * U_42_))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (l_m <= 3.35e+143) tmp = sqrt(((2.0 * n) * (U * (t - (2.0 * ((l_m ^ 2.0) / Om)))))); else tmp = sqrt((((2.0 * n) * U) * ((n * ((l_m / Om) ^ 2.0)) * U_42_))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[l$95$m, 3.35e+143], N[Sqrt[N[(N[(2.0 * n), $MachinePrecision] * N[(U * N[(t - N[(2.0 * N[(N[Power[l$95$m, 2.0], $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Sqrt[N[(N[(N[(2.0 * n), $MachinePrecision] * U), $MachinePrecision] * N[(N[(n * N[Power[N[(l$95$m / Om), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * U$42$), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;l\_m \leq 3.35 \cdot 10^{+143}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(t - 2 \cdot \frac{{l\_m}^{2}}{Om}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(n \cdot {\left(\frac{l\_m}{Om}\right)}^{2}\right) \cdot U*\right)}\\
\end{array}
\end{array}
if l < 3.3500000000000001e143Initial program 56.3%
Simplified58.6%
Taylor expanded in n around 0 49.6%
if 3.3500000000000001e143 < l Initial program 21.0%
associate-*r/44.0%
*-commutative44.0%
add-sqr-sqrt44.0%
associate-*r*44.0%
Applied egg-rr44.0%
Taylor expanded in U* around inf 22.7%
associate-/l*22.7%
*-commutative22.7%
associate-/l*22.7%
unpow222.7%
unpow222.7%
times-frac31.0%
unpow231.0%
Simplified31.0%
Final simplification46.5%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (if (<= l_m 3.35e+143) (sqrt (* (* 2.0 n) (* U (- t (* 2.0 (/ (pow l_m 2.0) Om)))))) (* l_m (* (sqrt (* U U*)) (* n (/ (sqrt 2.0) Om))))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 3.35e+143) {
tmp = sqrt(((2.0 * n) * (U * (t - (2.0 * (pow(l_m, 2.0) / Om))))));
} else {
tmp = l_m * (sqrt((U * U_42_)) * (n * (sqrt(2.0) / Om)));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (l_m <= 3.35d+143) then
tmp = sqrt(((2.0d0 * n) * (u * (t - (2.0d0 * ((l_m ** 2.0d0) / om))))))
else
tmp = l_m * (sqrt((u * u_42)) * (n * (sqrt(2.0d0) / om)))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 3.35e+143) {
tmp = Math.sqrt(((2.0 * n) * (U * (t - (2.0 * (Math.pow(l_m, 2.0) / Om))))));
} else {
tmp = l_m * (Math.sqrt((U * U_42_)) * (n * (Math.sqrt(2.0) / Om)));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if l_m <= 3.35e+143: tmp = math.sqrt(((2.0 * n) * (U * (t - (2.0 * (math.pow(l_m, 2.0) / Om)))))) else: tmp = l_m * (math.sqrt((U * U_42_)) * (n * (math.sqrt(2.0) / Om))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (l_m <= 3.35e+143) tmp = sqrt(Float64(Float64(2.0 * n) * Float64(U * Float64(t - Float64(2.0 * Float64((l_m ^ 2.0) / Om)))))); else tmp = Float64(l_m * Float64(sqrt(Float64(U * U_42_)) * Float64(n * Float64(sqrt(2.0) / Om)))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (l_m <= 3.35e+143) tmp = sqrt(((2.0 * n) * (U * (t - (2.0 * ((l_m ^ 2.0) / Om)))))); else tmp = l_m * (sqrt((U * U_42_)) * (n * (sqrt(2.0) / Om))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[l$95$m, 3.35e+143], N[Sqrt[N[(N[(2.0 * n), $MachinePrecision] * N[(U * N[(t - N[(2.0 * N[(N[Power[l$95$m, 2.0], $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[(l$95$m * N[(N[Sqrt[N[(U * U$42$), $MachinePrecision]], $MachinePrecision] * N[(n * N[(N[Sqrt[2.0], $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;l\_m \leq 3.35 \cdot 10^{+143}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(t - 2 \cdot \frac{{l\_m}^{2}}{Om}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;l\_m \cdot \left(\sqrt{U \cdot U*} \cdot \left(n \cdot \frac{\sqrt{2}}{Om}\right)\right)\\
\end{array}
\end{array}
if l < 3.3500000000000001e143Initial program 56.3%
Simplified58.6%
Taylor expanded in n around 0 49.6%
if 3.3500000000000001e143 < l Initial program 21.0%
Simplified44.0%
associate-*r*44.0%
fma-define46.7%
associate-*r*46.6%
Applied egg-rr46.6%
clear-num46.5%
inv-pow46.5%
Applied egg-rr46.5%
unpow-146.5%
Simplified46.5%
Taylor expanded in U* around inf 24.9%
associate-/l*25.0%
associate-*r/25.0%
associate-*r*29.3%
*-commutative29.3%
Simplified29.3%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (if (<= t 9.5e+196) (sqrt (* 2.0 (* U (* n (- t (* 2.0 (/ (pow l_m 2.0) Om))))))) (sqrt (fabs (* 2.0 (* n (* U t)))))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (t <= 9.5e+196) {
tmp = sqrt((2.0 * (U * (n * (t - (2.0 * (pow(l_m, 2.0) / Om)))))));
} else {
tmp = sqrt(fabs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (t <= 9.5d+196) then
tmp = sqrt((2.0d0 * (u * (n * (t - (2.0d0 * ((l_m ** 2.0d0) / om)))))))
else
tmp = sqrt(abs((2.0d0 * (n * (u * t)))))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (t <= 9.5e+196) {
tmp = Math.sqrt((2.0 * (U * (n * (t - (2.0 * (Math.pow(l_m, 2.0) / Om)))))));
} else {
tmp = Math.sqrt(Math.abs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if t <= 9.5e+196: tmp = math.sqrt((2.0 * (U * (n * (t - (2.0 * (math.pow(l_m, 2.0) / Om))))))) else: tmp = math.sqrt(math.fabs((2.0 * (n * (U * t))))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (t <= 9.5e+196) tmp = sqrt(Float64(2.0 * Float64(U * Float64(n * Float64(t - Float64(2.0 * Float64((l_m ^ 2.0) / Om))))))); else tmp = sqrt(abs(Float64(2.0 * Float64(n * Float64(U * t))))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (t <= 9.5e+196) tmp = sqrt((2.0 * (U * (n * (t - (2.0 * ((l_m ^ 2.0) / Om))))))); else tmp = sqrt(abs((2.0 * (n * (U * t))))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[t, 9.5e+196], N[Sqrt[N[(2.0 * N[(U * N[(n * N[(t - N[(2.0 * N[(N[Power[l$95$m, 2.0], $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Sqrt[N[Abs[N[(2.0 * N[(n * N[(U * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 9.5 \cdot 10^{+196}:\\
\;\;\;\;\sqrt{2 \cdot \left(U \cdot \left(n \cdot \left(t - 2 \cdot \frac{{l\_m}^{2}}{Om}\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left|2 \cdot \left(n \cdot \left(U \cdot t\right)\right)\right|}\\
\end{array}
\end{array}
if t < 9.5000000000000004e196Initial program 51.8%
Simplified56.2%
Taylor expanded in n around 0 45.9%
if 9.5000000000000004e196 < t Initial program 40.8%
Simplified50.6%
Taylor expanded in t around inf 53.9%
add-sqr-sqrt53.9%
pow1/253.9%
pow1/254.3%
pow-prod-down39.3%
pow239.3%
Applied egg-rr39.3%
unpow1/239.3%
associate-*l*39.3%
associate-*r*36.5%
*-commutative36.5%
associate-*r*36.7%
unpow236.7%
rem-sqrt-square49.2%
associate-*r*51.6%
*-commutative51.6%
associate-*r*54.6%
Simplified54.6%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (if (<= l_m 2.2e+84) (sqrt (fabs (* 2.0 (* n (* U t))))) (sqrt (* -4.0 (/ (* U (* n (pow l_m 2.0))) Om)))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 2.2e+84) {
tmp = sqrt(fabs((2.0 * (n * (U * t)))));
} else {
tmp = sqrt((-4.0 * ((U * (n * pow(l_m, 2.0))) / Om)));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (l_m <= 2.2d+84) then
tmp = sqrt(abs((2.0d0 * (n * (u * t)))))
else
tmp = sqrt(((-4.0d0) * ((u * (n * (l_m ** 2.0d0))) / om)))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (l_m <= 2.2e+84) {
tmp = Math.sqrt(Math.abs((2.0 * (n * (U * t)))));
} else {
tmp = Math.sqrt((-4.0 * ((U * (n * Math.pow(l_m, 2.0))) / Om)));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if l_m <= 2.2e+84: tmp = math.sqrt(math.fabs((2.0 * (n * (U * t))))) else: tmp = math.sqrt((-4.0 * ((U * (n * math.pow(l_m, 2.0))) / Om))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (l_m <= 2.2e+84) tmp = sqrt(abs(Float64(2.0 * Float64(n * Float64(U * t))))); else tmp = sqrt(Float64(-4.0 * Float64(Float64(U * Float64(n * (l_m ^ 2.0))) / Om))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (l_m <= 2.2e+84) tmp = sqrt(abs((2.0 * (n * (U * t))))); else tmp = sqrt((-4.0 * ((U * (n * (l_m ^ 2.0))) / Om))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[l$95$m, 2.2e+84], N[Sqrt[N[Abs[N[(2.0 * N[(n * N[(U * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision], N[Sqrt[N[(-4.0 * N[(N[(U * N[(n * N[Power[l$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / Om), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;l\_m \leq 2.2 \cdot 10^{+84}:\\
\;\;\;\;\sqrt{\left|2 \cdot \left(n \cdot \left(U \cdot t\right)\right)\right|}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{-4 \cdot \frac{U \cdot \left(n \cdot {l\_m}^{2}\right)}{Om}}\\
\end{array}
\end{array}
if l < 2.1999999999999998e84Initial program 56.6%
Simplified58.6%
Taylor expanded in t around inf 39.6%
add-sqr-sqrt39.6%
pow1/239.6%
pow1/240.2%
pow-prod-down31.2%
pow231.2%
Applied egg-rr31.2%
unpow1/231.2%
associate-*l*31.2%
associate-*r*31.2%
*-commutative31.2%
associate-*r*31.2%
unpow231.2%
rem-sqrt-square41.7%
associate-*r*41.3%
*-commutative41.3%
associate-*r*41.0%
Simplified41.0%
if 2.1999999999999998e84 < l Initial program 23.9%
Simplified42.1%
Taylor expanded in Om around inf 22.1%
Taylor expanded in l around inf 22.2%
Final simplification37.5%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (if (<= U 1.12e-237) (sqrt (fabs (* t (* 2.0 (* n U))))) (* (sqrt (* 2.0 U)) (sqrt (* n t)))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (U <= 1.12e-237) {
tmp = sqrt(fabs((t * (2.0 * (n * U)))));
} else {
tmp = sqrt((2.0 * U)) * sqrt((n * t));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (u <= 1.12d-237) then
tmp = sqrt(abs((t * (2.0d0 * (n * u)))))
else
tmp = sqrt((2.0d0 * u)) * sqrt((n * t))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (U <= 1.12e-237) {
tmp = Math.sqrt(Math.abs((t * (2.0 * (n * U)))));
} else {
tmp = Math.sqrt((2.0 * U)) * Math.sqrt((n * t));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if U <= 1.12e-237: tmp = math.sqrt(math.fabs((t * (2.0 * (n * U))))) else: tmp = math.sqrt((2.0 * U)) * math.sqrt((n * t)) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (U <= 1.12e-237) tmp = sqrt(abs(Float64(t * Float64(2.0 * Float64(n * U))))); else tmp = Float64(sqrt(Float64(2.0 * U)) * sqrt(Float64(n * t))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (U <= 1.12e-237) tmp = sqrt(abs((t * (2.0 * (n * U))))); else tmp = sqrt((2.0 * U)) * sqrt((n * t)); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[U, 1.12e-237], N[Sqrt[N[Abs[N[(t * N[(2.0 * N[(n * U), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision], N[(N[Sqrt[N[(2.0 * U), $MachinePrecision]], $MachinePrecision] * N[Sqrt[N[(n * t), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;U \leq 1.12 \cdot 10^{-237}:\\
\;\;\;\;\sqrt{\left|t \cdot \left(2 \cdot \left(n \cdot U\right)\right)\right|}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2 \cdot U} \cdot \sqrt{n \cdot t}\\
\end{array}
\end{array}
if U < 1.12000000000000002e-237Initial program 49.8%
Simplified53.5%
Taylor expanded in t around inf 32.3%
add-sqr-sqrt32.3%
pow1/232.3%
pow1/234.6%
pow-prod-down27.8%
pow227.8%
associate-*r*27.8%
Applied egg-rr27.8%
unpow1/227.8%
unpow227.8%
rem-sqrt-square35.6%
associate-*r*36.5%
associate-*r*36.5%
Simplified36.5%
if 1.12000000000000002e-237 < U Initial program 51.3%
Simplified57.7%
Taylor expanded in t around inf 37.7%
pow1/238.6%
associate-*r*38.6%
unpow-prod-down44.5%
pow1/243.6%
Applied egg-rr43.6%
unpow1/243.6%
Simplified43.6%
Final simplification39.9%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (if (<= t -9.8e-31) (pow (* (* n t) (* 2.0 U)) 0.5) (sqrt (fabs (* 2.0 (* n (* U t)))))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (t <= -9.8e-31) {
tmp = pow(((n * t) * (2.0 * U)), 0.5);
} else {
tmp = sqrt(fabs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (t <= (-9.8d-31)) then
tmp = ((n * t) * (2.0d0 * u)) ** 0.5d0
else
tmp = sqrt(abs((2.0d0 * (n * (u * t)))))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (t <= -9.8e-31) {
tmp = Math.pow(((n * t) * (2.0 * U)), 0.5);
} else {
tmp = Math.sqrt(Math.abs((2.0 * (n * (U * t)))));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if t <= -9.8e-31: tmp = math.pow(((n * t) * (2.0 * U)), 0.5) else: tmp = math.sqrt(math.fabs((2.0 * (n * (U * t))))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (t <= -9.8e-31) tmp = Float64(Float64(n * t) * Float64(2.0 * U)) ^ 0.5; else tmp = sqrt(abs(Float64(2.0 * Float64(n * Float64(U * t))))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (t <= -9.8e-31) tmp = ((n * t) * (2.0 * U)) ^ 0.5; else tmp = sqrt(abs((2.0 * (n * (U * t))))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[t, -9.8e-31], N[Power[N[(N[(n * t), $MachinePrecision] * N[(2.0 * U), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision], N[Sqrt[N[Abs[N[(2.0 * N[(n * N[(U * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq -9.8 \cdot 10^{-31}:\\
\;\;\;\;{\left(\left(n \cdot t\right) \cdot \left(2 \cdot U\right)\right)}^{0.5}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left|2 \cdot \left(n \cdot \left(U \cdot t\right)\right)\right|}\\
\end{array}
\end{array}
if t < -9.80000000000000047e-31Initial program 54.4%
Simplified52.8%
Taylor expanded in t around inf 47.3%
pow1/250.6%
associate-*r*50.6%
Applied egg-rr50.6%
if -9.80000000000000047e-31 < t Initial program 49.2%
Simplified56.3%
Taylor expanded in t around inf 33.2%
add-sqr-sqrt33.2%
pow1/233.2%
pow1/233.7%
pow-prod-down26.6%
pow226.6%
Applied egg-rr26.6%
unpow1/226.6%
associate-*l*26.6%
associate-*r*26.1%
*-commutative26.1%
associate-*r*25.7%
unpow225.7%
rem-sqrt-square33.1%
associate-*r*32.7%
*-commutative32.7%
associate-*r*34.5%
Simplified34.5%
Final simplification38.4%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (if (<= U -3.6e-69) (pow (* (* (* 2.0 n) U) t) 0.5) (sqrt (* (* 2.0 n) (* U t)))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (U <= -3.6e-69) {
tmp = pow((((2.0 * n) * U) * t), 0.5);
} else {
tmp = sqrt(((2.0 * n) * (U * t)));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (u <= (-3.6d-69)) then
tmp = (((2.0d0 * n) * u) * t) ** 0.5d0
else
tmp = sqrt(((2.0d0 * n) * (u * t)))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (U <= -3.6e-69) {
tmp = Math.pow((((2.0 * n) * U) * t), 0.5);
} else {
tmp = Math.sqrt(((2.0 * n) * (U * t)));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if U <= -3.6e-69: tmp = math.pow((((2.0 * n) * U) * t), 0.5) else: tmp = math.sqrt(((2.0 * n) * (U * t))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (U <= -3.6e-69) tmp = Float64(Float64(Float64(2.0 * n) * U) * t) ^ 0.5; else tmp = sqrt(Float64(Float64(2.0 * n) * Float64(U * t))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (U <= -3.6e-69) tmp = (((2.0 * n) * U) * t) ^ 0.5; else tmp = sqrt(((2.0 * n) * (U * t))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[U, -3.6e-69], N[Power[N[(N[(N[(2.0 * n), $MachinePrecision] * U), $MachinePrecision] * t), $MachinePrecision], 0.5], $MachinePrecision], N[Sqrt[N[(N[(2.0 * n), $MachinePrecision] * N[(U * t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;U \leq -3.6 \cdot 10^{-69}:\\
\;\;\;\;{\left(\left(\left(2 \cdot n\right) \cdot U\right) \cdot t\right)}^{0.5}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot t\right)}\\
\end{array}
\end{array}
if U < -3.60000000000000018e-69Initial program 54.5%
Simplified62.2%
associate-*r*62.2%
fma-define63.8%
associate-*r*59.0%
Applied egg-rr59.0%
clear-num59.1%
inv-pow59.1%
Applied egg-rr59.1%
unpow-159.1%
Simplified59.1%
Taylor expanded in t around inf 37.6%
pow1/240.9%
associate-*r*40.9%
*-commutative40.9%
*-commutative40.9%
Applied egg-rr40.9%
if -3.60000000000000018e-69 < U Initial program 49.2%
Simplified57.9%
Taylor expanded in t around inf 35.9%
Final simplification37.1%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (if (<= t -5e-152) (sqrt (* 2.0 (* U (* n t)))) (sqrt (* (* 2.0 n) (* U t)))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (t <= -5e-152) {
tmp = sqrt((2.0 * (U * (n * t))));
} else {
tmp = sqrt(((2.0 * n) * (U * t)));
}
return tmp;
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
real(8) :: tmp
if (t <= (-5d-152)) then
tmp = sqrt((2.0d0 * (u * (n * t))))
else
tmp = sqrt(((2.0d0 * n) * (u * t)))
end if
code = tmp
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
double tmp;
if (t <= -5e-152) {
tmp = Math.sqrt((2.0 * (U * (n * t))));
} else {
tmp = Math.sqrt(((2.0 * n) * (U * t)));
}
return tmp;
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): tmp = 0 if t <= -5e-152: tmp = math.sqrt((2.0 * (U * (n * t)))) else: tmp = math.sqrt(((2.0 * n) * (U * t))) return tmp
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) tmp = 0.0 if (t <= -5e-152) tmp = sqrt(Float64(2.0 * Float64(U * Float64(n * t)))); else tmp = sqrt(Float64(Float64(2.0 * n) * Float64(U * t))); end return tmp end
l_m = abs(l); function tmp_2 = code(n, U, t, l_m, Om, U_42_) tmp = 0.0; if (t <= -5e-152) tmp = sqrt((2.0 * (U * (n * t)))); else tmp = sqrt(((2.0 * n) * (U * t))); end tmp_2 = tmp; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := If[LessEqual[t, -5e-152], N[Sqrt[N[(2.0 * N[(U * N[(n * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Sqrt[N[(N[(2.0 * n), $MachinePrecision] * N[(U * t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq -5 \cdot 10^{-152}:\\
\;\;\;\;\sqrt{2 \cdot \left(U \cdot \left(n \cdot t\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot t\right)}\\
\end{array}
\end{array}
if t < -4.9999999999999997e-152Initial program 49.4%
Simplified53.1%
Taylor expanded in t around inf 43.7%
if -4.9999999999999997e-152 < t Initial program 50.9%
Simplified56.5%
Taylor expanded in t around inf 33.4%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (pow (* (* n t) (* 2.0 U)) 0.5))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
return pow(((n * t) * (2.0 * U)), 0.5);
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
code = ((n * t) * (2.0d0 * u)) ** 0.5d0
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
return Math.pow(((n * t) * (2.0 * U)), 0.5);
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): return math.pow(((n * t) * (2.0 * U)), 0.5)
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) return Float64(Float64(n * t) * Float64(2.0 * U)) ^ 0.5 end
l_m = abs(l); function tmp = code(n, U, t, l_m, Om, U_42_) tmp = ((n * t) * (2.0 * U)) ^ 0.5; end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := N[Power[N[(N[(n * t), $MachinePrecision] * N[(2.0 * U), $MachinePrecision]), $MachinePrecision], 0.5], $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
{\left(\left(n \cdot t\right) \cdot \left(2 \cdot U\right)\right)}^{0.5}
\end{array}
Initial program 50.5%
Simplified55.5%
Taylor expanded in t around inf 34.9%
pow1/236.5%
associate-*r*36.5%
Applied egg-rr36.5%
Final simplification36.5%
l_m = (fabs.f64 l) (FPCore (n U t l_m Om U*) :precision binary64 (sqrt (* 2.0 (* U (* n t)))))
l_m = fabs(l);
double code(double n, double U, double t, double l_m, double Om, double U_42_) {
return sqrt((2.0 * (U * (n * t))));
}
l_m = abs(l)
real(8) function code(n, u, t, l_m, om, u_42)
real(8), intent (in) :: n
real(8), intent (in) :: u
real(8), intent (in) :: t
real(8), intent (in) :: l_m
real(8), intent (in) :: om
real(8), intent (in) :: u_42
code = sqrt((2.0d0 * (u * (n * t))))
end function
l_m = Math.abs(l);
public static double code(double n, double U, double t, double l_m, double Om, double U_42_) {
return Math.sqrt((2.0 * (U * (n * t))));
}
l_m = math.fabs(l) def code(n, U, t, l_m, Om, U_42_): return math.sqrt((2.0 * (U * (n * t))))
l_m = abs(l) function code(n, U, t, l_m, Om, U_42_) return sqrt(Float64(2.0 * Float64(U * Float64(n * t)))) end
l_m = abs(l); function tmp = code(n, U, t, l_m, Om, U_42_) tmp = sqrt((2.0 * (U * (n * t)))); end
l_m = N[Abs[l], $MachinePrecision] code[n_, U_, t_, l$95$m_, Om_, U$42$_] := N[Sqrt[N[(2.0 * N[(U * N[(n * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
l_m = \left|\ell\right|
\\
\sqrt{2 \cdot \left(U \cdot \left(n \cdot t\right)\right)}
\end{array}
Initial program 50.5%
Simplified55.5%
Taylor expanded in t around inf 34.9%
herbie shell --seed 2024156
(FPCore (n U t l Om U*)
:name "Toniolo and Linder, Equation (13)"
:precision binary64
(sqrt (* (* (* 2.0 n) U) (- (- t (* 2.0 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2.0)) (- U U*))))))