
(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 27 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}
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (pow (sin k) 3.0)))
(*
t_s
(if (<= t_m 8.8e-134)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(cbrt (* t_2 t_2)))
(/
2.0
(pow
(*
(* (/ t_m (pow (cbrt l) 2.0)) (cbrt (sin k)))
(cbrt (* (tan k) (+ 2.0 (pow (/ k t_m) 2.0)))))
3.0))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = pow(sin(k), 3.0);
double tmp;
if (t_m <= 8.8e-134) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / cbrt((t_2 * t_2));
} else {
tmp = 2.0 / pow((((t_m / pow(cbrt(l), 2.0)) * cbrt(sin(k))) * cbrt((tan(k) * (2.0 + pow((k / t_m), 2.0))))), 3.0);
}
return t_s * tmp;
}
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) {
double t_2 = Math.pow(Math.sin(k), 3.0);
double tmp;
if (t_m <= 8.8e-134) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / Math.cbrt((t_2 * t_2));
} else {
tmp = 2.0 / Math.pow((((t_m / Math.pow(Math.cbrt(l), 2.0)) * Math.cbrt(Math.sin(k))) * Math.cbrt((Math.tan(k) * (2.0 + Math.pow((k / t_m), 2.0))))), 3.0);
}
return t_s * tmp;
}
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = sin(k) ^ 3.0 tmp = 0.0 if (t_m <= 8.8e-134) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / cbrt(Float64(t_2 * t_2))); else tmp = Float64(2.0 / (Float64(Float64(Float64(t_m / (cbrt(l) ^ 2.0)) * cbrt(sin(k))) * cbrt(Float64(tan(k) * Float64(2.0 + (Float64(k / t_m) ^ 2.0))))) ^ 3.0)); end return Float64(t_s * tmp) end
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_] := Block[{t$95$2 = N[Power[N[Sin[k], $MachinePrecision], 3.0], $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 8.8e-134], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[(t$95$2 * t$95$2), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[(N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] * N[Power[N[(N[Tan[k], $MachinePrecision] * N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := {\sin k}^{3}\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 8.8 \cdot 10^{-134}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{\sqrt[3]{t_2 \cdot t_2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\left(\frac{t_m}{{\left(\sqrt[3]{\ell}\right)}^{2}} \cdot \sqrt[3]{\sin k}\right) \cdot \sqrt[3]{\tan k \cdot \left(2 + {\left(\frac{k}{t_m}\right)}^{2}\right)}\right)}^{3}}\\
\end{array}
\end{array}
\end{array}
if t < 8.7999999999999999e-134Initial program 48.5%
Simplified48.5%
Taylor expanded in k around inf 55.5%
associate-*r/55.5%
associate-*r*55.5%
associate-/r*55.6%
Simplified55.6%
unpow255.6%
add-cbrt-cube55.6%
add-cbrt-cube55.6%
cbrt-unprod53.8%
pow353.8%
pow353.8%
Applied egg-rr53.8%
if 8.7999999999999999e-134 < t Initial program 58.8%
associate-/r*63.9%
add-cube-cbrt63.8%
*-un-lft-identity63.8%
times-frac63.8%
pow263.8%
cbrt-div63.9%
rem-cbrt-cube63.9%
cbrt-div63.8%
rem-cbrt-cube81.6%
Applied egg-rr81.6%
add-cube-cbrt81.4%
pow381.4%
cbrt-prod81.3%
frac-times78.1%
unpow278.1%
*-un-lft-identity78.1%
cbrt-div81.3%
add-cbrt-cube88.9%
Applied egg-rr88.9%
add-cube-cbrt88.7%
pow388.8%
Applied egg-rr93.0%
associate-+r+93.0%
metadata-eval93.0%
Simplified93.0%
Final simplification67.0%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (* (/ t_m (pow (cbrt l) 2.0)) (cbrt (sin k))))
(t_3 (hypot 1.0 (hypot 1.0 (/ k t_m)))))
(*
t_s
(if (<= t_m 1.45e-100)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(pow (sin k) 2.0))
(if (<= t_m 2.7e+103)
(*
(/ l t_3)
(/ (/ 2.0 (* (sin k) (* (tan k) (pow t_m 3.0)))) (/ t_3 l)))
(if (or (<= t_m 7e+135) (not (<= t_m 2e+251)))
(/ 2.0 (pow (* t_2 (cbrt (* 2.0 k))) 3.0))
(/
2.0
(* (* (tan k) (+ 2.0 (pow (/ k t_m) 2.0))) (pow t_2 3.0)))))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = (t_m / pow(cbrt(l), 2.0)) * cbrt(sin(k));
double t_3 = hypot(1.0, hypot(1.0, (k / t_m)));
double tmp;
if (t_m <= 1.45e-100) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / pow(sin(k), 2.0);
} else if (t_m <= 2.7e+103) {
tmp = (l / t_3) * ((2.0 / (sin(k) * (tan(k) * pow(t_m, 3.0)))) / (t_3 / l));
} else if ((t_m <= 7e+135) || !(t_m <= 2e+251)) {
tmp = 2.0 / pow((t_2 * cbrt((2.0 * k))), 3.0);
} else {
tmp = 2.0 / ((tan(k) * (2.0 + pow((k / t_m), 2.0))) * pow(t_2, 3.0));
}
return t_s * tmp;
}
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) {
double t_2 = (t_m / Math.pow(Math.cbrt(l), 2.0)) * Math.cbrt(Math.sin(k));
double t_3 = Math.hypot(1.0, Math.hypot(1.0, (k / t_m)));
double tmp;
if (t_m <= 1.45e-100) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / Math.pow(Math.sin(k), 2.0);
} else if (t_m <= 2.7e+103) {
tmp = (l / t_3) * ((2.0 / (Math.sin(k) * (Math.tan(k) * Math.pow(t_m, 3.0)))) / (t_3 / l));
} else if ((t_m <= 7e+135) || !(t_m <= 2e+251)) {
tmp = 2.0 / Math.pow((t_2 * Math.cbrt((2.0 * k))), 3.0);
} else {
tmp = 2.0 / ((Math.tan(k) * (2.0 + Math.pow((k / t_m), 2.0))) * Math.pow(t_2, 3.0));
}
return t_s * tmp;
}
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(Float64(t_m / (cbrt(l) ^ 2.0)) * cbrt(sin(k))) t_3 = hypot(1.0, hypot(1.0, Float64(k / t_m))) tmp = 0.0 if (t_m <= 1.45e-100) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / (sin(k) ^ 2.0)); elseif (t_m <= 2.7e+103) tmp = Float64(Float64(l / t_3) * Float64(Float64(2.0 / Float64(sin(k) * Float64(tan(k) * (t_m ^ 3.0)))) / Float64(t_3 / l))); elseif ((t_m <= 7e+135) || !(t_m <= 2e+251)) tmp = Float64(2.0 / (Float64(t_2 * cbrt(Float64(2.0 * k))) ^ 3.0)); else tmp = Float64(2.0 / Float64(Float64(tan(k) * Float64(2.0 + (Float64(k / t_m) ^ 2.0))) * (t_2 ^ 3.0))); end return Float64(t_s * tmp) end
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_] := Block[{t$95$2 = N[(N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[Sqrt[1.0 ^ 2 + N[Sqrt[1.0 ^ 2 + N[(k / t$95$m), $MachinePrecision] ^ 2], $MachinePrecision] ^ 2], $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 1.45e-100], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 2.7e+103], N[(N[(l / t$95$3), $MachinePrecision] * N[(N[(2.0 / N[(N[Sin[k], $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$3 / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t$95$m, 7e+135], N[Not[LessEqual[t$95$m, 2e+251]], $MachinePrecision]], N[(2.0 / N[Power[N[(t$95$2 * N[Power[N[(2.0 * k), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Power[t$95$2, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]), $MachinePrecision]]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \frac{t_m}{{\left(\sqrt[3]{\ell}\right)}^{2}} \cdot \sqrt[3]{\sin k}\\
t_3 := \mathsf{hypot}\left(1, \mathsf{hypot}\left(1, \frac{k}{t_m}\right)\right)\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 1.45 \cdot 10^{-100}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{{\sin k}^{2}}\\
\mathbf{elif}\;t_m \leq 2.7 \cdot 10^{+103}:\\
\;\;\;\;\frac{\ell}{t_3} \cdot \frac{\frac{2}{\sin k \cdot \left(\tan k \cdot {t_m}^{3}\right)}}{\frac{t_3}{\ell}}\\
\mathbf{elif}\;t_m \leq 7 \cdot 10^{+135} \lor \neg \left(t_m \leq 2 \cdot 10^{+251}\right):\\
\;\;\;\;\frac{2}{{\left(t_2 \cdot \sqrt[3]{2 \cdot k}\right)}^{3}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot \left(2 + {\left(\frac{k}{t_m}\right)}^{2}\right)\right) \cdot {t_2}^{3}}\\
\end{array}
\end{array}
\end{array}
if t < 1.44999999999999988e-100Initial program 48.9%
Simplified48.3%
Taylor expanded in k around inf 55.6%
associate-*r/55.6%
associate-*r*55.5%
associate-/r*56.7%
Simplified56.7%
if 1.44999999999999988e-100 < t < 2.69999999999999993e103Initial program 75.8%
Simplified80.1%
associate-*r*86.5%
add-sqr-sqrt86.4%
times-frac86.7%
metadata-eval86.7%
associate-+r+86.7%
add-sqr-sqrt86.7%
hypot-1-def86.7%
unpow286.7%
hypot-1-def86.7%
metadata-eval86.7%
Applied egg-rr93.1%
associate-/l*95.1%
associate-*r*95.1%
Simplified95.1%
if 2.69999999999999993e103 < t < 7.0000000000000005e135 or 2.0000000000000001e251 < t Initial program 41.1%
associate-/r*49.0%
add-cube-cbrt49.0%
*-un-lft-identity49.0%
times-frac49.0%
pow249.0%
cbrt-div49.0%
rem-cbrt-cube49.0%
cbrt-div49.0%
rem-cbrt-cube68.1%
Applied egg-rr68.1%
add-cube-cbrt68.1%
pow368.1%
cbrt-prod68.2%
frac-times68.2%
unpow268.2%
*-un-lft-identity68.2%
cbrt-div68.1%
add-cbrt-cube80.7%
Applied egg-rr80.7%
add-cube-cbrt80.4%
pow380.4%
Applied egg-rr92.4%
associate-+r+92.4%
metadata-eval92.4%
Simplified92.4%
Taylor expanded in k around 0 79.9%
*-commutative79.9%
Simplified79.9%
if 7.0000000000000005e135 < t < 2.0000000000000001e251Initial program 33.3%
associate-/r*40.8%
add-cube-cbrt40.8%
*-un-lft-identity40.8%
times-frac40.8%
pow240.8%
cbrt-div40.8%
rem-cbrt-cube40.8%
cbrt-div40.8%
rem-cbrt-cube93.9%
Applied egg-rr93.9%
add-cube-cbrt93.9%
pow393.9%
cbrt-prod93.4%
frac-times77.9%
unpow277.9%
*-un-lft-identity77.9%
cbrt-div77.8%
add-cbrt-cube98.6%
Applied egg-rr98.6%
add-cube-cbrt98.0%
pow398.0%
Applied egg-rr98.5%
*-commutative98.5%
cube-prod98.6%
rem-cube-cbrt98.8%
associate-+r+98.8%
metadata-eval98.8%
Simplified98.8%
Final simplification67.6%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (* (/ t_m (pow (cbrt l) 2.0)) (cbrt (sin k))))
(t_3 (hypot 1.0 (hypot 1.0 (/ k t_m)))))
(*
t_s
(if (<= t_m 1.55e-99)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(pow (sin k) 2.0))
(if (<= t_m 2.7e+103)
(*
(/ (/ (/ 2.0 (tan k)) (* (sin k) (pow t_m 3.0))) (/ t_3 l))
(/ l t_3))
(if (or (<= t_m 2.5e+134) (not (<= t_m 3.7e+251)))
(/ 2.0 (pow (* t_2 (cbrt (* 2.0 k))) 3.0))
(/
2.0
(* (* (tan k) (+ 2.0 (pow (/ k t_m) 2.0))) (pow t_2 3.0)))))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = (t_m / pow(cbrt(l), 2.0)) * cbrt(sin(k));
double t_3 = hypot(1.0, hypot(1.0, (k / t_m)));
double tmp;
if (t_m <= 1.55e-99) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / pow(sin(k), 2.0);
} else if (t_m <= 2.7e+103) {
tmp = (((2.0 / tan(k)) / (sin(k) * pow(t_m, 3.0))) / (t_3 / l)) * (l / t_3);
} else if ((t_m <= 2.5e+134) || !(t_m <= 3.7e+251)) {
tmp = 2.0 / pow((t_2 * cbrt((2.0 * k))), 3.0);
} else {
tmp = 2.0 / ((tan(k) * (2.0 + pow((k / t_m), 2.0))) * pow(t_2, 3.0));
}
return t_s * tmp;
}
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) {
double t_2 = (t_m / Math.pow(Math.cbrt(l), 2.0)) * Math.cbrt(Math.sin(k));
double t_3 = Math.hypot(1.0, Math.hypot(1.0, (k / t_m)));
double tmp;
if (t_m <= 1.55e-99) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / Math.pow(Math.sin(k), 2.0);
} else if (t_m <= 2.7e+103) {
tmp = (((2.0 / Math.tan(k)) / (Math.sin(k) * Math.pow(t_m, 3.0))) / (t_3 / l)) * (l / t_3);
} else if ((t_m <= 2.5e+134) || !(t_m <= 3.7e+251)) {
tmp = 2.0 / Math.pow((t_2 * Math.cbrt((2.0 * k))), 3.0);
} else {
tmp = 2.0 / ((Math.tan(k) * (2.0 + Math.pow((k / t_m), 2.0))) * Math.pow(t_2, 3.0));
}
return t_s * tmp;
}
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(Float64(t_m / (cbrt(l) ^ 2.0)) * cbrt(sin(k))) t_3 = hypot(1.0, hypot(1.0, Float64(k / t_m))) tmp = 0.0 if (t_m <= 1.55e-99) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / (sin(k) ^ 2.0)); elseif (t_m <= 2.7e+103) tmp = Float64(Float64(Float64(Float64(2.0 / tan(k)) / Float64(sin(k) * (t_m ^ 3.0))) / Float64(t_3 / l)) * Float64(l / t_3)); elseif ((t_m <= 2.5e+134) || !(t_m <= 3.7e+251)) tmp = Float64(2.0 / (Float64(t_2 * cbrt(Float64(2.0 * k))) ^ 3.0)); else tmp = Float64(2.0 / Float64(Float64(tan(k) * Float64(2.0 + (Float64(k / t_m) ^ 2.0))) * (t_2 ^ 3.0))); end return Float64(t_s * tmp) end
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_] := Block[{t$95$2 = N[(N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[Sqrt[1.0 ^ 2 + N[Sqrt[1.0 ^ 2 + N[(k / t$95$m), $MachinePrecision] ^ 2], $MachinePrecision] ^ 2], $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 1.55e-99], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 2.7e+103], N[(N[(N[(N[(2.0 / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$3 / l), $MachinePrecision]), $MachinePrecision] * N[(l / t$95$3), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t$95$m, 2.5e+134], N[Not[LessEqual[t$95$m, 3.7e+251]], $MachinePrecision]], N[(2.0 / N[Power[N[(t$95$2 * N[Power[N[(2.0 * k), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Power[t$95$2, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]), $MachinePrecision]]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \frac{t_m}{{\left(\sqrt[3]{\ell}\right)}^{2}} \cdot \sqrt[3]{\sin k}\\
t_3 := \mathsf{hypot}\left(1, \mathsf{hypot}\left(1, \frac{k}{t_m}\right)\right)\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 1.55 \cdot 10^{-99}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{{\sin k}^{2}}\\
\mathbf{elif}\;t_m \leq 2.7 \cdot 10^{+103}:\\
\;\;\;\;\frac{\frac{\frac{2}{\tan k}}{\sin k \cdot {t_m}^{3}}}{\frac{t_3}{\ell}} \cdot \frac{\ell}{t_3}\\
\mathbf{elif}\;t_m \leq 2.5 \cdot 10^{+134} \lor \neg \left(t_m \leq 3.7 \cdot 10^{+251}\right):\\
\;\;\;\;\frac{2}{{\left(t_2 \cdot \sqrt[3]{2 \cdot k}\right)}^{3}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot \left(2 + {\left(\frac{k}{t_m}\right)}^{2}\right)\right) \cdot {t_2}^{3}}\\
\end{array}
\end{array}
\end{array}
if t < 1.5499999999999999e-99Initial program 48.9%
Simplified48.3%
Taylor expanded in k around inf 55.6%
associate-*r/55.6%
associate-*r*55.5%
associate-/r*56.7%
Simplified56.7%
if 1.5499999999999999e-99 < t < 2.69999999999999993e103Initial program 75.8%
Simplified80.1%
associate-*r*86.5%
add-sqr-sqrt86.4%
times-frac86.7%
metadata-eval86.7%
associate-+r+86.7%
add-sqr-sqrt86.7%
hypot-1-def86.7%
unpow286.7%
hypot-1-def86.7%
metadata-eval86.7%
Applied egg-rr93.1%
associate-/l*95.1%
associate-/r*95.0%
Simplified95.0%
if 2.69999999999999993e103 < t < 2.4999999999999999e134 or 3.6999999999999999e251 < t Initial program 41.1%
associate-/r*49.0%
add-cube-cbrt49.0%
*-un-lft-identity49.0%
times-frac49.0%
pow249.0%
cbrt-div49.0%
rem-cbrt-cube49.0%
cbrt-div49.0%
rem-cbrt-cube68.1%
Applied egg-rr68.1%
add-cube-cbrt68.1%
pow368.1%
cbrt-prod68.2%
frac-times68.2%
unpow268.2%
*-un-lft-identity68.2%
cbrt-div68.1%
add-cbrt-cube80.7%
Applied egg-rr80.7%
add-cube-cbrt80.4%
pow380.4%
Applied egg-rr92.4%
associate-+r+92.4%
metadata-eval92.4%
Simplified92.4%
Taylor expanded in k around 0 79.9%
*-commutative79.9%
Simplified79.9%
if 2.4999999999999999e134 < t < 3.6999999999999999e251Initial program 33.3%
associate-/r*40.8%
add-cube-cbrt40.8%
*-un-lft-identity40.8%
times-frac40.8%
pow240.8%
cbrt-div40.8%
rem-cbrt-cube40.8%
cbrt-div40.8%
rem-cbrt-cube93.9%
Applied egg-rr93.9%
add-cube-cbrt93.9%
pow393.9%
cbrt-prod93.4%
frac-times77.9%
unpow277.9%
*-un-lft-identity77.9%
cbrt-div77.8%
add-cbrt-cube98.6%
Applied egg-rr98.6%
add-cube-cbrt98.0%
pow398.0%
Applied egg-rr98.5%
*-commutative98.5%
cube-prod98.6%
rem-cube-cbrt98.8%
associate-+r+98.8%
metadata-eval98.8%
Simplified98.8%
Final simplification67.6%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (cbrt (sin k))))
(*
t_s
(if (<= t_m 6.5e-55)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(pow (sin k) 2.0))
(if (<= t_m 4.2e+251)
(/
2.0
(*
(* (tan k) (pow (* t_2 (/ (/ t_m (cbrt l)) (cbrt l))) 3.0))
(+ 1.0 (+ (pow (/ k t_m) 2.0) 1.0))))
(/
2.0
(pow
(* (* (/ t_m (pow (cbrt l) 2.0)) t_2) (cbrt (* 2.0 k)))
3.0)))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = cbrt(sin(k));
double tmp;
if (t_m <= 6.5e-55) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / pow(sin(k), 2.0);
} else if (t_m <= 4.2e+251) {
tmp = 2.0 / ((tan(k) * pow((t_2 * ((t_m / cbrt(l)) / cbrt(l))), 3.0)) * (1.0 + (pow((k / t_m), 2.0) + 1.0)));
} else {
tmp = 2.0 / pow((((t_m / pow(cbrt(l), 2.0)) * t_2) * cbrt((2.0 * k))), 3.0);
}
return t_s * tmp;
}
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) {
double t_2 = Math.cbrt(Math.sin(k));
double tmp;
if (t_m <= 6.5e-55) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / Math.pow(Math.sin(k), 2.0);
} else if (t_m <= 4.2e+251) {
tmp = 2.0 / ((Math.tan(k) * Math.pow((t_2 * ((t_m / Math.cbrt(l)) / Math.cbrt(l))), 3.0)) * (1.0 + (Math.pow((k / t_m), 2.0) + 1.0)));
} else {
tmp = 2.0 / Math.pow((((t_m / Math.pow(Math.cbrt(l), 2.0)) * t_2) * Math.cbrt((2.0 * k))), 3.0);
}
return t_s * tmp;
}
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = cbrt(sin(k)) tmp = 0.0 if (t_m <= 6.5e-55) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / (sin(k) ^ 2.0)); elseif (t_m <= 4.2e+251) tmp = Float64(2.0 / Float64(Float64(tan(k) * (Float64(t_2 * Float64(Float64(t_m / cbrt(l)) / cbrt(l))) ^ 3.0)) * Float64(1.0 + Float64((Float64(k / t_m) ^ 2.0) + 1.0)))); else tmp = Float64(2.0 / (Float64(Float64(Float64(t_m / (cbrt(l) ^ 2.0)) * t_2) * cbrt(Float64(2.0 * k))) ^ 3.0)); end return Float64(t_s * tmp) end
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_] := Block[{t$95$2 = N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 6.5e-55], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 4.2e+251], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * N[Power[N[(t$95$2 * N[(N[(t$95$m / N[Power[l, 1/3], $MachinePrecision]), $MachinePrecision] / N[Power[l, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[(N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision] * N[Power[N[(2.0 * k), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sqrt[3]{\sin k}\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 6.5 \cdot 10^{-55}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{{\sin k}^{2}}\\
\mathbf{elif}\;t_m \leq 4.2 \cdot 10^{+251}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot {\left(t_2 \cdot \frac{\frac{t_m}{\sqrt[3]{\ell}}}{\sqrt[3]{\ell}}\right)}^{3}\right) \cdot \left(1 + \left({\left(\frac{k}{t_m}\right)}^{2} + 1\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\left(\frac{t_m}{{\left(\sqrt[3]{\ell}\right)}^{2}} \cdot t_2\right) \cdot \sqrt[3]{2 \cdot k}\right)}^{3}}\\
\end{array}
\end{array}
\end{array}
if t < 6.50000000000000006e-55Initial program 50.0%
Simplified49.5%
Taylor expanded in k around inf 57.0%
associate-*r/57.0%
associate-*r*57.0%
associate-/r*58.1%
Simplified58.1%
if 6.50000000000000006e-55 < t < 4.2000000000000001e251Initial program 57.7%
associate-/r*62.9%
add-cube-cbrt62.7%
*-un-lft-identity62.7%
times-frac62.7%
pow262.7%
cbrt-div62.9%
rem-cbrt-cube62.8%
cbrt-div62.7%
rem-cbrt-cube82.4%
Applied egg-rr82.4%
add-cube-cbrt82.3%
pow382.2%
cbrt-prod82.1%
frac-times77.7%
unpow277.7%
*-un-lft-identity77.7%
cbrt-div82.1%
add-cbrt-cube89.4%
Applied egg-rr89.4%
if 4.2000000000000001e251 < t Initial program 51.0%
associate-/r*64.8%
add-cube-cbrt64.8%
*-un-lft-identity64.8%
times-frac64.8%
pow264.8%
cbrt-div64.8%
rem-cbrt-cube64.8%
cbrt-div64.8%
rem-cbrt-cube64.8%
Applied egg-rr64.8%
add-cube-cbrt64.8%
pow364.8%
cbrt-prod64.8%
frac-times64.8%
unpow264.8%
*-un-lft-identity64.8%
cbrt-div64.8%
add-cbrt-cube88.3%
Applied egg-rr88.3%
add-cube-cbrt87.9%
pow387.8%
Applied egg-rr98.8%
associate-+r+98.8%
metadata-eval98.8%
Simplified98.8%
Taylor expanded in k around 0 87.2%
*-commutative87.2%
Simplified87.2%
Final simplification66.7%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (* (/ t_m (pow (cbrt l) 2.0)) (cbrt (sin k)))))
(*
t_s
(if (<= t_m 2.45e-67)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(pow (sin k) 2.0))
(if (<= t_m 3.6e+251)
(/ 2.0 (* (* (tan k) (+ 2.0 (pow (/ k t_m) 2.0))) (pow t_2 3.0)))
(/ 2.0 (pow (* t_2 (cbrt (* 2.0 k))) 3.0)))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = (t_m / pow(cbrt(l), 2.0)) * cbrt(sin(k));
double tmp;
if (t_m <= 2.45e-67) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / pow(sin(k), 2.0);
} else if (t_m <= 3.6e+251) {
tmp = 2.0 / ((tan(k) * (2.0 + pow((k / t_m), 2.0))) * pow(t_2, 3.0));
} else {
tmp = 2.0 / pow((t_2 * cbrt((2.0 * k))), 3.0);
}
return t_s * tmp;
}
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) {
double t_2 = (t_m / Math.pow(Math.cbrt(l), 2.0)) * Math.cbrt(Math.sin(k));
double tmp;
if (t_m <= 2.45e-67) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / Math.pow(Math.sin(k), 2.0);
} else if (t_m <= 3.6e+251) {
tmp = 2.0 / ((Math.tan(k) * (2.0 + Math.pow((k / t_m), 2.0))) * Math.pow(t_2, 3.0));
} else {
tmp = 2.0 / Math.pow((t_2 * Math.cbrt((2.0 * k))), 3.0);
}
return t_s * tmp;
}
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(Float64(t_m / (cbrt(l) ^ 2.0)) * cbrt(sin(k))) tmp = 0.0 if (t_m <= 2.45e-67) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / (sin(k) ^ 2.0)); elseif (t_m <= 3.6e+251) tmp = Float64(2.0 / Float64(Float64(tan(k) * Float64(2.0 + (Float64(k / t_m) ^ 2.0))) * (t_2 ^ 3.0))); else tmp = Float64(2.0 / (Float64(t_2 * cbrt(Float64(2.0 * k))) ^ 3.0)); end return Float64(t_s * tmp) end
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_] := Block[{t$95$2 = N[(N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 2.45e-67], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 3.6e+251], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Power[t$95$2, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(t$95$2 * N[Power[N[(2.0 * k), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \frac{t_m}{{\left(\sqrt[3]{\ell}\right)}^{2}} \cdot \sqrt[3]{\sin k}\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 2.45 \cdot 10^{-67}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{{\sin k}^{2}}\\
\mathbf{elif}\;t_m \leq 3.6 \cdot 10^{+251}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot \left(2 + {\left(\frac{k}{t_m}\right)}^{2}\right)\right) \cdot {t_2}^{3}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(t_2 \cdot \sqrt[3]{2 \cdot k}\right)}^{3}}\\
\end{array}
\end{array}
\end{array}
if t < 2.44999999999999997e-67Initial program 50.0%
Simplified49.5%
Taylor expanded in k around inf 57.0%
associate-*r/57.0%
associate-*r*57.0%
associate-/r*58.1%
Simplified58.1%
if 2.44999999999999997e-67 < t < 3.59999999999999997e251Initial program 57.7%
associate-/r*62.9%
add-cube-cbrt62.7%
*-un-lft-identity62.7%
times-frac62.7%
pow262.7%
cbrt-div62.9%
rem-cbrt-cube62.8%
cbrt-div62.7%
rem-cbrt-cube82.4%
Applied egg-rr82.4%
add-cube-cbrt82.3%
pow382.2%
cbrt-prod82.1%
frac-times77.7%
unpow277.7%
*-un-lft-identity77.7%
cbrt-div82.1%
add-cbrt-cube89.4%
Applied egg-rr89.4%
add-cube-cbrt89.3%
pow389.3%
Applied egg-rr92.3%
*-commutative92.3%
cube-prod89.3%
rem-cube-cbrt89.4%
associate-+r+89.4%
metadata-eval89.4%
Simplified89.4%
if 3.59999999999999997e251 < t Initial program 51.0%
associate-/r*64.8%
add-cube-cbrt64.8%
*-un-lft-identity64.8%
times-frac64.8%
pow264.8%
cbrt-div64.8%
rem-cbrt-cube64.8%
cbrt-div64.8%
rem-cbrt-cube64.8%
Applied egg-rr64.8%
add-cube-cbrt64.8%
pow364.8%
cbrt-prod64.8%
frac-times64.8%
unpow264.8%
*-un-lft-identity64.8%
cbrt-div64.8%
add-cbrt-cube88.3%
Applied egg-rr88.3%
add-cube-cbrt87.9%
pow387.8%
Applied egg-rr98.8%
associate-+r+98.8%
metadata-eval98.8%
Simplified98.8%
Taylor expanded in k around 0 87.2%
*-commutative87.2%
Simplified87.2%
Final simplification66.7%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (+ 2.0 (pow (/ k t_m) 2.0))) (t_3 (/ t_m (pow (cbrt l) 2.0))))
(*
t_s
(if (<= t_m 5.2e-68)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(pow (sin k) 2.0))
(if (<= t_m 7e+97)
(* (* l (/ (/ 2.0 (tan k)) (* (sin k) (pow t_m 3.0)))) (/ l t_2))
(if (<= t_m 1.65e+251)
(/ 2.0 (* (* (tan k) t_2) (* (sin k) (pow t_3 3.0))))
(/ 2.0 (pow (* (* t_3 (cbrt (sin k))) (cbrt (* 2.0 k))) 3.0))))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = 2.0 + pow((k / t_m), 2.0);
double t_3 = t_m / pow(cbrt(l), 2.0);
double tmp;
if (t_m <= 5.2e-68) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / pow(sin(k), 2.0);
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / tan(k)) / (sin(k) * pow(t_m, 3.0)))) * (l / t_2);
} else if (t_m <= 1.65e+251) {
tmp = 2.0 / ((tan(k) * t_2) * (sin(k) * pow(t_3, 3.0)));
} else {
tmp = 2.0 / pow(((t_3 * cbrt(sin(k))) * cbrt((2.0 * k))), 3.0);
}
return t_s * tmp;
}
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) {
double t_2 = 2.0 + Math.pow((k / t_m), 2.0);
double t_3 = t_m / Math.pow(Math.cbrt(l), 2.0);
double tmp;
if (t_m <= 5.2e-68) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / Math.pow(Math.sin(k), 2.0);
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / Math.tan(k)) / (Math.sin(k) * Math.pow(t_m, 3.0)))) * (l / t_2);
} else if (t_m <= 1.65e+251) {
tmp = 2.0 / ((Math.tan(k) * t_2) * (Math.sin(k) * Math.pow(t_3, 3.0)));
} else {
tmp = 2.0 / Math.pow(((t_3 * Math.cbrt(Math.sin(k))) * Math.cbrt((2.0 * k))), 3.0);
}
return t_s * tmp;
}
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(2.0 + (Float64(k / t_m) ^ 2.0)) t_3 = Float64(t_m / (cbrt(l) ^ 2.0)) tmp = 0.0 if (t_m <= 5.2e-68) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / (sin(k) ^ 2.0)); elseif (t_m <= 7e+97) tmp = Float64(Float64(l * Float64(Float64(2.0 / tan(k)) / Float64(sin(k) * (t_m ^ 3.0)))) * Float64(l / t_2)); elseif (t_m <= 1.65e+251) tmp = Float64(2.0 / Float64(Float64(tan(k) * t_2) * Float64(sin(k) * (t_3 ^ 3.0)))); else tmp = Float64(2.0 / (Float64(Float64(t_3 * cbrt(sin(k))) * cbrt(Float64(2.0 * k))) ^ 3.0)); end return Float64(t_s * tmp) end
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_] := Block[{t$95$2 = N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 5.2e-68], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 7e+97], N[(N[(l * N[(N[(2.0 / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / t$95$2), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 1.65e+251], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * t$95$2), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$3, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[(t$95$3 * N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] * N[Power[N[(2.0 * k), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]]]]), $MachinePrecision]]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := 2 + {\left(\frac{k}{t_m}\right)}^{2}\\
t_3 := \frac{t_m}{{\left(\sqrt[3]{\ell}\right)}^{2}}\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 5.2 \cdot 10^{-68}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{{\sin k}^{2}}\\
\mathbf{elif}\;t_m \leq 7 \cdot 10^{+97}:\\
\;\;\;\;\left(\ell \cdot \frac{\frac{2}{\tan k}}{\sin k \cdot {t_m}^{3}}\right) \cdot \frac{\ell}{t_2}\\
\mathbf{elif}\;t_m \leq 1.65 \cdot 10^{+251}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot t_2\right) \cdot \left(\sin k \cdot {t_3}^{3}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\left(t_3 \cdot \sqrt[3]{\sin k}\right) \cdot \sqrt[3]{2 \cdot k}\right)}^{3}}\\
\end{array}
\end{array}
\end{array}
if t < 5.1999999999999996e-68Initial program 49.8%
Simplified49.2%
Taylor expanded in k around inf 56.8%
associate-*r/56.8%
associate-*r*56.7%
associate-/r*57.9%
Simplified57.9%
if 5.1999999999999996e-68 < t < 7.0000000000000001e97Initial program 74.3%
Simplified79.3%
associate-*r*86.8%
*-un-lft-identity86.8%
times-frac86.9%
Applied egg-rr86.9%
/-rgt-identity86.9%
*-commutative86.9%
associate-/r*86.9%
Simplified86.9%
if 7.0000000000000001e97 < t < 1.65000000000000009e251Initial program 34.9%
associate-/r*40.4%
add-cube-cbrt40.4%
*-un-lft-identity40.4%
times-frac40.4%
pow240.4%
cbrt-div40.4%
rem-cbrt-cube40.4%
cbrt-div40.4%
rem-cbrt-cube88.2%
Applied egg-rr88.2%
add-cube-cbrt88.2%
pow388.2%
cbrt-prod87.9%
frac-times77.2%
unpow277.2%
*-un-lft-identity77.2%
cbrt-div77.1%
add-cbrt-cube91.5%
Applied egg-rr91.5%
distribute-lft-in91.5%
Applied egg-rr91.4%
distribute-lft-out91.4%
+-commutative91.4%
associate-*r*91.4%
associate-+r+91.4%
metadata-eval91.4%
Simplified91.4%
if 1.65000000000000009e251 < t Initial program 51.0%
associate-/r*64.8%
add-cube-cbrt64.8%
*-un-lft-identity64.8%
times-frac64.8%
pow264.8%
cbrt-div64.8%
rem-cbrt-cube64.8%
cbrt-div64.8%
rem-cbrt-cube64.8%
Applied egg-rr64.8%
add-cube-cbrt64.8%
pow364.8%
cbrt-prod64.8%
frac-times64.8%
unpow264.8%
*-un-lft-identity64.8%
cbrt-div64.8%
add-cbrt-cube88.3%
Applied egg-rr88.3%
add-cube-cbrt87.9%
pow387.8%
Applied egg-rr98.8%
associate-+r+98.8%
metadata-eval98.8%
Simplified98.8%
Taylor expanded in k around 0 87.2%
*-commutative87.2%
Simplified87.2%
Final simplification66.5%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.22e-67)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(pow (sin k) 2.0))
(if (<= t_m 6.4e+72)
(/
2.0
(*
(+ 1.0 (+ (pow (/ k t_m) 2.0) 1.0))
(* (tan k) (/ (* (sin k) (/ (pow t_m 3.0) l)) l))))
(/
2.0
(pow
(* (* (/ t_m (pow (cbrt l) 2.0)) (cbrt (sin k))) (cbrt (* 2.0 k)))
3.0))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.22e-67) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / pow(sin(k), 2.0);
} else if (t_m <= 6.4e+72) {
tmp = 2.0 / ((1.0 + (pow((k / t_m), 2.0) + 1.0)) * (tan(k) * ((sin(k) * (pow(t_m, 3.0) / l)) / l)));
} else {
tmp = 2.0 / pow((((t_m / pow(cbrt(l), 2.0)) * cbrt(sin(k))) * cbrt((2.0 * k))), 3.0);
}
return t_s * tmp;
}
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) {
double tmp;
if (t_m <= 1.22e-67) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / Math.pow(Math.sin(k), 2.0);
} else if (t_m <= 6.4e+72) {
tmp = 2.0 / ((1.0 + (Math.pow((k / t_m), 2.0) + 1.0)) * (Math.tan(k) * ((Math.sin(k) * (Math.pow(t_m, 3.0) / l)) / l)));
} else {
tmp = 2.0 / Math.pow((((t_m / Math.pow(Math.cbrt(l), 2.0)) * Math.cbrt(Math.sin(k))) * Math.cbrt((2.0 * k))), 3.0);
}
return t_s * tmp;
}
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.22e-67) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / (sin(k) ^ 2.0)); elseif (t_m <= 6.4e+72) tmp = Float64(2.0 / Float64(Float64(1.0 + Float64((Float64(k / t_m) ^ 2.0) + 1.0)) * Float64(tan(k) * Float64(Float64(sin(k) * Float64((t_m ^ 3.0) / l)) / l)))); else tmp = Float64(2.0 / (Float64(Float64(Float64(t_m / (cbrt(l) ^ 2.0)) * cbrt(sin(k))) * cbrt(Float64(2.0 * k))) ^ 3.0)); end return Float64(t_s * tmp) end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 1.22e-67], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 6.4e+72], N[(2.0 / N[(N[(1.0 + N[(N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(N[(N[Sin[k], $MachinePrecision] * N[(N[Power[t$95$m, 3.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[(N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] * N[Power[N[(2.0 * k), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 1.22 \cdot 10^{-67}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{{\sin k}^{2}}\\
\mathbf{elif}\;t_m \leq 6.4 \cdot 10^{+72}:\\
\;\;\;\;\frac{2}{\left(1 + \left({\left(\frac{k}{t_m}\right)}^{2} + 1\right)\right) \cdot \left(\tan k \cdot \frac{\sin k \cdot \frac{{t_m}^{3}}{\ell}}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\left(\frac{t_m}{{\left(\sqrt[3]{\ell}\right)}^{2}} \cdot \sqrt[3]{\sin k}\right) \cdot \sqrt[3]{2 \cdot k}\right)}^{3}}\\
\end{array}
\end{array}
if t < 1.22e-67Initial program 49.8%
Simplified49.2%
Taylor expanded in k around inf 56.8%
associate-*r/56.8%
associate-*r*56.7%
associate-/r*57.9%
Simplified57.9%
if 1.22e-67 < t < 6.4000000000000003e72Initial program 86.1%
associate-/r*86.0%
associate-*l/92.7%
Applied egg-rr92.7%
if 6.4000000000000003e72 < t Initial program 38.3%
associate-/r*48.5%
add-cube-cbrt48.5%
*-un-lft-identity48.5%
times-frac48.5%
pow248.5%
cbrt-div48.6%
rem-cbrt-cube48.6%
cbrt-div48.6%
rem-cbrt-cube77.5%
Applied egg-rr77.5%
add-cube-cbrt77.3%
pow377.3%
cbrt-prod77.1%
frac-times70.7%
unpow270.7%
*-un-lft-identity70.7%
cbrt-div72.9%
add-cbrt-cube88.0%
Applied egg-rr88.0%
add-cube-cbrt87.6%
pow387.6%
Applied egg-rr94.0%
associate-+r+94.0%
metadata-eval94.0%
Simplified94.0%
Taylor expanded in k around 0 76.0%
*-commutative76.0%
Simplified76.0%
Final simplification64.9%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (pow (/ k t_m) 2.0)))
(*
t_s
(if (<= t_m 9.5e-69)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(pow (sin k) 2.0))
(if (<= t_m 7e+97)
(*
(* l (/ (/ 2.0 (tan k)) (* (sin k) (pow t_m 3.0))))
(/ l (+ 2.0 t_2)))
(/
2.0
(*
(+ 1.0 (+ t_2 1.0))
(*
(tan k)
(* (sin k) (pow (* t_m (pow l -0.6666666666666666)) 3.0))))))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = pow((k / t_m), 2.0);
double tmp;
if (t_m <= 9.5e-69) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / pow(sin(k), 2.0);
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / tan(k)) / (sin(k) * pow(t_m, 3.0)))) * (l / (2.0 + t_2));
} else {
tmp = 2.0 / ((1.0 + (t_2 + 1.0)) * (tan(k) * (sin(k) * pow((t_m * pow(l, -0.6666666666666666)), 3.0))));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: tmp
t_2 = (k / t_m) ** 2.0d0
if (t_m <= 9.5d-69) then
tmp = ((2.0d0 * ((l ** 2.0d0) * cos(k))) / (t_m * (k ** 2.0d0))) / (sin(k) ** 2.0d0)
else if (t_m <= 7d+97) then
tmp = (l * ((2.0d0 / tan(k)) / (sin(k) * (t_m ** 3.0d0)))) * (l / (2.0d0 + t_2))
else
tmp = 2.0d0 / ((1.0d0 + (t_2 + 1.0d0)) * (tan(k) * (sin(k) * ((t_m * (l ** (-0.6666666666666666d0))) ** 3.0d0))))
end if
code = t_s * tmp
end function
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) {
double t_2 = Math.pow((k / t_m), 2.0);
double tmp;
if (t_m <= 9.5e-69) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / Math.pow(Math.sin(k), 2.0);
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / Math.tan(k)) / (Math.sin(k) * Math.pow(t_m, 3.0)))) * (l / (2.0 + t_2));
} else {
tmp = 2.0 / ((1.0 + (t_2 + 1.0)) * (Math.tan(k) * (Math.sin(k) * Math.pow((t_m * Math.pow(l, -0.6666666666666666)), 3.0))));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): t_2 = math.pow((k / t_m), 2.0) tmp = 0 if t_m <= 9.5e-69: tmp = ((2.0 * (math.pow(l, 2.0) * math.cos(k))) / (t_m * math.pow(k, 2.0))) / math.pow(math.sin(k), 2.0) elif t_m <= 7e+97: tmp = (l * ((2.0 / math.tan(k)) / (math.sin(k) * math.pow(t_m, 3.0)))) * (l / (2.0 + t_2)) else: tmp = 2.0 / ((1.0 + (t_2 + 1.0)) * (math.tan(k) * (math.sin(k) * math.pow((t_m * math.pow(l, -0.6666666666666666)), 3.0)))) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(k / t_m) ^ 2.0 tmp = 0.0 if (t_m <= 9.5e-69) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / (sin(k) ^ 2.0)); elseif (t_m <= 7e+97) tmp = Float64(Float64(l * Float64(Float64(2.0 / tan(k)) / Float64(sin(k) * (t_m ^ 3.0)))) * Float64(l / Float64(2.0 + t_2))); else tmp = Float64(2.0 / Float64(Float64(1.0 + Float64(t_2 + 1.0)) * Float64(tan(k) * Float64(sin(k) * (Float64(t_m * (l ^ -0.6666666666666666)) ^ 3.0))))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) t_2 = (k / t_m) ^ 2.0; tmp = 0.0; if (t_m <= 9.5e-69) tmp = ((2.0 * ((l ^ 2.0) * cos(k))) / (t_m * (k ^ 2.0))) / (sin(k) ^ 2.0); elseif (t_m <= 7e+97) tmp = (l * ((2.0 / tan(k)) / (sin(k) * (t_m ^ 3.0)))) * (l / (2.0 + t_2)); else tmp = 2.0 / ((1.0 + (t_2 + 1.0)) * (tan(k) * (sin(k) * ((t_m * (l ^ -0.6666666666666666)) ^ 3.0)))); end tmp_2 = t_s * tmp; end
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_] := Block[{t$95$2 = N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 9.5e-69], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 7e+97], N[(N[(l * N[(N[(2.0 / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[(2.0 + t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(1.0 + N[(t$95$2 + 1.0), $MachinePrecision]), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[Power[N[(t$95$m * N[Power[l, -0.6666666666666666], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := {\left(\frac{k}{t_m}\right)}^{2}\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 9.5 \cdot 10^{-69}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{{\sin k}^{2}}\\
\mathbf{elif}\;t_m \leq 7 \cdot 10^{+97}:\\
\;\;\;\;\left(\ell \cdot \frac{\frac{2}{\tan k}}{\sin k \cdot {t_m}^{3}}\right) \cdot \frac{\ell}{2 + t_2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(1 + \left(t_2 + 1\right)\right) \cdot \left(\tan k \cdot \left(\sin k \cdot {\left(t_m \cdot {\ell}^{-0.6666666666666666}\right)}^{3}\right)\right)}\\
\end{array}
\end{array}
\end{array}
if t < 9.50000000000000094e-69Initial program 49.8%
Simplified49.2%
Taylor expanded in k around inf 56.8%
associate-*r/56.8%
associate-*r*56.7%
associate-/r*57.9%
Simplified57.9%
if 9.50000000000000094e-69 < t < 7.0000000000000001e97Initial program 74.3%
Simplified79.3%
associate-*r*86.8%
*-un-lft-identity86.8%
times-frac86.9%
Applied egg-rr86.9%
/-rgt-identity86.9%
*-commutative86.9%
associate-/r*86.9%
Simplified86.9%
if 7.0000000000000001e97 < t Initial program 38.7%
add-cube-cbrt38.7%
pow338.7%
div-inv38.7%
cbrt-prod38.7%
rem-cbrt-cube51.8%
pow251.8%
pow-flip51.8%
metadata-eval51.8%
Applied egg-rr51.8%
pow1/351.5%
pow-pow41.2%
metadata-eval41.2%
Applied egg-rr41.2%
Final simplification60.0%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (+ 2.0 (pow (/ k t_m) 2.0))))
(*
t_s
(if (<= t_m 3e-68)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(pow (sin k) 2.0))
(if (<= t_m 7e+97)
(* (* l (/ (/ 2.0 (tan k)) (* (sin k) (pow t_m 3.0)))) (/ l t_2))
(/
2.0
(* (* (tan k) t_2) (* (sin k) (/ (pow (/ t_m (cbrt l)) 3.0) l)))))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = 2.0 + pow((k / t_m), 2.0);
double tmp;
if (t_m <= 3e-68) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / pow(sin(k), 2.0);
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / tan(k)) / (sin(k) * pow(t_m, 3.0)))) * (l / t_2);
} else {
tmp = 2.0 / ((tan(k) * t_2) * (sin(k) * (pow((t_m / cbrt(l)), 3.0) / l)));
}
return t_s * tmp;
}
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) {
double t_2 = 2.0 + Math.pow((k / t_m), 2.0);
double tmp;
if (t_m <= 3e-68) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / Math.pow(Math.sin(k), 2.0);
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / Math.tan(k)) / (Math.sin(k) * Math.pow(t_m, 3.0)))) * (l / t_2);
} else {
tmp = 2.0 / ((Math.tan(k) * t_2) * (Math.sin(k) * (Math.pow((t_m / Math.cbrt(l)), 3.0) / l)));
}
return t_s * tmp;
}
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(2.0 + (Float64(k / t_m) ^ 2.0)) tmp = 0.0 if (t_m <= 3e-68) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / (sin(k) ^ 2.0)); elseif (t_m <= 7e+97) tmp = Float64(Float64(l * Float64(Float64(2.0 / tan(k)) / Float64(sin(k) * (t_m ^ 3.0)))) * Float64(l / t_2)); else tmp = Float64(2.0 / Float64(Float64(tan(k) * t_2) * Float64(sin(k) * Float64((Float64(t_m / cbrt(l)) ^ 3.0) / l)))); end return Float64(t_s * tmp) end
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_] := Block[{t$95$2 = N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 3e-68], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 7e+97], N[(N[(l * N[(N[(2.0 / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / t$95$2), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * t$95$2), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[Power[N[(t$95$m / N[Power[l, 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := 2 + {\left(\frac{k}{t_m}\right)}^{2}\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 3 \cdot 10^{-68}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{{\sin k}^{2}}\\
\mathbf{elif}\;t_m \leq 7 \cdot 10^{+97}:\\
\;\;\;\;\left(\ell \cdot \frac{\frac{2}{\tan k}}{\sin k \cdot {t_m}^{3}}\right) \cdot \frac{\ell}{t_2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot t_2\right) \cdot \left(\sin k \cdot \frac{{\left(\frac{t_m}{\sqrt[3]{\ell}}\right)}^{3}}{\ell}\right)}\\
\end{array}
\end{array}
\end{array}
if t < 3e-68Initial program 49.8%
Simplified49.2%
Taylor expanded in k around inf 56.8%
associate-*r/56.8%
associate-*r*56.7%
associate-/r*57.9%
Simplified57.9%
if 3e-68 < t < 7.0000000000000001e97Initial program 74.3%
Simplified79.3%
associate-*r*86.8%
*-un-lft-identity86.8%
times-frac86.9%
Applied egg-rr86.9%
/-rgt-identity86.9%
*-commutative86.9%
associate-/r*86.9%
Simplified86.9%
if 7.0000000000000001e97 < t Initial program 38.7%
associate-*l*38.7%
*-commutative38.7%
*-commutative38.7%
associate-/r*46.1%
distribute-rgt-in46.1%
unpow246.1%
times-frac33.9%
sqr-neg33.9%
times-frac46.1%
unpow246.1%
distribute-rgt-in46.1%
+-commutative46.1%
Simplified46.1%
add-cube-cbrt46.1%
pow346.1%
cbrt-div46.1%
rem-cbrt-cube74.4%
Applied egg-rr74.4%
Final simplification64.4%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (+ 2.0 (pow (/ k t_m) 2.0))))
(*
t_s
(if (<= t_m 4.6e-68)
(*
2.0
(/ (* (pow l 2.0) (cos k)) (* (pow k 2.0) (* t_m (pow (sin k) 2.0)))))
(if (<= t_m 7e+97)
(* (* l (/ (/ 2.0 (tan k)) (* (sin k) (pow t_m 3.0)))) (/ l t_2))
(/
2.0
(*
(* (tan k) t_2)
(* (sin k) (/ (* (pow t_m 2.0) (/ t_m l)) l)))))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = 2.0 + pow((k / t_m), 2.0);
double tmp;
if (t_m <= 4.6e-68) {
tmp = 2.0 * ((pow(l, 2.0) * cos(k)) / (pow(k, 2.0) * (t_m * pow(sin(k), 2.0))));
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / tan(k)) / (sin(k) * pow(t_m, 3.0)))) * (l / t_2);
} else {
tmp = 2.0 / ((tan(k) * t_2) * (sin(k) * ((pow(t_m, 2.0) * (t_m / l)) / l)));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: tmp
t_2 = 2.0d0 + ((k / t_m) ** 2.0d0)
if (t_m <= 4.6d-68) then
tmp = 2.0d0 * (((l ** 2.0d0) * cos(k)) / ((k ** 2.0d0) * (t_m * (sin(k) ** 2.0d0))))
else if (t_m <= 7d+97) then
tmp = (l * ((2.0d0 / tan(k)) / (sin(k) * (t_m ** 3.0d0)))) * (l / t_2)
else
tmp = 2.0d0 / ((tan(k) * t_2) * (sin(k) * (((t_m ** 2.0d0) * (t_m / l)) / l)))
end if
code = t_s * tmp
end function
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) {
double t_2 = 2.0 + Math.pow((k / t_m), 2.0);
double tmp;
if (t_m <= 4.6e-68) {
tmp = 2.0 * ((Math.pow(l, 2.0) * Math.cos(k)) / (Math.pow(k, 2.0) * (t_m * Math.pow(Math.sin(k), 2.0))));
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / Math.tan(k)) / (Math.sin(k) * Math.pow(t_m, 3.0)))) * (l / t_2);
} else {
tmp = 2.0 / ((Math.tan(k) * t_2) * (Math.sin(k) * ((Math.pow(t_m, 2.0) * (t_m / l)) / l)));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): t_2 = 2.0 + math.pow((k / t_m), 2.0) tmp = 0 if t_m <= 4.6e-68: tmp = 2.0 * ((math.pow(l, 2.0) * math.cos(k)) / (math.pow(k, 2.0) * (t_m * math.pow(math.sin(k), 2.0)))) elif t_m <= 7e+97: tmp = (l * ((2.0 / math.tan(k)) / (math.sin(k) * math.pow(t_m, 3.0)))) * (l / t_2) else: tmp = 2.0 / ((math.tan(k) * t_2) * (math.sin(k) * ((math.pow(t_m, 2.0) * (t_m / l)) / l))) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(2.0 + (Float64(k / t_m) ^ 2.0)) tmp = 0.0 if (t_m <= 4.6e-68) tmp = Float64(2.0 * Float64(Float64((l ^ 2.0) * cos(k)) / Float64((k ^ 2.0) * Float64(t_m * (sin(k) ^ 2.0))))); elseif (t_m <= 7e+97) tmp = Float64(Float64(l * Float64(Float64(2.0 / tan(k)) / Float64(sin(k) * (t_m ^ 3.0)))) * Float64(l / t_2)); else tmp = Float64(2.0 / Float64(Float64(tan(k) * t_2) * Float64(sin(k) * Float64(Float64((t_m ^ 2.0) * Float64(t_m / l)) / l)))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) t_2 = 2.0 + ((k / t_m) ^ 2.0); tmp = 0.0; if (t_m <= 4.6e-68) tmp = 2.0 * (((l ^ 2.0) * cos(k)) / ((k ^ 2.0) * (t_m * (sin(k) ^ 2.0)))); elseif (t_m <= 7e+97) tmp = (l * ((2.0 / tan(k)) / (sin(k) * (t_m ^ 3.0)))) * (l / t_2); else tmp = 2.0 / ((tan(k) * t_2) * (sin(k) * (((t_m ^ 2.0) * (t_m / l)) / l))); end tmp_2 = t_s * tmp; end
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_] := Block[{t$95$2 = N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 4.6e-68], N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision] / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 7e+97], N[(N[(l * N[(N[(2.0 / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / t$95$2), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * t$95$2), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[(N[Power[t$95$m, 2.0], $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := 2 + {\left(\frac{k}{t_m}\right)}^{2}\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 4.6 \cdot 10^{-68}:\\
\;\;\;\;2 \cdot \frac{{\ell}^{2} \cdot \cos k}{{k}^{2} \cdot \left(t_m \cdot {\sin k}^{2}\right)}\\
\mathbf{elif}\;t_m \leq 7 \cdot 10^{+97}:\\
\;\;\;\;\left(\ell \cdot \frac{\frac{2}{\tan k}}{\sin k \cdot {t_m}^{3}}\right) \cdot \frac{\ell}{t_2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot t_2\right) \cdot \left(\sin k \cdot \frac{{t_m}^{2} \cdot \frac{t_m}{\ell}}{\ell}\right)}\\
\end{array}
\end{array}
\end{array}
if t < 4.59999999999999994e-68Initial program 49.8%
Simplified49.2%
Taylor expanded in k around inf 56.8%
if 4.59999999999999994e-68 < t < 7.0000000000000001e97Initial program 74.3%
Simplified79.3%
associate-*r*86.8%
*-un-lft-identity86.8%
times-frac86.9%
Applied egg-rr86.9%
/-rgt-identity86.9%
*-commutative86.9%
associate-/r*86.9%
Simplified86.9%
if 7.0000000000000001e97 < t Initial program 38.7%
associate-*l*38.7%
*-commutative38.7%
*-commutative38.7%
associate-/r*46.1%
distribute-rgt-in46.1%
unpow246.1%
times-frac33.9%
sqr-neg33.9%
times-frac46.1%
unpow246.1%
distribute-rgt-in46.1%
+-commutative46.1%
Simplified46.1%
unpow346.1%
*-un-lft-identity46.1%
times-frac68.9%
pow268.9%
Applied egg-rr68.9%
Final simplification62.8%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (+ 2.0 (pow (/ k t_m) 2.0))))
(*
t_s
(if (<= t_m 1.36e-67)
(*
(/ 2.0 (pow k 2.0))
(* (/ (pow l 2.0) (/ t_m (cos k))) (pow (sin k) -2.0)))
(if (<= t_m 7e+97)
(* (* l (/ (/ 2.0 (tan k)) (* (sin k) (pow t_m 3.0)))) (/ l t_2))
(/
2.0
(*
(* (tan k) t_2)
(* (sin k) (/ (* (pow t_m 2.0) (/ t_m l)) l)))))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = 2.0 + pow((k / t_m), 2.0);
double tmp;
if (t_m <= 1.36e-67) {
tmp = (2.0 / pow(k, 2.0)) * ((pow(l, 2.0) / (t_m / cos(k))) * pow(sin(k), -2.0));
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / tan(k)) / (sin(k) * pow(t_m, 3.0)))) * (l / t_2);
} else {
tmp = 2.0 / ((tan(k) * t_2) * (sin(k) * ((pow(t_m, 2.0) * (t_m / l)) / l)));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: tmp
t_2 = 2.0d0 + ((k / t_m) ** 2.0d0)
if (t_m <= 1.36d-67) then
tmp = (2.0d0 / (k ** 2.0d0)) * (((l ** 2.0d0) / (t_m / cos(k))) * (sin(k) ** (-2.0d0)))
else if (t_m <= 7d+97) then
tmp = (l * ((2.0d0 / tan(k)) / (sin(k) * (t_m ** 3.0d0)))) * (l / t_2)
else
tmp = 2.0d0 / ((tan(k) * t_2) * (sin(k) * (((t_m ** 2.0d0) * (t_m / l)) / l)))
end if
code = t_s * tmp
end function
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) {
double t_2 = 2.0 + Math.pow((k / t_m), 2.0);
double tmp;
if (t_m <= 1.36e-67) {
tmp = (2.0 / Math.pow(k, 2.0)) * ((Math.pow(l, 2.0) / (t_m / Math.cos(k))) * Math.pow(Math.sin(k), -2.0));
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / Math.tan(k)) / (Math.sin(k) * Math.pow(t_m, 3.0)))) * (l / t_2);
} else {
tmp = 2.0 / ((Math.tan(k) * t_2) * (Math.sin(k) * ((Math.pow(t_m, 2.0) * (t_m / l)) / l)));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): t_2 = 2.0 + math.pow((k / t_m), 2.0) tmp = 0 if t_m <= 1.36e-67: tmp = (2.0 / math.pow(k, 2.0)) * ((math.pow(l, 2.0) / (t_m / math.cos(k))) * math.pow(math.sin(k), -2.0)) elif t_m <= 7e+97: tmp = (l * ((2.0 / math.tan(k)) / (math.sin(k) * math.pow(t_m, 3.0)))) * (l / t_2) else: tmp = 2.0 / ((math.tan(k) * t_2) * (math.sin(k) * ((math.pow(t_m, 2.0) * (t_m / l)) / l))) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(2.0 + (Float64(k / t_m) ^ 2.0)) tmp = 0.0 if (t_m <= 1.36e-67) tmp = Float64(Float64(2.0 / (k ^ 2.0)) * Float64(Float64((l ^ 2.0) / Float64(t_m / cos(k))) * (sin(k) ^ -2.0))); elseif (t_m <= 7e+97) tmp = Float64(Float64(l * Float64(Float64(2.0 / tan(k)) / Float64(sin(k) * (t_m ^ 3.0)))) * Float64(l / t_2)); else tmp = Float64(2.0 / Float64(Float64(tan(k) * t_2) * Float64(sin(k) * Float64(Float64((t_m ^ 2.0) * Float64(t_m / l)) / l)))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) t_2 = 2.0 + ((k / t_m) ^ 2.0); tmp = 0.0; if (t_m <= 1.36e-67) tmp = (2.0 / (k ^ 2.0)) * (((l ^ 2.0) / (t_m / cos(k))) * (sin(k) ^ -2.0)); elseif (t_m <= 7e+97) tmp = (l * ((2.0 / tan(k)) / (sin(k) * (t_m ^ 3.0)))) * (l / t_2); else tmp = 2.0 / ((tan(k) * t_2) * (sin(k) * (((t_m ^ 2.0) * (t_m / l)) / l))); end tmp_2 = t_s * tmp; end
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_] := Block[{t$95$2 = N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 1.36e-67], N[(N[(2.0 / N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[(t$95$m / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Power[N[Sin[k], $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 7e+97], N[(N[(l * N[(N[(2.0 / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / t$95$2), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * t$95$2), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[(N[Power[t$95$m, 2.0], $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := 2 + {\left(\frac{k}{t_m}\right)}^{2}\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 1.36 \cdot 10^{-67}:\\
\;\;\;\;\frac{2}{{k}^{2}} \cdot \left(\frac{{\ell}^{2}}{\frac{t_m}{\cos k}} \cdot {\sin k}^{-2}\right)\\
\mathbf{elif}\;t_m \leq 7 \cdot 10^{+97}:\\
\;\;\;\;\left(\ell \cdot \frac{\frac{2}{\tan k}}{\sin k \cdot {t_m}^{3}}\right) \cdot \frac{\ell}{t_2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot t_2\right) \cdot \left(\sin k \cdot \frac{{t_m}^{2} \cdot \frac{t_m}{\ell}}{\ell}\right)}\\
\end{array}
\end{array}
\end{array}
if t < 1.36000000000000007e-67Initial program 49.8%
Simplified49.2%
Taylor expanded in k around inf 56.8%
associate-*r/56.8%
associate-*r*56.7%
associate-/r*57.9%
Simplified57.9%
expm1-log1p-u41.4%
expm1-udef38.1%
div-inv38.1%
times-frac38.3%
pow-flip38.3%
metadata-eval38.3%
Applied egg-rr38.3%
expm1-def41.2%
expm1-log1p57.2%
associate-*l*57.3%
associate-/l*57.3%
Simplified57.3%
if 1.36000000000000007e-67 < t < 7.0000000000000001e97Initial program 74.3%
Simplified79.3%
associate-*r*86.8%
*-un-lft-identity86.8%
times-frac86.9%
Applied egg-rr86.9%
/-rgt-identity86.9%
*-commutative86.9%
associate-/r*86.9%
Simplified86.9%
if 7.0000000000000001e97 < t Initial program 38.7%
associate-*l*38.7%
*-commutative38.7%
*-commutative38.7%
associate-/r*46.1%
distribute-rgt-in46.1%
unpow246.1%
times-frac33.9%
sqr-neg33.9%
times-frac46.1%
unpow246.1%
distribute-rgt-in46.1%
+-commutative46.1%
Simplified46.1%
unpow346.1%
*-un-lft-identity46.1%
times-frac68.9%
pow268.9%
Applied egg-rr68.9%
Final simplification63.2%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (+ 2.0 (pow (/ k t_m) 2.0))))
(*
t_s
(if (<= t_m 1.05e-67)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(pow (sin k) 2.0))
(if (<= t_m 7e+97)
(* (* l (/ (/ 2.0 (tan k)) (* (sin k) (pow t_m 3.0)))) (/ l t_2))
(/
2.0
(*
(* (tan k) t_2)
(* (sin k) (/ (* (pow t_m 2.0) (/ t_m l)) l)))))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double t_2 = 2.0 + pow((k / t_m), 2.0);
double tmp;
if (t_m <= 1.05e-67) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / pow(sin(k), 2.0);
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / tan(k)) / (sin(k) * pow(t_m, 3.0)))) * (l / t_2);
} else {
tmp = 2.0 / ((tan(k) * t_2) * (sin(k) * ((pow(t_m, 2.0) * (t_m / l)) / l)));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_2
real(8) :: tmp
t_2 = 2.0d0 + ((k / t_m) ** 2.0d0)
if (t_m <= 1.05d-67) then
tmp = ((2.0d0 * ((l ** 2.0d0) * cos(k))) / (t_m * (k ** 2.0d0))) / (sin(k) ** 2.0d0)
else if (t_m <= 7d+97) then
tmp = (l * ((2.0d0 / tan(k)) / (sin(k) * (t_m ** 3.0d0)))) * (l / t_2)
else
tmp = 2.0d0 / ((tan(k) * t_2) * (sin(k) * (((t_m ** 2.0d0) * (t_m / l)) / l)))
end if
code = t_s * tmp
end function
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) {
double t_2 = 2.0 + Math.pow((k / t_m), 2.0);
double tmp;
if (t_m <= 1.05e-67) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / Math.pow(Math.sin(k), 2.0);
} else if (t_m <= 7e+97) {
tmp = (l * ((2.0 / Math.tan(k)) / (Math.sin(k) * Math.pow(t_m, 3.0)))) * (l / t_2);
} else {
tmp = 2.0 / ((Math.tan(k) * t_2) * (Math.sin(k) * ((Math.pow(t_m, 2.0) * (t_m / l)) / l)));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): t_2 = 2.0 + math.pow((k / t_m), 2.0) tmp = 0 if t_m <= 1.05e-67: tmp = ((2.0 * (math.pow(l, 2.0) * math.cos(k))) / (t_m * math.pow(k, 2.0))) / math.pow(math.sin(k), 2.0) elif t_m <= 7e+97: tmp = (l * ((2.0 / math.tan(k)) / (math.sin(k) * math.pow(t_m, 3.0)))) * (l / t_2) else: tmp = 2.0 / ((math.tan(k) * t_2) * (math.sin(k) * ((math.pow(t_m, 2.0) * (t_m / l)) / l))) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) t_2 = Float64(2.0 + (Float64(k / t_m) ^ 2.0)) tmp = 0.0 if (t_m <= 1.05e-67) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / (sin(k) ^ 2.0)); elseif (t_m <= 7e+97) tmp = Float64(Float64(l * Float64(Float64(2.0 / tan(k)) / Float64(sin(k) * (t_m ^ 3.0)))) * Float64(l / t_2)); else tmp = Float64(2.0 / Float64(Float64(tan(k) * t_2) * Float64(sin(k) * Float64(Float64((t_m ^ 2.0) * Float64(t_m / l)) / l)))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) t_2 = 2.0 + ((k / t_m) ^ 2.0); tmp = 0.0; if (t_m <= 1.05e-67) tmp = ((2.0 * ((l ^ 2.0) * cos(k))) / (t_m * (k ^ 2.0))) / (sin(k) ^ 2.0); elseif (t_m <= 7e+97) tmp = (l * ((2.0 / tan(k)) / (sin(k) * (t_m ^ 3.0)))) * (l / t_2); else tmp = 2.0 / ((tan(k) * t_2) * (sin(k) * (((t_m ^ 2.0) * (t_m / l)) / l))); end tmp_2 = t_s * tmp; end
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_] := Block[{t$95$2 = N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 1.05e-67], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 7e+97], N[(N[(l * N[(N[(2.0 / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / t$95$2), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * t$95$2), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[(N[Power[t$95$m, 2.0], $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := 2 + {\left(\frac{k}{t_m}\right)}^{2}\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 1.05 \cdot 10^{-67}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{{\sin k}^{2}}\\
\mathbf{elif}\;t_m \leq 7 \cdot 10^{+97}:\\
\;\;\;\;\left(\ell \cdot \frac{\frac{2}{\tan k}}{\sin k \cdot {t_m}^{3}}\right) \cdot \frac{\ell}{t_2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot t_2\right) \cdot \left(\sin k \cdot \frac{{t_m}^{2} \cdot \frac{t_m}{\ell}}{\ell}\right)}\\
\end{array}
\end{array}
\end{array}
if t < 1.0500000000000001e-67Initial program 49.8%
Simplified49.2%
Taylor expanded in k around inf 56.8%
associate-*r/56.8%
associate-*r*56.7%
associate-/r*57.9%
Simplified57.9%
if 1.0500000000000001e-67 < t < 7.0000000000000001e97Initial program 74.3%
Simplified79.3%
associate-*r*86.8%
*-un-lft-identity86.8%
times-frac86.9%
Applied egg-rr86.9%
/-rgt-identity86.9%
*-commutative86.9%
associate-/r*86.9%
Simplified86.9%
if 7.0000000000000001e97 < t Initial program 38.7%
associate-*l*38.7%
*-commutative38.7%
*-commutative38.7%
associate-/r*46.1%
distribute-rgt-in46.1%
unpow246.1%
times-frac33.9%
sqr-neg33.9%
times-frac46.1%
unpow246.1%
distribute-rgt-in46.1%
+-commutative46.1%
Simplified46.1%
unpow346.1%
*-un-lft-identity46.1%
times-frac68.9%
pow268.9%
Applied egg-rr68.9%
Final simplification63.7%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.35e-133)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(- 0.5 (/ (cos (* 2.0 k)) 2.0)))
(/
2.0
(*
(+ 1.0 (+ (pow (/ k t_m) 2.0) 1.0))
(* (tan k) (* (sin k) (* (/ t_m l) (/ (pow t_m 2.0) l)))))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.35e-133) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / (0.5 - (cos((2.0 * k)) / 2.0));
} else {
tmp = 2.0 / ((1.0 + (pow((k / t_m), 2.0) + 1.0)) * (tan(k) * (sin(k) * ((t_m / l) * (pow(t_m, 2.0) / l)))));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1.35d-133) then
tmp = ((2.0d0 * ((l ** 2.0d0) * cos(k))) / (t_m * (k ** 2.0d0))) / (0.5d0 - (cos((2.0d0 * k)) / 2.0d0))
else
tmp = 2.0d0 / ((1.0d0 + (((k / t_m) ** 2.0d0) + 1.0d0)) * (tan(k) * (sin(k) * ((t_m / l) * ((t_m ** 2.0d0) / l)))))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 1.35e-133) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / (0.5 - (Math.cos((2.0 * k)) / 2.0));
} else {
tmp = 2.0 / ((1.0 + (Math.pow((k / t_m), 2.0) + 1.0)) * (Math.tan(k) * (Math.sin(k) * ((t_m / l) * (Math.pow(t_m, 2.0) / l)))));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1.35e-133: tmp = ((2.0 * (math.pow(l, 2.0) * math.cos(k))) / (t_m * math.pow(k, 2.0))) / (0.5 - (math.cos((2.0 * k)) / 2.0)) else: tmp = 2.0 / ((1.0 + (math.pow((k / t_m), 2.0) + 1.0)) * (math.tan(k) * (math.sin(k) * ((t_m / l) * (math.pow(t_m, 2.0) / l))))) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.35e-133) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / Float64(0.5 - Float64(cos(Float64(2.0 * k)) / 2.0))); else tmp = Float64(2.0 / Float64(Float64(1.0 + Float64((Float64(k / t_m) ^ 2.0) + 1.0)) * Float64(tan(k) * Float64(sin(k) * Float64(Float64(t_m / l) * Float64((t_m ^ 2.0) / l)))))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1.35e-133) tmp = ((2.0 * ((l ^ 2.0) * cos(k))) / (t_m * (k ^ 2.0))) / (0.5 - (cos((2.0 * k)) / 2.0)); else tmp = 2.0 / ((1.0 + (((k / t_m) ^ 2.0) + 1.0)) * (tan(k) * (sin(k) * ((t_m / l) * ((t_m ^ 2.0) / l))))); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 1.35e-133], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.5 - N[(N[Cos[N[(2.0 * k), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(1.0 + N[(N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[(t$95$m / l), $MachinePrecision] * N[(N[Power[t$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 1.35 \cdot 10^{-133}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{0.5 - \frac{\cos \left(2 \cdot k\right)}{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(1 + \left({\left(\frac{k}{t_m}\right)}^{2} + 1\right)\right) \cdot \left(\tan k \cdot \left(\sin k \cdot \left(\frac{t_m}{\ell} \cdot \frac{{t_m}^{2}}{\ell}\right)\right)\right)}\\
\end{array}
\end{array}
if t < 1.3499999999999999e-133Initial program 48.5%
Simplified48.5%
Taylor expanded in k around inf 55.5%
associate-*r/55.5%
associate-*r*55.5%
associate-/r*55.6%
Simplified55.6%
unpow255.6%
sin-mult52.3%
Applied egg-rr52.3%
div-sub52.3%
+-inverses52.3%
cos-052.3%
metadata-eval52.3%
count-252.3%
*-commutative52.3%
Simplified52.3%
if 1.3499999999999999e-133 < t Initial program 58.8%
unpow358.7%
times-frac76.3%
pow276.3%
Applied egg-rr76.3%
Final simplification60.3%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 6.8e-134)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(- 0.5 (/ (cos (* 2.0 k)) 2.0)))
(/
2.0
(*
(* (tan k) (+ 2.0 (pow (/ k t_m) 2.0)))
(* (sin k) (/ (* (pow t_m 2.0) (/ t_m l)) l)))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 6.8e-134) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / (0.5 - (cos((2.0 * k)) / 2.0));
} else {
tmp = 2.0 / ((tan(k) * (2.0 + pow((k / t_m), 2.0))) * (sin(k) * ((pow(t_m, 2.0) * (t_m / l)) / l)));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 6.8d-134) then
tmp = ((2.0d0 * ((l ** 2.0d0) * cos(k))) / (t_m * (k ** 2.0d0))) / (0.5d0 - (cos((2.0d0 * k)) / 2.0d0))
else
tmp = 2.0d0 / ((tan(k) * (2.0d0 + ((k / t_m) ** 2.0d0))) * (sin(k) * (((t_m ** 2.0d0) * (t_m / l)) / l)))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 6.8e-134) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / (0.5 - (Math.cos((2.0 * k)) / 2.0));
} else {
tmp = 2.0 / ((Math.tan(k) * (2.0 + Math.pow((k / t_m), 2.0))) * (Math.sin(k) * ((Math.pow(t_m, 2.0) * (t_m / l)) / l)));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 6.8e-134: tmp = ((2.0 * (math.pow(l, 2.0) * math.cos(k))) / (t_m * math.pow(k, 2.0))) / (0.5 - (math.cos((2.0 * k)) / 2.0)) else: tmp = 2.0 / ((math.tan(k) * (2.0 + math.pow((k / t_m), 2.0))) * (math.sin(k) * ((math.pow(t_m, 2.0) * (t_m / l)) / l))) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 6.8e-134) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / Float64(0.5 - Float64(cos(Float64(2.0 * k)) / 2.0))); else tmp = Float64(2.0 / Float64(Float64(tan(k) * Float64(2.0 + (Float64(k / t_m) ^ 2.0))) * Float64(sin(k) * Float64(Float64((t_m ^ 2.0) * Float64(t_m / l)) / l)))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 6.8e-134) tmp = ((2.0 * ((l ^ 2.0) * cos(k))) / (t_m * (k ^ 2.0))) / (0.5 - (cos((2.0 * k)) / 2.0)); else tmp = 2.0 / ((tan(k) * (2.0 + ((k / t_m) ^ 2.0))) * (sin(k) * (((t_m ^ 2.0) * (t_m / l)) / l))); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 6.8e-134], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.5 - N[(N[Cos[N[(2.0 * k), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[(N[Power[t$95$m, 2.0], $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
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^{-134}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{0.5 - \frac{\cos \left(2 \cdot k\right)}{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot \left(2 + {\left(\frac{k}{t_m}\right)}^{2}\right)\right) \cdot \left(\sin k \cdot \frac{{t_m}^{2} \cdot \frac{t_m}{\ell}}{\ell}\right)}\\
\end{array}
\end{array}
if t < 6.79999999999999954e-134Initial program 48.5%
Simplified48.5%
Taylor expanded in k around inf 55.5%
associate-*r/55.5%
associate-*r*55.5%
associate-/r*55.6%
Simplified55.6%
unpow255.6%
sin-mult52.3%
Applied egg-rr52.3%
div-sub52.3%
+-inverses52.3%
cos-052.3%
metadata-eval52.3%
count-252.3%
*-commutative52.3%
Simplified52.3%
if 6.79999999999999954e-134 < t Initial program 58.8%
associate-*l*58.8%
*-commutative58.8%
*-commutative58.8%
associate-/r*63.9%
distribute-rgt-in63.9%
unpow263.9%
times-frac59.1%
sqr-neg59.1%
times-frac63.9%
unpow263.9%
distribute-rgt-in63.9%
+-commutative63.9%
Simplified63.9%
unpow363.9%
*-un-lft-identity63.9%
times-frac76.3%
pow276.3%
Applied egg-rr76.3%
Final simplification60.3%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.95e-69)
(/ (* 2.0 (/ (/ (pow l 2.0) (pow k 2.0)) t_m)) (pow (sin k) 2.0))
(*
(* l (/ (/ 2.0 (tan k)) (* (sin k) (pow t_m 3.0))))
(/ l (+ 2.0 (pow (/ k t_m) 2.0)))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.95e-69) {
tmp = (2.0 * ((pow(l, 2.0) / pow(k, 2.0)) / t_m)) / pow(sin(k), 2.0);
} else {
tmp = (l * ((2.0 / tan(k)) / (sin(k) * pow(t_m, 3.0)))) * (l / (2.0 + pow((k / t_m), 2.0)));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1.95d-69) then
tmp = (2.0d0 * (((l ** 2.0d0) / (k ** 2.0d0)) / t_m)) / (sin(k) ** 2.0d0)
else
tmp = (l * ((2.0d0 / tan(k)) / (sin(k) * (t_m ** 3.0d0)))) * (l / (2.0d0 + ((k / t_m) ** 2.0d0)))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 1.95e-69) {
tmp = (2.0 * ((Math.pow(l, 2.0) / Math.pow(k, 2.0)) / t_m)) / Math.pow(Math.sin(k), 2.0);
} else {
tmp = (l * ((2.0 / Math.tan(k)) / (Math.sin(k) * Math.pow(t_m, 3.0)))) * (l / (2.0 + Math.pow((k / t_m), 2.0)));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1.95e-69: tmp = (2.0 * ((math.pow(l, 2.0) / math.pow(k, 2.0)) / t_m)) / math.pow(math.sin(k), 2.0) else: tmp = (l * ((2.0 / math.tan(k)) / (math.sin(k) * math.pow(t_m, 3.0)))) * (l / (2.0 + math.pow((k / t_m), 2.0))) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.95e-69) tmp = Float64(Float64(2.0 * Float64(Float64((l ^ 2.0) / (k ^ 2.0)) / t_m)) / (sin(k) ^ 2.0)); else tmp = Float64(Float64(l * Float64(Float64(2.0 / tan(k)) / Float64(sin(k) * (t_m ^ 3.0)))) * Float64(l / Float64(2.0 + (Float64(k / t_m) ^ 2.0)))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1.95e-69) tmp = (2.0 * (((l ^ 2.0) / (k ^ 2.0)) / t_m)) / (sin(k) ^ 2.0); else tmp = (l * ((2.0 / tan(k)) / (sin(k) * (t_m ^ 3.0)))) * (l / (2.0 + ((k / t_m) ^ 2.0))); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 1.95e-69], N[(N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(l * N[(N[(2.0 / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 1.95 \cdot 10^{-69}:\\
\;\;\;\;\frac{2 \cdot \frac{\frac{{\ell}^{2}}{{k}^{2}}}{t_m}}{{\sin k}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \frac{\frac{2}{\tan k}}{\sin k \cdot {t_m}^{3}}\right) \cdot \frac{\ell}{2 + {\left(\frac{k}{t_m}\right)}^{2}}\\
\end{array}
\end{array}
if t < 1.9499999999999999e-69Initial program 49.8%
Simplified49.2%
Taylor expanded in k around inf 56.8%
associate-*r/56.8%
associate-*r*56.7%
associate-/r*57.9%
Simplified57.9%
Taylor expanded in k around 0 47.8%
associate-/r*47.2%
Simplified47.2%
if 1.9499999999999999e-69 < t Initial program 57.5%
Simplified60.2%
associate-*r*67.6%
*-un-lft-identity67.6%
times-frac67.7%
Applied egg-rr67.7%
/-rgt-identity67.7%
*-commutative67.7%
associate-/r*67.7%
Simplified67.7%
Final simplification52.9%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 2.55e-70)
(/
(/ (* 2.0 (* (pow l 2.0) (cos k))) (* t_m (pow k 2.0)))
(- 0.5 (/ (cos (* 2.0 k)) 2.0)))
(*
(* l (/ (/ 2.0 (tan k)) (* (sin k) (pow t_m 3.0))))
(/ l (+ 2.0 (pow (/ k t_m) 2.0)))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 2.55e-70) {
tmp = ((2.0 * (pow(l, 2.0) * cos(k))) / (t_m * pow(k, 2.0))) / (0.5 - (cos((2.0 * k)) / 2.0));
} else {
tmp = (l * ((2.0 / tan(k)) / (sin(k) * pow(t_m, 3.0)))) * (l / (2.0 + pow((k / t_m), 2.0)));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 2.55d-70) then
tmp = ((2.0d0 * ((l ** 2.0d0) * cos(k))) / (t_m * (k ** 2.0d0))) / (0.5d0 - (cos((2.0d0 * k)) / 2.0d0))
else
tmp = (l * ((2.0d0 / tan(k)) / (sin(k) * (t_m ** 3.0d0)))) * (l / (2.0d0 + ((k / t_m) ** 2.0d0)))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 2.55e-70) {
tmp = ((2.0 * (Math.pow(l, 2.0) * Math.cos(k))) / (t_m * Math.pow(k, 2.0))) / (0.5 - (Math.cos((2.0 * k)) / 2.0));
} else {
tmp = (l * ((2.0 / Math.tan(k)) / (Math.sin(k) * Math.pow(t_m, 3.0)))) * (l / (2.0 + Math.pow((k / t_m), 2.0)));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 2.55e-70: tmp = ((2.0 * (math.pow(l, 2.0) * math.cos(k))) / (t_m * math.pow(k, 2.0))) / (0.5 - (math.cos((2.0 * k)) / 2.0)) else: tmp = (l * ((2.0 / math.tan(k)) / (math.sin(k) * math.pow(t_m, 3.0)))) * (l / (2.0 + math.pow((k / t_m), 2.0))) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 2.55e-70) tmp = Float64(Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k))) / Float64(t_m * (k ^ 2.0))) / Float64(0.5 - Float64(cos(Float64(2.0 * k)) / 2.0))); else tmp = Float64(Float64(l * Float64(Float64(2.0 / tan(k)) / Float64(sin(k) * (t_m ^ 3.0)))) * Float64(l / Float64(2.0 + (Float64(k / t_m) ^ 2.0)))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 2.55e-70) tmp = ((2.0 * ((l ^ 2.0) * cos(k))) / (t_m * (k ^ 2.0))) / (0.5 - (cos((2.0 * k)) / 2.0)); else tmp = (l * ((2.0 / tan(k)) / (sin(k) * (t_m ^ 3.0)))) * (l / (2.0 + ((k / t_m) ^ 2.0))); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 2.55e-70], N[(N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.5 - N[(N[Cos[N[(2.0 * k), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * N[(N[(2.0 / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
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.55 \cdot 10^{-70}:\\
\;\;\;\;\frac{\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\right)}{t_m \cdot {k}^{2}}}{0.5 - \frac{\cos \left(2 \cdot k\right)}{2}}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \frac{\frac{2}{\tan k}}{\sin k \cdot {t_m}^{3}}\right) \cdot \frac{\ell}{2 + {\left(\frac{k}{t_m}\right)}^{2}}\\
\end{array}
\end{array}
if t < 2.55000000000000013e-70Initial program 50.0%
Simplified49.4%
Taylor expanded in k around inf 56.5%
associate-*r/56.5%
associate-*r*56.5%
associate-/r*57.7%
Simplified57.7%
unpow257.7%
sin-mult53.5%
Applied egg-rr53.5%
div-sub53.5%
+-inverses53.5%
cos-053.5%
metadata-eval53.5%
count-253.5%
*-commutative53.5%
Simplified53.5%
if 2.55000000000000013e-70 < t Initial program 56.9%
Simplified59.5%
associate-*r*66.8%
*-un-lft-identity66.8%
times-frac66.9%
Applied egg-rr66.9%
/-rgt-identity66.9%
*-commutative66.9%
associate-/r*66.9%
Simplified66.9%
Final simplification57.3%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 4.8e-19)
(/ (* 2.0 (/ (/ (pow l 2.0) (pow k 2.0)) t_m)) (pow (sin k) 2.0))
(/
2.0
(* 2.0 (/ k (/ (* (pow l 2.0) (cos k)) (* (sin k) (pow t_m 3.0)))))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 4.8e-19) {
tmp = (2.0 * ((pow(l, 2.0) / pow(k, 2.0)) / t_m)) / pow(sin(k), 2.0);
} else {
tmp = 2.0 / (2.0 * (k / ((pow(l, 2.0) * cos(k)) / (sin(k) * pow(t_m, 3.0)))));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 4.8d-19) then
tmp = (2.0d0 * (((l ** 2.0d0) / (k ** 2.0d0)) / t_m)) / (sin(k) ** 2.0d0)
else
tmp = 2.0d0 / (2.0d0 * (k / (((l ** 2.0d0) * cos(k)) / (sin(k) * (t_m ** 3.0d0)))))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 4.8e-19) {
tmp = (2.0 * ((Math.pow(l, 2.0) / Math.pow(k, 2.0)) / t_m)) / Math.pow(Math.sin(k), 2.0);
} else {
tmp = 2.0 / (2.0 * (k / ((Math.pow(l, 2.0) * Math.cos(k)) / (Math.sin(k) * Math.pow(t_m, 3.0)))));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 4.8e-19: tmp = (2.0 * ((math.pow(l, 2.0) / math.pow(k, 2.0)) / t_m)) / math.pow(math.sin(k), 2.0) else: tmp = 2.0 / (2.0 * (k / ((math.pow(l, 2.0) * math.cos(k)) / (math.sin(k) * math.pow(t_m, 3.0))))) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 4.8e-19) tmp = Float64(Float64(2.0 * Float64(Float64((l ^ 2.0) / (k ^ 2.0)) / t_m)) / (sin(k) ^ 2.0)); else tmp = Float64(2.0 / Float64(2.0 * Float64(k / Float64(Float64((l ^ 2.0) * cos(k)) / Float64(sin(k) * (t_m ^ 3.0)))))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 4.8e-19) tmp = (2.0 * (((l ^ 2.0) / (k ^ 2.0)) / t_m)) / (sin(k) ^ 2.0); else tmp = 2.0 / (2.0 * (k / (((l ^ 2.0) * cos(k)) / (sin(k) * (t_m ^ 3.0))))); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 4.8e-19], N[(N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(2.0 * N[(k / N[(N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k], $MachinePrecision]), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
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.8 \cdot 10^{-19}:\\
\;\;\;\;\frac{2 \cdot \frac{\frac{{\ell}^{2}}{{k}^{2}}}{t_m}}{{\sin k}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{k}{\frac{{\ell}^{2} \cdot \cos k}{\sin k \cdot {t_m}^{3}}}}\\
\end{array}
\end{array}
if t < 4.80000000000000046e-19Initial program 52.3%
Simplified51.8%
Taylor expanded in k around inf 58.5%
associate-*r/58.5%
associate-*r*58.5%
associate-/r*59.6%
Simplified59.6%
Taylor expanded in k around 0 49.2%
associate-/r*48.6%
Simplified48.6%
if 4.80000000000000046e-19 < t Initial program 50.7%
Taylor expanded in k around 0 51.4%
Taylor expanded in t around inf 49.2%
associate-/l*51.4%
*-commutative51.4%
Simplified51.4%
Final simplification49.3%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.42e-27)
(/ (* 2.0 (/ (/ (pow l 2.0) (pow k 2.0)) t_m)) (pow (sin k) 2.0))
(/
(* (/ 2.0 (* (tan k) (* k (pow t_m 3.0)))) (* l l))
(+ 2.0 (pow (/ k t_m) 2.0))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.42e-27) {
tmp = (2.0 * ((pow(l, 2.0) / pow(k, 2.0)) / t_m)) / pow(sin(k), 2.0);
} else {
tmp = ((2.0 / (tan(k) * (k * pow(t_m, 3.0)))) * (l * l)) / (2.0 + pow((k / t_m), 2.0));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1.42d-27) then
tmp = (2.0d0 * (((l ** 2.0d0) / (k ** 2.0d0)) / t_m)) / (sin(k) ** 2.0d0)
else
tmp = ((2.0d0 / (tan(k) * (k * (t_m ** 3.0d0)))) * (l * l)) / (2.0d0 + ((k / t_m) ** 2.0d0))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 1.42e-27) {
tmp = (2.0 * ((Math.pow(l, 2.0) / Math.pow(k, 2.0)) / t_m)) / Math.pow(Math.sin(k), 2.0);
} else {
tmp = ((2.0 / (Math.tan(k) * (k * Math.pow(t_m, 3.0)))) * (l * l)) / (2.0 + Math.pow((k / t_m), 2.0));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1.42e-27: tmp = (2.0 * ((math.pow(l, 2.0) / math.pow(k, 2.0)) / t_m)) / math.pow(math.sin(k), 2.0) else: tmp = ((2.0 / (math.tan(k) * (k * math.pow(t_m, 3.0)))) * (l * l)) / (2.0 + math.pow((k / t_m), 2.0)) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.42e-27) tmp = Float64(Float64(2.0 * Float64(Float64((l ^ 2.0) / (k ^ 2.0)) / t_m)) / (sin(k) ^ 2.0)); else tmp = Float64(Float64(Float64(2.0 / Float64(tan(k) * Float64(k * (t_m ^ 3.0)))) * Float64(l * l)) / Float64(2.0 + (Float64(k / t_m) ^ 2.0))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1.42e-27) tmp = (2.0 * (((l ^ 2.0) / (k ^ 2.0)) / t_m)) / (sin(k) ^ 2.0); else tmp = ((2.0 / (tan(k) * (k * (t_m ^ 3.0)))) * (l * l)) / (2.0 + ((k / t_m) ^ 2.0)); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 1.42e-27], N[(N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / N[(N[Tan[k], $MachinePrecision] * N[(k * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l * l), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 1.42 \cdot 10^{-27}:\\
\;\;\;\;\frac{2 \cdot \frac{\frac{{\ell}^{2}}{{k}^{2}}}{t_m}}{{\sin k}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\tan k \cdot \left(k \cdot {t_m}^{3}\right)} \cdot \left(\ell \cdot \ell\right)}{2 + {\left(\frac{k}{t_m}\right)}^{2}}\\
\end{array}
\end{array}
if t < 1.41999999999999996e-27Initial program 52.1%
Simplified51.6%
Taylor expanded in k around inf 58.3%
associate-*r/58.3%
associate-*r*58.3%
associate-/r*59.4%
Simplified59.4%
Taylor expanded in k around 0 48.9%
associate-/r*48.3%
Simplified48.3%
if 1.41999999999999996e-27 < t Initial program 51.5%
Simplified54.5%
Taylor expanded in k around 0 50.0%
*-commutative50.0%
Simplified50.0%
Final simplification48.7%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 4.8e+18)
(* 2.0 (* (/ (pow l 2.0) (pow k 4.0)) (/ (cos k) t_m)))
(/
(* (/ 2.0 (* (tan k) (* k (pow t_m 3.0)))) (* l l))
(+ 2.0 (pow (/ k t_m) 2.0))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 4.8e+18) {
tmp = 2.0 * ((pow(l, 2.0) / pow(k, 4.0)) * (cos(k) / t_m));
} else {
tmp = ((2.0 / (tan(k) * (k * pow(t_m, 3.0)))) * (l * l)) / (2.0 + pow((k / t_m), 2.0));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 4.8d+18) then
tmp = 2.0d0 * (((l ** 2.0d0) / (k ** 4.0d0)) * (cos(k) / t_m))
else
tmp = ((2.0d0 / (tan(k) * (k * (t_m ** 3.0d0)))) * (l * l)) / (2.0d0 + ((k / t_m) ** 2.0d0))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 4.8e+18) {
tmp = 2.0 * ((Math.pow(l, 2.0) / Math.pow(k, 4.0)) * (Math.cos(k) / t_m));
} else {
tmp = ((2.0 / (Math.tan(k) * (k * Math.pow(t_m, 3.0)))) * (l * l)) / (2.0 + Math.pow((k / t_m), 2.0));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 4.8e+18: tmp = 2.0 * ((math.pow(l, 2.0) / math.pow(k, 4.0)) * (math.cos(k) / t_m)) else: tmp = ((2.0 / (math.tan(k) * (k * math.pow(t_m, 3.0)))) * (l * l)) / (2.0 + math.pow((k / t_m), 2.0)) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 4.8e+18) tmp = Float64(2.0 * Float64(Float64((l ^ 2.0) / (k ^ 4.0)) * Float64(cos(k) / t_m))); else tmp = Float64(Float64(Float64(2.0 / Float64(tan(k) * Float64(k * (t_m ^ 3.0)))) * Float64(l * l)) / Float64(2.0 + (Float64(k / t_m) ^ 2.0))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 4.8e+18) tmp = 2.0 * (((l ^ 2.0) / (k ^ 4.0)) * (cos(k) / t_m)); else tmp = ((2.0 / (tan(k) * (k * (t_m ^ 3.0)))) * (l * l)) / (2.0 + ((k / t_m) ^ 2.0)); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 4.8e+18], N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k], $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / N[(N[Tan[k], $MachinePrecision] * N[(k * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l * l), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
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.8 \cdot 10^{+18}:\\
\;\;\;\;2 \cdot \left(\frac{{\ell}^{2}}{{k}^{4}} \cdot \frac{\cos k}{t_m}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\tan k \cdot \left(k \cdot {t_m}^{3}\right)} \cdot \left(\ell \cdot \ell\right)}{2 + {\left(\frac{k}{t_m}\right)}^{2}}\\
\end{array}
\end{array}
if t < 4.8e18Initial program 53.3%
Simplified52.3%
Taylor expanded in k around inf 58.8%
associate-*r/58.8%
associate-*r*58.7%
associate-/r*59.8%
Simplified59.8%
Taylor expanded in k around 0 50.7%
Taylor expanded in l around 0 49.6%
times-frac51.0%
Simplified51.0%
if 4.8e18 < t Initial program 47.3%
Simplified52.3%
Taylor expanded in k around 0 47.2%
*-commutative47.2%
Simplified47.2%
Final simplification50.2%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.7e+38)
(* 2.0 (* (/ (pow l 2.0) t_m) (/ (cos k) (pow k 4.0))))
(* (pow l 2.0) (* (pow t_m -3.0) (pow k -2.0))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.7e+38) {
tmp = 2.0 * ((pow(l, 2.0) / t_m) * (cos(k) / pow(k, 4.0)));
} else {
tmp = pow(l, 2.0) * (pow(t_m, -3.0) * pow(k, -2.0));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1.7d+38) then
tmp = 2.0d0 * (((l ** 2.0d0) / t_m) * (cos(k) / (k ** 4.0d0)))
else
tmp = (l ** 2.0d0) * ((t_m ** (-3.0d0)) * (k ** (-2.0d0)))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 1.7e+38) {
tmp = 2.0 * ((Math.pow(l, 2.0) / t_m) * (Math.cos(k) / Math.pow(k, 4.0)));
} else {
tmp = Math.pow(l, 2.0) * (Math.pow(t_m, -3.0) * Math.pow(k, -2.0));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1.7e+38: tmp = 2.0 * ((math.pow(l, 2.0) / t_m) * (math.cos(k) / math.pow(k, 4.0))) else: tmp = math.pow(l, 2.0) * (math.pow(t_m, -3.0) * math.pow(k, -2.0)) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.7e+38) tmp = Float64(2.0 * Float64(Float64((l ^ 2.0) / t_m) * Float64(cos(k) / (k ^ 4.0)))); else tmp = Float64((l ^ 2.0) * Float64((t_m ^ -3.0) * (k ^ -2.0))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1.7e+38) tmp = 2.0 * (((l ^ 2.0) / t_m) * (cos(k) / (k ^ 4.0))); else tmp = (l ^ 2.0) * ((t_m ^ -3.0) * (k ^ -2.0)); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 1.7e+38], N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / t$95$m), $MachinePrecision] * N[(N[Cos[k], $MachinePrecision] / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Power[l, 2.0], $MachinePrecision] * N[(N[Power[t$95$m, -3.0], $MachinePrecision] * N[Power[k, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 1.7 \cdot 10^{+38}:\\
\;\;\;\;2 \cdot \left(\frac{{\ell}^{2}}{t_m} \cdot \frac{\cos k}{{k}^{4}}\right)\\
\mathbf{else}:\\
\;\;\;\;{\ell}^{2} \cdot \left({t_m}^{-3} \cdot {k}^{-2}\right)\\
\end{array}
\end{array}
if t < 1.69999999999999998e38Initial program 53.5%
Simplified52.5%
Taylor expanded in k around inf 59.0%
associate-*r/59.0%
associate-*r*59.0%
associate-/r*60.0%
Simplified60.0%
Taylor expanded in k around 0 51.0%
Taylor expanded in l around 0 49.8%
*-commutative49.8%
times-frac51.6%
Simplified51.6%
if 1.69999999999999998e38 < t Initial program 46.3%
Simplified51.5%
Taylor expanded in k around 0 32.8%
*-commutative32.8%
associate-/r*30.9%
Simplified30.9%
expm1-log1p-u30.4%
expm1-udef30.5%
div-inv29.2%
div-inv29.2%
pow-flip29.2%
metadata-eval29.2%
pow-flip29.2%
metadata-eval29.2%
Applied egg-rr29.2%
expm1-def29.2%
expm1-log1p29.6%
associate-*l*31.6%
Simplified31.6%
Final simplification47.3%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 1.85e+40)
(* 2.0 (* (/ (pow l 2.0) (pow k 4.0)) (/ (cos k) t_m)))
(* (pow l 2.0) (* (pow t_m -3.0) (pow k -2.0))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 1.85e+40) {
tmp = 2.0 * ((pow(l, 2.0) / pow(k, 4.0)) * (cos(k) / t_m));
} else {
tmp = pow(l, 2.0) * (pow(t_m, -3.0) * pow(k, -2.0));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 1.85d+40) then
tmp = 2.0d0 * (((l ** 2.0d0) / (k ** 4.0d0)) * (cos(k) / t_m))
else
tmp = (l ** 2.0d0) * ((t_m ** (-3.0d0)) * (k ** (-2.0d0)))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 1.85e+40) {
tmp = 2.0 * ((Math.pow(l, 2.0) / Math.pow(k, 4.0)) * (Math.cos(k) / t_m));
} else {
tmp = Math.pow(l, 2.0) * (Math.pow(t_m, -3.0) * Math.pow(k, -2.0));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 1.85e+40: tmp = 2.0 * ((math.pow(l, 2.0) / math.pow(k, 4.0)) * (math.cos(k) / t_m)) else: tmp = math.pow(l, 2.0) * (math.pow(t_m, -3.0) * math.pow(k, -2.0)) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 1.85e+40) tmp = Float64(2.0 * Float64(Float64((l ^ 2.0) / (k ^ 4.0)) * Float64(cos(k) / t_m))); else tmp = Float64((l ^ 2.0) * Float64((t_m ^ -3.0) * (k ^ -2.0))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 1.85e+40) tmp = 2.0 * (((l ^ 2.0) / (k ^ 4.0)) * (cos(k) / t_m)); else tmp = (l ^ 2.0) * ((t_m ^ -3.0) * (k ^ -2.0)); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 1.85e+40], N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k], $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Power[l, 2.0], $MachinePrecision] * N[(N[Power[t$95$m, -3.0], $MachinePrecision] * N[Power[k, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 1.85 \cdot 10^{+40}:\\
\;\;\;\;2 \cdot \left(\frac{{\ell}^{2}}{{k}^{4}} \cdot \frac{\cos k}{t_m}\right)\\
\mathbf{else}:\\
\;\;\;\;{\ell}^{2} \cdot \left({t_m}^{-3} \cdot {k}^{-2}\right)\\
\end{array}
\end{array}
if t < 1.85e40Initial program 53.5%
Simplified52.5%
Taylor expanded in k around inf 59.0%
associate-*r/59.0%
associate-*r*59.0%
associate-/r*60.0%
Simplified60.0%
Taylor expanded in k around 0 51.0%
Taylor expanded in l around 0 49.8%
times-frac51.2%
Simplified51.2%
if 1.85e40 < t Initial program 46.3%
Simplified51.5%
Taylor expanded in k around 0 32.8%
*-commutative32.8%
associate-/r*30.9%
Simplified30.9%
expm1-log1p-u30.4%
expm1-udef30.5%
div-inv29.2%
div-inv29.2%
pow-flip29.2%
metadata-eval29.2%
pow-flip29.2%
metadata-eval29.2%
Applied egg-rr29.2%
expm1-def29.2%
expm1-log1p29.6%
associate-*l*31.6%
Simplified31.6%
Final simplification47.0%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 7e-51)
(* (/ (pow l 2.0) (pow k 4.0)) (/ 2.0 t_m))
(/ 1.0 (/ (/ (pow k 2.0) (pow l 2.0)) (pow t_m -3.0))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 7e-51) {
tmp = (pow(l, 2.0) / pow(k, 4.0)) * (2.0 / t_m);
} else {
tmp = 1.0 / ((pow(k, 2.0) / pow(l, 2.0)) / pow(t_m, -3.0));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 7d-51) then
tmp = ((l ** 2.0d0) / (k ** 4.0d0)) * (2.0d0 / t_m)
else
tmp = 1.0d0 / (((k ** 2.0d0) / (l ** 2.0d0)) / (t_m ** (-3.0d0)))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 7e-51) {
tmp = (Math.pow(l, 2.0) / Math.pow(k, 4.0)) * (2.0 / t_m);
} else {
tmp = 1.0 / ((Math.pow(k, 2.0) / Math.pow(l, 2.0)) / Math.pow(t_m, -3.0));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 7e-51: tmp = (math.pow(l, 2.0) / math.pow(k, 4.0)) * (2.0 / t_m) else: tmp = 1.0 / ((math.pow(k, 2.0) / math.pow(l, 2.0)) / math.pow(t_m, -3.0)) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 7e-51) tmp = Float64(Float64((l ^ 2.0) / (k ^ 4.0)) * Float64(2.0 / t_m)); else tmp = Float64(1.0 / Float64(Float64((k ^ 2.0) / (l ^ 2.0)) / (t_m ^ -3.0))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 7e-51) tmp = ((l ^ 2.0) / (k ^ 4.0)) * (2.0 / t_m); else tmp = 1.0 / (((k ^ 2.0) / (l ^ 2.0)) / (t_m ^ -3.0)); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 7e-51], N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision] * N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(N[(N[Power[k, 2.0], $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] / N[Power[t$95$m, -3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t_s \cdot \begin{array}{l}
\mathbf{if}\;t_m \leq 7 \cdot 10^{-51}:\\
\;\;\;\;\frac{{\ell}^{2}}{{k}^{4}} \cdot \frac{2}{t_m}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{{k}^{2}}{{\ell}^{2}}}{{t_m}^{-3}}}\\
\end{array}
\end{array}
if t < 6.9999999999999995e-51Initial program 50.3%
Simplified49.8%
Taylor expanded in k around inf 57.2%
associate-*r/57.2%
associate-*r*57.2%
associate-/r*58.3%
Simplified58.3%
Taylor expanded in k around 0 46.5%
associate-*r/46.5%
*-commutative46.5%
times-frac48.1%
Simplified48.1%
if 6.9999999999999995e-51 < t Initial program 56.3%
Simplified59.0%
Taylor expanded in k around 0 43.1%
*-commutative43.1%
associate-/r*41.6%
Simplified41.6%
clear-num41.6%
inv-pow41.6%
div-inv41.6%
pow-flip41.6%
metadata-eval41.6%
Applied egg-rr41.6%
unpow-141.6%
associate-/r*41.2%
Simplified41.2%
Final simplification46.2%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 3.4e-35)
(* (/ (pow l 2.0) (pow k 4.0)) (/ 2.0 t_m))
(* (pow k -2.0) (* (pow l 2.0) (pow t_m -3.0))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 3.4e-35) {
tmp = (pow(l, 2.0) / pow(k, 4.0)) * (2.0 / t_m);
} else {
tmp = pow(k, -2.0) * (pow(l, 2.0) * pow(t_m, -3.0));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 3.4d-35) then
tmp = ((l ** 2.0d0) / (k ** 4.0d0)) * (2.0d0 / t_m)
else
tmp = (k ** (-2.0d0)) * ((l ** 2.0d0) * (t_m ** (-3.0d0)))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 3.4e-35) {
tmp = (Math.pow(l, 2.0) / Math.pow(k, 4.0)) * (2.0 / t_m);
} else {
tmp = Math.pow(k, -2.0) * (Math.pow(l, 2.0) * Math.pow(t_m, -3.0));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 3.4e-35: tmp = (math.pow(l, 2.0) / math.pow(k, 4.0)) * (2.0 / t_m) else: tmp = math.pow(k, -2.0) * (math.pow(l, 2.0) * math.pow(t_m, -3.0)) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 3.4e-35) tmp = Float64(Float64((l ^ 2.0) / (k ^ 4.0)) * Float64(2.0 / t_m)); else tmp = Float64((k ^ -2.0) * Float64((l ^ 2.0) * (t_m ^ -3.0))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 3.4e-35) tmp = ((l ^ 2.0) / (k ^ 4.0)) * (2.0 / t_m); else tmp = (k ^ -2.0) * ((l ^ 2.0) * (t_m ^ -3.0)); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 3.4e-35], N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision] * N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision], N[(N[Power[k, -2.0], $MachinePrecision] * N[(N[Power[l, 2.0], $MachinePrecision] * N[Power[t$95$m, -3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
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.4 \cdot 10^{-35}:\\
\;\;\;\;\frac{{\ell}^{2}}{{k}^{4}} \cdot \frac{2}{t_m}\\
\mathbf{else}:\\
\;\;\;\;{k}^{-2} \cdot \left({\ell}^{2} \cdot {t_m}^{-3}\right)\\
\end{array}
\end{array}
if t < 3.4000000000000003e-35Initial program 51.6%
Simplified51.1%
Taylor expanded in k around inf 58.3%
associate-*r/58.3%
associate-*r*58.3%
associate-/r*59.4%
Simplified59.4%
Taylor expanded in k around 0 47.4%
associate-*r/47.4%
*-commutative47.4%
times-frac48.9%
Simplified48.9%
if 3.4000000000000003e-35 < t Initial program 53.0%
Simplified55.9%
Taylor expanded in k around 0 40.1%
*-commutative40.1%
associate-/r*38.5%
Simplified38.5%
div-inv37.4%
div-inv37.5%
pow-flip37.5%
metadata-eval37.5%
pow-flip37.5%
metadata-eval37.5%
Applied egg-rr37.5%
Final simplification46.0%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= t_m 6.8e-50)
(* (/ (pow l 2.0) (pow k 4.0)) (/ 2.0 t_m))
(* (pow l 2.0) (* (pow t_m -3.0) (pow k -2.0))))))t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
double tmp;
if (t_m <= 6.8e-50) {
tmp = (pow(l, 2.0) / pow(k, 4.0)) * (2.0 / t_m);
} else {
tmp = pow(l, 2.0) * (pow(t_m, -3.0) * pow(k, -2.0));
}
return t_s * tmp;
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t_m <= 6.8d-50) then
tmp = ((l ** 2.0d0) / (k ** 4.0d0)) * (2.0d0 / t_m)
else
tmp = (l ** 2.0d0) * ((t_m ** (-3.0d0)) * (k ** (-2.0d0)))
end if
code = t_s * tmp
end function
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) {
double tmp;
if (t_m <= 6.8e-50) {
tmp = (Math.pow(l, 2.0) / Math.pow(k, 4.0)) * (2.0 / t_m);
} else {
tmp = Math.pow(l, 2.0) * (Math.pow(t_m, -3.0) * Math.pow(k, -2.0));
}
return t_s * tmp;
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): tmp = 0 if t_m <= 6.8e-50: tmp = (math.pow(l, 2.0) / math.pow(k, 4.0)) * (2.0 / t_m) else: tmp = math.pow(l, 2.0) * (math.pow(t_m, -3.0) * math.pow(k, -2.0)) return t_s * tmp
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) tmp = 0.0 if (t_m <= 6.8e-50) tmp = Float64(Float64((l ^ 2.0) / (k ^ 4.0)) * Float64(2.0 / t_m)); else tmp = Float64((l ^ 2.0) * Float64((t_m ^ -3.0) * (k ^ -2.0))); end return Float64(t_s * tmp) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k) tmp = 0.0; if (t_m <= 6.8e-50) tmp = ((l ^ 2.0) / (k ^ 4.0)) * (2.0 / t_m); else tmp = (l ^ 2.0) * ((t_m ^ -3.0) * (k ^ -2.0)); end tmp_2 = t_s * tmp; end
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_] := N[(t$95$s * If[LessEqual[t$95$m, 6.8e-50], N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision] * N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision], N[(N[Power[l, 2.0], $MachinePrecision] * N[(N[Power[t$95$m, -3.0], $MachinePrecision] * N[Power[k, -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
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^{-50}:\\
\;\;\;\;\frac{{\ell}^{2}}{{k}^{4}} \cdot \frac{2}{t_m}\\
\mathbf{else}:\\
\;\;\;\;{\ell}^{2} \cdot \left({t_m}^{-3} \cdot {k}^{-2}\right)\\
\end{array}
\end{array}
if t < 6.80000000000000029e-50Initial program 50.3%
Simplified49.8%
Taylor expanded in k around inf 57.2%
associate-*r/57.2%
associate-*r*57.2%
associate-/r*58.3%
Simplified58.3%
Taylor expanded in k around 0 46.5%
associate-*r/46.5%
*-commutative46.5%
times-frac48.1%
Simplified48.1%
if 6.80000000000000029e-50 < t Initial program 56.3%
Simplified59.0%
Taylor expanded in k around 0 43.1%
*-commutative43.1%
associate-/r*41.6%
Simplified41.6%
expm1-log1p-u41.0%
expm1-udef41.1%
div-inv40.1%
div-inv40.1%
pow-flip40.1%
metadata-eval40.1%
pow-flip40.1%
metadata-eval40.1%
Applied egg-rr40.1%
expm1-def40.0%
expm1-log1p40.6%
associate-*l*42.1%
Simplified42.1%
Final simplification46.5%
t_m = (fabs.f64 t) t_s = (copysign.f64 1 t) (FPCore (t_s t_m l k) :precision binary64 (* t_s (* 2.0 (/ (pow l 2.0) (* t_m (pow k 4.0))))))
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
return t_s * (2.0 * (pow(l, 2.0) / (t_m * pow(k, 4.0))));
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
code = t_s * (2.0d0 * ((l ** 2.0d0) / (t_m * (k ** 4.0d0))))
end function
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) {
return t_s * (2.0 * (Math.pow(l, 2.0) / (t_m * Math.pow(k, 4.0))));
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): return t_s * (2.0 * (math.pow(l, 2.0) / (t_m * math.pow(k, 4.0))))
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) return Float64(t_s * Float64(2.0 * Float64((l ^ 2.0) / Float64(t_m * (k ^ 4.0))))) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k) tmp = t_s * (2.0 * ((l ^ 2.0) / (t_m * (k ^ 4.0)))); end
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_] := N[(t$95$s * N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] / N[(t$95$m * N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t_s \cdot \left(2 \cdot \frac{{\ell}^{2}}{t_m \cdot {k}^{4}}\right)
\end{array}
Initial program 51.9%
Simplified52.3%
Taylor expanded in k around inf 52.5%
associate-*r/52.5%
associate-*r*52.4%
associate-/r*53.3%
Simplified53.3%
Taylor expanded in k around 0 44.0%
Final simplification44.0%
t_m = (fabs.f64 t) t_s = (copysign.f64 1 t) (FPCore (t_s t_m l k) :precision binary64 (* t_s (* (/ (pow l 2.0) (pow k 4.0)) (/ 2.0 t_m))))
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
return t_s * ((pow(l, 2.0) / pow(k, 4.0)) * (2.0 / t_m));
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
code = t_s * (((l ** 2.0d0) / (k ** 4.0d0)) * (2.0d0 / t_m))
end function
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) {
return t_s * ((Math.pow(l, 2.0) / Math.pow(k, 4.0)) * (2.0 / t_m));
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): return t_s * ((math.pow(l, 2.0) / math.pow(k, 4.0)) * (2.0 / t_m))
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) return Float64(t_s * Float64(Float64((l ^ 2.0) / (k ^ 4.0)) * Float64(2.0 / t_m))) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k) tmp = t_s * (((l ^ 2.0) / (k ^ 4.0)) * (2.0 / t_m)); end
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_] := N[(t$95$s * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision] * N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t_s \cdot \left(\frac{{\ell}^{2}}{{k}^{4}} \cdot \frac{2}{t_m}\right)
\end{array}
Initial program 51.9%
Simplified52.3%
Taylor expanded in k around inf 52.5%
associate-*r/52.5%
associate-*r*52.4%
associate-/r*53.3%
Simplified53.3%
Taylor expanded in k around 0 44.0%
associate-*r/44.0%
*-commutative44.0%
times-frac45.2%
Simplified45.2%
Final simplification45.2%
t_m = (fabs.f64 t) t_s = (copysign.f64 1 t) (FPCore (t_s t_m l k) :precision binary64 (* t_s (/ (* (pow l 2.0) (/ 2.0 t_m)) (pow k 4.0))))
t_m = fabs(t);
t_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k) {
return t_s * ((pow(l, 2.0) * (2.0 / t_m)) / pow(k, 4.0));
}
t_m = abs(t)
t_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k
code = t_s * (((l ** 2.0d0) * (2.0d0 / t_m)) / (k ** 4.0d0))
end function
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) {
return t_s * ((Math.pow(l, 2.0) * (2.0 / t_m)) / Math.pow(k, 4.0));
}
t_m = math.fabs(t) t_s = math.copysign(1.0, t) def code(t_s, t_m, l, k): return t_s * ((math.pow(l, 2.0) * (2.0 / t_m)) / math.pow(k, 4.0))
t_m = abs(t) t_s = copysign(1.0, t) function code(t_s, t_m, l, k) return Float64(t_s * Float64(Float64((l ^ 2.0) * Float64(2.0 / t_m)) / (k ^ 4.0))) end
t_m = abs(t); t_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k) tmp = t_s * (((l ^ 2.0) * (2.0 / t_m)) / (k ^ 4.0)); end
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_] := N[(t$95$s * N[(N[(N[Power[l, 2.0], $MachinePrecision] * N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision] / N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
t_s \cdot \frac{{\ell}^{2} \cdot \frac{2}{t_m}}{{k}^{4}}
\end{array}
Initial program 51.9%
Simplified52.3%
Taylor expanded in k around inf 52.5%
associate-*r/52.5%
associate-*r*52.4%
associate-/r*53.3%
Simplified53.3%
Taylor expanded in k around 0 44.0%
associate-*r/44.0%
*-commutative44.0%
times-frac45.2%
Simplified45.2%
associate-*r/45.4%
Applied egg-rr45.4%
Final simplification45.4%
herbie shell --seed 2023322
(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))))