
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (- (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) - 1.0));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t) ** 2.0d0)) - 1.0d0))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t), 2.0)) - 1.0));
}
def code(t, l, k): return 2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * ((1.0 + math.pow((k / t), 2.0)) - 1.0))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) - 1.0))) end
function tmp = code(t, l, k) tmp = 2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t) ^ 2.0)) - 1.0)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 23 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (- (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) - 1.0));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t) ** 2.0d0)) - 1.0d0))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t), 2.0)) - 1.0));
}
def code(t, l, k): return 2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * ((1.0 + math.pow((k / t), 2.0)) - 1.0))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) - 1.0))) end
function tmp = code(t, l, k) tmp = 2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t) ^ 2.0)) - 1.0)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}
\end{array}
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (* (/ t_m (pow (cbrt l) 2.0)) (cbrt (* (sin k_m) (tan k_m))))))
(*
t_s
(if (<= k_m 2e-24)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0))
(if (<= k_m 1.52e+147)
(/
2.0
(/
(* (pow k_m 2.0) (* t_m (pow (sin k_m) 2.0)))
(* (cos k_m) (pow l 2.0))))
(*
(/ (/ (* t_m (sqrt 2.0)) k_m) (pow t_2 2.0))
(/ (exp (* (log 2.0) 0.5)) (* t_2 (/ k_m t_m)))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = (t_m / pow(cbrt(l), 2.0)) * cbrt((sin(k_m) * tan(k_m)));
double tmp;
if (k_m <= 2e-24) {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
} else if (k_m <= 1.52e+147) {
tmp = 2.0 / ((pow(k_m, 2.0) * (t_m * pow(sin(k_m), 2.0))) / (cos(k_m) * pow(l, 2.0)));
} else {
tmp = (((t_m * sqrt(2.0)) / k_m) / pow(t_2, 2.0)) * (exp((log(2.0) * 0.5)) / (t_2 * (k_m / t_m)));
}
return t_s * tmp;
}
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = (t_m / Math.pow(Math.cbrt(l), 2.0)) * Math.cbrt((Math.sin(k_m) * Math.tan(k_m)));
double tmp;
if (k_m <= 2e-24) {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else if (k_m <= 1.52e+147) {
tmp = 2.0 / ((Math.pow(k_m, 2.0) * (t_m * Math.pow(Math.sin(k_m), 2.0))) / (Math.cos(k_m) * Math.pow(l, 2.0)));
} else {
tmp = (((t_m * Math.sqrt(2.0)) / k_m) / Math.pow(t_2, 2.0)) * (Math.exp((Math.log(2.0) * 0.5)) / (t_2 * (k_m / t_m)));
}
return t_s * tmp;
}
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = Float64(Float64(t_m / (cbrt(l) ^ 2.0)) * cbrt(Float64(sin(k_m) * tan(k_m)))) tmp = 0.0 if (k_m <= 2e-24) tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); elseif (k_m <= 1.52e+147) tmp = Float64(2.0 / Float64(Float64((k_m ^ 2.0) * Float64(t_m * (sin(k_m) ^ 2.0))) / Float64(cos(k_m) * (l ^ 2.0)))); else tmp = Float64(Float64(Float64(Float64(t_m * sqrt(2.0)) / k_m) / (t_2 ^ 2.0)) * Float64(exp(Float64(log(2.0) * 0.5)) / Float64(t_2 * Float64(k_m / t_m)))); end return Float64(t_s * tmp) end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[(N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 2e-24], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 1.52e+147], N[(2.0 / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(t$95$m * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision] / N[Power[t$95$2, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[N[(N[Log[2.0], $MachinePrecision] * 0.5), $MachinePrecision]], $MachinePrecision] / N[(t$95$2 * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \frac{t\_m}{{\left(\sqrt[3]{\ell}\right)}^{2}} \cdot \sqrt[3]{\sin k\_m \cdot \tan k\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 2 \cdot 10^{-24}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 1.52 \cdot 10^{+147}:\\
\;\;\;\;\frac{2}{\frac{{k\_m}^{2} \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}{\cos k\_m \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{t\_m \cdot \sqrt{2}}{k\_m}}{{t\_2}^{2}} \cdot \frac{e^{\log 2 \cdot 0.5}}{t\_2 \cdot \frac{k\_m}{t\_m}}\\
\end{array}
\end{array}
\end{array}
if k < 1.99999999999999985e-24Initial program 29.0%
*-commutative29.0%
associate-/r*29.0%
Simplified34.2%
Applied egg-rr28.3%
Taylor expanded in t around 0 48.6%
associate-/l*49.0%
Simplified49.0%
div-inv49.0%
unpow-prod-down49.0%
pow249.0%
pow1/249.0%
pow1/249.0%
pow-prod-up49.0%
metadata-eval49.0%
metadata-eval49.0%
Applied egg-rr49.0%
if 1.99999999999999985e-24 < k < 1.51999999999999992e147Initial program 22.6%
Simplified22.6%
Taylor expanded in t around 0 83.6%
if 1.51999999999999992e147 < k Initial program 38.1%
*-commutative38.1%
associate-/r*38.1%
Simplified48.0%
add-sqr-sqrt48.0%
add-cube-cbrt48.0%
times-frac48.0%
Applied egg-rr90.1%
associate-/r/90.0%
associate-/l/90.0%
Simplified90.0%
associate-*l/90.1%
Applied egg-rr90.1%
pow1/290.1%
pow-to-exp90.1%
Applied egg-rr90.1%
Final simplification57.3%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (* (/ t_m (pow (cbrt l) 2.0)) (cbrt (* (sin k_m) (tan k_m))))))
(*
t_s
(if (<= k_m 3.7e-25)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0))
(if (<= k_m 1.8e+147)
(/
2.0
(/
(* (pow k_m 2.0) (* t_m (pow (sin k_m) 2.0)))
(* (cos k_m) (pow l 2.0))))
(*
(/ (/ (* t_m (sqrt 2.0)) k_m) (pow t_2 2.0))
(/ (sqrt 2.0) (* t_2 (/ k_m t_m)))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = (t_m / pow(cbrt(l), 2.0)) * cbrt((sin(k_m) * tan(k_m)));
double tmp;
if (k_m <= 3.7e-25) {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
} else if (k_m <= 1.8e+147) {
tmp = 2.0 / ((pow(k_m, 2.0) * (t_m * pow(sin(k_m), 2.0))) / (cos(k_m) * pow(l, 2.0)));
} else {
tmp = (((t_m * sqrt(2.0)) / k_m) / pow(t_2, 2.0)) * (sqrt(2.0) / (t_2 * (k_m / t_m)));
}
return t_s * tmp;
}
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = (t_m / Math.pow(Math.cbrt(l), 2.0)) * Math.cbrt((Math.sin(k_m) * Math.tan(k_m)));
double tmp;
if (k_m <= 3.7e-25) {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else if (k_m <= 1.8e+147) {
tmp = 2.0 / ((Math.pow(k_m, 2.0) * (t_m * Math.pow(Math.sin(k_m), 2.0))) / (Math.cos(k_m) * Math.pow(l, 2.0)));
} else {
tmp = (((t_m * Math.sqrt(2.0)) / k_m) / Math.pow(t_2, 2.0)) * (Math.sqrt(2.0) / (t_2 * (k_m / t_m)));
}
return t_s * tmp;
}
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = Float64(Float64(t_m / (cbrt(l) ^ 2.0)) * cbrt(Float64(sin(k_m) * tan(k_m)))) tmp = 0.0 if (k_m <= 3.7e-25) tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); elseif (k_m <= 1.8e+147) tmp = Float64(2.0 / Float64(Float64((k_m ^ 2.0) * Float64(t_m * (sin(k_m) ^ 2.0))) / Float64(cos(k_m) * (l ^ 2.0)))); else tmp = Float64(Float64(Float64(Float64(t_m * sqrt(2.0)) / k_m) / (t_2 ^ 2.0)) * Float64(sqrt(2.0) / Float64(t_2 * Float64(k_m / t_m)))); end return Float64(t_s * tmp) end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[(N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 3.7e-25], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 1.8e+147], N[(2.0 / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(t$95$m * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision] / k$95$m), $MachinePrecision] / N[Power[t$95$2, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Sqrt[2.0], $MachinePrecision] / N[(t$95$2 * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \frac{t\_m}{{\left(\sqrt[3]{\ell}\right)}^{2}} \cdot \sqrt[3]{\sin k\_m \cdot \tan k\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 3.7 \cdot 10^{-25}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 1.8 \cdot 10^{+147}:\\
\;\;\;\;\frac{2}{\frac{{k\_m}^{2} \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}{\cos k\_m \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{t\_m \cdot \sqrt{2}}{k\_m}}{{t\_2}^{2}} \cdot \frac{\sqrt{2}}{t\_2 \cdot \frac{k\_m}{t\_m}}\\
\end{array}
\end{array}
\end{array}
if k < 3.70000000000000009e-25Initial program 29.0%
*-commutative29.0%
associate-/r*29.0%
Simplified34.2%
Applied egg-rr28.3%
Taylor expanded in t around 0 48.6%
associate-/l*49.0%
Simplified49.0%
div-inv49.0%
unpow-prod-down49.0%
pow249.0%
pow1/249.0%
pow1/249.0%
pow-prod-up49.0%
metadata-eval49.0%
metadata-eval49.0%
Applied egg-rr49.0%
if 3.70000000000000009e-25 < k < 1.8000000000000001e147Initial program 22.6%
Simplified22.6%
Taylor expanded in t around 0 83.6%
if 1.8000000000000001e147 < k Initial program 38.1%
*-commutative38.1%
associate-/r*38.1%
Simplified48.0%
add-sqr-sqrt48.0%
add-cube-cbrt48.0%
times-frac48.0%
Applied egg-rr90.1%
associate-/r/90.0%
associate-/l/90.0%
Simplified90.0%
associate-*l/90.1%
Applied egg-rr90.1%
Final simplification57.3%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (* (/ t_m (pow (cbrt l) 2.0)) (cbrt (* (sin k_m) (tan k_m))))))
(*
t_s
(if (<= k_m 1.58e-24)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0))
(if (<= k_m 1.35e+147)
(/
2.0
(/
(* (pow k_m 2.0) (* t_m (pow (sin k_m) 2.0)))
(* (cos k_m) (pow l 2.0))))
(*
(/ (sqrt 2.0) (* t_2 (/ k_m t_m)))
(/ (* t_m (/ (sqrt 2.0) k_m)) (pow t_2 2.0))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = (t_m / pow(cbrt(l), 2.0)) * cbrt((sin(k_m) * tan(k_m)));
double tmp;
if (k_m <= 1.58e-24) {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
} else if (k_m <= 1.35e+147) {
tmp = 2.0 / ((pow(k_m, 2.0) * (t_m * pow(sin(k_m), 2.0))) / (cos(k_m) * pow(l, 2.0)));
} else {
tmp = (sqrt(2.0) / (t_2 * (k_m / t_m))) * ((t_m * (sqrt(2.0) / k_m)) / pow(t_2, 2.0));
}
return t_s * tmp;
}
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = (t_m / Math.pow(Math.cbrt(l), 2.0)) * Math.cbrt((Math.sin(k_m) * Math.tan(k_m)));
double tmp;
if (k_m <= 1.58e-24) {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else if (k_m <= 1.35e+147) {
tmp = 2.0 / ((Math.pow(k_m, 2.0) * (t_m * Math.pow(Math.sin(k_m), 2.0))) / (Math.cos(k_m) * Math.pow(l, 2.0)));
} else {
tmp = (Math.sqrt(2.0) / (t_2 * (k_m / t_m))) * ((t_m * (Math.sqrt(2.0) / k_m)) / Math.pow(t_2, 2.0));
}
return t_s * tmp;
}
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = Float64(Float64(t_m / (cbrt(l) ^ 2.0)) * cbrt(Float64(sin(k_m) * tan(k_m)))) tmp = 0.0 if (k_m <= 1.58e-24) tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); elseif (k_m <= 1.35e+147) tmp = Float64(2.0 / Float64(Float64((k_m ^ 2.0) * Float64(t_m * (sin(k_m) ^ 2.0))) / Float64(cos(k_m) * (l ^ 2.0)))); else tmp = Float64(Float64(sqrt(2.0) / Float64(t_2 * Float64(k_m / t_m))) * Float64(Float64(t_m * Float64(sqrt(2.0) / k_m)) / (t_2 ^ 2.0))); end return Float64(t_s * tmp) end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[(N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 1.58e-24], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 1.35e+147], N[(2.0 / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[Sqrt[2.0], $MachinePrecision] / N[(t$95$2 * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(t$95$m * N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision] / N[Power[t$95$2, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \frac{t\_m}{{\left(\sqrt[3]{\ell}\right)}^{2}} \cdot \sqrt[3]{\sin k\_m \cdot \tan k\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 1.58 \cdot 10^{-24}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 1.35 \cdot 10^{+147}:\\
\;\;\;\;\frac{2}{\frac{{k\_m}^{2} \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}{\cos k\_m \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{2}}{t\_2 \cdot \frac{k\_m}{t\_m}} \cdot \frac{t\_m \cdot \frac{\sqrt{2}}{k\_m}}{{t\_2}^{2}}\\
\end{array}
\end{array}
\end{array}
if k < 1.5799999999999999e-24Initial program 29.0%
*-commutative29.0%
associate-/r*29.0%
Simplified34.2%
Applied egg-rr28.3%
Taylor expanded in t around 0 48.6%
associate-/l*49.0%
Simplified49.0%
div-inv49.0%
unpow-prod-down49.0%
pow249.0%
pow1/249.0%
pow1/249.0%
pow-prod-up49.0%
metadata-eval49.0%
metadata-eval49.0%
Applied egg-rr49.0%
if 1.5799999999999999e-24 < k < 1.34999999999999999e147Initial program 22.6%
Simplified22.6%
Taylor expanded in t around 0 83.6%
if 1.34999999999999999e147 < k Initial program 38.1%
*-commutative38.1%
associate-/r*38.1%
Simplified48.0%
add-sqr-sqrt48.0%
add-cube-cbrt48.0%
times-frac48.0%
Applied egg-rr90.1%
associate-/r/90.0%
associate-/l/90.0%
Simplified90.0%
Final simplification57.3%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (cbrt (* (sin k_m) (tan k_m)))) (t_3 (pow (cbrt l) -2.0)))
(*
t_s
(if (<= k_m 3.7e-25)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0))
(if (<= k_m 4.9e+148)
(/
2.0
(/
(* (pow k_m 2.0) (* t_m (pow (sin k_m) 2.0)))
(* (cos k_m) (pow l 2.0))))
(*
(sqrt 2.0)
(/
(* t_m (/ (sqrt 2.0) k_m))
(*
(pow (* t_m (* t_2 t_3)) 2.0)
(/ (* (* k_m t_2) (* t_m t_3)) t_m)))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = cbrt((sin(k_m) * tan(k_m)));
double t_3 = pow(cbrt(l), -2.0);
double tmp;
if (k_m <= 3.7e-25) {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
} else if (k_m <= 4.9e+148) {
tmp = 2.0 / ((pow(k_m, 2.0) * (t_m * pow(sin(k_m), 2.0))) / (cos(k_m) * pow(l, 2.0)));
} else {
tmp = sqrt(2.0) * ((t_m * (sqrt(2.0) / k_m)) / (pow((t_m * (t_2 * t_3)), 2.0) * (((k_m * t_2) * (t_m * t_3)) / t_m)));
}
return t_s * tmp;
}
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = Math.cbrt((Math.sin(k_m) * Math.tan(k_m)));
double t_3 = Math.pow(Math.cbrt(l), -2.0);
double tmp;
if (k_m <= 3.7e-25) {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else if (k_m <= 4.9e+148) {
tmp = 2.0 / ((Math.pow(k_m, 2.0) * (t_m * Math.pow(Math.sin(k_m), 2.0))) / (Math.cos(k_m) * Math.pow(l, 2.0)));
} else {
tmp = Math.sqrt(2.0) * ((t_m * (Math.sqrt(2.0) / k_m)) / (Math.pow((t_m * (t_2 * t_3)), 2.0) * (((k_m * t_2) * (t_m * t_3)) / t_m)));
}
return t_s * tmp;
}
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = cbrt(Float64(sin(k_m) * tan(k_m))) t_3 = cbrt(l) ^ -2.0 tmp = 0.0 if (k_m <= 3.7e-25) tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); elseif (k_m <= 4.9e+148) tmp = Float64(2.0 / Float64(Float64((k_m ^ 2.0) * Float64(t_m * (sin(k_m) ^ 2.0))) / Float64(cos(k_m) * (l ^ 2.0)))); else tmp = Float64(sqrt(2.0) * Float64(Float64(t_m * Float64(sqrt(2.0) / k_m)) / Float64((Float64(t_m * Float64(t_2 * t_3)) ^ 2.0) * Float64(Float64(Float64(k_m * t_2) * Float64(t_m * t_3)) / t_m)))); end return Float64(t_s * tmp) end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]}, Block[{t$95$3 = N[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 3.7e-25], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 4.9e+148], N[(2.0 / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Sqrt[2.0], $MachinePrecision] * N[(N[(t$95$m * N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision] / N[(N[Power[N[(t$95$m * N[(t$95$2 * t$95$3), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] * N[(N[(N[(k$95$m * t$95$2), $MachinePrecision] * N[(t$95$m * t$95$3), $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sqrt[3]{\sin k\_m \cdot \tan k\_m}\\
t_3 := {\left(\sqrt[3]{\ell}\right)}^{-2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 3.7 \cdot 10^{-25}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 4.9 \cdot 10^{+148}:\\
\;\;\;\;\frac{2}{\frac{{k\_m}^{2} \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}{\cos k\_m \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{2} \cdot \frac{t\_m \cdot \frac{\sqrt{2}}{k\_m}}{{\left(t\_m \cdot \left(t\_2 \cdot t\_3\right)\right)}^{2} \cdot \frac{\left(k\_m \cdot t\_2\right) \cdot \left(t\_m \cdot t\_3\right)}{t\_m}}\\
\end{array}
\end{array}
\end{array}
if k < 3.70000000000000009e-25Initial program 29.0%
*-commutative29.0%
associate-/r*29.0%
Simplified34.2%
Applied egg-rr28.3%
Taylor expanded in t around 0 48.6%
associate-/l*49.0%
Simplified49.0%
div-inv49.0%
unpow-prod-down49.0%
pow249.0%
pow1/249.0%
pow1/249.0%
pow-prod-up49.0%
metadata-eval49.0%
metadata-eval49.0%
Applied egg-rr49.0%
if 3.70000000000000009e-25 < k < 4.9e148Initial program 21.4%
Simplified21.4%
Taylor expanded in t around 0 82.0%
if 4.9e148 < k Initial program 42.1%
*-commutative42.1%
associate-/r*42.1%
Simplified47.4%
add-sqr-sqrt47.4%
add-cube-cbrt47.4%
times-frac47.4%
Applied egg-rr89.2%
associate-/r/89.1%
associate-/l/89.1%
Simplified89.1%
associate-*l/89.2%
Applied egg-rr89.2%
frac-times89.0%
associate-/l*89.0%
div-inv89.0%
pow-flip89.0%
metadata-eval89.0%
associate-*l*88.9%
Applied egg-rr88.9%
*-commutative88.9%
associate-/l*89.1%
associate-*r/89.0%
*-commutative89.0%
associate-/l*89.0%
associate-*l*89.2%
Simplified89.3%
Final simplification56.9%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 2e-24)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0))
(if (<= k_m 1.52e+147)
(/
2.0
(/
(* (pow k_m 2.0) (* t_m (pow (sin k_m) 2.0)))
(* (cos k_m) (pow l 2.0))))
(/
2.0
(pow
(*
(/ t_m (pow (cbrt l) 2.0))
(cbrt (* (sin k_m) (* (tan k_m) (pow (/ k_m t_m) 2.0)))))
3.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2e-24) {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
} else if (k_m <= 1.52e+147) {
tmp = 2.0 / ((pow(k_m, 2.0) * (t_m * pow(sin(k_m), 2.0))) / (cos(k_m) * pow(l, 2.0)));
} else {
tmp = 2.0 / pow(((t_m / pow(cbrt(l), 2.0)) * cbrt((sin(k_m) * (tan(k_m) * pow((k_m / t_m), 2.0))))), 3.0);
}
return t_s * tmp;
}
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2e-24) {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else if (k_m <= 1.52e+147) {
tmp = 2.0 / ((Math.pow(k_m, 2.0) * (t_m * Math.pow(Math.sin(k_m), 2.0))) / (Math.cos(k_m) * Math.pow(l, 2.0)));
} else {
tmp = 2.0 / Math.pow(((t_m / Math.pow(Math.cbrt(l), 2.0)) * Math.cbrt((Math.sin(k_m) * (Math.tan(k_m) * Math.pow((k_m / t_m), 2.0))))), 3.0);
}
return t_s * tmp;
}
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 2e-24) tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); elseif (k_m <= 1.52e+147) tmp = Float64(2.0 / Float64(Float64((k_m ^ 2.0) * Float64(t_m * (sin(k_m) ^ 2.0))) / Float64(cos(k_m) * (l ^ 2.0)))); else tmp = Float64(2.0 / (Float64(Float64(t_m / (cbrt(l) ^ 2.0)) * cbrt(Float64(sin(k_m) * Float64(tan(k_m) * (Float64(k_m / t_m) ^ 2.0))))) ^ 3.0)); end return Float64(t_s * tmp) end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 2e-24], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 1.52e+147], N[(2.0 / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Tan[k$95$m], $MachinePrecision] * N[Power[N[(k$95$m / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 2 \cdot 10^{-24}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 1.52 \cdot 10^{+147}:\\
\;\;\;\;\frac{2}{\frac{{k\_m}^{2} \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}{\cos k\_m \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{t\_m}{{\left(\sqrt[3]{\ell}\right)}^{2}} \cdot \sqrt[3]{\sin k\_m \cdot \left(\tan k\_m \cdot {\left(\frac{k\_m}{t\_m}\right)}^{2}\right)}\right)}^{3}}\\
\end{array}
\end{array}
if k < 1.99999999999999985e-24Initial program 29.0%
*-commutative29.0%
associate-/r*29.0%
Simplified34.2%
Applied egg-rr28.3%
Taylor expanded in t around 0 48.6%
associate-/l*49.0%
Simplified49.0%
div-inv49.0%
unpow-prod-down49.0%
pow249.0%
pow1/249.0%
pow1/249.0%
pow-prod-up49.0%
metadata-eval49.0%
metadata-eval49.0%
Applied egg-rr49.0%
if 1.99999999999999985e-24 < k < 1.51999999999999992e147Initial program 22.6%
Simplified22.6%
Taylor expanded in t around 0 83.6%
if 1.51999999999999992e147 < k Initial program 38.1%
Simplified38.1%
add-cube-cbrt38.1%
pow338.1%
Applied egg-rr81.3%
Final simplification56.5%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 2e-24)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0))
(/
2.0
(/
(* (pow k_m 2.0) (* t_m (pow (sin k_m) 2.0)))
(* (cos k_m) (pow l 2.0)))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2e-24) {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
} else {
tmp = 2.0 / ((pow(k_m, 2.0) * (t_m * pow(sin(k_m), 2.0))) / (cos(k_m) * pow(l, 2.0)));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 2d-24) then
tmp = 2.0d0 * ((1.0d0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ** 2.0d0)
else
tmp = 2.0d0 / (((k_m ** 2.0d0) * (t_m * (sin(k_m) ** 2.0d0))) / (cos(k_m) * (l ** 2.0d0)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2e-24) {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else {
tmp = 2.0 / ((Math.pow(k_m, 2.0) * (t_m * Math.pow(Math.sin(k_m), 2.0))) / (Math.cos(k_m) * Math.pow(l, 2.0)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 2e-24: tmp = 2.0 * math.pow((1.0 / ((k_m * (math.sin(k_m) / l)) * math.sqrt((t_m / math.cos(k_m))))), 2.0) else: tmp = 2.0 / ((math.pow(k_m, 2.0) * (t_m * math.pow(math.sin(k_m), 2.0))) / (math.cos(k_m) * math.pow(l, 2.0))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 2e-24) tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); else tmp = Float64(2.0 / Float64(Float64((k_m ^ 2.0) * Float64(t_m * (sin(k_m) ^ 2.0))) / Float64(cos(k_m) * (l ^ 2.0)))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 2e-24) tmp = 2.0 * ((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ^ 2.0); else tmp = 2.0 / (((k_m ^ 2.0) * (t_m * (sin(k_m) ^ 2.0))) / (cos(k_m) * (l ^ 2.0))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 2e-24], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 2 \cdot 10^{-24}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{{k\_m}^{2} \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}{\cos k\_m \cdot {\ell}^{2}}}\\
\end{array}
\end{array}
if k < 1.99999999999999985e-24Initial program 29.0%
*-commutative29.0%
associate-/r*29.0%
Simplified34.2%
Applied egg-rr28.3%
Taylor expanded in t around 0 48.6%
associate-/l*49.0%
Simplified49.0%
div-inv49.0%
unpow-prod-down49.0%
pow249.0%
pow1/249.0%
pow1/249.0%
pow-prod-up49.0%
metadata-eval49.0%
metadata-eval49.0%
Applied egg-rr49.0%
if 1.99999999999999985e-24 < k Initial program 28.3%
Simplified28.3%
Taylor expanded in t around 0 77.7%
Final simplification55.4%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 8.2e-25)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0))
(*
(/ 2.0 (* t_m (pow k_m 2.0)))
(/ (* (cos k_m) (pow l 2.0)) (pow (sin k_m) 2.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 8.2e-25) {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
} else {
tmp = (2.0 / (t_m * pow(k_m, 2.0))) * ((cos(k_m) * pow(l, 2.0)) / pow(sin(k_m), 2.0));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 8.2d-25) then
tmp = 2.0d0 * ((1.0d0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ** 2.0d0)
else
tmp = (2.0d0 / (t_m * (k_m ** 2.0d0))) * ((cos(k_m) * (l ** 2.0d0)) / (sin(k_m) ** 2.0d0))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 8.2e-25) {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else {
tmp = (2.0 / (t_m * Math.pow(k_m, 2.0))) * ((Math.cos(k_m) * Math.pow(l, 2.0)) / Math.pow(Math.sin(k_m), 2.0));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 8.2e-25: tmp = 2.0 * math.pow((1.0 / ((k_m * (math.sin(k_m) / l)) * math.sqrt((t_m / math.cos(k_m))))), 2.0) else: tmp = (2.0 / (t_m * math.pow(k_m, 2.0))) * ((math.cos(k_m) * math.pow(l, 2.0)) / math.pow(math.sin(k_m), 2.0)) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 8.2e-25) tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); else tmp = Float64(Float64(2.0 / Float64(t_m * (k_m ^ 2.0))) * Float64(Float64(cos(k_m) * (l ^ 2.0)) / (sin(k_m) ^ 2.0))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 8.2e-25) tmp = 2.0 * ((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ^ 2.0); else tmp = (2.0 / (t_m * (k_m ^ 2.0))) * ((cos(k_m) * (l ^ 2.0)) / (sin(k_m) ^ 2.0)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 8.2e-25], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 8.2 \cdot 10^{-25}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t\_m \cdot {k\_m}^{2}} \cdot \frac{\cos k\_m \cdot {\ell}^{2}}{{\sin k\_m}^{2}}\\
\end{array}
\end{array}
if k < 8.19999999999999974e-25Initial program 29.0%
*-commutative29.0%
associate-/r*29.0%
Simplified34.2%
Applied egg-rr28.3%
Taylor expanded in t around 0 48.6%
associate-/l*49.0%
Simplified49.0%
div-inv49.0%
unpow-prod-down49.0%
pow249.0%
pow1/249.0%
pow1/249.0%
pow-prod-up49.0%
metadata-eval49.0%
metadata-eval49.0%
Applied egg-rr49.0%
if 8.19999999999999974e-25 < k Initial program 28.3%
Simplified40.7%
Taylor expanded in t around 0 77.7%
associate-*r/77.7%
associate-*r*77.7%
times-frac77.7%
*-commutative77.7%
Simplified77.7%
Final simplification55.4%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 1e-204)
(pow (/ (sqrt 2.0) (* (/ (pow k_m 2.0) l) (sqrt t_m))) 2.0)
(*
2.0
(pow (* l (/ 1.0 (* (sqrt (/ t_m (cos k_m))) (* k_m (sin k_m))))) 2.0)))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 1e-204) {
tmp = pow((sqrt(2.0) / ((pow(k_m, 2.0) / l) * sqrt(t_m))), 2.0);
} else {
tmp = 2.0 * pow((l * (1.0 / (sqrt((t_m / cos(k_m))) * (k_m * sin(k_m))))), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 1d-204) then
tmp = (sqrt(2.0d0) / (((k_m ** 2.0d0) / l) * sqrt(t_m))) ** 2.0d0
else
tmp = 2.0d0 * ((l * (1.0d0 / (sqrt((t_m / cos(k_m))) * (k_m * sin(k_m))))) ** 2.0d0)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 1e-204) {
tmp = Math.pow((Math.sqrt(2.0) / ((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m))), 2.0);
} else {
tmp = 2.0 * Math.pow((l * (1.0 / (Math.sqrt((t_m / Math.cos(k_m))) * (k_m * Math.sin(k_m))))), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 1e-204: tmp = math.pow((math.sqrt(2.0) / ((math.pow(k_m, 2.0) / l) * math.sqrt(t_m))), 2.0) else: tmp = 2.0 * math.pow((l * (1.0 / (math.sqrt((t_m / math.cos(k_m))) * (k_m * math.sin(k_m))))), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 1e-204) tmp = Float64(sqrt(2.0) / Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m))) ^ 2.0; else tmp = Float64(2.0 * (Float64(l * Float64(1.0 / Float64(sqrt(Float64(t_m / cos(k_m))) * Float64(k_m * sin(k_m))))) ^ 2.0)); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 1e-204) tmp = (sqrt(2.0) / (((k_m ^ 2.0) / l) * sqrt(t_m))) ^ 2.0; else tmp = 2.0 * ((l * (1.0 / (sqrt((t_m / cos(k_m))) * (k_m * sin(k_m))))) ^ 2.0); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 1e-204], N[Power[N[(N[Sqrt[2.0], $MachinePrecision] / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(2.0 * N[Power[N[(l * N[(1.0 / N[(N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 10^{-204}:\\
\;\;\;\;{\left(\frac{\sqrt{2}}{\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot {\left(\ell \cdot \frac{1}{\sqrt{\frac{t\_m}{\cos k\_m}} \cdot \left(k\_m \cdot \sin k\_m\right)}\right)}^{2}\\
\end{array}
\end{array}
if (*.f64 l l) < 1e-204Initial program 22.4%
*-commutative22.4%
associate-/r*22.4%
Simplified29.4%
Applied egg-rr23.3%
Taylor expanded in k around 0 42.5%
if 1e-204 < (*.f64 l l) Initial program 32.4%
*-commutative32.4%
associate-/r*32.4%
Simplified38.4%
Applied egg-rr29.1%
Taylor expanded in t around 0 52.1%
associate-/l*52.1%
Simplified52.1%
div-inv52.1%
unpow-prod-down52.0%
pow252.0%
pow1/252.0%
pow1/252.0%
pow-prod-up52.1%
metadata-eval52.1%
metadata-eval52.1%
Applied egg-rr52.1%
associate-*r/52.1%
associate-*l/51.6%
associate-/r/51.6%
Simplified51.6%
Final simplification48.4%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 3.7e-25)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0))
(*
(/ (* 2.0 (cos k_m)) (* (pow (sin k_m) 2.0) (* t_m (pow k_m 2.0))))
(* l l)))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 3.7e-25) {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
} else {
tmp = ((2.0 * cos(k_m)) / (pow(sin(k_m), 2.0) * (t_m * pow(k_m, 2.0)))) * (l * l);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 3.7d-25) then
tmp = 2.0d0 * ((1.0d0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ** 2.0d0)
else
tmp = ((2.0d0 * cos(k_m)) / ((sin(k_m) ** 2.0d0) * (t_m * (k_m ** 2.0d0)))) * (l * l)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 3.7e-25) {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else {
tmp = ((2.0 * Math.cos(k_m)) / (Math.pow(Math.sin(k_m), 2.0) * (t_m * Math.pow(k_m, 2.0)))) * (l * l);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 3.7e-25: tmp = 2.0 * math.pow((1.0 / ((k_m * (math.sin(k_m) / l)) * math.sqrt((t_m / math.cos(k_m))))), 2.0) else: tmp = ((2.0 * math.cos(k_m)) / (math.pow(math.sin(k_m), 2.0) * (t_m * math.pow(k_m, 2.0)))) * (l * l) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 3.7e-25) tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); else tmp = Float64(Float64(Float64(2.0 * cos(k_m)) / Float64((sin(k_m) ^ 2.0) * Float64(t_m * (k_m ^ 2.0)))) * Float64(l * l)); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 3.7e-25) tmp = 2.0 * ((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ^ 2.0); else tmp = ((2.0 * cos(k_m)) / ((sin(k_m) ^ 2.0) * (t_m * (k_m ^ 2.0)))) * (l * l); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 3.7e-25], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] * N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l * l), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 3.7 \cdot 10^{-25}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \cos k\_m}{{\sin k\_m}^{2} \cdot \left(t\_m \cdot {k\_m}^{2}\right)} \cdot \left(\ell \cdot \ell\right)\\
\end{array}
\end{array}
if k < 3.70000000000000009e-25Initial program 29.0%
*-commutative29.0%
associate-/r*29.0%
Simplified34.2%
Applied egg-rr28.3%
Taylor expanded in t around 0 48.6%
associate-/l*49.0%
Simplified49.0%
div-inv49.0%
unpow-prod-down49.0%
pow249.0%
pow1/249.0%
pow1/249.0%
pow-prod-up49.0%
metadata-eval49.0%
metadata-eval49.0%
Applied egg-rr49.0%
if 3.70000000000000009e-25 < k Initial program 28.3%
Simplified40.7%
Taylor expanded in t around 0 77.6%
associate-*r/77.6%
associate-*r*77.7%
Simplified77.7%
Final simplification55.4%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 1e-25)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0))
(*
(* l l)
(/ 2.0 (/ (* (pow k_m 2.0) (* t_m (pow (sin k_m) 2.0))) (cos k_m)))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1e-25) {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
} else {
tmp = (l * l) * (2.0 / ((pow(k_m, 2.0) * (t_m * pow(sin(k_m), 2.0))) / cos(k_m)));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 1d-25) then
tmp = 2.0d0 * ((1.0d0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ** 2.0d0)
else
tmp = (l * l) * (2.0d0 / (((k_m ** 2.0d0) * (t_m * (sin(k_m) ** 2.0d0))) / cos(k_m)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1e-25) {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else {
tmp = (l * l) * (2.0 / ((Math.pow(k_m, 2.0) * (t_m * Math.pow(Math.sin(k_m), 2.0))) / Math.cos(k_m)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 1e-25: tmp = 2.0 * math.pow((1.0 / ((k_m * (math.sin(k_m) / l)) * math.sqrt((t_m / math.cos(k_m))))), 2.0) else: tmp = (l * l) * (2.0 / ((math.pow(k_m, 2.0) * (t_m * math.pow(math.sin(k_m), 2.0))) / math.cos(k_m))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 1e-25) tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); else tmp = Float64(Float64(l * l) * Float64(2.0 / Float64(Float64((k_m ^ 2.0) * Float64(t_m * (sin(k_m) ^ 2.0))) / cos(k_m)))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 1e-25) tmp = 2.0 * ((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ^ 2.0); else tmp = (l * l) * (2.0 / (((k_m ^ 2.0) * (t_m * (sin(k_m) ^ 2.0))) / cos(k_m))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 1e-25], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 10^{-25}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{\frac{{k\_m}^{2} \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}{\cos k\_m}}\\
\end{array}
\end{array}
if k < 1.00000000000000004e-25Initial program 29.0%
*-commutative29.0%
associate-/r*29.0%
Simplified34.2%
Applied egg-rr28.3%
Taylor expanded in t around 0 48.6%
associate-/l*49.0%
Simplified49.0%
div-inv49.0%
unpow-prod-down49.0%
pow249.0%
pow1/249.0%
pow1/249.0%
pow-prod-up49.0%
metadata-eval49.0%
metadata-eval49.0%
Applied egg-rr49.0%
if 1.00000000000000004e-25 < k Initial program 28.3%
Simplified40.7%
Taylor expanded in t around 0 77.6%
Final simplification55.4%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 2e-24)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0))
(*
(* l l)
(* 2.0 (/ (/ (cos k_m) (pow k_m 2.0)) (* t_m (pow (sin k_m) 2.0))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2e-24) {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
} else {
tmp = (l * l) * (2.0 * ((cos(k_m) / pow(k_m, 2.0)) / (t_m * pow(sin(k_m), 2.0))));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 2d-24) then
tmp = 2.0d0 * ((1.0d0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ** 2.0d0)
else
tmp = (l * l) * (2.0d0 * ((cos(k_m) / (k_m ** 2.0d0)) / (t_m * (sin(k_m) ** 2.0d0))))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 2e-24) {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else {
tmp = (l * l) * (2.0 * ((Math.cos(k_m) / Math.pow(k_m, 2.0)) / (t_m * Math.pow(Math.sin(k_m), 2.0))));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 2e-24: tmp = 2.0 * math.pow((1.0 / ((k_m * (math.sin(k_m) / l)) * math.sqrt((t_m / math.cos(k_m))))), 2.0) else: tmp = (l * l) * (2.0 * ((math.cos(k_m) / math.pow(k_m, 2.0)) / (t_m * math.pow(math.sin(k_m), 2.0)))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 2e-24) tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); else tmp = Float64(Float64(l * l) * Float64(2.0 * Float64(Float64(cos(k_m) / (k_m ^ 2.0)) / Float64(t_m * (sin(k_m) ^ 2.0))))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 2e-24) tmp = 2.0 * ((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ^ 2.0); else tmp = (l * l) * (2.0 * ((cos(k_m) / (k_m ^ 2.0)) / (t_m * (sin(k_m) ^ 2.0)))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 2e-24], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 2 \cdot 10^{-24}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \left(2 \cdot \frac{\frac{\cos k\_m}{{k\_m}^{2}}}{t\_m \cdot {\sin k\_m}^{2}}\right)\\
\end{array}
\end{array}
if k < 1.99999999999999985e-24Initial program 29.0%
*-commutative29.0%
associate-/r*29.0%
Simplified34.2%
Applied egg-rr28.3%
Taylor expanded in t around 0 48.6%
associate-/l*49.0%
Simplified49.0%
div-inv49.0%
unpow-prod-down49.0%
pow249.0%
pow1/249.0%
pow1/249.0%
pow-prod-up49.0%
metadata-eval49.0%
metadata-eval49.0%
Applied egg-rr49.0%
if 1.99999999999999985e-24 < k Initial program 28.3%
Simplified40.7%
Taylor expanded in t around 0 77.6%
associate-/r*77.6%
Simplified77.6%
Final simplification55.4%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 5e-11)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0))
(*
2.0
(*
(/ (cos k_m) (pow k_m 2.0))
(/ (* l l) (* t_m (pow (sin k_m) 2.0))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 5e-11) {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
} else {
tmp = 2.0 * ((cos(k_m) / pow(k_m, 2.0)) * ((l * l) / (t_m * pow(sin(k_m), 2.0))));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 5d-11) then
tmp = 2.0d0 * ((1.0d0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ** 2.0d0)
else
tmp = 2.0d0 * ((cos(k_m) / (k_m ** 2.0d0)) * ((l * l) / (t_m * (sin(k_m) ** 2.0d0))))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 5e-11) {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
} else {
tmp = 2.0 * ((Math.cos(k_m) / Math.pow(k_m, 2.0)) * ((l * l) / (t_m * Math.pow(Math.sin(k_m), 2.0))));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 5e-11: tmp = 2.0 * math.pow((1.0 / ((k_m * (math.sin(k_m) / l)) * math.sqrt((t_m / math.cos(k_m))))), 2.0) else: tmp = 2.0 * ((math.cos(k_m) / math.pow(k_m, 2.0)) * ((l * l) / (t_m * math.pow(math.sin(k_m), 2.0)))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 5e-11) tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); else tmp = Float64(2.0 * Float64(Float64(cos(k_m) / (k_m ^ 2.0)) * Float64(Float64(l * l) / Float64(t_m * (sin(k_m) ^ 2.0))))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 5e-11) tmp = 2.0 * ((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ^ 2.0); else tmp = 2.0 * ((cos(k_m) / (k_m ^ 2.0)) * ((l * l) / (t_m * (sin(k_m) ^ 2.0)))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 5e-11], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] / N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 5 \cdot 10^{-11}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{\cos k\_m}{{k\_m}^{2}} \cdot \frac{\ell \cdot \ell}{t\_m \cdot {\sin k\_m}^{2}}\right)\\
\end{array}
\end{array}
if k < 5.00000000000000018e-11Initial program 29.0%
*-commutative29.0%
associate-/r*29.0%
Simplified34.5%
Applied egg-rr28.7%
Taylor expanded in t around 0 48.7%
associate-/l*49.0%
Simplified49.0%
div-inv49.0%
unpow-prod-down49.0%
pow249.0%
pow1/249.0%
pow1/249.0%
pow-prod-up49.1%
metadata-eval49.1%
metadata-eval49.1%
Applied egg-rr49.1%
if 5.00000000000000018e-11 < k Initial program 28.5%
*-commutative28.5%
associate-/r*28.5%
Simplified38.1%
add-sqr-sqrt38.1%
add-cube-cbrt38.1%
times-frac38.1%
Applied egg-rr82.1%
associate-/r/82.1%
associate-/l/82.1%
Simplified82.1%
Taylor expanded in k around inf 75.6%
associate-*r*75.7%
unpow275.7%
rem-square-sqrt76.0%
associate-*l/76.0%
*-commutative76.0%
*-commutative76.0%
times-frac72.5%
*-commutative72.5%
Simplified72.5%
pow272.5%
Applied egg-rr72.5%
Final simplification53.9%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= l 1.9e-102)
(pow (/ (sqrt 2.0) (* (/ (pow k_m 2.0) l) (sqrt t_m))) 2.0)
(*
2.0
(pow (/ 1.0 (* (* k_m (/ (sin k_m) l)) (sqrt (/ t_m (cos k_m))))) 2.0)))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (l <= 1.9e-102) {
tmp = pow((sqrt(2.0) / ((pow(k_m, 2.0) / l) * sqrt(t_m))), 2.0);
} else {
tmp = 2.0 * pow((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (l <= 1.9d-102) then
tmp = (sqrt(2.0d0) / (((k_m ** 2.0d0) / l) * sqrt(t_m))) ** 2.0d0
else
tmp = 2.0d0 * ((1.0d0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ** 2.0d0)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (l <= 1.9e-102) {
tmp = Math.pow((Math.sqrt(2.0) / ((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m))), 2.0);
} else {
tmp = 2.0 * Math.pow((1.0 / ((k_m * (Math.sin(k_m) / l)) * Math.sqrt((t_m / Math.cos(k_m))))), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if l <= 1.9e-102: tmp = math.pow((math.sqrt(2.0) / ((math.pow(k_m, 2.0) / l) * math.sqrt(t_m))), 2.0) else: tmp = 2.0 * math.pow((1.0 / ((k_m * (math.sin(k_m) / l)) * math.sqrt((t_m / math.cos(k_m))))), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (l <= 1.9e-102) tmp = Float64(sqrt(2.0) / Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m))) ^ 2.0; else tmp = Float64(2.0 * (Float64(1.0 / Float64(Float64(k_m * Float64(sin(k_m) / l)) * sqrt(Float64(t_m / cos(k_m))))) ^ 2.0)); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (l <= 1.9e-102) tmp = (sqrt(2.0) / (((k_m ^ 2.0) / l) * sqrt(t_m))) ^ 2.0; else tmp = 2.0 * ((1.0 / ((k_m * (sin(k_m) / l)) * sqrt((t_m / cos(k_m))))) ^ 2.0); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[l, 1.9e-102], N[Power[N[(N[Sqrt[2.0], $MachinePrecision] / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(2.0 * N[Power[N[(1.0 / N[(N[(k$95$m * N[(N[Sin[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \leq 1.9 \cdot 10^{-102}:\\
\;\;\;\;{\left(\frac{\sqrt{2}}{\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot {\left(\frac{1}{\left(k\_m \cdot \frac{\sin k\_m}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k\_m}}}\right)}^{2}\\
\end{array}
\end{array}
if l < 1.90000000000000013e-102Initial program 28.3%
*-commutative28.3%
associate-/r*28.3%
Simplified36.1%
Applied egg-rr28.7%
Taylor expanded in k around 0 43.3%
if 1.90000000000000013e-102 < l Initial program 29.9%
*-commutative29.9%
associate-/r*29.9%
Simplified33.5%
Applied egg-rr24.0%
Taylor expanded in t around 0 45.2%
associate-/l*45.2%
Simplified45.2%
div-inv45.2%
unpow-prod-down45.1%
pow245.1%
pow1/245.1%
pow1/245.1%
pow-prod-up45.2%
metadata-eval45.2%
metadata-eval45.2%
Applied egg-rr45.2%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= l 6e+215)
(pow (/ (sqrt 2.0) (* (/ (pow k_m 2.0) l) (sqrt t_m))) 2.0)
(*
2.0
(* (/ (cos k_m) (pow k_m 2.0)) (/ (/ (pow l 2.0) t_m) (pow k_m 2.0)))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (l <= 6e+215) {
tmp = pow((sqrt(2.0) / ((pow(k_m, 2.0) / l) * sqrt(t_m))), 2.0);
} else {
tmp = 2.0 * ((cos(k_m) / pow(k_m, 2.0)) * ((pow(l, 2.0) / t_m) / pow(k_m, 2.0)));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (l <= 6d+215) then
tmp = (sqrt(2.0d0) / (((k_m ** 2.0d0) / l) * sqrt(t_m))) ** 2.0d0
else
tmp = 2.0d0 * ((cos(k_m) / (k_m ** 2.0d0)) * (((l ** 2.0d0) / t_m) / (k_m ** 2.0d0)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (l <= 6e+215) {
tmp = Math.pow((Math.sqrt(2.0) / ((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m))), 2.0);
} else {
tmp = 2.0 * ((Math.cos(k_m) / Math.pow(k_m, 2.0)) * ((Math.pow(l, 2.0) / t_m) / Math.pow(k_m, 2.0)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if l <= 6e+215: tmp = math.pow((math.sqrt(2.0) / ((math.pow(k_m, 2.0) / l) * math.sqrt(t_m))), 2.0) else: tmp = 2.0 * ((math.cos(k_m) / math.pow(k_m, 2.0)) * ((math.pow(l, 2.0) / t_m) / math.pow(k_m, 2.0))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (l <= 6e+215) tmp = Float64(sqrt(2.0) / Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m))) ^ 2.0; else tmp = Float64(2.0 * Float64(Float64(cos(k_m) / (k_m ^ 2.0)) * Float64(Float64((l ^ 2.0) / t_m) / (k_m ^ 2.0)))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (l <= 6e+215) tmp = (sqrt(2.0) / (((k_m ^ 2.0) / l) * sqrt(t_m))) ^ 2.0; else tmp = 2.0 * ((cos(k_m) / (k_m ^ 2.0)) * (((l ^ 2.0) / t_m) / (k_m ^ 2.0))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[l, 6e+215], N[Power[N[(N[Sqrt[2.0], $MachinePrecision] / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(N[Power[l, 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \leq 6 \cdot 10^{+215}:\\
\;\;\;\;{\left(\frac{\sqrt{2}}{\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{\cos k\_m}{{k\_m}^{2}} \cdot \frac{\frac{{\ell}^{2}}{t\_m}}{{k\_m}^{2}}\right)\\
\end{array}
\end{array}
if l < 5.9999999999999998e215Initial program 27.9%
*-commutative27.9%
associate-/r*27.9%
Simplified34.8%
Applied egg-rr26.4%
Taylor expanded in k around 0 37.3%
if 5.9999999999999998e215 < l Initial program 40.0%
*-commutative40.0%
associate-/r*40.0%
Simplified40.0%
add-sqr-sqrt40.0%
add-cube-cbrt40.0%
times-frac40.0%
Applied egg-rr80.1%
associate-/r/80.1%
associate-/l/80.1%
Simplified80.1%
Taylor expanded in k around inf 90.5%
associate-*r*90.5%
unpow290.5%
rem-square-sqrt90.5%
associate-*l/90.5%
*-commutative90.5%
*-commutative90.5%
times-frac90.5%
*-commutative90.5%
Simplified90.5%
Taylor expanded in k around 0 90.5%
*-commutative90.5%
associate-/r*90.5%
Simplified90.5%
Final simplification41.5%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 19500.0)
(pow (/ (sqrt 2.0) (* (/ (pow k_m 2.0) l) (sqrt t_m))) 2.0)
(if (<= k_m 8e+150)
(/
2.0
(*
(* (* (sin k_m) (tan k_m)) (* (/ (pow t_m 2.0) l) (/ t_m l)))
(/ (/ k_m t_m) (/ t_m k_m))))
(pow (* (pow k_m -2.0) (* l (sqrt (/ 2.0 t_m)))) 2.0)))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 19500.0) {
tmp = pow((sqrt(2.0) / ((pow(k_m, 2.0) / l) * sqrt(t_m))), 2.0);
} else if (k_m <= 8e+150) {
tmp = 2.0 / (((sin(k_m) * tan(k_m)) * ((pow(t_m, 2.0) / l) * (t_m / l))) * ((k_m / t_m) / (t_m / k_m)));
} else {
tmp = pow((pow(k_m, -2.0) * (l * sqrt((2.0 / t_m)))), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 19500.0d0) then
tmp = (sqrt(2.0d0) / (((k_m ** 2.0d0) / l) * sqrt(t_m))) ** 2.0d0
else if (k_m <= 8d+150) then
tmp = 2.0d0 / (((sin(k_m) * tan(k_m)) * (((t_m ** 2.0d0) / l) * (t_m / l))) * ((k_m / t_m) / (t_m / k_m)))
else
tmp = ((k_m ** (-2.0d0)) * (l * sqrt((2.0d0 / t_m)))) ** 2.0d0
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 19500.0) {
tmp = Math.pow((Math.sqrt(2.0) / ((Math.pow(k_m, 2.0) / l) * Math.sqrt(t_m))), 2.0);
} else if (k_m <= 8e+150) {
tmp = 2.0 / (((Math.sin(k_m) * Math.tan(k_m)) * ((Math.pow(t_m, 2.0) / l) * (t_m / l))) * ((k_m / t_m) / (t_m / k_m)));
} else {
tmp = Math.pow((Math.pow(k_m, -2.0) * (l * Math.sqrt((2.0 / t_m)))), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 19500.0: tmp = math.pow((math.sqrt(2.0) / ((math.pow(k_m, 2.0) / l) * math.sqrt(t_m))), 2.0) elif k_m <= 8e+150: tmp = 2.0 / (((math.sin(k_m) * math.tan(k_m)) * ((math.pow(t_m, 2.0) / l) * (t_m / l))) * ((k_m / t_m) / (t_m / k_m))) else: tmp = math.pow((math.pow(k_m, -2.0) * (l * math.sqrt((2.0 / t_m)))), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 19500.0) tmp = Float64(sqrt(2.0) / Float64(Float64((k_m ^ 2.0) / l) * sqrt(t_m))) ^ 2.0; elseif (k_m <= 8e+150) tmp = Float64(2.0 / Float64(Float64(Float64(sin(k_m) * tan(k_m)) * Float64(Float64((t_m ^ 2.0) / l) * Float64(t_m / l))) * Float64(Float64(k_m / t_m) / Float64(t_m / k_m)))); else tmp = Float64((k_m ^ -2.0) * Float64(l * sqrt(Float64(2.0 / t_m)))) ^ 2.0; end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 19500.0) tmp = (sqrt(2.0) / (((k_m ^ 2.0) / l) * sqrt(t_m))) ^ 2.0; elseif (k_m <= 8e+150) tmp = 2.0 / (((sin(k_m) * tan(k_m)) * (((t_m ^ 2.0) / l) * (t_m / l))) * ((k_m / t_m) / (t_m / k_m))); else tmp = ((k_m ^ -2.0) * (l * sqrt((2.0 / t_m)))) ^ 2.0; end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 19500.0], N[Power[N[(N[Sqrt[2.0], $MachinePrecision] / N[(N[(N[Power[k$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[t$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[k$95$m, 8e+150], N[(2.0 / N[(N[(N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[(N[Power[t$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(k$95$m / t$95$m), $MachinePrecision] / N[(t$95$m / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[Power[k$95$m, -2.0], $MachinePrecision] * N[(l * N[Sqrt[N[(2.0 / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 19500:\\
\;\;\;\;{\left(\frac{\sqrt{2}}{\frac{{k\_m}^{2}}{\ell} \cdot \sqrt{t\_m}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 8 \cdot 10^{+150}:\\
\;\;\;\;\frac{2}{\left(\left(\sin k\_m \cdot \tan k\_m\right) \cdot \left(\frac{{t\_m}^{2}}{\ell} \cdot \frac{t\_m}{\ell}\right)\right) \cdot \frac{\frac{k\_m}{t\_m}}{\frac{t\_m}{k\_m}}}\\
\mathbf{else}:\\
\;\;\;\;{\left({k\_m}^{-2} \cdot \left(\ell \cdot \sqrt{\frac{2}{t\_m}}\right)\right)}^{2}\\
\end{array}
\end{array}
if k < 19500Initial program 29.1%
*-commutative29.1%
associate-/r*29.1%
Simplified34.4%
Applied egg-rr29.3%
Taylor expanded in k around 0 39.5%
if 19500 < k < 7.99999999999999985e150Initial program 17.6%
Simplified17.7%
unpow317.7%
times-frac25.1%
pow225.1%
Applied egg-rr25.1%
+-commutative25.1%
associate-+l-48.9%
metadata-eval48.9%
--rgt-identity48.9%
unpow248.9%
clear-num48.8%
un-div-inv48.9%
Applied egg-rr48.9%
if 7.99999999999999985e150 < k Initial program 44.4%
*-commutative44.4%
associate-/r*44.4%
Simplified50.0%
Applied egg-rr22.4%
Taylor expanded in k around 0 50.8%
associate-*l/50.8%
associate-*l*50.8%
associate-/l*50.8%
Simplified50.8%
pow150.8%
div-inv50.8%
sqrt-unprod50.8%
div-inv50.8%
pow-flip50.8%
metadata-eval50.8%
Applied egg-rr50.8%
unpow150.8%
associate-*r*50.8%
*-commutative50.8%
Simplified50.8%
Final simplification41.4%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (sqrt (/ 2.0 t_m))))
(*
t_s
(if (<= k_m 25000.0)
(pow (* l (/ t_2 (pow k_m 2.0))) 2.0)
(if (<= k_m 4.5e+150)
(/
2.0
(*
(* (* (sin k_m) (tan k_m)) (* (/ (pow t_m 2.0) l) (/ t_m l)))
(/ (/ k_m t_m) (/ t_m k_m))))
(pow (* (pow k_m -2.0) (* l t_2)) 2.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = sqrt((2.0 / t_m));
double tmp;
if (k_m <= 25000.0) {
tmp = pow((l * (t_2 / pow(k_m, 2.0))), 2.0);
} else if (k_m <= 4.5e+150) {
tmp = 2.0 / (((sin(k_m) * tan(k_m)) * ((pow(t_m, 2.0) / l) * (t_m / l))) * ((k_m / t_m) / (t_m / k_m)));
} else {
tmp = pow((pow(k_m, -2.0) * (l * t_2)), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: t_2
real(8) :: tmp
t_2 = sqrt((2.0d0 / t_m))
if (k_m <= 25000.0d0) then
tmp = (l * (t_2 / (k_m ** 2.0d0))) ** 2.0d0
else if (k_m <= 4.5d+150) then
tmp = 2.0d0 / (((sin(k_m) * tan(k_m)) * (((t_m ** 2.0d0) / l) * (t_m / l))) * ((k_m / t_m) / (t_m / k_m)))
else
tmp = ((k_m ** (-2.0d0)) * (l * t_2)) ** 2.0d0
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = Math.sqrt((2.0 / t_m));
double tmp;
if (k_m <= 25000.0) {
tmp = Math.pow((l * (t_2 / Math.pow(k_m, 2.0))), 2.0);
} else if (k_m <= 4.5e+150) {
tmp = 2.0 / (((Math.sin(k_m) * Math.tan(k_m)) * ((Math.pow(t_m, 2.0) / l) * (t_m / l))) * ((k_m / t_m) / (t_m / k_m)));
} else {
tmp = Math.pow((Math.pow(k_m, -2.0) * (l * t_2)), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): t_2 = math.sqrt((2.0 / t_m)) tmp = 0 if k_m <= 25000.0: tmp = math.pow((l * (t_2 / math.pow(k_m, 2.0))), 2.0) elif k_m <= 4.5e+150: tmp = 2.0 / (((math.sin(k_m) * math.tan(k_m)) * ((math.pow(t_m, 2.0) / l) * (t_m / l))) * ((k_m / t_m) / (t_m / k_m))) else: tmp = math.pow((math.pow(k_m, -2.0) * (l * t_2)), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = sqrt(Float64(2.0 / t_m)) tmp = 0.0 if (k_m <= 25000.0) tmp = Float64(l * Float64(t_2 / (k_m ^ 2.0))) ^ 2.0; elseif (k_m <= 4.5e+150) tmp = Float64(2.0 / Float64(Float64(Float64(sin(k_m) * tan(k_m)) * Float64(Float64((t_m ^ 2.0) / l) * Float64(t_m / l))) * Float64(Float64(k_m / t_m) / Float64(t_m / k_m)))); else tmp = Float64((k_m ^ -2.0) * Float64(l * t_2)) ^ 2.0; end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) t_2 = sqrt((2.0 / t_m)); tmp = 0.0; if (k_m <= 25000.0) tmp = (l * (t_2 / (k_m ^ 2.0))) ^ 2.0; elseif (k_m <= 4.5e+150) tmp = 2.0 / (((sin(k_m) * tan(k_m)) * (((t_m ^ 2.0) / l) * (t_m / l))) * ((k_m / t_m) / (t_m / k_m))); else tmp = ((k_m ^ -2.0) * (l * t_2)) ^ 2.0; end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[Sqrt[N[(2.0 / t$95$m), $MachinePrecision]], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 25000.0], N[Power[N[(l * N[(t$95$2 / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[k$95$m, 4.5e+150], N[(2.0 / N[(N[(N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[(N[Power[t$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(k$95$m / t$95$m), $MachinePrecision] / N[(t$95$m / k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[Power[k$95$m, -2.0], $MachinePrecision] * N[(l * t$95$2), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sqrt{\frac{2}{t\_m}}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 25000:\\
\;\;\;\;{\left(\ell \cdot \frac{t\_2}{{k\_m}^{2}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 4.5 \cdot 10^{+150}:\\
\;\;\;\;\frac{2}{\left(\left(\sin k\_m \cdot \tan k\_m\right) \cdot \left(\frac{{t\_m}^{2}}{\ell} \cdot \frac{t\_m}{\ell}\right)\right) \cdot \frac{\frac{k\_m}{t\_m}}{\frac{t\_m}{k\_m}}}\\
\mathbf{else}:\\
\;\;\;\;{\left({k\_m}^{-2} \cdot \left(\ell \cdot t\_2\right)\right)}^{2}\\
\end{array}
\end{array}
\end{array}
if k < 25000Initial program 29.1%
*-commutative29.1%
associate-/r*29.1%
Simplified34.4%
Applied egg-rr29.3%
Taylor expanded in k around 0 39.5%
associate-*l/39.1%
associate-*l*39.1%
associate-/l*39.5%
Simplified39.5%
associate-*r/39.1%
sqrt-unprod39.1%
div-inv39.1%
Applied egg-rr39.1%
associate-/l*39.5%
Simplified39.5%
if 25000 < k < 4.5e150Initial program 17.6%
Simplified17.7%
unpow317.7%
times-frac25.1%
pow225.1%
Applied egg-rr25.1%
+-commutative25.1%
associate-+l-48.9%
metadata-eval48.9%
--rgt-identity48.9%
unpow248.9%
clear-num48.8%
un-div-inv48.9%
Applied egg-rr48.9%
if 4.5e150 < k Initial program 44.4%
*-commutative44.4%
associate-/r*44.4%
Simplified50.0%
Applied egg-rr22.4%
Taylor expanded in k around 0 50.8%
associate-*l/50.8%
associate-*l*50.8%
associate-/l*50.8%
Simplified50.8%
pow150.8%
div-inv50.8%
sqrt-unprod50.8%
div-inv50.8%
pow-flip50.8%
metadata-eval50.8%
Applied egg-rr50.8%
unpow150.8%
associate-*r*50.8%
*-commutative50.8%
Simplified50.8%
Final simplification41.4%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (sqrt (/ 2.0 t_m))))
(*
t_s
(if (<= k_m 25000.0)
(pow (* l (/ t_2 (pow k_m 2.0))) 2.0)
(if (<= k_m 5.2e+150)
(/
2.0
(*
(* (* (sin k_m) (tan k_m)) (* (/ (pow t_m 2.0) l) (/ t_m l)))
(* (/ k_m t_m) (/ k_m t_m))))
(pow (* (pow k_m -2.0) (* l t_2)) 2.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = sqrt((2.0 / t_m));
double tmp;
if (k_m <= 25000.0) {
tmp = pow((l * (t_2 / pow(k_m, 2.0))), 2.0);
} else if (k_m <= 5.2e+150) {
tmp = 2.0 / (((sin(k_m) * tan(k_m)) * ((pow(t_m, 2.0) / l) * (t_m / l))) * ((k_m / t_m) * (k_m / t_m)));
} else {
tmp = pow((pow(k_m, -2.0) * (l * t_2)), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: t_2
real(8) :: tmp
t_2 = sqrt((2.0d0 / t_m))
if (k_m <= 25000.0d0) then
tmp = (l * (t_2 / (k_m ** 2.0d0))) ** 2.0d0
else if (k_m <= 5.2d+150) then
tmp = 2.0d0 / (((sin(k_m) * tan(k_m)) * (((t_m ** 2.0d0) / l) * (t_m / l))) * ((k_m / t_m) * (k_m / t_m)))
else
tmp = ((k_m ** (-2.0d0)) * (l * t_2)) ** 2.0d0
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = Math.sqrt((2.0 / t_m));
double tmp;
if (k_m <= 25000.0) {
tmp = Math.pow((l * (t_2 / Math.pow(k_m, 2.0))), 2.0);
} else if (k_m <= 5.2e+150) {
tmp = 2.0 / (((Math.sin(k_m) * Math.tan(k_m)) * ((Math.pow(t_m, 2.0) / l) * (t_m / l))) * ((k_m / t_m) * (k_m / t_m)));
} else {
tmp = Math.pow((Math.pow(k_m, -2.0) * (l * t_2)), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): t_2 = math.sqrt((2.0 / t_m)) tmp = 0 if k_m <= 25000.0: tmp = math.pow((l * (t_2 / math.pow(k_m, 2.0))), 2.0) elif k_m <= 5.2e+150: tmp = 2.0 / (((math.sin(k_m) * math.tan(k_m)) * ((math.pow(t_m, 2.0) / l) * (t_m / l))) * ((k_m / t_m) * (k_m / t_m))) else: tmp = math.pow((math.pow(k_m, -2.0) * (l * t_2)), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = sqrt(Float64(2.0 / t_m)) tmp = 0.0 if (k_m <= 25000.0) tmp = Float64(l * Float64(t_2 / (k_m ^ 2.0))) ^ 2.0; elseif (k_m <= 5.2e+150) tmp = Float64(2.0 / Float64(Float64(Float64(sin(k_m) * tan(k_m)) * Float64(Float64((t_m ^ 2.0) / l) * Float64(t_m / l))) * Float64(Float64(k_m / t_m) * Float64(k_m / t_m)))); else tmp = Float64((k_m ^ -2.0) * Float64(l * t_2)) ^ 2.0; end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) t_2 = sqrt((2.0 / t_m)); tmp = 0.0; if (k_m <= 25000.0) tmp = (l * (t_2 / (k_m ^ 2.0))) ^ 2.0; elseif (k_m <= 5.2e+150) tmp = 2.0 / (((sin(k_m) * tan(k_m)) * (((t_m ^ 2.0) / l) * (t_m / l))) * ((k_m / t_m) * (k_m / t_m))); else tmp = ((k_m ^ -2.0) * (l * t_2)) ^ 2.0; end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[Sqrt[N[(2.0 / t$95$m), $MachinePrecision]], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 25000.0], N[Power[N[(l * N[(t$95$2 / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[k$95$m, 5.2e+150], N[(2.0 / N[(N[(N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(N[(N[Power[t$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[Power[k$95$m, -2.0], $MachinePrecision] * N[(l * t$95$2), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sqrt{\frac{2}{t\_m}}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 25000:\\
\;\;\;\;{\left(\ell \cdot \frac{t\_2}{{k\_m}^{2}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 5.2 \cdot 10^{+150}:\\
\;\;\;\;\frac{2}{\left(\left(\sin k\_m \cdot \tan k\_m\right) \cdot \left(\frac{{t\_m}^{2}}{\ell} \cdot \frac{t\_m}{\ell}\right)\right) \cdot \left(\frac{k\_m}{t\_m} \cdot \frac{k\_m}{t\_m}\right)}\\
\mathbf{else}:\\
\;\;\;\;{\left({k\_m}^{-2} \cdot \left(\ell \cdot t\_2\right)\right)}^{2}\\
\end{array}
\end{array}
\end{array}
if k < 25000Initial program 29.1%
*-commutative29.1%
associate-/r*29.1%
Simplified34.4%
Applied egg-rr29.3%
Taylor expanded in k around 0 39.5%
associate-*l/39.1%
associate-*l*39.1%
associate-/l*39.5%
Simplified39.5%
associate-*r/39.1%
sqrt-unprod39.1%
div-inv39.1%
Applied egg-rr39.1%
associate-/l*39.5%
Simplified39.5%
if 25000 < k < 5.20000000000000012e150Initial program 17.6%
Simplified17.7%
unpow317.7%
times-frac25.1%
pow225.1%
Applied egg-rr25.1%
+-commutative25.1%
associate-+l-48.9%
metadata-eval48.9%
--rgt-identity48.9%
unpow248.9%
Applied egg-rr48.9%
if 5.20000000000000012e150 < k Initial program 44.4%
*-commutative44.4%
associate-/r*44.4%
Simplified50.0%
Applied egg-rr22.4%
Taylor expanded in k around 0 50.8%
associate-*l/50.8%
associate-*l*50.8%
associate-/l*50.8%
Simplified50.8%
pow150.8%
div-inv50.8%
sqrt-unprod50.8%
div-inv50.8%
pow-flip50.8%
metadata-eval50.8%
Applied egg-rr50.8%
unpow150.8%
associate-*r*50.8%
*-commutative50.8%
Simplified50.8%
Final simplification41.4%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (pow (* l (/ (sqrt (/ 2.0 t_m)) (pow k_m 2.0))) 2.0)))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * pow((l * (sqrt((2.0 / t_m)) / pow(k_m, 2.0))), 2.0);
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((l * (sqrt((2.0d0 / t_m)) / (k_m ** 2.0d0))) ** 2.0d0)
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * Math.pow((l * (Math.sqrt((2.0 / t_m)) / Math.pow(k_m, 2.0))), 2.0);
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * math.pow((l * (math.sqrt((2.0 / t_m)) / math.pow(k_m, 2.0))), 2.0)
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * (Float64(l * Float64(sqrt(Float64(2.0 / t_m)) / (k_m ^ 2.0))) ^ 2.0)) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((l * (sqrt((2.0 / t_m)) / (k_m ^ 2.0))) ^ 2.0); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[Power[N[(l * N[(N[Sqrt[N[(2.0 / t$95$m), $MachinePrecision]], $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot {\left(\ell \cdot \frac{\sqrt{\frac{2}{t\_m}}}{{k\_m}^{2}}\right)}^{2}
\end{array}
Initial program 28.9%
*-commutative28.9%
associate-/r*28.9%
Simplified35.2%
Applied egg-rr27.1%
Taylor expanded in k around 0 37.2%
associate-*l/36.9%
associate-*l*36.9%
associate-/l*37.2%
Simplified37.2%
associate-*r/36.9%
sqrt-unprod36.9%
div-inv36.9%
Applied egg-rr36.9%
associate-/l*37.2%
Simplified37.2%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (pow (* l (* (pow k_m -2.0) (sqrt (/ 2.0 t_m)))) 2.0)))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * pow((l * (pow(k_m, -2.0) * sqrt((2.0 / t_m)))), 2.0);
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((l * ((k_m ** (-2.0d0)) * sqrt((2.0d0 / t_m)))) ** 2.0d0)
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * Math.pow((l * (Math.pow(k_m, -2.0) * Math.sqrt((2.0 / t_m)))), 2.0);
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * math.pow((l * (math.pow(k_m, -2.0) * math.sqrt((2.0 / t_m)))), 2.0)
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * (Float64(l * Float64((k_m ^ -2.0) * sqrt(Float64(2.0 / t_m)))) ^ 2.0)) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((l * ((k_m ^ -2.0) * sqrt((2.0 / t_m)))) ^ 2.0); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[Power[N[(l * N[(N[Power[k$95$m, -2.0], $MachinePrecision] * N[Sqrt[N[(2.0 / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot {\left(\ell \cdot \left({k\_m}^{-2} \cdot \sqrt{\frac{2}{t\_m}}\right)\right)}^{2}
\end{array}
Initial program 28.9%
*-commutative28.9%
associate-/r*28.9%
Simplified35.2%
Applied egg-rr27.1%
Taylor expanded in k around 0 37.2%
associate-*l/36.9%
associate-*l*36.9%
associate-/l*37.2%
Simplified37.2%
pow137.2%
div-inv37.1%
sqrt-unprod37.1%
div-inv37.1%
pow-flip37.1%
metadata-eval37.1%
Applied egg-rr37.1%
unpow137.1%
Simplified37.1%
Final simplification37.1%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(*
(* l l)
(/
2.0
(* (pow k_m 4.0) (+ t_m (* (* t_m (pow k_m 2.0)) 0.16666666666666666)))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (2.0 / (pow(k_m, 4.0) * (t_m + ((t_m * pow(k_m, 2.0)) * 0.16666666666666666)))));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((l * l) * (2.0d0 / ((k_m ** 4.0d0) * (t_m + ((t_m * (k_m ** 2.0d0)) * 0.16666666666666666d0)))))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (2.0 / (Math.pow(k_m, 4.0) * (t_m + ((t_m * Math.pow(k_m, 2.0)) * 0.16666666666666666)))));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * ((l * l) * (2.0 / (math.pow(k_m, 4.0) * (t_m + ((t_m * math.pow(k_m, 2.0)) * 0.16666666666666666)))))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(Float64(l * l) * Float64(2.0 / Float64((k_m ^ 4.0) * Float64(t_m + Float64(Float64(t_m * (k_m ^ 2.0)) * 0.16666666666666666)))))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((l * l) * (2.0 / ((k_m ^ 4.0) * (t_m + ((t_m * (k_m ^ 2.0)) * 0.16666666666666666))))); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[Power[k$95$m, 4.0], $MachinePrecision] * N[(t$95$m + N[(N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(\ell \cdot \ell\right) \cdot \frac{2}{{k\_m}^{4} \cdot \left(t\_m + \left(t\_m \cdot {k\_m}^{2}\right) \cdot 0.16666666666666666\right)}\right)
\end{array}
Initial program 28.9%
Simplified37.6%
Taylor expanded in k around 0 60.6%
Final simplification60.6%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* 2.0 (/ (pow l 2.0) (* t_m (pow k_m 4.0))))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * (pow(l, 2.0) / (t_m * pow(k_m, 4.0))));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (2.0d0 * ((l ** 2.0d0) / (t_m * (k_m ** 4.0d0))))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (2.0 * (Math.pow(l, 2.0) / (t_m * Math.pow(k_m, 4.0))));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (2.0 * (math.pow(l, 2.0) / (t_m * math.pow(k_m, 4.0))))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(2.0 * Float64((l ^ 2.0) / Float64(t_m * (k_m ^ 4.0))))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (2.0 * ((l ^ 2.0) / (t_m * (k_m ^ 4.0)))); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] / N[(t$95$m * N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(2 \cdot \frac{{\ell}^{2}}{t\_m \cdot {k\_m}^{4}}\right)
\end{array}
Initial program 28.9%
Simplified37.6%
Taylor expanded in k around 0 60.6%
Final simplification60.6%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* (* l l) (/ 2.0 (* t_m (pow k_m 4.0))))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (2.0 / (t_m * pow(k_m, 4.0))));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((l * l) * (2.0d0 / (t_m * (k_m ** 4.0d0))))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (2.0 / (t_m * Math.pow(k_m, 4.0))));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * ((l * l) * (2.0 / (t_m * math.pow(k_m, 4.0))))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(Float64(l * l) * Float64(2.0 / Float64(t_m * (k_m ^ 4.0))))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((l * l) * (2.0 / (t_m * (k_m ^ 4.0)))); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(t$95$m * N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(\ell \cdot \ell\right) \cdot \frac{2}{t\_m \cdot {k\_m}^{4}}\right)
\end{array}
Initial program 28.9%
Simplified37.6%
Taylor expanded in k around 0 60.6%
Final simplification60.6%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* (* l l) (* 2.0 (/ (pow k_m -4.0) t_m)))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (2.0 * (pow(k_m, -4.0) / t_m)));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((l * l) * (2.0d0 * ((k_m ** (-4.0d0)) / t_m)))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (2.0 * (Math.pow(k_m, -4.0) / t_m)));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * ((l * l) * (2.0 * (math.pow(k_m, -4.0) / t_m)))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(Float64(l * l) * Float64(2.0 * Float64((k_m ^ -4.0) / t_m)))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((l * l) * (2.0 * ((k_m ^ -4.0) / t_m))); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(N[(l * l), $MachinePrecision] * N[(2.0 * N[(N[Power[k$95$m, -4.0], $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(\ell \cdot \ell\right) \cdot \left(2 \cdot \frac{{k\_m}^{-4}}{t\_m}\right)\right)
\end{array}
Initial program 28.9%
Simplified37.6%
Taylor expanded in k around 0 60.6%
div-inv60.6%
*-commutative60.6%
Applied egg-rr60.6%
associate-*r/60.6%
metadata-eval60.6%
associate-/r*60.6%
Simplified60.6%
div-inv60.6%
pow-flip60.6%
metadata-eval60.6%
Applied egg-rr60.6%
associate-*l/60.6%
metadata-eval60.6%
pow-sqr60.6%
unpow260.6%
associate-/l*60.6%
unpow260.6%
pow-sqr60.6%
metadata-eval60.6%
Simplified60.6%
Final simplification60.6%
herbie shell --seed 2024137
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10-)"
:precision binary64
(/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (- (+ 1.0 (pow (/ k t) 2.0)) 1.0))))