
(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 22 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 (cbrt (* (sin k_m) (tan k_m)))) (t_3 (pow (cbrt l) -2.0)))
(*
t_s
(if (<= k_m 11200000000000.0)
(pow
(* (* (/ l k_m) (/ (sqrt 2.0) (sin k_m))) (sqrt (/ (cos k_m) t_m)))
2.0)
(/
(pow (/ (/ (/ (sqrt 2.0) k_m) t_3) (log (exp t_2))) 2.0)
(* t_m (* t_3 t_2)))))))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 <= 11200000000000.0) {
tmp = pow((((l / k_m) * (sqrt(2.0) / sin(k_m))) * sqrt((cos(k_m) / t_m))), 2.0);
} else {
tmp = pow((((sqrt(2.0) / k_m) / t_3) / log(exp(t_2))), 2.0) / (t_m * (t_3 * t_2));
}
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 <= 11200000000000.0) {
tmp = Math.pow((((l / k_m) * (Math.sqrt(2.0) / Math.sin(k_m))) * Math.sqrt((Math.cos(k_m) / t_m))), 2.0);
} else {
tmp = Math.pow((((Math.sqrt(2.0) / k_m) / t_3) / Math.log(Math.exp(t_2))), 2.0) / (t_m * (t_3 * t_2));
}
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 <= 11200000000000.0) tmp = Float64(Float64(Float64(l / k_m) * Float64(sqrt(2.0) / sin(k_m))) * sqrt(Float64(cos(k_m) / t_m))) ^ 2.0; else tmp = Float64((Float64(Float64(Float64(sqrt(2.0) / k_m) / t_3) / log(exp(t_2))) ^ 2.0) / Float64(t_m * Float64(t_3 * t_2))); 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, 11200000000000.0], N[Power[N[(N[(N[(l / k$95$m), $MachinePrecision] * N[(N[Sqrt[2.0], $MachinePrecision] / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / t$95$3), $MachinePrecision] / N[Log[N[Exp[t$95$2], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] / N[(t$95$m * N[(t$95$3 * t$95$2), $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 11200000000000:\\
\;\;\;\;{\left(\left(\frac{\ell}{k\_m} \cdot \frac{\sqrt{2}}{\sin k\_m}\right) \cdot \sqrt{\frac{\cos k\_m}{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\left(\frac{\frac{\frac{\sqrt{2}}{k\_m}}{t\_3}}{\log \left(e^{t\_2}\right)}\right)}^{2}}{t\_m \cdot \left(t\_3 \cdot t\_2\right)}\\
\end{array}
\end{array}
\end{array}
if k < 1.12e13Initial program 37.8%
*-commutative37.8%
associate-/r*37.8%
Simplified43.9%
add-sqr-sqrt28.3%
pow228.3%
Applied egg-rr31.3%
Taylor expanded in k around inf 48.9%
times-frac50.8%
Simplified50.8%
if 1.12e13 < k Initial program 28.2%
*-commutative28.2%
associate-/r*28.6%
Simplified30.3%
add-sqr-sqrt30.3%
add-cube-cbrt30.2%
times-frac30.2%
Applied egg-rr76.2%
associate-/r/76.2%
associate-/r*76.2%
associate-/r/76.2%
Simplified76.2%
frac-times73.2%
associate-/r/73.2%
associate-/l*73.3%
div-inv73.3%
pow-flip73.3%
metadata-eval73.3%
Applied egg-rr73.2%
times-frac76.2%
associate-*r/76.2%
associate-/r/76.2%
associate-/r*76.2%
associate-*l/73.2%
Simplified90.0%
add-log-exp90.0%
Applied egg-rr90.0%
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 (pow (cbrt l) -2.0)))
(*
t_s
(if (<= k_m 11200000000000.0)
(pow
(* (* (/ l k_m) (/ (sqrt 2.0) (sin k_m))) (sqrt (/ (cos k_m) t_m)))
2.0)
(/
(pow (/ (/ (sqrt 2.0) k_m) (* t_2 (cbrt (* (sin k_m) (tan k_m))))) 2.0)
(* t_m (* t_2 (* (cbrt (tan k_m)) (cbrt (sin 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 t_2 = pow(cbrt(l), -2.0);
double tmp;
if (k_m <= 11200000000000.0) {
tmp = pow((((l / k_m) * (sqrt(2.0) / sin(k_m))) * sqrt((cos(k_m) / t_m))), 2.0);
} else {
tmp = pow(((sqrt(2.0) / k_m) / (t_2 * cbrt((sin(k_m) * tan(k_m))))), 2.0) / (t_m * (t_2 * (cbrt(tan(k_m)) * cbrt(sin(k_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.pow(Math.cbrt(l), -2.0);
double tmp;
if (k_m <= 11200000000000.0) {
tmp = Math.pow((((l / k_m) * (Math.sqrt(2.0) / Math.sin(k_m))) * Math.sqrt((Math.cos(k_m) / t_m))), 2.0);
} else {
tmp = Math.pow(((Math.sqrt(2.0) / k_m) / (t_2 * Math.cbrt((Math.sin(k_m) * Math.tan(k_m))))), 2.0) / (t_m * (t_2 * (Math.cbrt(Math.tan(k_m)) * Math.cbrt(Math.sin(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) t_2 = cbrt(l) ^ -2.0 tmp = 0.0 if (k_m <= 11200000000000.0) tmp = Float64(Float64(Float64(l / k_m) * Float64(sqrt(2.0) / sin(k_m))) * sqrt(Float64(cos(k_m) / t_m))) ^ 2.0; else tmp = Float64((Float64(Float64(sqrt(2.0) / k_m) / Float64(t_2 * cbrt(Float64(sin(k_m) * tan(k_m))))) ^ 2.0) / Float64(t_m * Float64(t_2 * Float64(cbrt(tan(k_m)) * cbrt(sin(k_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[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 11200000000000.0], N[Power[N[(N[(N[(l / k$95$m), $MachinePrecision] * N[(N[Sqrt[2.0], $MachinePrecision] / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[Power[N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / N[(t$95$2 * N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] / N[(t$95$m * N[(t$95$2 * N[(N[Power[N[Tan[k$95$m], $MachinePrecision], 1/3], $MachinePrecision] * N[Power[N[Sin[k$95$m], $MachinePrecision], 1/3], $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)
\\
\begin{array}{l}
t_2 := {\left(\sqrt[3]{\ell}\right)}^{-2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 11200000000000:\\
\;\;\;\;{\left(\left(\frac{\ell}{k\_m} \cdot \frac{\sqrt{2}}{\sin k\_m}\right) \cdot \sqrt{\frac{\cos k\_m}{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\left(\frac{\frac{\sqrt{2}}{k\_m}}{t\_2 \cdot \sqrt[3]{\sin k\_m \cdot \tan k\_m}}\right)}^{2}}{t\_m \cdot \left(t\_2 \cdot \left(\sqrt[3]{\tan k\_m} \cdot \sqrt[3]{\sin k\_m}\right)\right)}\\
\end{array}
\end{array}
\end{array}
if k < 1.12e13Initial program 37.8%
*-commutative37.8%
associate-/r*37.8%
Simplified43.9%
add-sqr-sqrt28.3%
pow228.3%
Applied egg-rr31.3%
Taylor expanded in k around inf 48.9%
times-frac50.8%
Simplified50.8%
if 1.12e13 < k Initial program 28.2%
*-commutative28.2%
associate-/r*28.6%
Simplified30.3%
add-sqr-sqrt30.3%
add-cube-cbrt30.2%
times-frac30.2%
Applied egg-rr76.2%
associate-/r/76.2%
associate-/r*76.2%
associate-/r/76.2%
Simplified76.2%
frac-times73.2%
associate-/r/73.2%
associate-/l*73.3%
div-inv73.3%
pow-flip73.3%
metadata-eval73.3%
Applied egg-rr73.2%
times-frac76.2%
associate-*r/76.2%
associate-/r/76.2%
associate-/r*76.2%
associate-*l/73.2%
Simplified90.0%
*-commutative90.0%
cbrt-prod90.0%
Applied egg-rr90.0%
unpow290.0%
associate-/l/90.0%
*-commutative90.0%
associate-/l/90.0%
*-commutative90.0%
Applied egg-rr90.0%
unpow290.0%
Simplified90.0%
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 (pow (cbrt l) -2.0)) (t_3 (cbrt (* (sin k_m) (tan k_m)))))
(*
t_s
(if (<= k_m 105000000000.0)
(pow
(* (* (/ l k_m) (/ (sqrt 2.0) (sin k_m))) (sqrt (/ (cos k_m) t_m)))
2.0)
(/ (pow (/ (/ (/ (sqrt 2.0) k_m) t_2) t_3) 2.0) (* t_m (* t_2 t_3)))))))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 = pow(cbrt(l), -2.0);
double t_3 = cbrt((sin(k_m) * tan(k_m)));
double tmp;
if (k_m <= 105000000000.0) {
tmp = pow((((l / k_m) * (sqrt(2.0) / sin(k_m))) * sqrt((cos(k_m) / t_m))), 2.0);
} else {
tmp = pow((((sqrt(2.0) / k_m) / t_2) / t_3), 2.0) / (t_m * (t_2 * t_3));
}
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.pow(Math.cbrt(l), -2.0);
double t_3 = Math.cbrt((Math.sin(k_m) * Math.tan(k_m)));
double tmp;
if (k_m <= 105000000000.0) {
tmp = Math.pow((((l / k_m) * (Math.sqrt(2.0) / Math.sin(k_m))) * Math.sqrt((Math.cos(k_m) / t_m))), 2.0);
} else {
tmp = Math.pow((((Math.sqrt(2.0) / k_m) / t_2) / t_3), 2.0) / (t_m * (t_2 * t_3));
}
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(l) ^ -2.0 t_3 = cbrt(Float64(sin(k_m) * tan(k_m))) tmp = 0.0 if (k_m <= 105000000000.0) tmp = Float64(Float64(Float64(l / k_m) * Float64(sqrt(2.0) / sin(k_m))) * sqrt(Float64(cos(k_m) / t_m))) ^ 2.0; else tmp = Float64((Float64(Float64(Float64(sqrt(2.0) / k_m) / t_2) / t_3) ^ 2.0) / Float64(t_m * Float64(t_2 * t_3))); 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[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]}, Block[{t$95$3 = N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 105000000000.0], N[Power[N[(N[(N[(l / k$95$m), $MachinePrecision] * N[(N[Sqrt[2.0], $MachinePrecision] / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / t$95$2), $MachinePrecision] / t$95$3), $MachinePrecision], 2.0], $MachinePrecision] / N[(t$95$m * N[(t$95$2 * t$95$3), $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 := {\left(\sqrt[3]{\ell}\right)}^{-2}\\
t_3 := \sqrt[3]{\sin k\_m \cdot \tan k\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 105000000000:\\
\;\;\;\;{\left(\left(\frac{\ell}{k\_m} \cdot \frac{\sqrt{2}}{\sin k\_m}\right) \cdot \sqrt{\frac{\cos k\_m}{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\left(\frac{\frac{\frac{\sqrt{2}}{k\_m}}{t\_2}}{t\_3}\right)}^{2}}{t\_m \cdot \left(t\_2 \cdot t\_3\right)}\\
\end{array}
\end{array}
\end{array}
if k < 1.05e11Initial program 37.8%
*-commutative37.8%
associate-/r*37.8%
Simplified43.9%
add-sqr-sqrt28.3%
pow228.3%
Applied egg-rr31.3%
Taylor expanded in k around inf 48.9%
times-frac50.8%
Simplified50.8%
if 1.05e11 < k Initial program 28.2%
*-commutative28.2%
associate-/r*28.6%
Simplified30.3%
add-sqr-sqrt30.3%
add-cube-cbrt30.2%
times-frac30.2%
Applied egg-rr76.2%
associate-/r/76.2%
associate-/r*76.2%
associate-/r/76.2%
Simplified76.2%
frac-times73.2%
associate-/r/73.2%
associate-/l*73.3%
div-inv73.3%
pow-flip73.3%
metadata-eval73.3%
Applied egg-rr73.2%
times-frac76.2%
associate-*r/76.2%
associate-/r/76.2%
associate-/r*76.2%
associate-*l/73.2%
Simplified90.0%
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 (* (pow (cbrt l) -2.0) (cbrt (* (sin k_m) (tan k_m))))))
(*
t_s
(if (<= k_m 11200000000000.0)
(pow
(* (* (/ l k_m) (/ (sqrt 2.0) (sin k_m))) (sqrt (/ (cos k_m) t_m)))
2.0)
(/ (pow (/ (/ (sqrt 2.0) k_m) t_2) 2.0) (* t_m t_2))))))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 = pow(cbrt(l), -2.0) * cbrt((sin(k_m) * tan(k_m)));
double tmp;
if (k_m <= 11200000000000.0) {
tmp = pow((((l / k_m) * (sqrt(2.0) / sin(k_m))) * sqrt((cos(k_m) / t_m))), 2.0);
} else {
tmp = pow(((sqrt(2.0) / k_m) / t_2), 2.0) / (t_m * t_2);
}
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.pow(Math.cbrt(l), -2.0) * Math.cbrt((Math.sin(k_m) * Math.tan(k_m)));
double tmp;
if (k_m <= 11200000000000.0) {
tmp = Math.pow((((l / k_m) * (Math.sqrt(2.0) / Math.sin(k_m))) * Math.sqrt((Math.cos(k_m) / t_m))), 2.0);
} else {
tmp = Math.pow(((Math.sqrt(2.0) / k_m) / t_2), 2.0) / (t_m * t_2);
}
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((cbrt(l) ^ -2.0) * cbrt(Float64(sin(k_m) * tan(k_m)))) tmp = 0.0 if (k_m <= 11200000000000.0) tmp = Float64(Float64(Float64(l / k_m) * Float64(sqrt(2.0) / sin(k_m))) * sqrt(Float64(cos(k_m) / t_m))) ^ 2.0; else tmp = Float64((Float64(Float64(sqrt(2.0) / k_m) / t_2) ^ 2.0) / Float64(t_m * t_2)); 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[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $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, 11200000000000.0], N[Power[N[(N[(N[(l / k$95$m), $MachinePrecision] * N[(N[Sqrt[2.0], $MachinePrecision] / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[Power[N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / t$95$2), $MachinePrecision], 2.0], $MachinePrecision] / N[(t$95$m * t$95$2), $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 := {\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 11200000000000:\\
\;\;\;\;{\left(\left(\frac{\ell}{k\_m} \cdot \frac{\sqrt{2}}{\sin k\_m}\right) \cdot \sqrt{\frac{\cos k\_m}{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\left(\frac{\frac{\sqrt{2}}{k\_m}}{t\_2}\right)}^{2}}{t\_m \cdot t\_2}\\
\end{array}
\end{array}
\end{array}
if k < 1.12e13Initial program 37.8%
*-commutative37.8%
associate-/r*37.8%
Simplified43.9%
add-sqr-sqrt28.3%
pow228.3%
Applied egg-rr31.3%
Taylor expanded in k around inf 48.9%
times-frac50.8%
Simplified50.8%
if 1.12e13 < k Initial program 28.2%
*-commutative28.2%
associate-/r*28.6%
Simplified30.3%
add-sqr-sqrt30.3%
add-cube-cbrt30.2%
times-frac30.2%
Applied egg-rr76.2%
associate-/r/76.2%
associate-/r*76.2%
associate-/r/76.2%
Simplified76.2%
frac-times73.2%
associate-/r/73.2%
associate-/l*73.3%
div-inv73.3%
pow-flip73.3%
metadata-eval73.3%
Applied egg-rr73.2%
times-frac76.2%
associate-*r/76.2%
associate-/r/76.2%
associate-/r*76.2%
associate-*l/73.2%
Simplified90.0%
div-inv88.6%
associate-/l/88.6%
*-commutative88.6%
Applied egg-rr88.6%
associate-*r/89.9%
*-rgt-identity89.9%
Simplified89.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
(let* ((t_2 (pow (cbrt l) -2.0)) (t_3 (cbrt (* (sin k_m) (tan k_m)))))
(*
t_s
(if (<= k_m 11200000000000.0)
(pow
(* (* (/ l k_m) (/ (sqrt 2.0) (sin k_m))) (sqrt (/ (cos k_m) t_m)))
2.0)
(/ (/ (pow (/ (/ (/ (sqrt 2.0) k_m) t_3) t_2) 2.0) t_m) (* t_2 t_3))))))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 = pow(cbrt(l), -2.0);
double t_3 = cbrt((sin(k_m) * tan(k_m)));
double tmp;
if (k_m <= 11200000000000.0) {
tmp = pow((((l / k_m) * (sqrt(2.0) / sin(k_m))) * sqrt((cos(k_m) / t_m))), 2.0);
} else {
tmp = (pow((((sqrt(2.0) / k_m) / t_3) / t_2), 2.0) / t_m) / (t_2 * t_3);
}
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.pow(Math.cbrt(l), -2.0);
double t_3 = Math.cbrt((Math.sin(k_m) * Math.tan(k_m)));
double tmp;
if (k_m <= 11200000000000.0) {
tmp = Math.pow((((l / k_m) * (Math.sqrt(2.0) / Math.sin(k_m))) * Math.sqrt((Math.cos(k_m) / t_m))), 2.0);
} else {
tmp = (Math.pow((((Math.sqrt(2.0) / k_m) / t_3) / t_2), 2.0) / t_m) / (t_2 * t_3);
}
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(l) ^ -2.0 t_3 = cbrt(Float64(sin(k_m) * tan(k_m))) tmp = 0.0 if (k_m <= 11200000000000.0) tmp = Float64(Float64(Float64(l / k_m) * Float64(sqrt(2.0) / sin(k_m))) * sqrt(Float64(cos(k_m) / t_m))) ^ 2.0; else tmp = Float64(Float64((Float64(Float64(Float64(sqrt(2.0) / k_m) / t_3) / t_2) ^ 2.0) / t_m) / Float64(t_2 * t_3)); 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[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]}, Block[{t$95$3 = N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 11200000000000.0], N[Power[N[(N[(N[(l / k$95$m), $MachinePrecision] * N[(N[Sqrt[2.0], $MachinePrecision] / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[(N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] / t$95$3), $MachinePrecision] / t$95$2), $MachinePrecision], 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] / N[(t$95$2 * t$95$3), $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 := {\left(\sqrt[3]{\ell}\right)}^{-2}\\
t_3 := \sqrt[3]{\sin k\_m \cdot \tan k\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 11200000000000:\\
\;\;\;\;{\left(\left(\frac{\ell}{k\_m} \cdot \frac{\sqrt{2}}{\sin k\_m}\right) \cdot \sqrt{\frac{\cos k\_m}{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{{\left(\frac{\frac{\frac{\sqrt{2}}{k\_m}}{t\_3}}{t\_2}\right)}^{2}}{t\_m}}{t\_2 \cdot t\_3}\\
\end{array}
\end{array}
\end{array}
if k < 1.12e13Initial program 37.8%
*-commutative37.8%
associate-/r*37.8%
Simplified43.9%
add-sqr-sqrt28.3%
pow228.3%
Applied egg-rr31.3%
Taylor expanded in k around inf 48.9%
times-frac50.8%
Simplified50.8%
if 1.12e13 < k Initial program 28.2%
*-commutative28.2%
associate-/r*28.6%
Simplified30.3%
add-sqr-sqrt30.3%
add-cube-cbrt30.2%
times-frac30.2%
Applied egg-rr76.2%
associate-/r/76.2%
associate-/r*76.2%
associate-/r/76.2%
Simplified76.2%
frac-times73.2%
associate-/r/73.2%
associate-/l*73.3%
div-inv73.3%
pow-flip73.3%
metadata-eval73.3%
Applied egg-rr73.2%
times-frac76.2%
associate-*r/76.2%
associate-/r/76.2%
associate-/r*76.2%
associate-*l/73.2%
Simplified90.0%
*-commutative90.0%
cbrt-prod90.0%
Applied egg-rr90.0%
div-inv88.7%
associate-/l/88.7%
cbrt-prod88.6%
*-commutative88.6%
*-commutative88.6%
*-commutative88.6%
cbrt-prod88.7%
Applied egg-rr88.6%
associate-*r/89.9%
*-rgt-identity89.9%
associate-/r*88.6%
*-commutative88.6%
associate-/r*88.7%
Simplified88.7%
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 (/ (cos k_m) t_m)))
(*
t_s
(if (<= k_m 114000000000.0)
(pow (* (* (/ l k_m) (/ (sqrt 2.0) (sin k_m))) (sqrt t_2)) 2.0)
(if (<= k_m 6e+160)
(* (* 2.0 (* t_2 (pow (* k_m (sin k_m)) -2.0))) (* l l))
(pow
(/
(cbrt (* 2.0 (pow (/ k_m t_m) -2.0)))
(* (cbrt (* (sin k_m) (tan k_m))) (* t_m (pow (cbrt l) -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 t_2 = cos(k_m) / t_m;
double tmp;
if (k_m <= 114000000000.0) {
tmp = pow((((l / k_m) * (sqrt(2.0) / sin(k_m))) * sqrt(t_2)), 2.0);
} else if (k_m <= 6e+160) {
tmp = (2.0 * (t_2 * pow((k_m * sin(k_m)), -2.0))) * (l * l);
} else {
tmp = pow((cbrt((2.0 * pow((k_m / t_m), -2.0))) / (cbrt((sin(k_m) * tan(k_m))) * (t_m * pow(cbrt(l), -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 t_2 = Math.cos(k_m) / t_m;
double tmp;
if (k_m <= 114000000000.0) {
tmp = Math.pow((((l / k_m) * (Math.sqrt(2.0) / Math.sin(k_m))) * Math.sqrt(t_2)), 2.0);
} else if (k_m <= 6e+160) {
tmp = (2.0 * (t_2 * Math.pow((k_m * Math.sin(k_m)), -2.0))) * (l * l);
} else {
tmp = Math.pow((Math.cbrt((2.0 * Math.pow((k_m / t_m), -2.0))) / (Math.cbrt((Math.sin(k_m) * Math.tan(k_m))) * (t_m * Math.pow(Math.cbrt(l), -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) t_2 = Float64(cos(k_m) / t_m) tmp = 0.0 if (k_m <= 114000000000.0) tmp = Float64(Float64(Float64(l / k_m) * Float64(sqrt(2.0) / sin(k_m))) * sqrt(t_2)) ^ 2.0; elseif (k_m <= 6e+160) tmp = Float64(Float64(2.0 * Float64(t_2 * (Float64(k_m * sin(k_m)) ^ -2.0))) * Float64(l * l)); else tmp = Float64(cbrt(Float64(2.0 * (Float64(k_m / t_m) ^ -2.0))) / Float64(cbrt(Float64(sin(k_m) * tan(k_m))) * Float64(t_m * (cbrt(l) ^ -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_] := Block[{t$95$2 = N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 114000000000.0], N[Power[N[(N[(N[(l / k$95$m), $MachinePrecision] * N[(N[Sqrt[2.0], $MachinePrecision] / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[t$95$2], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[k$95$m, 6e+160], N[(N[(2.0 * N[(t$95$2 * N[Power[N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l * l), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[Power[N[(2.0 * N[Power[N[(k$95$m / t$95$m), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] / N[(N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * N[(t$95$m * N[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.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 := \frac{\cos k\_m}{t\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 114000000000:\\
\;\;\;\;{\left(\left(\frac{\ell}{k\_m} \cdot \frac{\sqrt{2}}{\sin k\_m}\right) \cdot \sqrt{t\_2}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 6 \cdot 10^{+160}:\\
\;\;\;\;\left(2 \cdot \left(t\_2 \cdot {\left(k\_m \cdot \sin k\_m\right)}^{-2}\right)\right) \cdot \left(\ell \cdot \ell\right)\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{\sqrt[3]{2 \cdot {\left(\frac{k\_m}{t\_m}\right)}^{-2}}}{\sqrt[3]{\sin k\_m \cdot \tan k\_m} \cdot \left(t\_m \cdot {\left(\sqrt[3]{\ell}\right)}^{-2}\right)}\right)}^{3}\\
\end{array}
\end{array}
\end{array}
if k < 1.14e11Initial program 37.8%
*-commutative37.8%
associate-/r*37.8%
Simplified43.9%
add-sqr-sqrt28.3%
pow228.3%
Applied egg-rr31.3%
Taylor expanded in k around inf 48.9%
times-frac50.8%
Simplified50.8%
if 1.14e11 < k < 5.9999999999999997e160Initial program 27.6%
Simplified35.3%
Taylor expanded in t around 0 85.7%
associate-*r/85.7%
associate-*r*85.8%
*-commutative85.8%
Simplified85.8%
associate-/l*85.8%
associate-*r*85.8%
pow-prod-down85.7%
Applied egg-rr85.7%
*-commutative85.7%
associate-/r*86.0%
*-commutative86.0%
Simplified86.0%
div-inv86.0%
pow-flip86.0%
metadata-eval86.0%
Applied egg-rr86.0%
if 5.9999999999999997e160 < k Initial program 29.3%
*-commutative29.3%
associate-/r*29.3%
Simplified29.3%
add-cube-cbrt29.3%
pow329.3%
cbrt-prod29.3%
cbrt-div29.3%
rem-cbrt-cube46.2%
cbrt-prod58.2%
pow258.2%
Applied egg-rr58.2%
add-cube-cbrt58.0%
pow358.0%
Applied egg-rr63.6%
Final simplification57.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
(let* ((t_2 (/ (cos k_m) t_m)))
(*
t_s
(if (<= k_m 11200000000000.0)
(pow (* (* (/ l k_m) (/ (sqrt 2.0) (sin k_m))) (sqrt t_2)) 2.0)
(* (* 2.0 (* t_2 (pow (* k_m (sin 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 t_2 = cos(k_m) / t_m;
double tmp;
if (k_m <= 11200000000000.0) {
tmp = pow((((l / k_m) * (sqrt(2.0) / sin(k_m))) * sqrt(t_2)), 2.0);
} else {
tmp = (2.0 * (t_2 * pow((k_m * sin(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) :: t_2
real(8) :: tmp
t_2 = cos(k_m) / t_m
if (k_m <= 11200000000000.0d0) then
tmp = (((l / k_m) * (sqrt(2.0d0) / sin(k_m))) * sqrt(t_2)) ** 2.0d0
else
tmp = (2.0d0 * (t_2 * ((k_m * sin(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 t_2 = Math.cos(k_m) / t_m;
double tmp;
if (k_m <= 11200000000000.0) {
tmp = Math.pow((((l / k_m) * (Math.sqrt(2.0) / Math.sin(k_m))) * Math.sqrt(t_2)), 2.0);
} else {
tmp = (2.0 * (t_2 * Math.pow((k_m * Math.sin(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): t_2 = math.cos(k_m) / t_m tmp = 0 if k_m <= 11200000000000.0: tmp = math.pow((((l / k_m) * (math.sqrt(2.0) / math.sin(k_m))) * math.sqrt(t_2)), 2.0) else: tmp = (2.0 * (t_2 * math.pow((k_m * math.sin(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) t_2 = Float64(cos(k_m) / t_m) tmp = 0.0 if (k_m <= 11200000000000.0) tmp = Float64(Float64(Float64(l / k_m) * Float64(sqrt(2.0) / sin(k_m))) * sqrt(t_2)) ^ 2.0; else tmp = Float64(Float64(2.0 * Float64(t_2 * (Float64(k_m * sin(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) t_2 = cos(k_m) / t_m; tmp = 0.0; if (k_m <= 11200000000000.0) tmp = (((l / k_m) * (sqrt(2.0) / sin(k_m))) * sqrt(t_2)) ^ 2.0; else tmp = (2.0 * (t_2 * ((k_m * sin(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_] := Block[{t$95$2 = N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 11200000000000.0], N[Power[N[(N[(N[(l / k$95$m), $MachinePrecision] * N[(N[Sqrt[2.0], $MachinePrecision] / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[t$95$2], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[(2.0 * N[(t$95$2 * N[Power[N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], -2.0], $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)
\\
\begin{array}{l}
t_2 := \frac{\cos k\_m}{t\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 11200000000000:\\
\;\;\;\;{\left(\left(\frac{\ell}{k\_m} \cdot \frac{\sqrt{2}}{\sin k\_m}\right) \cdot \sqrt{t\_2}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \left(t\_2 \cdot {\left(k\_m \cdot \sin k\_m\right)}^{-2}\right)\right) \cdot \left(\ell \cdot \ell\right)\\
\end{array}
\end{array}
\end{array}
if k < 1.12e13Initial program 37.8%
*-commutative37.8%
associate-/r*37.8%
Simplified43.9%
add-sqr-sqrt28.3%
pow228.3%
Applied egg-rr31.3%
Taylor expanded in k around inf 48.9%
times-frac50.8%
Simplified50.8%
if 1.12e13 < k Initial program 28.2%
Simplified33.0%
Taylor expanded in t around 0 71.1%
associate-*r/71.1%
associate-*r*71.1%
*-commutative71.1%
Simplified71.1%
associate-/l*71.1%
associate-*r*71.1%
pow-prod-down71.1%
Applied egg-rr71.1%
*-commutative71.1%
associate-/r*71.3%
*-commutative71.3%
Simplified71.3%
div-inv71.2%
pow-flip71.3%
metadata-eval71.3%
Applied egg-rr71.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 (* k_m (sin k_m))) (t_3 (/ (cos k_m) t_m)))
(*
t_s
(if (<= k_m 105000.0)
(pow (* (sqrt t_3) (* l (/ (sqrt 2.0) t_2))) 2.0)
(* (* 2.0 (* t_3 (pow t_2 -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 t_2 = k_m * sin(k_m);
double t_3 = cos(k_m) / t_m;
double tmp;
if (k_m <= 105000.0) {
tmp = pow((sqrt(t_3) * (l * (sqrt(2.0) / t_2))), 2.0);
} else {
tmp = (2.0 * (t_3 * pow(t_2, -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) :: t_2
real(8) :: t_3
real(8) :: tmp
t_2 = k_m * sin(k_m)
t_3 = cos(k_m) / t_m
if (k_m <= 105000.0d0) then
tmp = (sqrt(t_3) * (l * (sqrt(2.0d0) / t_2))) ** 2.0d0
else
tmp = (2.0d0 * (t_3 * (t_2 ** (-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 t_2 = k_m * Math.sin(k_m);
double t_3 = Math.cos(k_m) / t_m;
double tmp;
if (k_m <= 105000.0) {
tmp = Math.pow((Math.sqrt(t_3) * (l * (Math.sqrt(2.0) / t_2))), 2.0);
} else {
tmp = (2.0 * (t_3 * Math.pow(t_2, -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): t_2 = k_m * math.sin(k_m) t_3 = math.cos(k_m) / t_m tmp = 0 if k_m <= 105000.0: tmp = math.pow((math.sqrt(t_3) * (l * (math.sqrt(2.0) / t_2))), 2.0) else: tmp = (2.0 * (t_3 * math.pow(t_2, -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) t_2 = Float64(k_m * sin(k_m)) t_3 = Float64(cos(k_m) / t_m) tmp = 0.0 if (k_m <= 105000.0) tmp = Float64(sqrt(t_3) * Float64(l * Float64(sqrt(2.0) / t_2))) ^ 2.0; else tmp = Float64(Float64(2.0 * Float64(t_3 * (t_2 ^ -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) t_2 = k_m * sin(k_m); t_3 = cos(k_m) / t_m; tmp = 0.0; if (k_m <= 105000.0) tmp = (sqrt(t_3) * (l * (sqrt(2.0) / t_2))) ^ 2.0; else tmp = (2.0 * (t_3 * (t_2 ^ -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_] := Block[{t$95$2 = N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 105000.0], N[Power[N[(N[Sqrt[t$95$3], $MachinePrecision] * N[(l * N[(N[Sqrt[2.0], $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[(2.0 * N[(t$95$3 * N[Power[t$95$2, -2.0], $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)
\\
\begin{array}{l}
t_2 := k\_m \cdot \sin k\_m\\
t_3 := \frac{\cos k\_m}{t\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 105000:\\
\;\;\;\;{\left(\sqrt{t\_3} \cdot \left(\ell \cdot \frac{\sqrt{2}}{t\_2}\right)\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \left(t\_3 \cdot {t\_2}^{-2}\right)\right) \cdot \left(\ell \cdot \ell\right)\\
\end{array}
\end{array}
\end{array}
if k < 105000Initial program 38.0%
*-commutative38.0%
associate-/r*38.0%
Simplified43.6%
add-sqr-sqrt28.4%
pow228.4%
Applied egg-rr31.5%
Taylor expanded in k around inf 49.1%
associate-/l*49.1%
Simplified49.1%
if 105000 < k Initial program 27.8%
Simplified34.0%
Taylor expanded in t around 0 71.5%
associate-*r/71.5%
associate-*r*71.5%
*-commutative71.5%
Simplified71.5%
associate-/l*71.5%
associate-*r*71.5%
pow-prod-down71.5%
Applied egg-rr71.5%
*-commutative71.5%
associate-/r*71.7%
*-commutative71.7%
Simplified71.7%
div-inv71.6%
pow-flip71.7%
metadata-eval71.7%
Applied egg-rr71.7%
Final simplification54.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 7.2e-6)
(pow (* (/ (* l (sqrt 2.0)) (pow k_m 2.0)) (sqrt (/ 1.0 t_m))) 2.0)
(* (* 2.0 (* (/ (cos k_m) t_m) (pow (* k_m (sin 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 <= 7.2e-6) {
tmp = pow((((l * sqrt(2.0)) / pow(k_m, 2.0)) * sqrt((1.0 / t_m))), 2.0);
} else {
tmp = (2.0 * ((cos(k_m) / t_m) * pow((k_m * sin(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 <= 7.2d-6) then
tmp = (((l * sqrt(2.0d0)) / (k_m ** 2.0d0)) * sqrt((1.0d0 / t_m))) ** 2.0d0
else
tmp = (2.0d0 * ((cos(k_m) / t_m) * ((k_m * sin(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 <= 7.2e-6) {
tmp = Math.pow((((l * Math.sqrt(2.0)) / Math.pow(k_m, 2.0)) * Math.sqrt((1.0 / t_m))), 2.0);
} else {
tmp = (2.0 * ((Math.cos(k_m) / t_m) * Math.pow((k_m * Math.sin(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 <= 7.2e-6: tmp = math.pow((((l * math.sqrt(2.0)) / math.pow(k_m, 2.0)) * math.sqrt((1.0 / t_m))), 2.0) else: tmp = (2.0 * ((math.cos(k_m) / t_m) * math.pow((k_m * math.sin(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 <= 7.2e-6) tmp = Float64(Float64(Float64(l * sqrt(2.0)) / (k_m ^ 2.0)) * sqrt(Float64(1.0 / t_m))) ^ 2.0; else tmp = Float64(Float64(2.0 * Float64(Float64(cos(k_m) / t_m) * (Float64(k_m * sin(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 <= 7.2e-6) tmp = (((l * sqrt(2.0)) / (k_m ^ 2.0)) * sqrt((1.0 / t_m))) ^ 2.0; else tmp = (2.0 * ((cos(k_m) / t_m) * ((k_m * sin(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, 7.2e-6], N[Power[N[(N[(N[(l * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision] * N[Power[N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], -2.0], $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 7.2 \cdot 10^{-6}:\\
\;\;\;\;{\left(\frac{\ell \cdot \sqrt{2}}{{k\_m}^{2}} \cdot \sqrt{\frac{1}{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \left(\frac{\cos k\_m}{t\_m} \cdot {\left(k\_m \cdot \sin k\_m\right)}^{-2}\right)\right) \cdot \left(\ell \cdot \ell\right)\\
\end{array}
\end{array}
if k < 7.19999999999999967e-6Initial program 38.4%
*-commutative38.4%
associate-/r*38.4%
Simplified43.5%
add-sqr-sqrt28.2%
pow228.2%
Applied egg-rr30.7%
Taylor expanded in k around 0 41.2%
if 7.19999999999999967e-6 < k Initial program 27.0%
Simplified34.5%
Taylor expanded in t around 0 72.4%
associate-*r/72.4%
associate-*r*72.4%
*-commutative72.4%
Simplified72.4%
associate-/l*72.4%
associate-*r*72.4%
pow-prod-down72.3%
Applied egg-rr72.3%
*-commutative72.3%
associate-/r*72.5%
*-commutative72.5%
Simplified72.5%
div-inv72.5%
pow-flip72.5%
metadata-eval72.5%
Applied egg-rr72.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 6.8e-6)
(pow (* (sqrt (/ 1.0 t_m)) (* l (/ (sqrt 2.0) (pow k_m 2.0)))) 2.0)
(* (* 2.0 (* (/ (cos k_m) t_m) (pow (* k_m (sin 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 <= 6.8e-6) {
tmp = pow((sqrt((1.0 / t_m)) * (l * (sqrt(2.0) / pow(k_m, 2.0)))), 2.0);
} else {
tmp = (2.0 * ((cos(k_m) / t_m) * pow((k_m * sin(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 <= 6.8d-6) then
tmp = (sqrt((1.0d0 / t_m)) * (l * (sqrt(2.0d0) / (k_m ** 2.0d0)))) ** 2.0d0
else
tmp = (2.0d0 * ((cos(k_m) / t_m) * ((k_m * sin(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 <= 6.8e-6) {
tmp = Math.pow((Math.sqrt((1.0 / t_m)) * (l * (Math.sqrt(2.0) / Math.pow(k_m, 2.0)))), 2.0);
} else {
tmp = (2.0 * ((Math.cos(k_m) / t_m) * Math.pow((k_m * Math.sin(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 <= 6.8e-6: tmp = math.pow((math.sqrt((1.0 / t_m)) * (l * (math.sqrt(2.0) / math.pow(k_m, 2.0)))), 2.0) else: tmp = (2.0 * ((math.cos(k_m) / t_m) * math.pow((k_m * math.sin(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 <= 6.8e-6) tmp = Float64(sqrt(Float64(1.0 / t_m)) * Float64(l * Float64(sqrt(2.0) / (k_m ^ 2.0)))) ^ 2.0; else tmp = Float64(Float64(2.0 * Float64(Float64(cos(k_m) / t_m) * (Float64(k_m * sin(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 <= 6.8e-6) tmp = (sqrt((1.0 / t_m)) * (l * (sqrt(2.0) / (k_m ^ 2.0)))) ^ 2.0; else tmp = (2.0 * ((cos(k_m) / t_m) * ((k_m * sin(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, 6.8e-6], N[Power[N[(N[Sqrt[N[(1.0 / t$95$m), $MachinePrecision]], $MachinePrecision] * N[(l * N[(N[Sqrt[2.0], $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision] * N[Power[N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], -2.0], $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 6.8 \cdot 10^{-6}:\\
\;\;\;\;{\left(\sqrt{\frac{1}{t\_m}} \cdot \left(\ell \cdot \frac{\sqrt{2}}{{k\_m}^{2}}\right)\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \left(\frac{\cos k\_m}{t\_m} \cdot {\left(k\_m \cdot \sin k\_m\right)}^{-2}\right)\right) \cdot \left(\ell \cdot \ell\right)\\
\end{array}
\end{array}
if k < 6.80000000000000012e-6Initial program 38.4%
*-commutative38.4%
associate-/r*38.4%
Simplified43.5%
add-sqr-sqrt28.2%
pow228.2%
Applied egg-rr30.7%
Taylor expanded in k around 0 41.2%
associate-/l*41.2%
Simplified41.2%
if 6.80000000000000012e-6 < k Initial program 27.0%
Simplified34.5%
Taylor expanded in t around 0 72.4%
associate-*r/72.4%
associate-*r*72.4%
*-commutative72.4%
Simplified72.4%
associate-/l*72.4%
associate-*r*72.4%
pow-prod-down72.3%
Applied egg-rr72.3%
*-commutative72.3%
associate-/r*72.5%
*-commutative72.5%
Simplified72.5%
div-inv72.5%
pow-flip72.5%
metadata-eval72.5%
Applied egg-rr72.5%
Final simplification49.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 6.5e-6)
(pow (* l (/ (* (sqrt 2.0) (sqrt (/ 1.0 t_m))) (pow k_m 2.0))) 2.0)
(* (* 2.0 (* (/ (cos k_m) t_m) (pow (* k_m (sin 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 <= 6.5e-6) {
tmp = pow((l * ((sqrt(2.0) * sqrt((1.0 / t_m))) / pow(k_m, 2.0))), 2.0);
} else {
tmp = (2.0 * ((cos(k_m) / t_m) * pow((k_m * sin(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 <= 6.5d-6) then
tmp = (l * ((sqrt(2.0d0) * sqrt((1.0d0 / t_m))) / (k_m ** 2.0d0))) ** 2.0d0
else
tmp = (2.0d0 * ((cos(k_m) / t_m) * ((k_m * sin(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 <= 6.5e-6) {
tmp = Math.pow((l * ((Math.sqrt(2.0) * Math.sqrt((1.0 / t_m))) / Math.pow(k_m, 2.0))), 2.0);
} else {
tmp = (2.0 * ((Math.cos(k_m) / t_m) * Math.pow((k_m * Math.sin(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 <= 6.5e-6: tmp = math.pow((l * ((math.sqrt(2.0) * math.sqrt((1.0 / t_m))) / math.pow(k_m, 2.0))), 2.0) else: tmp = (2.0 * ((math.cos(k_m) / t_m) * math.pow((k_m * math.sin(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 <= 6.5e-6) tmp = Float64(l * Float64(Float64(sqrt(2.0) * sqrt(Float64(1.0 / t_m))) / (k_m ^ 2.0))) ^ 2.0; else tmp = Float64(Float64(2.0 * Float64(Float64(cos(k_m) / t_m) * (Float64(k_m * sin(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 <= 6.5e-6) tmp = (l * ((sqrt(2.0) * sqrt((1.0 / t_m))) / (k_m ^ 2.0))) ^ 2.0; else tmp = (2.0 * ((cos(k_m) / t_m) * ((k_m * sin(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, 6.5e-6], N[Power[N[(l * N[(N[(N[Sqrt[2.0], $MachinePrecision] * N[Sqrt[N[(1.0 / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision] * N[Power[N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], -2.0], $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 6.5 \cdot 10^{-6}:\\
\;\;\;\;{\left(\ell \cdot \frac{\sqrt{2} \cdot \sqrt{\frac{1}{t\_m}}}{{k\_m}^{2}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \left(\frac{\cos k\_m}{t\_m} \cdot {\left(k\_m \cdot \sin k\_m\right)}^{-2}\right)\right) \cdot \left(\ell \cdot \ell\right)\\
\end{array}
\end{array}
if k < 6.4999999999999996e-6Initial program 38.4%
*-commutative38.4%
associate-/r*38.4%
Simplified43.5%
add-sqr-sqrt28.2%
pow228.2%
Applied egg-rr30.7%
Taylor expanded in k around 0 31.2%
*-commutative31.2%
Simplified31.2%
Taylor expanded in k around 0 41.2%
associate-*l/39.6%
associate-*l*39.6%
associate-/l*41.1%
Simplified41.1%
if 6.4999999999999996e-6 < k Initial program 27.0%
Simplified34.5%
Taylor expanded in t around 0 72.4%
associate-*r/72.4%
associate-*r*72.4%
*-commutative72.4%
Simplified72.4%
associate-/l*72.4%
associate-*r*72.4%
pow-prod-down72.3%
Applied egg-rr72.3%
*-commutative72.3%
associate-/r*72.5%
*-commutative72.5%
Simplified72.5%
div-inv72.5%
pow-flip72.5%
metadata-eval72.5%
Applied egg-rr72.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
(let* ((t_2 (* t_m (pow k_m 2.0))))
(*
t_s
(if (<= t_m 4.4e-183)
(* (* l l) (/ 2.0 (* (pow (sin k_m) 2.0) t_2)))
(if (<= t_m 6e+132)
(pow (/ (/ (sqrt 2.0) (/ k_m t_m)) (* k_m (/ (pow t_m 1.5) l))) 2.0)
(* (* l l) (/ 2.0 (* (pow k_m 2.0) t_2))))))))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(k_m, 2.0);
double tmp;
if (t_m <= 4.4e-183) {
tmp = (l * l) * (2.0 / (pow(sin(k_m), 2.0) * t_2));
} else if (t_m <= 6e+132) {
tmp = pow(((sqrt(2.0) / (k_m / t_m)) / (k_m * (pow(t_m, 1.5) / l))), 2.0);
} else {
tmp = (l * l) * (2.0 / (pow(k_m, 2.0) * t_2));
}
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 = t_m * (k_m ** 2.0d0)
if (t_m <= 4.4d-183) then
tmp = (l * l) * (2.0d0 / ((sin(k_m) ** 2.0d0) * t_2))
else if (t_m <= 6d+132) then
tmp = ((sqrt(2.0d0) / (k_m / t_m)) / (k_m * ((t_m ** 1.5d0) / l))) ** 2.0d0
else
tmp = (l * l) * (2.0d0 / ((k_m ** 2.0d0) * t_2))
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 = t_m * Math.pow(k_m, 2.0);
double tmp;
if (t_m <= 4.4e-183) {
tmp = (l * l) * (2.0 / (Math.pow(Math.sin(k_m), 2.0) * t_2));
} else if (t_m <= 6e+132) {
tmp = Math.pow(((Math.sqrt(2.0) / (k_m / t_m)) / (k_m * (Math.pow(t_m, 1.5) / l))), 2.0);
} else {
tmp = (l * l) * (2.0 / (Math.pow(k_m, 2.0) * t_2));
}
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 = t_m * math.pow(k_m, 2.0) tmp = 0 if t_m <= 4.4e-183: tmp = (l * l) * (2.0 / (math.pow(math.sin(k_m), 2.0) * t_2)) elif t_m <= 6e+132: tmp = math.pow(((math.sqrt(2.0) / (k_m / t_m)) / (k_m * (math.pow(t_m, 1.5) / l))), 2.0) else: tmp = (l * l) * (2.0 / (math.pow(k_m, 2.0) * t_2)) 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(t_m * (k_m ^ 2.0)) tmp = 0.0 if (t_m <= 4.4e-183) tmp = Float64(Float64(l * l) * Float64(2.0 / Float64((sin(k_m) ^ 2.0) * t_2))); elseif (t_m <= 6e+132) tmp = Float64(Float64(sqrt(2.0) / Float64(k_m / t_m)) / Float64(k_m * Float64((t_m ^ 1.5) / l))) ^ 2.0; else tmp = Float64(Float64(l * l) * Float64(2.0 / Float64((k_m ^ 2.0) * t_2))); 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 = t_m * (k_m ^ 2.0); tmp = 0.0; if (t_m <= 4.4e-183) tmp = (l * l) * (2.0 / ((sin(k_m) ^ 2.0) * t_2)); elseif (t_m <= 6e+132) tmp = ((sqrt(2.0) / (k_m / t_m)) / (k_m * ((t_m ^ 1.5) / l))) ^ 2.0; else tmp = (l * l) * (2.0 / ((k_m ^ 2.0) * t_2)); 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[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 4.4e-183], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 6e+132], N[Power[N[(N[(N[Sqrt[2.0], $MachinePrecision] / N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * t$95$2), $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 := t\_m \cdot {k\_m}^{2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 4.4 \cdot 10^{-183}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{{\sin k\_m}^{2} \cdot t\_2}\\
\mathbf{elif}\;t\_m \leq 6 \cdot 10^{+132}:\\
\;\;\;\;{\left(\frac{\frac{\sqrt{2}}{\frac{k\_m}{t\_m}}}{k\_m \cdot \frac{{t\_m}^{1.5}}{\ell}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{{k\_m}^{2} \cdot t\_2}\\
\end{array}
\end{array}
\end{array}
if t < 4.3999999999999999e-183Initial program 33.6%
Simplified38.3%
Taylor expanded in t around 0 70.3%
associate-*r/70.3%
associate-*r*70.3%
*-commutative70.3%
Simplified70.3%
Taylor expanded in k around 0 61.2%
if 4.3999999999999999e-183 < t < 5.9999999999999996e132Initial program 53.5%
*-commutative53.5%
associate-/r*53.9%
Simplified57.8%
add-sqr-sqrt53.3%
pow253.3%
Applied egg-rr65.0%
Taylor expanded in k around 0 69.4%
if 5.9999999999999996e132 < t Initial program 10.8%
Simplified24.6%
Taylor expanded in t around 0 76.0%
associate-*r/76.0%
associate-*r*76.0%
*-commutative76.0%
Simplified76.0%
Taylor expanded in k around 0 73.5%
Taylor expanded in k around 0 73.5%
Final simplification65.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
(if (<= (* l l) 0.0)
(pow (/ (/ (sqrt 2.0) (/ k_m t_m)) (* k_m (/ (pow t_m 1.5) l))) 2.0)
(* (* 2.0 (* (/ (cos k_m) t_m) (pow (* k_m (sin 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 ((l * l) <= 0.0) {
tmp = pow(((sqrt(2.0) / (k_m / t_m)) / (k_m * (pow(t_m, 1.5) / l))), 2.0);
} else {
tmp = (2.0 * ((cos(k_m) / t_m) * pow((k_m * sin(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 ((l * l) <= 0.0d0) then
tmp = ((sqrt(2.0d0) / (k_m / t_m)) / (k_m * ((t_m ** 1.5d0) / l))) ** 2.0d0
else
tmp = (2.0d0 * ((cos(k_m) / t_m) * ((k_m * sin(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 ((l * l) <= 0.0) {
tmp = Math.pow(((Math.sqrt(2.0) / (k_m / t_m)) / (k_m * (Math.pow(t_m, 1.5) / l))), 2.0);
} else {
tmp = (2.0 * ((Math.cos(k_m) / t_m) * Math.pow((k_m * Math.sin(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 (l * l) <= 0.0: tmp = math.pow(((math.sqrt(2.0) / (k_m / t_m)) / (k_m * (math.pow(t_m, 1.5) / l))), 2.0) else: tmp = (2.0 * ((math.cos(k_m) / t_m) * math.pow((k_m * math.sin(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 (Float64(l * l) <= 0.0) tmp = Float64(Float64(sqrt(2.0) / Float64(k_m / t_m)) / Float64(k_m * Float64((t_m ^ 1.5) / l))) ^ 2.0; else tmp = Float64(Float64(2.0 * Float64(Float64(cos(k_m) / t_m) * (Float64(k_m * sin(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 ((l * l) <= 0.0) tmp = ((sqrt(2.0) / (k_m / t_m)) / (k_m * ((t_m ^ 1.5) / l))) ^ 2.0; else tmp = (2.0 * ((cos(k_m) / t_m) * ((k_m * sin(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[N[(l * l), $MachinePrecision], 0.0], N[Power[N[(N[(N[Sqrt[2.0], $MachinePrecision] / N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision] / N[(k$95$m * N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision] * N[Power[N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision], -2.0], $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}\;\ell \cdot \ell \leq 0:\\
\;\;\;\;{\left(\frac{\frac{\sqrt{2}}{\frac{k\_m}{t\_m}}}{k\_m \cdot \frac{{t\_m}^{1.5}}{\ell}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \left(\frac{\cos k\_m}{t\_m} \cdot {\left(k\_m \cdot \sin k\_m\right)}^{-2}\right)\right) \cdot \left(\ell \cdot \ell\right)\\
\end{array}
\end{array}
if (*.f64 l l) < 0.0Initial program 15.8%
*-commutative15.8%
associate-/r*15.8%
Simplified26.0%
add-sqr-sqrt26.0%
pow226.0%
Applied egg-rr23.1%
Taylor expanded in k around 0 34.8%
if 0.0 < (*.f64 l l) Initial program 42.8%
Simplified46.6%
Taylor expanded in t around 0 79.1%
associate-*r/79.1%
associate-*r*79.1%
*-commutative79.1%
Simplified79.1%
associate-/l*79.1%
associate-*r*76.0%
pow-prod-down76.0%
Applied egg-rr76.0%
*-commutative76.0%
associate-/r*76.1%
*-commutative76.1%
Simplified76.1%
div-inv76.1%
pow-flip76.3%
metadata-eval76.3%
Applied egg-rr76.3%
Final simplification65.0%
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 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) {
return t_s * ((l * l) * ((2.0 / pow(k_m, 2.0)) / (t_m * pow(sin(k_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 * l) * ((2.0d0 / (k_m ** 2.0d0)) / (t_m * (sin(k_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 * ((l * l) * ((2.0 / Math.pow(k_m, 2.0)) / (t_m * Math.pow(Math.sin(k_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 * ((l * l) * ((2.0 / math.pow(k_m, 2.0)) / (t_m * math.pow(math.sin(k_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(Float64(l * l) * Float64(Float64(2.0 / (k_m ^ 2.0)) / Float64(t_m * (sin(k_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 * l) * ((2.0 / (k_m ^ 2.0)) / (t_m * (sin(k_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[(N[(l * l), $MachinePrecision] * N[(N[(2.0 / 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]
\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{\frac{2}{{k\_m}^{2}}}{t\_m \cdot {\sin k\_m}^{2}}\right)
\end{array}
Initial program 35.4%
Simplified41.3%
Taylor expanded in t around 0 71.8%
associate-*r/71.8%
associate-*r*71.8%
*-commutative71.8%
Simplified71.8%
Taylor expanded in k around 0 63.0%
Taylor expanded in k around inf 63.0%
associate-/r*63.0%
Simplified63.0%
Final simplification63.0%
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 (sin k_m) 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) {
return t_s * ((l * l) * (2.0 / (pow(sin(k_m), 2.0) * (t_m * pow(k_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 * l) * (2.0d0 / ((sin(k_m) ** 2.0d0) * (t_m * (k_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 * ((l * l) * (2.0 / (Math.pow(Math.sin(k_m), 2.0) * (t_m * Math.pow(k_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 * ((l * l) * (2.0 / (math.pow(math.sin(k_m), 2.0) * (t_m * math.pow(k_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(Float64(l * l) * Float64(2.0 / Float64((sin(k_m) ^ 2.0) * Float64(t_m * (k_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 * l) * (2.0 / ((sin(k_m) ^ 2.0) * (t_m * (k_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[(N[(l * l), $MachinePrecision] * N[(2.0 / 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]), $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}{{\sin k\_m}^{2} \cdot \left(t\_m \cdot {k\_m}^{2}\right)}\right)
\end{array}
Initial program 35.4%
Simplified41.3%
Taylor expanded in t around 0 71.8%
associate-*r/71.8%
associate-*r*71.8%
*-commutative71.8%
Simplified71.8%
Taylor expanded in k around 0 63.0%
Final simplification63.0%
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 0.019)
(* (* l l) (/ 2.0 (* (pow k_m 2.0) (* t_m (pow k_m 2.0)))))
(* (* l l) (/ (/ 2.0 t_m) (pow (* 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 (k_m <= 0.019) {
tmp = (l * l) * (2.0 / (pow(k_m, 2.0) * (t_m * pow(k_m, 2.0))));
} else {
tmp = (l * l) * ((2.0 / t_m) / pow((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 (k_m <= 0.019d0) then
tmp = (l * l) * (2.0d0 / ((k_m ** 2.0d0) * (t_m * (k_m ** 2.0d0))))
else
tmp = (l * l) * ((2.0d0 / t_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 (k_m <= 0.019) {
tmp = (l * l) * (2.0 / (Math.pow(k_m, 2.0) * (t_m * Math.pow(k_m, 2.0))));
} else {
tmp = (l * l) * ((2.0 / t_m) / Math.pow((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 k_m <= 0.019: tmp = (l * l) * (2.0 / (math.pow(k_m, 2.0) * (t_m * math.pow(k_m, 2.0)))) else: tmp = (l * l) * ((2.0 / t_m) / math.pow((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 (k_m <= 0.019) tmp = Float64(Float64(l * l) * Float64(2.0 / Float64((k_m ^ 2.0) * Float64(t_m * (k_m ^ 2.0))))); else tmp = Float64(Float64(l * l) * Float64(Float64(2.0 / t_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 (k_m <= 0.019) tmp = (l * l) * (2.0 / ((k_m ^ 2.0) * (t_m * (k_m ^ 2.0)))); else tmp = (l * l) * ((2.0 / t_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[k$95$m, 0.019], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(N[(2.0 / t$95$m), $MachinePrecision] / N[Power[N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $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 0.019:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{{k\_m}^{2} \cdot \left(t\_m \cdot {k\_m}^{2}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{\frac{2}{t\_m}}{{\left(k\_m \cdot \sin k\_m\right)}^{2}}\\
\end{array}
\end{array}
if k < 0.0189999999999999995Initial program 38.2%
Simplified44.0%
Taylor expanded in t around 0 71.8%
associate-*r/71.8%
associate-*r*71.8%
*-commutative71.8%
Simplified71.8%
Taylor expanded in k around 0 66.4%
Taylor expanded in k around 0 65.0%
if 0.0189999999999999995 < k Initial program 27.4%
Simplified33.5%
Taylor expanded in t around 0 72.0%
associate-*r/72.0%
associate-*r*72.0%
*-commutative72.0%
Simplified72.0%
Taylor expanded in k around 0 53.1%
div-inv53.1%
associate-*r*53.1%
pow-prod-down53.1%
Applied egg-rr53.1%
associate-*r/53.1%
metadata-eval53.1%
*-commutative53.1%
associate-/r*53.1%
*-commutative53.1%
Simplified53.1%
Final simplification61.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 (* (* l l) (/ (/ 2.0 t_m) (pow (* 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) {
return t_s * ((l * l) * ((2.0 / t_m) / pow((k_m * sin(k_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 * l) * ((2.0d0 / t_m) / ((k_m * sin(k_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 * ((l * l) * ((2.0 / t_m) / Math.pow((k_m * Math.sin(k_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 * ((l * l) * ((2.0 / t_m) / math.pow((k_m * math.sin(k_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(Float64(l * l) * Float64(Float64(2.0 / t_m) / (Float64(k_m * sin(k_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 * l) * ((2.0 / t_m) / ((k_m * sin(k_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[(N[(l * l), $MachinePrecision] * N[(N[(2.0 / t$95$m), $MachinePrecision] / N[Power[N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $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 \left(\left(\ell \cdot \ell\right) \cdot \frac{\frac{2}{t\_m}}{{\left(k\_m \cdot \sin k\_m\right)}^{2}}\right)
\end{array}
Initial program 35.4%
Simplified41.3%
Taylor expanded in t around 0 71.8%
associate-*r/71.8%
associate-*r*71.8%
*-commutative71.8%
Simplified71.8%
Taylor expanded in k around 0 63.0%
div-inv63.0%
associate-*r*59.9%
pow-prod-down59.9%
Applied egg-rr59.9%
associate-*r/59.9%
metadata-eval59.9%
*-commutative59.9%
associate-/r*59.9%
*-commutative59.9%
Simplified59.9%
Final simplification59.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 (* (* l l) (* 2.0 (/ (/ (cos k_m) 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 * ((cos(k_m) / 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 * ((cos(k_m) / 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 * ((Math.cos(k_m) / 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 * ((math.cos(k_m) / 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(Float64(cos(k_m) / 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 * ((cos(k_m) / 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[(N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision] / 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 \left(2 \cdot \frac{\frac{\cos k\_m}{t\_m}}{{k\_m}^{4}}\right)\right)
\end{array}
Initial program 35.4%
Simplified41.3%
Taylor expanded in t around 0 71.8%
associate-*r/71.8%
associate-*r*71.8%
*-commutative71.8%
Simplified71.8%
associate-/l*71.8%
associate-*r*68.7%
pow-prod-down68.8%
Applied egg-rr68.8%
*-commutative68.8%
associate-/r*68.8%
*-commutative68.8%
Simplified68.8%
Taylor expanded in k around 0 57.9%
Final simplification57.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 l) 5e+303)
(* (* l l) (/ (* 2.0 (pow k_m -4.0)) t_m))
(* (* 2.0 (/ (* l l) t_m)) 0.006944444444444444))))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) <= 5e+303) {
tmp = (l * l) * ((2.0 * pow(k_m, -4.0)) / t_m);
} else {
tmp = (2.0 * ((l * l) / t_m)) * 0.006944444444444444;
}
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) <= 5d+303) then
tmp = (l * l) * ((2.0d0 * (k_m ** (-4.0d0))) / t_m)
else
tmp = (2.0d0 * ((l * l) / t_m)) * 0.006944444444444444d0
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) <= 5e+303) {
tmp = (l * l) * ((2.0 * Math.pow(k_m, -4.0)) / t_m);
} else {
tmp = (2.0 * ((l * l) / t_m)) * 0.006944444444444444;
}
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) <= 5e+303: tmp = (l * l) * ((2.0 * math.pow(k_m, -4.0)) / t_m) else: tmp = (2.0 * ((l * l) / t_m)) * 0.006944444444444444 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) <= 5e+303) tmp = Float64(Float64(l * l) * Float64(Float64(2.0 * (k_m ^ -4.0)) / t_m)); else tmp = Float64(Float64(2.0 * Float64(Float64(l * l) / t_m)) * 0.006944444444444444); 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) <= 5e+303) tmp = (l * l) * ((2.0 * (k_m ^ -4.0)) / t_m); else tmp = (2.0 * ((l * l) / t_m)) * 0.006944444444444444; 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], 5e+303], N[(N[(l * l), $MachinePrecision] * N[(N[(2.0 * N[Power[k$95$m, -4.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(N[(l * l), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision] * 0.006944444444444444), $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 5 \cdot 10^{+303}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2 \cdot {k\_m}^{-4}}{t\_m}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \frac{\ell \cdot \ell}{t\_m}\right) \cdot 0.006944444444444444\\
\end{array}
\end{array}
if (*.f64 l l) < 4.9999999999999997e303Initial program 35.3%
Simplified43.5%
Taylor expanded in k around 0 60.0%
*-commutative60.0%
associate-/r*60.0%
Simplified60.0%
div-inv60.0%
pow-flip60.0%
metadata-eval60.0%
Applied egg-rr60.0%
associate-*l/60.0%
Simplified60.0%
if 4.9999999999999997e303 < (*.f64 l l) Initial program 35.6%
*-commutative35.6%
associate-/r*35.6%
Simplified35.6%
add-sqr-sqrt19.9%
pow219.9%
Applied egg-rr36.6%
Taylor expanded in k around 0 31.1%
associate-*r*31.1%
distribute-rgt-out33.9%
Simplified33.9%
Taylor expanded in k around inf 57.3%
*-commutative57.3%
*-commutative57.3%
unpow257.3%
rem-square-sqrt57.3%
associate-*r/57.3%
Simplified57.3%
pow257.3%
Applied egg-rr57.3%
Final simplification59.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
(*
t_s
(if (<= (* l l) 5e+303)
(* (* l l) (/ 2.0 (* t_m (pow k_m 4.0))))
(* (* 2.0 (/ (* l l) t_m)) 0.006944444444444444))))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) <= 5e+303) {
tmp = (l * l) * (2.0 / (t_m * pow(k_m, 4.0)));
} else {
tmp = (2.0 * ((l * l) / t_m)) * 0.006944444444444444;
}
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) <= 5d+303) then
tmp = (l * l) * (2.0d0 / (t_m * (k_m ** 4.0d0)))
else
tmp = (2.0d0 * ((l * l) / t_m)) * 0.006944444444444444d0
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) <= 5e+303) {
tmp = (l * l) * (2.0 / (t_m * Math.pow(k_m, 4.0)));
} else {
tmp = (2.0 * ((l * l) / t_m)) * 0.006944444444444444;
}
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) <= 5e+303: tmp = (l * l) * (2.0 / (t_m * math.pow(k_m, 4.0))) else: tmp = (2.0 * ((l * l) / t_m)) * 0.006944444444444444 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) <= 5e+303) tmp = Float64(Float64(l * l) * Float64(2.0 / Float64(t_m * (k_m ^ 4.0)))); else tmp = Float64(Float64(2.0 * Float64(Float64(l * l) / t_m)) * 0.006944444444444444); 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) <= 5e+303) tmp = (l * l) * (2.0 / (t_m * (k_m ^ 4.0))); else tmp = (2.0 * ((l * l) / t_m)) * 0.006944444444444444; 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], 5e+303], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(t$95$m * N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(N[(l * l), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision] * 0.006944444444444444), $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 5 \cdot 10^{+303}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{t\_m \cdot {k\_m}^{4}}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \frac{\ell \cdot \ell}{t\_m}\right) \cdot 0.006944444444444444\\
\end{array}
\end{array}
if (*.f64 l l) < 4.9999999999999997e303Initial program 35.3%
Simplified43.5%
Taylor expanded in k around 0 60.0%
if 4.9999999999999997e303 < (*.f64 l l) Initial program 35.6%
*-commutative35.6%
associate-/r*35.6%
Simplified35.6%
add-sqr-sqrt19.9%
pow219.9%
Applied egg-rr36.6%
Taylor expanded in k around 0 31.1%
associate-*r*31.1%
distribute-rgt-out33.9%
Simplified33.9%
Taylor expanded in k around inf 57.3%
*-commutative57.3%
*-commutative57.3%
unpow257.3%
rem-square-sqrt57.3%
associate-*r/57.3%
Simplified57.3%
pow257.3%
Applied egg-rr57.3%
Final simplification59.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 (* t_s (/ (* 0.013888888888888888 (pow l 2.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 * ((0.013888888888888888 * pow(l, 2.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 * ((0.013888888888888888d0 * (l ** 2.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 * ((0.013888888888888888 * Math.pow(l, 2.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 * ((0.013888888888888888 * math.pow(l, 2.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(0.013888888888888888 * (l ^ 2.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 * ((0.013888888888888888 * (l ^ 2.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[(0.013888888888888888 * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $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 \frac{0.013888888888888888 \cdot {\ell}^{2}}{t\_m}
\end{array}
Initial program 35.4%
*-commutative35.4%
associate-/r*35.5%
Simplified40.5%
add-sqr-sqrt27.9%
pow227.9%
Applied egg-rr27.7%
Taylor expanded in k around 0 29.6%
associate-*r*29.6%
distribute-rgt-out30.4%
Simplified30.4%
Taylor expanded in k around inf 36.8%
*-commutative36.8%
*-commutative36.8%
unpow236.8%
rem-square-sqrt36.8%
associate-*r/36.8%
Simplified36.8%
Taylor expanded in l around 0 36.8%
associate-*r/36.8%
Simplified36.8%
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 (/ (* l l) t_m)) 0.006944444444444444)))
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 * ((l * l) / t_m)) * 0.006944444444444444);
}
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 * l) / t_m)) * 0.006944444444444444d0)
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 * ((l * l) / t_m)) * 0.006944444444444444);
}
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 * ((l * l) / t_m)) * 0.006944444444444444)
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(2.0 * Float64(Float64(l * l) / t_m)) * 0.006944444444444444)) 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 * l) / t_m)) * 0.006944444444444444); 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[(2.0 * N[(N[(l * l), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision] * 0.006944444444444444), $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(2 \cdot \frac{\ell \cdot \ell}{t\_m}\right) \cdot 0.006944444444444444\right)
\end{array}
Initial program 35.4%
*-commutative35.4%
associate-/r*35.5%
Simplified40.5%
add-sqr-sqrt27.9%
pow227.9%
Applied egg-rr27.7%
Taylor expanded in k around 0 29.6%
associate-*r*29.6%
distribute-rgt-out30.4%
Simplified30.4%
Taylor expanded in k around inf 36.8%
*-commutative36.8%
*-commutative36.8%
unpow236.8%
rem-square-sqrt36.8%
associate-*r/36.8%
Simplified36.8%
pow236.8%
Applied egg-rr36.8%
herbie shell --seed 2024139
(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))))