
(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 19 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 (pow (sin k_m) 2.0)))
(*
t_s
(if (<= k_m 6.2e-133)
(/
2.0
(*
(* (pow (* (/ (pow t_m 1.5) l) (sqrt (sin k_m))) 2.0) (tan k_m))
(+ 1.0 (+ 1.0 (* (/ k_m t_m) (/ k_m t_m))))))
(/
2.0
(*
t_m
(fma
2.0
(/ (* (pow (/ t_m l) 2.0) t_2) (cos k_m))
(* (/ t_2 (pow l 2.0)) (/ (pow 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 t_2 = pow(sin(k_m), 2.0);
double tmp;
if (k_m <= 6.2e-133) {
tmp = 2.0 / ((pow(((pow(t_m, 1.5) / l) * sqrt(sin(k_m))), 2.0) * tan(k_m)) * (1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))));
} else {
tmp = 2.0 / (t_m * fma(2.0, ((pow((t_m / l), 2.0) * t_2) / cos(k_m)), ((t_2 / pow(l, 2.0)) * (pow(k_m, 2.0) / 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) t_2 = sin(k_m) ^ 2.0 tmp = 0.0 if (k_m <= 6.2e-133) tmp = Float64(2.0 / Float64(Float64((Float64(Float64((t_m ^ 1.5) / l) * sqrt(sin(k_m))) ^ 2.0) * tan(k_m)) * Float64(1.0 + Float64(1.0 + Float64(Float64(k_m / t_m) * Float64(k_m / t_m)))))); else tmp = Float64(2.0 / Float64(t_m * fma(2.0, Float64(Float64((Float64(t_m / l) ^ 2.0) * t_2) / cos(k_m)), Float64(Float64(t_2 / (l ^ 2.0)) * Float64((k_m ^ 2.0) / cos(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[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 6.2e-133], N[(2.0 / N[(N[(N[Power[N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision] * N[Sqrt[N[Sin[k$95$m], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(1.0 + N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(t$95$m * N[(2.0 * N[(N[(N[Power[N[(t$95$m / l), $MachinePrecision], 2.0], $MachinePrecision] * t$95$2), $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] + N[(N[(t$95$2 / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Power[k$95$m, 2.0], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $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 := {\sin k\_m}^{2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 6.2 \cdot 10^{-133}:\\
\;\;\;\;\frac{2}{\left({\left(\frac{{t\_m}^{1.5}}{\ell} \cdot \sqrt{\sin k\_m}\right)}^{2} \cdot \tan k\_m\right) \cdot \left(1 + \left(1 + \frac{k\_m}{t\_m} \cdot \frac{k\_m}{t\_m}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t\_m \cdot \mathsf{fma}\left(2, \frac{{\left(\frac{t\_m}{\ell}\right)}^{2} \cdot t\_2}{\cos k\_m}, \frac{t\_2}{{\ell}^{2}} \cdot \frac{{k\_m}^{2}}{\cos k\_m}\right)}\\
\end{array}
\end{array}
\end{array}
if k < 6.20000000000000032e-133Initial program 61.2%
add-sqr-sqrt32.4%
pow232.4%
associate-/r*37.1%
*-commutative37.1%
sqrt-prod17.4%
associate-/r*15.0%
sqrt-div15.0%
sqrt-pow116.2%
metadata-eval16.2%
sqrt-prod10.7%
add-sqr-sqrt19.2%
Applied egg-rr19.2%
*-commutative19.2%
Simplified19.2%
unpow219.2%
Applied egg-rr19.2%
if 6.20000000000000032e-133 < k Initial program 57.4%
add-sqr-sqrt31.1%
pow231.1%
associate-/r*37.1%
*-commutative37.1%
sqrt-prod23.4%
associate-/r*19.6%
sqrt-div19.6%
sqrt-pow121.8%
metadata-eval21.8%
sqrt-prod17.0%
add-sqr-sqrt26.2%
Applied egg-rr26.2%
*-commutative26.2%
Simplified26.2%
div-inv26.2%
associate-*l*26.2%
*-commutative26.2%
unpow-prod-down26.1%
pow226.1%
add-sqr-sqrt33.1%
add-sqr-sqrt33.1%
pow233.1%
Applied egg-rr33.1%
associate-*r/33.1%
metadata-eval33.1%
associate-*l*32.0%
unpow232.0%
hypot-undefine32.0%
hypot-undefine32.0%
rem-square-sqrt32.0%
metadata-eval32.0%
hypot-undefine32.0%
metadata-eval32.0%
unpow232.0%
+-commutative32.0%
Simplified32.0%
Taylor expanded in t around 0 73.1%
fma-define73.1%
times-frac73.1%
associate-*r/73.1%
unpow273.1%
unpow273.1%
times-frac80.9%
unpow180.9%
pow-plus80.9%
metadata-eval80.9%
*-commutative80.9%
times-frac82.1%
Simplified82.1%
Final simplification40.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
(let* ((t_2 (/ (pow t_m 1.5) l)))
(*
t_s
(if (<= k_m 2.4e-111)
(/
2.0
(*
(* (pow (* t_2 (sqrt (sin k_m))) 2.0) (tan k_m))
(+ 1.0 (+ 1.0 (* (/ k_m t_m) (/ k_m t_m))))))
(if (<= k_m 1.95e-5)
(pow
(/
(sqrt 2.0)
(*
(hypot 1.0 (hypot 1.0 (/ k_m t_m)))
(* t_2 (sqrt (* (sin k_m) (tan k_m))))))
2.0)
(/
2.0
(*
(/ (* t_m (pow k_m 2.0)) (pow l 2.0))
(/ (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 t_2 = pow(t_m, 1.5) / l;
double tmp;
if (k_m <= 2.4e-111) {
tmp = 2.0 / ((pow((t_2 * sqrt(sin(k_m))), 2.0) * tan(k_m)) * (1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))));
} else if (k_m <= 1.95e-5) {
tmp = pow((sqrt(2.0) / (hypot(1.0, hypot(1.0, (k_m / t_m))) * (t_2 * sqrt((sin(k_m) * tan(k_m)))))), 2.0);
} else {
tmp = 2.0 / (((t_m * pow(k_m, 2.0)) / pow(l, 2.0)) * (pow(sin(k_m), 2.0) / cos(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(t_m, 1.5) / l;
double tmp;
if (k_m <= 2.4e-111) {
tmp = 2.0 / ((Math.pow((t_2 * Math.sqrt(Math.sin(k_m))), 2.0) * Math.tan(k_m)) * (1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))));
} else if (k_m <= 1.95e-5) {
tmp = Math.pow((Math.sqrt(2.0) / (Math.hypot(1.0, Math.hypot(1.0, (k_m / t_m))) * (t_2 * Math.sqrt((Math.sin(k_m) * Math.tan(k_m)))))), 2.0);
} else {
tmp = 2.0 / (((t_m * Math.pow(k_m, 2.0)) / Math.pow(l, 2.0)) * (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): t_2 = math.pow(t_m, 1.5) / l tmp = 0 if k_m <= 2.4e-111: tmp = 2.0 / ((math.pow((t_2 * math.sqrt(math.sin(k_m))), 2.0) * math.tan(k_m)) * (1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m))))) elif k_m <= 1.95e-5: tmp = math.pow((math.sqrt(2.0) / (math.hypot(1.0, math.hypot(1.0, (k_m / t_m))) * (t_2 * math.sqrt((math.sin(k_m) * math.tan(k_m)))))), 2.0) else: tmp = 2.0 / (((t_m * math.pow(k_m, 2.0)) / math.pow(l, 2.0)) * (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) t_2 = Float64((t_m ^ 1.5) / l) tmp = 0.0 if (k_m <= 2.4e-111) tmp = Float64(2.0 / Float64(Float64((Float64(t_2 * sqrt(sin(k_m))) ^ 2.0) * tan(k_m)) * Float64(1.0 + Float64(1.0 + Float64(Float64(k_m / t_m) * Float64(k_m / t_m)))))); elseif (k_m <= 1.95e-5) tmp = Float64(sqrt(2.0) / Float64(hypot(1.0, hypot(1.0, Float64(k_m / t_m))) * Float64(t_2 * sqrt(Float64(sin(k_m) * tan(k_m)))))) ^ 2.0; else tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k_m ^ 2.0)) / (l ^ 2.0)) * Float64((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) t_2 = (t_m ^ 1.5) / l; tmp = 0.0; if (k_m <= 2.4e-111) tmp = 2.0 / ((((t_2 * sqrt(sin(k_m))) ^ 2.0) * tan(k_m)) * (1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m))))); elseif (k_m <= 1.95e-5) tmp = (sqrt(2.0) / (hypot(1.0, hypot(1.0, (k_m / t_m))) * (t_2 * sqrt((sin(k_m) * tan(k_m)))))) ^ 2.0; else tmp = 2.0 / (((t_m * (k_m ^ 2.0)) / (l ^ 2.0)) * ((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_] := Block[{t$95$2 = N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 2.4e-111], N[(2.0 / N[(N[(N[Power[N[(t$95$2 * N[Sqrt[N[Sin[k$95$m], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(1.0 + N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 1.95e-5], N[Power[N[(N[Sqrt[2.0], $MachinePrecision] / N[(N[Sqrt[1.0 ^ 2 + N[Sqrt[1.0 ^ 2 + N[(k$95$m / t$95$m), $MachinePrecision] ^ 2], $MachinePrecision] ^ 2], $MachinePrecision] * N[(t$95$2 * N[Sqrt[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $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)
\\
\begin{array}{l}
t_2 := \frac{{t\_m}^{1.5}}{\ell}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 2.4 \cdot 10^{-111}:\\
\;\;\;\;\frac{2}{\left({\left(t\_2 \cdot \sqrt{\sin k\_m}\right)}^{2} \cdot \tan k\_m\right) \cdot \left(1 + \left(1 + \frac{k\_m}{t\_m} \cdot \frac{k\_m}{t\_m}\right)\right)}\\
\mathbf{elif}\;k\_m \leq 1.95 \cdot 10^{-5}:\\
\;\;\;\;{\left(\frac{\sqrt{2}}{\mathsf{hypot}\left(1, \mathsf{hypot}\left(1, \frac{k\_m}{t\_m}\right)\right) \cdot \left(t\_2 \cdot \sqrt{\sin k\_m \cdot \tan k\_m}\right)}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k\_m}^{2}}{{\ell}^{2}} \cdot \frac{{\sin k\_m}^{2}}{\cos k\_m}}\\
\end{array}
\end{array}
\end{array}
if k < 2.4000000000000001e-111Initial program 61.5%
add-sqr-sqrt32.8%
pow232.8%
associate-/r*37.4%
*-commutative37.4%
sqrt-prod18.2%
associate-/r*15.9%
sqrt-div15.8%
sqrt-pow117.0%
metadata-eval17.0%
sqrt-prod11.1%
add-sqr-sqrt19.9%
Applied egg-rr19.9%
*-commutative19.9%
Simplified19.9%
unpow219.9%
Applied egg-rr19.9%
if 2.4000000000000001e-111 < k < 1.95e-5Initial program 63.6%
add-sqr-sqrt42.9%
Applied egg-rr37.3%
unpow237.3%
Simplified37.3%
if 1.95e-5 < k Initial program 53.8%
Simplified58.8%
add-cube-cbrt58.7%
*-un-lft-identity58.7%
times-frac58.7%
pow258.7%
cbrt-div58.7%
rem-cbrt-cube58.7%
cbrt-div58.7%
rem-cbrt-cube65.2%
Applied egg-rr65.2%
add-cube-cbrt65.1%
pow365.1%
Applied egg-rr73.6%
*-commutative73.6%
cube-prod66.8%
rem-cube-cbrt66.8%
associate-*l*66.8%
Simplified66.8%
Taylor expanded in k around inf 76.8%
associate-*r*76.8%
times-frac76.9%
Simplified76.9%
Final simplification34.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 t_m 1.5) l)))
(*
t_s
(if (<= k_m 1.3e-101)
(/
2.0
(*
(* (pow (* t_2 (sqrt (sin k_m))) 2.0) (tan k_m))
(+ 1.0 (+ 1.0 (* (/ k_m t_m) (/ k_m t_m))))))
(if (<= k_m 1.65e-5)
(/
2.0
(pow
(*
(hypot 1.0 (hypot 1.0 (/ k_m t_m)))
(* t_2 (sqrt (* (sin k_m) (tan k_m)))))
2.0))
(/
2.0
(*
(/ (* t_m (pow k_m 2.0)) (pow l 2.0))
(/ (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 t_2 = pow(t_m, 1.5) / l;
double tmp;
if (k_m <= 1.3e-101) {
tmp = 2.0 / ((pow((t_2 * sqrt(sin(k_m))), 2.0) * tan(k_m)) * (1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))));
} else if (k_m <= 1.65e-5) {
tmp = 2.0 / pow((hypot(1.0, hypot(1.0, (k_m / t_m))) * (t_2 * sqrt((sin(k_m) * tan(k_m))))), 2.0);
} else {
tmp = 2.0 / (((t_m * pow(k_m, 2.0)) / pow(l, 2.0)) * (pow(sin(k_m), 2.0) / cos(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(t_m, 1.5) / l;
double tmp;
if (k_m <= 1.3e-101) {
tmp = 2.0 / ((Math.pow((t_2 * Math.sqrt(Math.sin(k_m))), 2.0) * Math.tan(k_m)) * (1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))));
} else if (k_m <= 1.65e-5) {
tmp = 2.0 / Math.pow((Math.hypot(1.0, Math.hypot(1.0, (k_m / t_m))) * (t_2 * Math.sqrt((Math.sin(k_m) * Math.tan(k_m))))), 2.0);
} else {
tmp = 2.0 / (((t_m * Math.pow(k_m, 2.0)) / Math.pow(l, 2.0)) * (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): t_2 = math.pow(t_m, 1.5) / l tmp = 0 if k_m <= 1.3e-101: tmp = 2.0 / ((math.pow((t_2 * math.sqrt(math.sin(k_m))), 2.0) * math.tan(k_m)) * (1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m))))) elif k_m <= 1.65e-5: tmp = 2.0 / math.pow((math.hypot(1.0, math.hypot(1.0, (k_m / t_m))) * (t_2 * math.sqrt((math.sin(k_m) * math.tan(k_m))))), 2.0) else: tmp = 2.0 / (((t_m * math.pow(k_m, 2.0)) / math.pow(l, 2.0)) * (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) t_2 = Float64((t_m ^ 1.5) / l) tmp = 0.0 if (k_m <= 1.3e-101) tmp = Float64(2.0 / Float64(Float64((Float64(t_2 * sqrt(sin(k_m))) ^ 2.0) * tan(k_m)) * Float64(1.0 + Float64(1.0 + Float64(Float64(k_m / t_m) * Float64(k_m / t_m)))))); elseif (k_m <= 1.65e-5) tmp = Float64(2.0 / (Float64(hypot(1.0, hypot(1.0, Float64(k_m / t_m))) * Float64(t_2 * sqrt(Float64(sin(k_m) * tan(k_m))))) ^ 2.0)); else tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k_m ^ 2.0)) / (l ^ 2.0)) * Float64((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) t_2 = (t_m ^ 1.5) / l; tmp = 0.0; if (k_m <= 1.3e-101) tmp = 2.0 / ((((t_2 * sqrt(sin(k_m))) ^ 2.0) * tan(k_m)) * (1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m))))); elseif (k_m <= 1.65e-5) tmp = 2.0 / ((hypot(1.0, hypot(1.0, (k_m / t_m))) * (t_2 * sqrt((sin(k_m) * tan(k_m))))) ^ 2.0); else tmp = 2.0 / (((t_m * (k_m ^ 2.0)) / (l ^ 2.0)) * ((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_] := Block[{t$95$2 = N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 1.3e-101], N[(2.0 / N[(N[(N[Power[N[(t$95$2 * N[Sqrt[N[Sin[k$95$m], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(1.0 + N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 1.65e-5], N[(2.0 / N[Power[N[(N[Sqrt[1.0 ^ 2 + N[Sqrt[1.0 ^ 2 + N[(k$95$m / t$95$m), $MachinePrecision] ^ 2], $MachinePrecision] ^ 2], $MachinePrecision] * N[(t$95$2 * N[Sqrt[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $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)
\\
\begin{array}{l}
t_2 := \frac{{t\_m}^{1.5}}{\ell}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 1.3 \cdot 10^{-101}:\\
\;\;\;\;\frac{2}{\left({\left(t\_2 \cdot \sqrt{\sin k\_m}\right)}^{2} \cdot \tan k\_m\right) \cdot \left(1 + \left(1 + \frac{k\_m}{t\_m} \cdot \frac{k\_m}{t\_m}\right)\right)}\\
\mathbf{elif}\;k\_m \leq 1.65 \cdot 10^{-5}:\\
\;\;\;\;\frac{2}{{\left(\mathsf{hypot}\left(1, \mathsf{hypot}\left(1, \frac{k\_m}{t\_m}\right)\right) \cdot \left(t\_2 \cdot \sqrt{\sin k\_m \cdot \tan k\_m}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k\_m}^{2}}{{\ell}^{2}} \cdot \frac{{\sin k\_m}^{2}}{\cos k\_m}}\\
\end{array}
\end{array}
\end{array}
if k < 1.3000000000000001e-101Initial program 61.5%
add-sqr-sqrt32.8%
pow232.8%
associate-/r*37.4%
*-commutative37.4%
sqrt-prod18.2%
associate-/r*15.9%
sqrt-div15.8%
sqrt-pow117.0%
metadata-eval17.0%
sqrt-prod11.1%
add-sqr-sqrt19.9%
Applied egg-rr19.9%
*-commutative19.9%
Simplified19.9%
unpow219.9%
Applied egg-rr19.9%
if 1.3000000000000001e-101 < k < 1.6500000000000001e-5Initial program 63.6%
div-inv63.6%
add-sqr-sqrt21.8%
pow221.8%
Applied egg-rr37.3%
associate-*r/37.3%
metadata-eval37.3%
Simplified37.3%
if 1.6500000000000001e-5 < k Initial program 53.8%
Simplified58.8%
add-cube-cbrt58.7%
*-un-lft-identity58.7%
times-frac58.7%
pow258.7%
cbrt-div58.7%
rem-cbrt-cube58.7%
cbrt-div58.7%
rem-cbrt-cube65.2%
Applied egg-rr65.2%
add-cube-cbrt65.1%
pow365.1%
Applied egg-rr73.6%
*-commutative73.6%
cube-prod66.8%
rem-cube-cbrt66.8%
associate-*l*66.8%
Simplified66.8%
Taylor expanded in k around inf 76.8%
associate-*r*76.8%
times-frac76.9%
Simplified76.9%
Final simplification34.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 (<= t_m 4.7e-98)
(/
2.0
(*
(/ (* t_m (pow k_m 2.0)) (pow l 2.0))
(/ (pow (sin k_m) 2.0) (cos k_m))))
(/
(/ (/ 2.0 (sin k_m)) (pow (/ (pow t_m 1.5) l) 2.0))
(* (tan k_m) (+ 2.0 (pow (/ k_m 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 (t_m <= 4.7e-98) {
tmp = 2.0 / (((t_m * pow(k_m, 2.0)) / pow(l, 2.0)) * (pow(sin(k_m), 2.0) / cos(k_m)));
} else {
tmp = ((2.0 / sin(k_m)) / pow((pow(t_m, 1.5) / l), 2.0)) / (tan(k_m) * (2.0 + pow((k_m / 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 (t_m <= 4.7d-98) then
tmp = 2.0d0 / (((t_m * (k_m ** 2.0d0)) / (l ** 2.0d0)) * ((sin(k_m) ** 2.0d0) / cos(k_m)))
else
tmp = ((2.0d0 / sin(k_m)) / (((t_m ** 1.5d0) / l) ** 2.0d0)) / (tan(k_m) * (2.0d0 + ((k_m / 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 (t_m <= 4.7e-98) {
tmp = 2.0 / (((t_m * Math.pow(k_m, 2.0)) / Math.pow(l, 2.0)) * (Math.pow(Math.sin(k_m), 2.0) / Math.cos(k_m)));
} else {
tmp = ((2.0 / Math.sin(k_m)) / Math.pow((Math.pow(t_m, 1.5) / l), 2.0)) / (Math.tan(k_m) * (2.0 + Math.pow((k_m / 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 t_m <= 4.7e-98: tmp = 2.0 / (((t_m * math.pow(k_m, 2.0)) / math.pow(l, 2.0)) * (math.pow(math.sin(k_m), 2.0) / math.cos(k_m))) else: tmp = ((2.0 / math.sin(k_m)) / math.pow((math.pow(t_m, 1.5) / l), 2.0)) / (math.tan(k_m) * (2.0 + math.pow((k_m / 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 (t_m <= 4.7e-98) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k_m ^ 2.0)) / (l ^ 2.0)) * Float64((sin(k_m) ^ 2.0) / cos(k_m)))); else tmp = Float64(Float64(Float64(2.0 / sin(k_m)) / (Float64((t_m ^ 1.5) / l) ^ 2.0)) / Float64(tan(k_m) * Float64(2.0 + (Float64(k_m / 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 (t_m <= 4.7e-98) tmp = 2.0 / (((t_m * (k_m ^ 2.0)) / (l ^ 2.0)) * ((sin(k_m) ^ 2.0) / cos(k_m))); else tmp = ((2.0 / sin(k_m)) / (((t_m ^ 1.5) / l) ^ 2.0)) / (tan(k_m) * (2.0 + ((k_m / 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[t$95$m, 4.7e-98], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision] / N[Power[N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[(N[Tan[k$95$m], $MachinePrecision] * N[(2.0 + N[Power[N[(k$95$m / t$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 \begin{array}{l}
\mathbf{if}\;t\_m \leq 4.7 \cdot 10^{-98}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k\_m}^{2}}{{\ell}^{2}} \cdot \frac{{\sin k\_m}^{2}}{\cos k\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\frac{2}{\sin k\_m}}{{\left(\frac{{t\_m}^{1.5}}{\ell}\right)}^{2}}}{\tan k\_m \cdot \left(2 + {\left(\frac{k\_m}{t\_m}\right)}^{2}\right)}\\
\end{array}
\end{array}
if t < 4.70000000000000005e-98Initial program 54.1%
Simplified56.2%
add-cube-cbrt56.1%
*-un-lft-identity56.1%
times-frac56.1%
pow256.1%
cbrt-div56.2%
rem-cbrt-cube56.2%
cbrt-div56.1%
rem-cbrt-cube65.3%
Applied egg-rr65.3%
add-cube-cbrt65.2%
pow365.2%
Applied egg-rr72.5%
*-commutative72.5%
cube-prod66.6%
rem-cube-cbrt66.6%
associate-*l*66.9%
Simplified66.9%
Taylor expanded in k around inf 65.9%
associate-*r*65.9%
times-frac67.8%
Simplified67.8%
if 4.70000000000000005e-98 < t Initial program 70.2%
add-sqr-sqrt39.5%
pow239.5%
associate-/r*44.9%
*-commutative44.9%
sqrt-prod45.0%
associate-/r*39.5%
sqrt-div39.5%
sqrt-pow139.5%
metadata-eval39.5%
sqrt-prod26.1%
add-sqr-sqrt46.0%
Applied egg-rr46.0%
*-commutative46.0%
Simplified46.0%
div-inv46.0%
associate-*l*46.0%
*-commutative46.0%
unpow-prod-down44.9%
pow244.9%
add-sqr-sqrt80.0%
add-sqr-sqrt80.0%
pow280.0%
Applied egg-rr80.0%
associate-*r/80.0%
metadata-eval80.0%
associate-*l*79.9%
unpow279.9%
hypot-undefine79.9%
hypot-undefine79.9%
rem-square-sqrt80.0%
metadata-eval80.0%
hypot-undefine80.0%
metadata-eval80.0%
unpow280.0%
+-commutative80.0%
Simplified80.0%
div-inv80.0%
associate-*r*80.0%
pow280.0%
associate-+r+80.0%
metadata-eval80.0%
pow280.0%
Applied egg-rr80.0%
associate-*r/80.0%
metadata-eval80.0%
associate-/r*80.0%
associate-/r*80.0%
*-commutative80.0%
Simplified80.0%
Final simplification72.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 (<= t_m 6.8e-100)
(/
2.0
(*
(/ (* t_m (pow k_m 2.0)) (pow l 2.0))
(/ (pow (sin k_m) 2.0) (cos k_m))))
(/
2.0
(*
(sin k_m)
(*
(pow (/ (pow t_m 1.5) l) 2.0)
(* (tan k_m) (+ 2.0 (pow (/ k_m 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 (t_m <= 6.8e-100) {
tmp = 2.0 / (((t_m * pow(k_m, 2.0)) / pow(l, 2.0)) * (pow(sin(k_m), 2.0) / cos(k_m)));
} else {
tmp = 2.0 / (sin(k_m) * (pow((pow(t_m, 1.5) / l), 2.0) * (tan(k_m) * (2.0 + pow((k_m / 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 (t_m <= 6.8d-100) then
tmp = 2.0d0 / (((t_m * (k_m ** 2.0d0)) / (l ** 2.0d0)) * ((sin(k_m) ** 2.0d0) / cos(k_m)))
else
tmp = 2.0d0 / (sin(k_m) * ((((t_m ** 1.5d0) / l) ** 2.0d0) * (tan(k_m) * (2.0d0 + ((k_m / 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 (t_m <= 6.8e-100) {
tmp = 2.0 / (((t_m * Math.pow(k_m, 2.0)) / Math.pow(l, 2.0)) * (Math.pow(Math.sin(k_m), 2.0) / Math.cos(k_m)));
} else {
tmp = 2.0 / (Math.sin(k_m) * (Math.pow((Math.pow(t_m, 1.5) / l), 2.0) * (Math.tan(k_m) * (2.0 + Math.pow((k_m / 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 t_m <= 6.8e-100: tmp = 2.0 / (((t_m * math.pow(k_m, 2.0)) / math.pow(l, 2.0)) * (math.pow(math.sin(k_m), 2.0) / math.cos(k_m))) else: tmp = 2.0 / (math.sin(k_m) * (math.pow((math.pow(t_m, 1.5) / l), 2.0) * (math.tan(k_m) * (2.0 + math.pow((k_m / 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 (t_m <= 6.8e-100) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k_m ^ 2.0)) / (l ^ 2.0)) * Float64((sin(k_m) ^ 2.0) / cos(k_m)))); else tmp = Float64(2.0 / Float64(sin(k_m) * Float64((Float64((t_m ^ 1.5) / l) ^ 2.0) * Float64(tan(k_m) * Float64(2.0 + (Float64(k_m / 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 (t_m <= 6.8e-100) tmp = 2.0 / (((t_m * (k_m ^ 2.0)) / (l ^ 2.0)) * ((sin(k_m) ^ 2.0) / cos(k_m))); else tmp = 2.0 / (sin(k_m) * ((((t_m ^ 1.5) / l) ^ 2.0) * (tan(k_m) * (2.0 + ((k_m / 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[t$95$m, 6.8e-100], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Power[N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision], 2.0], $MachinePrecision] * N[(N[Tan[k$95$m], $MachinePrecision] * N[(2.0 + N[Power[N[(k$95$m / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $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}\;t\_m \leq 6.8 \cdot 10^{-100}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k\_m}^{2}}{{\ell}^{2}} \cdot \frac{{\sin k\_m}^{2}}{\cos k\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\sin k\_m \cdot \left({\left(\frac{{t\_m}^{1.5}}{\ell}\right)}^{2} \cdot \left(\tan k\_m \cdot \left(2 + {\left(\frac{k\_m}{t\_m}\right)}^{2}\right)\right)\right)}\\
\end{array}
\end{array}
if t < 6.79999999999999953e-100Initial program 54.1%
Simplified56.2%
add-cube-cbrt56.1%
*-un-lft-identity56.1%
times-frac56.1%
pow256.1%
cbrt-div56.2%
rem-cbrt-cube56.2%
cbrt-div56.1%
rem-cbrt-cube65.3%
Applied egg-rr65.3%
add-cube-cbrt65.2%
pow365.2%
Applied egg-rr72.5%
*-commutative72.5%
cube-prod66.6%
rem-cube-cbrt66.6%
associate-*l*66.9%
Simplified66.9%
Taylor expanded in k around inf 65.9%
associate-*r*65.9%
times-frac67.8%
Simplified67.8%
if 6.79999999999999953e-100 < t Initial program 70.2%
add-sqr-sqrt39.5%
pow239.5%
associate-/r*44.9%
*-commutative44.9%
sqrt-prod45.0%
associate-/r*39.5%
sqrt-div39.5%
sqrt-pow139.5%
metadata-eval39.5%
sqrt-prod26.1%
add-sqr-sqrt46.0%
Applied egg-rr46.0%
*-commutative46.0%
Simplified46.0%
div-inv46.0%
associate-*l*46.0%
*-commutative46.0%
unpow-prod-down44.9%
pow244.9%
add-sqr-sqrt80.0%
add-sqr-sqrt80.0%
pow280.0%
Applied egg-rr80.0%
associate-*r/80.0%
metadata-eval80.0%
associate-*l*79.9%
unpow279.9%
hypot-undefine79.9%
hypot-undefine79.9%
rem-square-sqrt80.0%
metadata-eval80.0%
hypot-undefine80.0%
metadata-eval80.0%
unpow280.0%
+-commutative80.0%
Simplified80.0%
associate-*r*80.0%
add-sqr-sqrt44.9%
pow244.9%
unpow-prod-down46.0%
*-commutative46.0%
+-commutative46.0%
pow246.0%
associate-*l*46.0%
distribute-lft-in46.0%
Applied egg-rr80.0%
distribute-lft-out80.0%
*-commutative80.0%
+-commutative80.0%
unpow280.0%
times-frac69.0%
unpow269.0%
unpow269.0%
associate-+l+69.0%
metadata-eval69.0%
+-commutative69.0%
unpow269.0%
unpow269.0%
times-frac80.0%
unpow280.0%
associate-*r*80.0%
associate-*l*80.0%
Simplified80.0%
Final simplification72.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 (<= t_m 3.2e-99)
(/
2.0
(*
(/ (* t_m (pow k_m 2.0)) (pow l 2.0))
(/ (pow (sin k_m) 2.0) (cos k_m))))
(/
2.0
(*
(+ 1.0 (+ 1.0 (* (/ k_m t_m) (/ k_m t_m))))
(* (tan k_m) (* (sin k_m) (pow (/ (pow t_m 1.5) 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 (t_m <= 3.2e-99) {
tmp = 2.0 / (((t_m * pow(k_m, 2.0)) / pow(l, 2.0)) * (pow(sin(k_m), 2.0) / cos(k_m)));
} else {
tmp = 2.0 / ((1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))) * (tan(k_m) * (sin(k_m) * pow((pow(t_m, 1.5) / 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 (t_m <= 3.2d-99) then
tmp = 2.0d0 / (((t_m * (k_m ** 2.0d0)) / (l ** 2.0d0)) * ((sin(k_m) ** 2.0d0) / cos(k_m)))
else
tmp = 2.0d0 / ((1.0d0 + (1.0d0 + ((k_m / t_m) * (k_m / t_m)))) * (tan(k_m) * (sin(k_m) * (((t_m ** 1.5d0) / 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 (t_m <= 3.2e-99) {
tmp = 2.0 / (((t_m * Math.pow(k_m, 2.0)) / Math.pow(l, 2.0)) * (Math.pow(Math.sin(k_m), 2.0) / Math.cos(k_m)));
} else {
tmp = 2.0 / ((1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))) * (Math.tan(k_m) * (Math.sin(k_m) * Math.pow((Math.pow(t_m, 1.5) / 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 t_m <= 3.2e-99: tmp = 2.0 / (((t_m * math.pow(k_m, 2.0)) / math.pow(l, 2.0)) * (math.pow(math.sin(k_m), 2.0) / math.cos(k_m))) else: tmp = 2.0 / ((1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))) * (math.tan(k_m) * (math.sin(k_m) * math.pow((math.pow(t_m, 1.5) / 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 (t_m <= 3.2e-99) tmp = Float64(2.0 / Float64(Float64(Float64(t_m * (k_m ^ 2.0)) / (l ^ 2.0)) * Float64((sin(k_m) ^ 2.0) / cos(k_m)))); else tmp = Float64(2.0 / Float64(Float64(1.0 + Float64(1.0 + Float64(Float64(k_m / t_m) * Float64(k_m / t_m)))) * Float64(tan(k_m) * Float64(sin(k_m) * (Float64((t_m ^ 1.5) / 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 (t_m <= 3.2e-99) tmp = 2.0 / (((t_m * (k_m ^ 2.0)) / (l ^ 2.0)) * ((sin(k_m) ^ 2.0) / cos(k_m))); else tmp = 2.0 / ((1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))) * (tan(k_m) * (sin(k_m) * (((t_m ^ 1.5) / 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[t$95$m, 3.2e-99], N[(2.0 / N[(N[(N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(1.0 + N[(1.0 + N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Tan[k$95$m], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Power[N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $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}\;t\_m \leq 3.2 \cdot 10^{-99}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k\_m}^{2}}{{\ell}^{2}} \cdot \frac{{\sin k\_m}^{2}}{\cos k\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(1 + \left(1 + \frac{k\_m}{t\_m} \cdot \frac{k\_m}{t\_m}\right)\right) \cdot \left(\tan k\_m \cdot \left(\sin k\_m \cdot {\left(\frac{{t\_m}^{1.5}}{\ell}\right)}^{2}\right)\right)}\\
\end{array}
\end{array}
if t < 3.2000000000000001e-99Initial program 54.1%
Simplified56.2%
add-cube-cbrt56.1%
*-un-lft-identity56.1%
times-frac56.1%
pow256.1%
cbrt-div56.2%
rem-cbrt-cube56.2%
cbrt-div56.1%
rem-cbrt-cube65.3%
Applied egg-rr65.3%
add-cube-cbrt65.2%
pow365.2%
Applied egg-rr72.5%
*-commutative72.5%
cube-prod66.6%
rem-cube-cbrt66.6%
associate-*l*66.9%
Simplified66.9%
Taylor expanded in k around inf 65.9%
associate-*r*65.9%
times-frac67.8%
Simplified67.8%
if 3.2000000000000001e-99 < t Initial program 70.2%
add-sqr-sqrt39.5%
pow239.5%
associate-/r*44.9%
*-commutative44.9%
sqrt-prod45.0%
associate-/r*39.5%
sqrt-div39.5%
sqrt-pow139.5%
metadata-eval39.5%
sqrt-prod26.1%
add-sqr-sqrt46.0%
Applied egg-rr46.0%
*-commutative46.0%
Simplified46.0%
unpow246.0%
Applied egg-rr46.0%
unpow-prod-down44.9%
pow244.9%
add-sqr-sqrt80.0%
Applied egg-rr80.0%
Final simplification72.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 (<= t_m 2.6e-194)
(*
2.0
(/
(* (cos k_m) (* l l))
(* (pow k_m 2.0) (* t_m (- 0.5 (/ (cos (* k_m 2.0)) 2.0))))))
(/
2.0
(*
(+ 1.0 (+ 1.0 (* (/ k_m t_m) (/ k_m t_m))))
(* (tan k_m) (* (sin k_m) (pow (/ (pow t_m 1.5) 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 (t_m <= 2.6e-194) {
tmp = 2.0 * ((cos(k_m) * (l * l)) / (pow(k_m, 2.0) * (t_m * (0.5 - (cos((k_m * 2.0)) / 2.0)))));
} else {
tmp = 2.0 / ((1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))) * (tan(k_m) * (sin(k_m) * pow((pow(t_m, 1.5) / 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 (t_m <= 2.6d-194) then
tmp = 2.0d0 * ((cos(k_m) * (l * l)) / ((k_m ** 2.0d0) * (t_m * (0.5d0 - (cos((k_m * 2.0d0)) / 2.0d0)))))
else
tmp = 2.0d0 / ((1.0d0 + (1.0d0 + ((k_m / t_m) * (k_m / t_m)))) * (tan(k_m) * (sin(k_m) * (((t_m ** 1.5d0) / 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 (t_m <= 2.6e-194) {
tmp = 2.0 * ((Math.cos(k_m) * (l * l)) / (Math.pow(k_m, 2.0) * (t_m * (0.5 - (Math.cos((k_m * 2.0)) / 2.0)))));
} else {
tmp = 2.0 / ((1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))) * (Math.tan(k_m) * (Math.sin(k_m) * Math.pow((Math.pow(t_m, 1.5) / 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 t_m <= 2.6e-194: tmp = 2.0 * ((math.cos(k_m) * (l * l)) / (math.pow(k_m, 2.0) * (t_m * (0.5 - (math.cos((k_m * 2.0)) / 2.0))))) else: tmp = 2.0 / ((1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))) * (math.tan(k_m) * (math.sin(k_m) * math.pow((math.pow(t_m, 1.5) / 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 (t_m <= 2.6e-194) tmp = Float64(2.0 * Float64(Float64(cos(k_m) * Float64(l * l)) / Float64((k_m ^ 2.0) * Float64(t_m * Float64(0.5 - Float64(cos(Float64(k_m * 2.0)) / 2.0)))))); else tmp = Float64(2.0 / Float64(Float64(1.0 + Float64(1.0 + Float64(Float64(k_m / t_m) * Float64(k_m / t_m)))) * Float64(tan(k_m) * Float64(sin(k_m) * (Float64((t_m ^ 1.5) / 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 (t_m <= 2.6e-194) tmp = 2.0 * ((cos(k_m) * (l * l)) / ((k_m ^ 2.0) * (t_m * (0.5 - (cos((k_m * 2.0)) / 2.0))))); else tmp = 2.0 / ((1.0 + (1.0 + ((k_m / t_m) * (k_m / t_m)))) * (tan(k_m) * (sin(k_m) * (((t_m ^ 1.5) / 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[t$95$m, 2.6e-194], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] * N[(l * l), $MachinePrecision]), $MachinePrecision] / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[(0.5 - N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(1.0 + N[(1.0 + N[(N[(k$95$m / t$95$m), $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Tan[k$95$m], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Power[N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $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}\;t\_m \leq 2.6 \cdot 10^{-194}:\\
\;\;\;\;2 \cdot \frac{\cos k\_m \cdot \left(\ell \cdot \ell\right)}{{k\_m}^{2} \cdot \left(t\_m \cdot \left(0.5 - \frac{\cos \left(k\_m \cdot 2\right)}{2}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(1 + \left(1 + \frac{k\_m}{t\_m} \cdot \frac{k\_m}{t\_m}\right)\right) \cdot \left(\tan k\_m \cdot \left(\sin k\_m \cdot {\left(\frac{{t\_m}^{1.5}}{\ell}\right)}^{2}\right)\right)}\\
\end{array}
\end{array}
if t < 2.60000000000000002e-194Initial program 57.2%
Taylor expanded in t around 0 65.8%
unpow265.8%
sin-mult60.2%
Applied egg-rr60.2%
div-sub60.2%
+-inverses60.2%
cos-060.2%
metadata-eval60.2%
count-260.2%
Simplified60.2%
pow260.2%
Applied egg-rr60.2%
if 2.60000000000000002e-194 < t Initial program 63.4%
add-sqr-sqrt35.0%
pow235.0%
associate-/r*41.7%
*-commutative41.7%
sqrt-prod41.7%
associate-/r*35.0%
sqrt-div35.0%
sqrt-pow138.6%
metadata-eval38.6%
sqrt-prod27.3%
add-sqr-sqrt46.6%
Applied egg-rr46.6%
*-commutative46.6%
Simplified46.6%
unpow246.6%
Applied egg-rr46.6%
unpow-prod-down45.7%
pow245.7%
add-sqr-sqrt79.6%
Applied egg-rr79.6%
Final simplification68.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 (<= t_m 2.6e-55)
(*
2.0
(/ (* (cos k_m) (pow l 2.0)) (* (* k_m k_m) (* t_m (pow (sin k_m) 2.0)))))
(*
(* l (/ (/ 2.0 (tan k_m)) (* (sin k_m) (pow t_m 3.0))))
(/ l (+ 2.0 (pow (/ k_m 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 (t_m <= 2.6e-55) {
tmp = 2.0 * ((cos(k_m) * pow(l, 2.0)) / ((k_m * k_m) * (t_m * pow(sin(k_m), 2.0))));
} else {
tmp = (l * ((2.0 / tan(k_m)) / (sin(k_m) * pow(t_m, 3.0)))) * (l / (2.0 + pow((k_m / 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 (t_m <= 2.6d-55) then
tmp = 2.0d0 * ((cos(k_m) * (l ** 2.0d0)) / ((k_m * k_m) * (t_m * (sin(k_m) ** 2.0d0))))
else
tmp = (l * ((2.0d0 / tan(k_m)) / (sin(k_m) * (t_m ** 3.0d0)))) * (l / (2.0d0 + ((k_m / 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 (t_m <= 2.6e-55) {
tmp = 2.0 * ((Math.cos(k_m) * Math.pow(l, 2.0)) / ((k_m * k_m) * (t_m * Math.pow(Math.sin(k_m), 2.0))));
} else {
tmp = (l * ((2.0 / Math.tan(k_m)) / (Math.sin(k_m) * Math.pow(t_m, 3.0)))) * (l / (2.0 + Math.pow((k_m / 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 t_m <= 2.6e-55: tmp = 2.0 * ((math.cos(k_m) * math.pow(l, 2.0)) / ((k_m * k_m) * (t_m * math.pow(math.sin(k_m), 2.0)))) else: tmp = (l * ((2.0 / math.tan(k_m)) / (math.sin(k_m) * math.pow(t_m, 3.0)))) * (l / (2.0 + math.pow((k_m / 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 (t_m <= 2.6e-55) tmp = Float64(2.0 * Float64(Float64(cos(k_m) * (l ^ 2.0)) / Float64(Float64(k_m * k_m) * Float64(t_m * (sin(k_m) ^ 2.0))))); else tmp = Float64(Float64(l * Float64(Float64(2.0 / tan(k_m)) / Float64(sin(k_m) * (t_m ^ 3.0)))) * Float64(l / Float64(2.0 + (Float64(k_m / 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 (t_m <= 2.6e-55) tmp = 2.0 * ((cos(k_m) * (l ^ 2.0)) / ((k_m * k_m) * (t_m * (sin(k_m) ^ 2.0)))); else tmp = (l * ((2.0 / tan(k_m)) / (sin(k_m) * (t_m ^ 3.0)))) * (l / (2.0 + ((k_m / 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[t$95$m, 2.6e-55], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * N[(N[(2.0 / N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[k$95$m], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[(2.0 + N[Power[N[(k$95$m / t$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 \begin{array}{l}
\mathbf{if}\;t\_m \leq 2.6 \cdot 10^{-55}:\\
\;\;\;\;2 \cdot \frac{\cos k\_m \cdot {\ell}^{2}}{\left(k\_m \cdot k\_m\right) \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \frac{\frac{2}{\tan k\_m}}{\sin k\_m \cdot {t\_m}^{3}}\right) \cdot \frac{\ell}{2 + {\left(\frac{k\_m}{t\_m}\right)}^{2}}\\
\end{array}
\end{array}
if t < 2.5999999999999999e-55Initial program 54.9%
Taylor expanded in t around 0 66.4%
unpow266.4%
Applied egg-rr66.4%
if 2.5999999999999999e-55 < t Initial program 69.6%
Simplified67.4%
associate-*r*75.8%
*-un-lft-identity75.8%
times-frac78.0%
div-inv78.0%
frac-times77.9%
metadata-eval77.9%
Applied egg-rr77.9%
/-rgt-identity77.9%
*-commutative77.9%
associate-/r*78.0%
*-commutative78.0%
Simplified78.0%
Final simplification70.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 (<= t_m 9.6e-55)
(*
2.0
(/ (* (cos k_m) (pow l 2.0)) (* (* k_m k_m) (* t_m (pow (sin k_m) 2.0)))))
(*
(/ l (+ 2.0 (pow (/ k_m t_m) 2.0)))
(* l (/ 2.0 (* (tan k_m) (* (sin k_m) (pow t_m 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 (t_m <= 9.6e-55) {
tmp = 2.0 * ((cos(k_m) * pow(l, 2.0)) / ((k_m * k_m) * (t_m * pow(sin(k_m), 2.0))));
} else {
tmp = (l / (2.0 + pow((k_m / t_m), 2.0))) * (l * (2.0 / (tan(k_m) * (sin(k_m) * pow(t_m, 3.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 (t_m <= 9.6d-55) then
tmp = 2.0d0 * ((cos(k_m) * (l ** 2.0d0)) / ((k_m * k_m) * (t_m * (sin(k_m) ** 2.0d0))))
else
tmp = (l / (2.0d0 + ((k_m / t_m) ** 2.0d0))) * (l * (2.0d0 / (tan(k_m) * (sin(k_m) * (t_m ** 3.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 (t_m <= 9.6e-55) {
tmp = 2.0 * ((Math.cos(k_m) * Math.pow(l, 2.0)) / ((k_m * k_m) * (t_m * Math.pow(Math.sin(k_m), 2.0))));
} else {
tmp = (l / (2.0 + Math.pow((k_m / t_m), 2.0))) * (l * (2.0 / (Math.tan(k_m) * (Math.sin(k_m) * Math.pow(t_m, 3.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 t_m <= 9.6e-55: tmp = 2.0 * ((math.cos(k_m) * math.pow(l, 2.0)) / ((k_m * k_m) * (t_m * math.pow(math.sin(k_m), 2.0)))) else: tmp = (l / (2.0 + math.pow((k_m / t_m), 2.0))) * (l * (2.0 / (math.tan(k_m) * (math.sin(k_m) * math.pow(t_m, 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 (t_m <= 9.6e-55) tmp = Float64(2.0 * Float64(Float64(cos(k_m) * (l ^ 2.0)) / Float64(Float64(k_m * k_m) * Float64(t_m * (sin(k_m) ^ 2.0))))); else tmp = Float64(Float64(l / Float64(2.0 + (Float64(k_m / t_m) ^ 2.0))) * Float64(l * Float64(2.0 / Float64(tan(k_m) * Float64(sin(k_m) * (t_m ^ 3.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 (t_m <= 9.6e-55) tmp = 2.0 * ((cos(k_m) * (l ^ 2.0)) / ((k_m * k_m) * (t_m * (sin(k_m) ^ 2.0)))); else tmp = (l / (2.0 + ((k_m / t_m) ^ 2.0))) * (l * (2.0 / (tan(k_m) * (sin(k_m) * (t_m ^ 3.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[t$95$m, 9.6e-55], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / N[(2.0 + N[Power[N[(k$95$m / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l * N[(2.0 / N[(N[Tan[k$95$m], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $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}\;t\_m \leq 9.6 \cdot 10^{-55}:\\
\;\;\;\;2 \cdot \frac{\cos k\_m \cdot {\ell}^{2}}{\left(k\_m \cdot k\_m\right) \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{2 + {\left(\frac{k\_m}{t\_m}\right)}^{2}} \cdot \left(\ell \cdot \frac{2}{\tan k\_m \cdot \left(\sin k\_m \cdot {t\_m}^{3}\right)}\right)\\
\end{array}
\end{array}
if t < 9.59999999999999966e-55Initial program 54.9%
Taylor expanded in t around 0 66.4%
unpow266.4%
Applied egg-rr66.4%
if 9.59999999999999966e-55 < t Initial program 69.6%
Simplified67.4%
associate-*r*75.8%
*-un-lft-identity75.8%
times-frac78.0%
div-inv78.0%
frac-times77.9%
metadata-eval77.9%
Applied egg-rr77.9%
/-rgt-identity77.9%
*-commutative77.9%
*-commutative77.9%
*-commutative77.9%
Simplified77.9%
Final simplification70.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 (<= k_m 2.45e+20)
(/ 2.0 (* (sin k_m) (* (pow (/ (pow t_m 1.5) l) 2.0) (* k_m 2.0))))
(*
2.0
(/
(* (cos k_m) (pow l 2.0))
(* (* k_m k_m) (* 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 <= 2.45e+20) {
tmp = 2.0 / (sin(k_m) * (pow((pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0)));
} else {
tmp = 2.0 * ((cos(k_m) * pow(l, 2.0)) / ((k_m * k_m) * (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 <= 2.45d+20) then
tmp = 2.0d0 / (sin(k_m) * ((((t_m ** 1.5d0) / l) ** 2.0d0) * (k_m * 2.0d0)))
else
tmp = 2.0d0 * ((cos(k_m) * (l ** 2.0d0)) / ((k_m * k_m) * (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 <= 2.45e+20) {
tmp = 2.0 / (Math.sin(k_m) * (Math.pow((Math.pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0)));
} else {
tmp = 2.0 * ((Math.cos(k_m) * Math.pow(l, 2.0)) / ((k_m * k_m) * (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 <= 2.45e+20: tmp = 2.0 / (math.sin(k_m) * (math.pow((math.pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0))) else: tmp = 2.0 * ((math.cos(k_m) * math.pow(l, 2.0)) / ((k_m * k_m) * (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 <= 2.45e+20) tmp = Float64(2.0 / Float64(sin(k_m) * Float64((Float64((t_m ^ 1.5) / l) ^ 2.0) * Float64(k_m * 2.0)))); else tmp = Float64(2.0 * Float64(Float64(cos(k_m) * (l ^ 2.0)) / Float64(Float64(k_m * k_m) * 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 <= 2.45e+20) tmp = 2.0 / (sin(k_m) * ((((t_m ^ 1.5) / l) ^ 2.0) * (k_m * 2.0))); else tmp = 2.0 * ((cos(k_m) * (l ^ 2.0)) / ((k_m * k_m) * (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, 2.45e+20], N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Power[N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision], 2.0], $MachinePrecision] * N[(k$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] / N[(N[(k$95$m * k$95$m), $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.45 \cdot 10^{+20}:\\
\;\;\;\;\frac{2}{\sin k\_m \cdot \left({\left(\frac{{t\_m}^{1.5}}{\ell}\right)}^{2} \cdot \left(k\_m \cdot 2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{\cos k\_m \cdot {\ell}^{2}}{\left(k\_m \cdot k\_m\right) \cdot \left(t\_m \cdot {\sin k\_m}^{2}\right)}\\
\end{array}
\end{array}
if k < 2.45e20Initial program 61.8%
add-sqr-sqrt32.0%
pow232.0%
associate-/r*37.1%
*-commutative37.1%
sqrt-prod20.0%
associate-/r*16.8%
sqrt-div16.8%
sqrt-pow118.3%
metadata-eval18.3%
sqrt-prod13.0%
add-sqr-sqrt22.2%
Applied egg-rr22.2%
*-commutative22.2%
Simplified22.2%
div-inv22.2%
associate-*l*22.2%
*-commutative22.2%
unpow-prod-down21.7%
pow221.7%
add-sqr-sqrt39.0%
add-sqr-sqrt39.0%
pow239.0%
Applied egg-rr39.0%
associate-*r/39.0%
metadata-eval39.0%
associate-*l*38.6%
unpow238.6%
hypot-undefine38.6%
hypot-undefine38.6%
rem-square-sqrt38.6%
metadata-eval38.6%
hypot-undefine38.6%
metadata-eval38.6%
unpow238.6%
+-commutative38.6%
Simplified38.6%
Taylor expanded in k around 0 36.3%
if 2.45e20 < k Initial program 53.1%
Taylor expanded in t around 0 75.6%
unpow275.6%
Applied egg-rr75.6%
Final simplification45.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 (<= k_m 2.45e+20)
(/ 2.0 (* (sin k_m) (* (pow (/ (pow t_m 1.5) l) 2.0) (* k_m 2.0))))
(*
2.0
(/
(* (cos k_m) (* l l))
(* (pow k_m 2.0) (* t_m (- 0.5 (/ (cos (* 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) {
double tmp;
if (k_m <= 2.45e+20) {
tmp = 2.0 / (sin(k_m) * (pow((pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0)));
} else {
tmp = 2.0 * ((cos(k_m) * (l * l)) / (pow(k_m, 2.0) * (t_m * (0.5 - (cos((k_m * 2.0)) / 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 <= 2.45d+20) then
tmp = 2.0d0 / (sin(k_m) * ((((t_m ** 1.5d0) / l) ** 2.0d0) * (k_m * 2.0d0)))
else
tmp = 2.0d0 * ((cos(k_m) * (l * l)) / ((k_m ** 2.0d0) * (t_m * (0.5d0 - (cos((k_m * 2.0d0)) / 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 <= 2.45e+20) {
tmp = 2.0 / (Math.sin(k_m) * (Math.pow((Math.pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0)));
} else {
tmp = 2.0 * ((Math.cos(k_m) * (l * l)) / (Math.pow(k_m, 2.0) * (t_m * (0.5 - (Math.cos((k_m * 2.0)) / 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 <= 2.45e+20: tmp = 2.0 / (math.sin(k_m) * (math.pow((math.pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0))) else: tmp = 2.0 * ((math.cos(k_m) * (l * l)) / (math.pow(k_m, 2.0) * (t_m * (0.5 - (math.cos((k_m * 2.0)) / 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 <= 2.45e+20) tmp = Float64(2.0 / Float64(sin(k_m) * Float64((Float64((t_m ^ 1.5) / l) ^ 2.0) * Float64(k_m * 2.0)))); else tmp = Float64(2.0 * Float64(Float64(cos(k_m) * Float64(l * l)) / Float64((k_m ^ 2.0) * Float64(t_m * Float64(0.5 - Float64(cos(Float64(k_m * 2.0)) / 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 <= 2.45e+20) tmp = 2.0 / (sin(k_m) * ((((t_m ^ 1.5) / l) ^ 2.0) * (k_m * 2.0))); else tmp = 2.0 * ((cos(k_m) * (l * l)) / ((k_m ^ 2.0) * (t_m * (0.5 - (cos((k_m * 2.0)) / 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, 2.45e+20], N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Power[N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision], 2.0], $MachinePrecision] * N[(k$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] * N[(l * l), $MachinePrecision]), $MachinePrecision] / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[(0.5 - N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $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.45 \cdot 10^{+20}:\\
\;\;\;\;\frac{2}{\sin k\_m \cdot \left({\left(\frac{{t\_m}^{1.5}}{\ell}\right)}^{2} \cdot \left(k\_m \cdot 2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{\cos k\_m \cdot \left(\ell \cdot \ell\right)}{{k\_m}^{2} \cdot \left(t\_m \cdot \left(0.5 - \frac{\cos \left(k\_m \cdot 2\right)}{2}\right)\right)}\\
\end{array}
\end{array}
if k < 2.45e20Initial program 61.8%
add-sqr-sqrt32.0%
pow232.0%
associate-/r*37.1%
*-commutative37.1%
sqrt-prod20.0%
associate-/r*16.8%
sqrt-div16.8%
sqrt-pow118.3%
metadata-eval18.3%
sqrt-prod13.0%
add-sqr-sqrt22.2%
Applied egg-rr22.2%
*-commutative22.2%
Simplified22.2%
div-inv22.2%
associate-*l*22.2%
*-commutative22.2%
unpow-prod-down21.7%
pow221.7%
add-sqr-sqrt39.0%
add-sqr-sqrt39.0%
pow239.0%
Applied egg-rr39.0%
associate-*r/39.0%
metadata-eval39.0%
associate-*l*38.6%
unpow238.6%
hypot-undefine38.6%
hypot-undefine38.6%
rem-square-sqrt38.6%
metadata-eval38.6%
hypot-undefine38.6%
metadata-eval38.6%
unpow238.6%
+-commutative38.6%
Simplified38.6%
Taylor expanded in k around 0 36.3%
if 2.45e20 < k Initial program 53.1%
Taylor expanded in t around 0 75.6%
unpow275.6%
sin-mult75.5%
Applied egg-rr75.5%
div-sub75.5%
+-inverses75.5%
cos-075.5%
metadata-eval75.5%
count-275.5%
Simplified75.5%
pow275.5%
Applied egg-rr75.5%
Final simplification45.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 (<= k_m 2.5e+99)
(/ 2.0 (* (sin k_m) (* (pow (/ (pow t_m 1.5) l) 2.0) (* k_m 2.0))))
(*
2.0
(/
(pow l 2.0)
(* (pow k_m 2.0) (* t_m (- 0.5 (/ (cos (* 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) {
double tmp;
if (k_m <= 2.5e+99) {
tmp = 2.0 / (sin(k_m) * (pow((pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0)));
} else {
tmp = 2.0 * (pow(l, 2.0) / (pow(k_m, 2.0) * (t_m * (0.5 - (cos((k_m * 2.0)) / 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 <= 2.5d+99) then
tmp = 2.0d0 / (sin(k_m) * ((((t_m ** 1.5d0) / l) ** 2.0d0) * (k_m * 2.0d0)))
else
tmp = 2.0d0 * ((l ** 2.0d0) / ((k_m ** 2.0d0) * (t_m * (0.5d0 - (cos((k_m * 2.0d0)) / 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 <= 2.5e+99) {
tmp = 2.0 / (Math.sin(k_m) * (Math.pow((Math.pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0)));
} else {
tmp = 2.0 * (Math.pow(l, 2.0) / (Math.pow(k_m, 2.0) * (t_m * (0.5 - (Math.cos((k_m * 2.0)) / 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 <= 2.5e+99: tmp = 2.0 / (math.sin(k_m) * (math.pow((math.pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0))) else: tmp = 2.0 * (math.pow(l, 2.0) / (math.pow(k_m, 2.0) * (t_m * (0.5 - (math.cos((k_m * 2.0)) / 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 <= 2.5e+99) tmp = Float64(2.0 / Float64(sin(k_m) * Float64((Float64((t_m ^ 1.5) / l) ^ 2.0) * Float64(k_m * 2.0)))); else tmp = Float64(2.0 * Float64((l ^ 2.0) / Float64((k_m ^ 2.0) * Float64(t_m * Float64(0.5 - Float64(cos(Float64(k_m * 2.0)) / 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 <= 2.5e+99) tmp = 2.0 / (sin(k_m) * ((((t_m ^ 1.5) / l) ^ 2.0) * (k_m * 2.0))); else tmp = 2.0 * ((l ^ 2.0) / ((k_m ^ 2.0) * (t_m * (0.5 - (cos((k_m * 2.0)) / 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, 2.5e+99], N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Power[N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision], 2.0], $MachinePrecision] * N[(k$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[(0.5 - N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $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.5 \cdot 10^{+99}:\\
\;\;\;\;\frac{2}{\sin k\_m \cdot \left({\left(\frac{{t\_m}^{1.5}}{\ell}\right)}^{2} \cdot \left(k\_m \cdot 2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{{\ell}^{2}}{{k\_m}^{2} \cdot \left(t\_m \cdot \left(0.5 - \frac{\cos \left(k\_m \cdot 2\right)}{2}\right)\right)}\\
\end{array}
\end{array}
if k < 2.50000000000000004e99Initial program 62.3%
add-sqr-sqrt32.0%
pow232.0%
associate-/r*37.2%
*-commutative37.2%
sqrt-prod20.2%
associate-/r*17.2%
sqrt-div17.2%
sqrt-pow118.6%
metadata-eval18.6%
sqrt-prod13.2%
add-sqr-sqrt22.3%
Applied egg-rr22.3%
*-commutative22.3%
Simplified22.3%
div-inv22.3%
associate-*l*22.3%
*-commutative22.3%
unpow-prod-down21.8%
pow221.8%
add-sqr-sqrt39.4%
add-sqr-sqrt39.4%
pow239.4%
Applied egg-rr39.4%
associate-*r/39.4%
metadata-eval39.4%
associate-*l*39.1%
unpow239.1%
hypot-undefine39.1%
hypot-undefine39.1%
rem-square-sqrt39.1%
metadata-eval39.1%
hypot-undefine39.1%
metadata-eval39.1%
unpow239.1%
+-commutative39.1%
Simplified39.1%
Taylor expanded in k around 0 37.0%
if 2.50000000000000004e99 < k Initial program 48.1%
Taylor expanded in t around 0 72.8%
unpow272.8%
sin-mult72.6%
Applied egg-rr72.6%
div-sub72.6%
+-inverses72.6%
cos-072.6%
metadata-eval72.6%
count-272.6%
Simplified72.6%
Taylor expanded in k around 0 59.5%
Final simplification40.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
(if (<= k_m 2.6e-162)
(/
(* (* l l) (/ (/ 2.0 k_m) (* k_m (pow t_m 3.0))))
(+ 2.0 (pow (/ k_m t_m) 2.0)))
(if (<= k_m 2.9e+20)
(/ 2.0 (* (pow (/ (pow t_m 1.5) l) 2.0) (* 2.0 (* k_m k_m))))
(* 2.0 (/ (* (cos k_m) (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) {
double tmp;
if (k_m <= 2.6e-162) {
tmp = ((l * l) * ((2.0 / k_m) / (k_m * pow(t_m, 3.0)))) / (2.0 + pow((k_m / t_m), 2.0));
} else if (k_m <= 2.9e+20) {
tmp = 2.0 / (pow((pow(t_m, 1.5) / l), 2.0) * (2.0 * (k_m * k_m)));
} else {
tmp = 2.0 * ((cos(k_m) * pow(l, 2.0)) / (t_m * pow(k_m, 4.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 <= 2.6d-162) then
tmp = ((l * l) * ((2.0d0 / k_m) / (k_m * (t_m ** 3.0d0)))) / (2.0d0 + ((k_m / t_m) ** 2.0d0))
else if (k_m <= 2.9d+20) then
tmp = 2.0d0 / ((((t_m ** 1.5d0) / l) ** 2.0d0) * (2.0d0 * (k_m * k_m)))
else
tmp = 2.0d0 * ((cos(k_m) * (l ** 2.0d0)) / (t_m * (k_m ** 4.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 <= 2.6e-162) {
tmp = ((l * l) * ((2.0 / k_m) / (k_m * Math.pow(t_m, 3.0)))) / (2.0 + Math.pow((k_m / t_m), 2.0));
} else if (k_m <= 2.9e+20) {
tmp = 2.0 / (Math.pow((Math.pow(t_m, 1.5) / l), 2.0) * (2.0 * (k_m * k_m)));
} else {
tmp = 2.0 * ((Math.cos(k_m) * Math.pow(l, 2.0)) / (t_m * Math.pow(k_m, 4.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 <= 2.6e-162: tmp = ((l * l) * ((2.0 / k_m) / (k_m * math.pow(t_m, 3.0)))) / (2.0 + math.pow((k_m / t_m), 2.0)) elif k_m <= 2.9e+20: tmp = 2.0 / (math.pow((math.pow(t_m, 1.5) / l), 2.0) * (2.0 * (k_m * k_m))) else: tmp = 2.0 * ((math.cos(k_m) * math.pow(l, 2.0)) / (t_m * math.pow(k_m, 4.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 <= 2.6e-162) tmp = Float64(Float64(Float64(l * l) * Float64(Float64(2.0 / k_m) / Float64(k_m * (t_m ^ 3.0)))) / Float64(2.0 + (Float64(k_m / t_m) ^ 2.0))); elseif (k_m <= 2.9e+20) tmp = Float64(2.0 / Float64((Float64((t_m ^ 1.5) / l) ^ 2.0) * Float64(2.0 * Float64(k_m * k_m)))); else tmp = Float64(2.0 * Float64(Float64(cos(k_m) * (l ^ 2.0)) / Float64(t_m * (k_m ^ 4.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 <= 2.6e-162) tmp = ((l * l) * ((2.0 / k_m) / (k_m * (t_m ^ 3.0)))) / (2.0 + ((k_m / t_m) ^ 2.0)); elseif (k_m <= 2.9e+20) tmp = 2.0 / ((((t_m ^ 1.5) / l) ^ 2.0) * (2.0 * (k_m * k_m))); else tmp = 2.0 * ((cos(k_m) * (l ^ 2.0)) / (t_m * (k_m ^ 4.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, 2.6e-162], N[(N[(N[(l * l), $MachinePrecision] * N[(N[(2.0 / k$95$m), $MachinePrecision] / N[(k$95$m * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[Power[N[(k$95$m / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2.9e+20], N[(2.0 / N[(N[Power[N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision], 2.0], $MachinePrecision] * N[(2.0 * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $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 \begin{array}{l}
\mathbf{if}\;k\_m \leq 2.6 \cdot 10^{-162}:\\
\;\;\;\;\frac{\left(\ell \cdot \ell\right) \cdot \frac{\frac{2}{k\_m}}{k\_m \cdot {t\_m}^{3}}}{2 + {\left(\frac{k\_m}{t\_m}\right)}^{2}}\\
\mathbf{elif}\;k\_m \leq 2.9 \cdot 10^{+20}:\\
\;\;\;\;\frac{2}{{\left(\frac{{t\_m}^{1.5}}{\ell}\right)}^{2} \cdot \left(2 \cdot \left(k\_m \cdot k\_m\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{\cos k\_m \cdot {\ell}^{2}}{t\_m \cdot {k\_m}^{4}}\\
\end{array}
\end{array}
if k < 2.6e-162Initial program 60.4%
Simplified59.1%
Taylor expanded in k around 0 56.7%
Taylor expanded in k around 0 56.7%
if 2.6e-162 < k < 2.9e20Initial program 68.3%
Simplified77.3%
Taylor expanded in k around 0 81.8%
unpow258.5%
Applied egg-rr81.8%
add-sqr-sqrt38.1%
pow238.1%
associate-/r*33.0%
sqrt-div30.3%
sqrt-pow130.4%
metadata-eval30.4%
sqrt-prod24.6%
add-sqr-sqrt35.5%
Applied egg-rr35.5%
if 2.9e20 < k Initial program 53.1%
Taylor expanded in t around 0 75.6%
Taylor expanded in k around 0 60.6%
Final simplification54.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 2.8e+20)
(/ 2.0 (* (sin k_m) (* (pow (/ (pow t_m 1.5) l) 2.0) (* k_m 2.0))))
(* 2.0 (/ (* (cos k_m) (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) {
double tmp;
if (k_m <= 2.8e+20) {
tmp = 2.0 / (sin(k_m) * (pow((pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0)));
} else {
tmp = 2.0 * ((cos(k_m) * pow(l, 2.0)) / (t_m * pow(k_m, 4.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 <= 2.8d+20) then
tmp = 2.0d0 / (sin(k_m) * ((((t_m ** 1.5d0) / l) ** 2.0d0) * (k_m * 2.0d0)))
else
tmp = 2.0d0 * ((cos(k_m) * (l ** 2.0d0)) / (t_m * (k_m ** 4.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 <= 2.8e+20) {
tmp = 2.0 / (Math.sin(k_m) * (Math.pow((Math.pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0)));
} else {
tmp = 2.0 * ((Math.cos(k_m) * Math.pow(l, 2.0)) / (t_m * Math.pow(k_m, 4.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 <= 2.8e+20: tmp = 2.0 / (math.sin(k_m) * (math.pow((math.pow(t_m, 1.5) / l), 2.0) * (k_m * 2.0))) else: tmp = 2.0 * ((math.cos(k_m) * math.pow(l, 2.0)) / (t_m * math.pow(k_m, 4.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 <= 2.8e+20) tmp = Float64(2.0 / Float64(sin(k_m) * Float64((Float64((t_m ^ 1.5) / l) ^ 2.0) * Float64(k_m * 2.0)))); else tmp = Float64(2.0 * Float64(Float64(cos(k_m) * (l ^ 2.0)) / Float64(t_m * (k_m ^ 4.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 <= 2.8e+20) tmp = 2.0 / (sin(k_m) * ((((t_m ^ 1.5) / l) ^ 2.0) * (k_m * 2.0))); else tmp = 2.0 * ((cos(k_m) * (l ^ 2.0)) / (t_m * (k_m ^ 4.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, 2.8e+20], N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Power[N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision], 2.0], $MachinePrecision] * N[(k$95$m * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k$95$m], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $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 \begin{array}{l}
\mathbf{if}\;k\_m \leq 2.8 \cdot 10^{+20}:\\
\;\;\;\;\frac{2}{\sin k\_m \cdot \left({\left(\frac{{t\_m}^{1.5}}{\ell}\right)}^{2} \cdot \left(k\_m \cdot 2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{\cos k\_m \cdot {\ell}^{2}}{t\_m \cdot {k\_m}^{4}}\\
\end{array}
\end{array}
if k < 2.8e20Initial program 61.8%
add-sqr-sqrt32.0%
pow232.0%
associate-/r*37.1%
*-commutative37.1%
sqrt-prod20.0%
associate-/r*16.8%
sqrt-div16.8%
sqrt-pow118.3%
metadata-eval18.3%
sqrt-prod13.0%
add-sqr-sqrt22.2%
Applied egg-rr22.2%
*-commutative22.2%
Simplified22.2%
div-inv22.2%
associate-*l*22.2%
*-commutative22.2%
unpow-prod-down21.7%
pow221.7%
add-sqr-sqrt39.0%
add-sqr-sqrt39.0%
pow239.0%
Applied egg-rr39.0%
associate-*r/39.0%
metadata-eval39.0%
associate-*l*38.6%
unpow238.6%
hypot-undefine38.6%
hypot-undefine38.6%
rem-square-sqrt38.6%
metadata-eval38.6%
hypot-undefine38.6%
metadata-eval38.6%
unpow238.6%
+-commutative38.6%
Simplified38.6%
Taylor expanded in k around 0 36.3%
if 2.8e20 < k Initial program 53.1%
Taylor expanded in t around 0 75.6%
Taylor expanded in k around 0 60.6%
Final simplification41.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
(*
t_s
(if (<= k_m 2.1e-162)
(/
(* (* l l) (/ (/ 2.0 k_m) (* k_m (pow t_m 3.0))))
(+ 2.0 (pow (/ k_m t_m) 2.0)))
(if (<= k_m 5e+69)
(/ 2.0 (* (pow (/ (pow t_m 1.5) l) 2.0) (* 2.0 (* k_m k_m))))
(* (/ 2.0 t_m) (/ (pow l 2.0) (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) {
double tmp;
if (k_m <= 2.1e-162) {
tmp = ((l * l) * ((2.0 / k_m) / (k_m * pow(t_m, 3.0)))) / (2.0 + pow((k_m / t_m), 2.0));
} else if (k_m <= 5e+69) {
tmp = 2.0 / (pow((pow(t_m, 1.5) / l), 2.0) * (2.0 * (k_m * k_m)));
} else {
tmp = (2.0 / t_m) * (pow(l, 2.0) / pow(k_m, 4.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 <= 2.1d-162) then
tmp = ((l * l) * ((2.0d0 / k_m) / (k_m * (t_m ** 3.0d0)))) / (2.0d0 + ((k_m / t_m) ** 2.0d0))
else if (k_m <= 5d+69) then
tmp = 2.0d0 / ((((t_m ** 1.5d0) / l) ** 2.0d0) * (2.0d0 * (k_m * k_m)))
else
tmp = (2.0d0 / t_m) * ((l ** 2.0d0) / (k_m ** 4.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 <= 2.1e-162) {
tmp = ((l * l) * ((2.0 / k_m) / (k_m * Math.pow(t_m, 3.0)))) / (2.0 + Math.pow((k_m / t_m), 2.0));
} else if (k_m <= 5e+69) {
tmp = 2.0 / (Math.pow((Math.pow(t_m, 1.5) / l), 2.0) * (2.0 * (k_m * k_m)));
} else {
tmp = (2.0 / t_m) * (Math.pow(l, 2.0) / Math.pow(k_m, 4.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 <= 2.1e-162: tmp = ((l * l) * ((2.0 / k_m) / (k_m * math.pow(t_m, 3.0)))) / (2.0 + math.pow((k_m / t_m), 2.0)) elif k_m <= 5e+69: tmp = 2.0 / (math.pow((math.pow(t_m, 1.5) / l), 2.0) * (2.0 * (k_m * k_m))) else: tmp = (2.0 / t_m) * (math.pow(l, 2.0) / math.pow(k_m, 4.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 <= 2.1e-162) tmp = Float64(Float64(Float64(l * l) * Float64(Float64(2.0 / k_m) / Float64(k_m * (t_m ^ 3.0)))) / Float64(2.0 + (Float64(k_m / t_m) ^ 2.0))); elseif (k_m <= 5e+69) tmp = Float64(2.0 / Float64((Float64((t_m ^ 1.5) / l) ^ 2.0) * Float64(2.0 * Float64(k_m * k_m)))); else tmp = Float64(Float64(2.0 / t_m) * Float64((l ^ 2.0) / (k_m ^ 4.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 <= 2.1e-162) tmp = ((l * l) * ((2.0 / k_m) / (k_m * (t_m ^ 3.0)))) / (2.0 + ((k_m / t_m) ^ 2.0)); elseif (k_m <= 5e+69) tmp = 2.0 / ((((t_m ^ 1.5) / l) ^ 2.0) * (2.0 * (k_m * k_m))); else tmp = (2.0 / t_m) * ((l ^ 2.0) / (k_m ^ 4.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, 2.1e-162], N[(N[(N[(l * l), $MachinePrecision] * N[(N[(2.0 / k$95$m), $MachinePrecision] / N[(k$95$m * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[Power[N[(k$95$m / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 5e+69], N[(2.0 / N[(N[Power[N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision], 2.0], $MachinePrecision] * N[(2.0 * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / t$95$m), $MachinePrecision] * N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k$95$m, 4.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 2.1 \cdot 10^{-162}:\\
\;\;\;\;\frac{\left(\ell \cdot \ell\right) \cdot \frac{\frac{2}{k\_m}}{k\_m \cdot {t\_m}^{3}}}{2 + {\left(\frac{k\_m}{t\_m}\right)}^{2}}\\
\mathbf{elif}\;k\_m \leq 5 \cdot 10^{+69}:\\
\;\;\;\;\frac{2}{{\left(\frac{{t\_m}^{1.5}}{\ell}\right)}^{2} \cdot \left(2 \cdot \left(k\_m \cdot k\_m\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t\_m} \cdot \frac{{\ell}^{2}}{{k\_m}^{4}}\\
\end{array}
\end{array}
if k < 2.1e-162Initial program 60.4%
Simplified59.1%
Taylor expanded in k around 0 56.7%
Taylor expanded in k around 0 56.7%
if 2.1e-162 < k < 5.00000000000000036e69Initial program 68.7%
Simplified76.4%
Taylor expanded in k around 0 78.0%
unpow262.9%
Applied egg-rr78.0%
add-sqr-sqrt38.9%
pow238.9%
associate-/r*34.6%
sqrt-div32.3%
sqrt-pow132.3%
metadata-eval32.3%
sqrt-prod25.2%
add-sqr-sqrt36.7%
Applied egg-rr36.7%
if 5.00000000000000036e69 < k Initial program 50.6%
Taylor expanded in t around 0 74.2%
unpow274.2%
sin-mult74.0%
Applied egg-rr74.0%
div-sub74.0%
+-inverses74.0%
cos-074.0%
metadata-eval74.0%
count-274.0%
Simplified74.0%
Taylor expanded in k around 0 58.8%
associate-*r/58.8%
*-commutative58.8%
times-frac58.8%
Simplified58.8%
Final simplification53.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
(*
t_s
(if (<= t_m 2.35e-80)
(* (/ 2.0 t_m) (/ (pow l 2.0) (pow k_m 4.0)))
(/ 2.0 (* (pow (/ (pow t_m 1.5) l) 2.0) (* 2.0 (* k_m 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 (t_m <= 2.35e-80) {
tmp = (2.0 / t_m) * (pow(l, 2.0) / pow(k_m, 4.0));
} else {
tmp = 2.0 / (pow((pow(t_m, 1.5) / l), 2.0) * (2.0 * (k_m * 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 (t_m <= 2.35d-80) then
tmp = (2.0d0 / t_m) * ((l ** 2.0d0) / (k_m ** 4.0d0))
else
tmp = 2.0d0 / ((((t_m ** 1.5d0) / l) ** 2.0d0) * (2.0d0 * (k_m * 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 (t_m <= 2.35e-80) {
tmp = (2.0 / t_m) * (Math.pow(l, 2.0) / Math.pow(k_m, 4.0));
} else {
tmp = 2.0 / (Math.pow((Math.pow(t_m, 1.5) / l), 2.0) * (2.0 * (k_m * 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 t_m <= 2.35e-80: tmp = (2.0 / t_m) * (math.pow(l, 2.0) / math.pow(k_m, 4.0)) else: tmp = 2.0 / (math.pow((math.pow(t_m, 1.5) / l), 2.0) * (2.0 * (k_m * 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 (t_m <= 2.35e-80) tmp = Float64(Float64(2.0 / t_m) * Float64((l ^ 2.0) / (k_m ^ 4.0))); else tmp = Float64(2.0 / Float64((Float64((t_m ^ 1.5) / l) ^ 2.0) * Float64(2.0 * Float64(k_m * 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 (t_m <= 2.35e-80) tmp = (2.0 / t_m) * ((l ^ 2.0) / (k_m ^ 4.0)); else tmp = 2.0 / ((((t_m ^ 1.5) / l) ^ 2.0) * (2.0 * (k_m * 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[t$95$m, 2.35e-80], N[(N[(2.0 / t$95$m), $MachinePrecision] * N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Power[N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision], 2.0], $MachinePrecision] * N[(2.0 * N[(k$95$m * 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}\;t\_m \leq 2.35 \cdot 10^{-80}:\\
\;\;\;\;\frac{2}{t\_m} \cdot \frac{{\ell}^{2}}{{k\_m}^{4}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{{t\_m}^{1.5}}{\ell}\right)}^{2} \cdot \left(2 \cdot \left(k\_m \cdot k\_m\right)\right)}\\
\end{array}
\end{array}
if t < 2.34999999999999986e-80Initial program 53.8%
Taylor expanded in t around 0 65.5%
unpow265.5%
sin-mult60.5%
Applied egg-rr60.5%
div-sub60.5%
+-inverses60.5%
cos-060.5%
metadata-eval60.5%
count-260.5%
Simplified60.5%
Taylor expanded in k around 0 56.3%
associate-*r/55.7%
*-commutative55.7%
times-frac57.4%
Simplified57.4%
if 2.34999999999999986e-80 < t Initial program 71.0%
Simplified66.5%
Taylor expanded in k around 0 62.2%
unpow252.7%
Applied egg-rr62.2%
add-sqr-sqrt62.2%
pow262.2%
associate-/r*56.5%
sqrt-div56.4%
sqrt-pow156.5%
metadata-eval56.5%
sqrt-prod31.0%
add-sqr-sqrt64.3%
Applied egg-rr64.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 (<= t_m 2.8e-83)
(* (/ 2.0 t_m) (/ (pow l 2.0) (pow k_m 4.0)))
(/ 2.0 (* (* 2.0 (* k_m k_m)) (* (/ t_m l) (/ (* t_m t_m) 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 (t_m <= 2.8e-83) {
tmp = (2.0 / t_m) * (pow(l, 2.0) / pow(k_m, 4.0));
} else {
tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / 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 (t_m <= 2.8d-83) then
tmp = (2.0d0 / t_m) * ((l ** 2.0d0) / (k_m ** 4.0d0))
else
tmp = 2.0d0 / ((2.0d0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / 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 (t_m <= 2.8e-83) {
tmp = (2.0 / t_m) * (Math.pow(l, 2.0) / Math.pow(k_m, 4.0));
} else {
tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / 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 t_m <= 2.8e-83: tmp = (2.0 / t_m) * (math.pow(l, 2.0) / math.pow(k_m, 4.0)) else: tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / 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 (t_m <= 2.8e-83) tmp = Float64(Float64(2.0 / t_m) * Float64((l ^ 2.0) / (k_m ^ 4.0))); else tmp = Float64(2.0 / Float64(Float64(2.0 * Float64(k_m * k_m)) * Float64(Float64(t_m / l) * Float64(Float64(t_m * t_m) / 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 (t_m <= 2.8e-83) tmp = (2.0 / t_m) * ((l ^ 2.0) / (k_m ^ 4.0)); else tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / 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[t$95$m, 2.8e-83], N[(N[(2.0 / t$95$m), $MachinePrecision] * N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(2.0 * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(t$95$m / l), $MachinePrecision] * N[(N[(t$95$m * t$95$m), $MachinePrecision] / l), $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}\;t\_m \leq 2.8 \cdot 10^{-83}:\\
\;\;\;\;\frac{2}{t\_m} \cdot \frac{{\ell}^{2}}{{k\_m}^{4}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(2 \cdot \left(k\_m \cdot k\_m\right)\right) \cdot \left(\frac{t\_m}{\ell} \cdot \frac{t\_m \cdot t\_m}{\ell}\right)}\\
\end{array}
\end{array}
if t < 2.8000000000000001e-83Initial program 53.8%
Taylor expanded in t around 0 65.5%
unpow265.5%
sin-mult60.5%
Applied egg-rr60.5%
div-sub60.5%
+-inverses60.5%
cos-060.5%
metadata-eval60.5%
count-260.5%
Simplified60.5%
Taylor expanded in k around 0 56.3%
associate-*r/55.7%
*-commutative55.7%
times-frac57.4%
Simplified57.4%
if 2.8000000000000001e-83 < t Initial program 71.0%
Simplified66.5%
Taylor expanded in k around 0 62.2%
unpow252.7%
Applied egg-rr62.2%
associate-/r*56.4%
unpow356.4%
times-frac62.2%
pow262.2%
Applied egg-rr62.2%
unpow262.2%
Applied egg-rr62.2%
Final simplification59.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 (<= k_m 2.1e+58)
(/ 2.0 (* (* 2.0 (* k_m k_m)) (* (/ t_m l) (/ (* t_m t_m) l))))
(* 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) {
double tmp;
if (k_m <= 2.1e+58) {
tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / l)));
} else {
tmp = 2.0 * (pow(l, 2.0) / (t_m * pow(k_m, 4.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 <= 2.1d+58) then
tmp = 2.0d0 / ((2.0d0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / l)))
else
tmp = 2.0d0 * ((l ** 2.0d0) / (t_m * (k_m ** 4.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 <= 2.1e+58) {
tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / l)));
} else {
tmp = 2.0 * (Math.pow(l, 2.0) / (t_m * Math.pow(k_m, 4.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 <= 2.1e+58: tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / l))) else: tmp = 2.0 * (math.pow(l, 2.0) / (t_m * math.pow(k_m, 4.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 <= 2.1e+58) tmp = Float64(2.0 / Float64(Float64(2.0 * Float64(k_m * k_m)) * Float64(Float64(t_m / l) * Float64(Float64(t_m * t_m) / l)))); else tmp = Float64(2.0 * Float64((l ^ 2.0) / Float64(t_m * (k_m ^ 4.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 <= 2.1e+58) tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / l))); else tmp = 2.0 * ((l ^ 2.0) / (t_m * (k_m ^ 4.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, 2.1e+58], N[(2.0 / N[(N[(2.0 * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(t$95$m / l), $MachinePrecision] * N[(N[(t$95$m * t$95$m), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 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 \begin{array}{l}
\mathbf{if}\;k\_m \leq 2.1 \cdot 10^{+58}:\\
\;\;\;\;\frac{2}{\left(2 \cdot \left(k\_m \cdot k\_m\right)\right) \cdot \left(\frac{t\_m}{\ell} \cdot \frac{t\_m \cdot t\_m}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{{\ell}^{2}}{t\_m \cdot {k\_m}^{4}}\\
\end{array}
\end{array}
if k < 2.10000000000000012e58Initial program 61.8%
Simplified60.0%
Taylor expanded in k around 0 59.5%
unpow257.4%
Applied egg-rr59.5%
associate-/r*54.3%
unpow354.3%
times-frac61.5%
pow261.5%
Applied egg-rr61.5%
unpow261.5%
Applied egg-rr61.5%
if 2.10000000000000012e58 < k Initial program 52.5%
Taylor expanded in t around 0 75.1%
Taylor expanded in k around 0 58.5%
Final simplification60.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 (* (* 2.0 (* k_m k_m)) (* (/ t_m l) (/ (* t_m t_m) 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) {
return t_s * (2.0 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / l))));
}
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 / ((2.0d0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / l))))
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 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / l))));
}
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 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / l))))
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(Float64(2.0 * Float64(k_m * k_m)) * Float64(Float64(t_m / l) * Float64(Float64(t_m * t_m) / l))))) 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 / ((2.0 * (k_m * k_m)) * ((t_m / l) * ((t_m * t_m) / l)))); 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[(2.0 * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(t$95$m / l), $MachinePrecision] * N[(N[(t$95$m * t$95$m), $MachinePrecision] / l), $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 \frac{2}{\left(2 \cdot \left(k\_m \cdot k\_m\right)\right) \cdot \left(\frac{t\_m}{\ell} \cdot \frac{t\_m \cdot t\_m}{\ell}\right)}
\end{array}
Initial program 59.9%
Simplified59.7%
Taylor expanded in k around 0 57.4%
unpow261.0%
Applied egg-rr57.4%
associate-/r*52.8%
unpow352.8%
times-frac59.0%
pow259.0%
Applied egg-rr59.0%
unpow259.0%
Applied egg-rr59.0%
Final simplification59.0%
herbie shell --seed 2024156
(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))))