
(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 21 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
(*
t_s
(if (<= t_m 2.05e-54)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) (/ l (tan k)))
(pow
(/
(/ (cbrt (/ 2.0 (sin k))) (* t_m (pow (cbrt l) -2.0)))
(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 tmp;
if (t_m <= 2.05e-54) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * (l / tan(k));
} else {
tmp = pow(((cbrt((2.0 / sin(k))) / (t_m * pow(cbrt(l), -2.0))) / 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 tmp;
if (t_m <= 2.05e-54) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * (l / Math.tan(k));
} else {
tmp = Math.pow(((Math.cbrt((2.0 / Math.sin(k))) / (t_m * Math.pow(Math.cbrt(l), -2.0))) / 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) tmp = 0.0 if (t_m <= 2.05e-54) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * Float64(l / tan(k))); else tmp = Float64(Float64(cbrt(Float64(2.0 / sin(k))) / Float64(t_m * (cbrt(l) ^ -2.0))) / 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_] := N[(t$95$s * If[LessEqual[t$95$m, 2.05e-54], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(N[Power[N[(2.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] / N[(t$95$m * N[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]), $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]
\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.05 \cdot 10^{-54}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot \frac{\ell}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{\frac{\sqrt[3]{\frac{2}{\sin k}}}{t\_m \cdot {\left(\sqrt[3]{\ell}\right)}^{-2}}}{\sqrt[3]{\tan k \cdot \left(2 + {\left(\frac{k}{t\_m}\right)}^{2}\right)}}\right)}^{3}\\
\end{array}
\end{array}
if t < 2.05e-54Initial program 47.0%
Simplified52.9%
associate-/l/52.9%
associate-/r/54.5%
times-frac54.4%
div-inv54.3%
clear-num54.4%
Applied egg-rr54.4%
Taylor expanded in k around inf 67.6%
if 2.05e-54 < t Initial program 78.6%
Simplified82.6%
associate-/l/82.6%
add-cube-cbrt82.6%
times-frac82.5%
Applied egg-rr94.9%
add-cube-cbrt94.8%
pow394.8%
Applied egg-rr96.1%
Final simplification75.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 4.8e-54)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) (/ l (tan k)))
(/
(/
(pow (/ (cbrt (/ 2.0 (sin k))) (/ t_m (pow (cbrt l) 2.0))) 3.0)
(tan k))
(+ 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-54) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * (l / tan(k));
} else {
tmp = (pow((cbrt((2.0 / sin(k))) / (t_m / pow(cbrt(l), 2.0))), 3.0) / tan(k)) / (2.0 + pow((k / t_m), 2.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 <= 4.8e-54) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * (l / Math.tan(k));
} else {
tmp = (Math.pow((Math.cbrt((2.0 / Math.sin(k))) / (t_m / Math.pow(Math.cbrt(l), 2.0))), 3.0) / Math.tan(k)) / (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-54) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * Float64(l / tan(k))); else tmp = Float64(Float64((Float64(cbrt(Float64(2.0 / sin(k))) / Float64(t_m / (cbrt(l) ^ 2.0))) ^ 3.0) / tan(k)) / Float64(2.0 + (Float64(k / t_m) ^ 2.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, 4.8e-54], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[Power[N[(N[Power[N[(2.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] / N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] / N[Tan[k], $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^{-54}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot \frac{\ell}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{{\left(\frac{\sqrt[3]{\frac{2}{\sin k}}}{\frac{t\_m}{{\left(\sqrt[3]{\ell}\right)}^{2}}}\right)}^{3}}{\tan k}}{2 + {\left(\frac{k}{t\_m}\right)}^{2}}\\
\end{array}
\end{array}
if t < 4.80000000000000026e-54Initial program 47.0%
Simplified52.9%
associate-/l/52.9%
associate-/r/54.5%
times-frac54.4%
div-inv54.3%
clear-num54.4%
Applied egg-rr54.4%
Taylor expanded in k around inf 67.6%
if 4.80000000000000026e-54 < t Initial program 78.6%
Simplified82.6%
add-cube-cbrt82.6%
pow382.5%
cbrt-div82.5%
associate-/r*78.6%
cbrt-div78.5%
rem-cbrt-cube80.2%
cbrt-prod92.1%
pow292.1%
Applied egg-rr92.1%
Final simplification73.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_3 (/ l (tan k))))
(*
t_s
(if (<= t_m 2.15e-55)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) t_3)
(if (<= t_m 7.4e+136)
(* t_3 (/ (pow (/ (cbrt (* l (/ 2.0 (sin k)))) t_m) 3.0) (+ 2.0 t_2)))
(/
2.0
(*
(* (tan k) (* (sin k) (pow (/ t_m (pow (cbrt l) 2.0)) 3.0)))
(+ 1.0 (+ t_2 1.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 t_3 = l / tan(k);
double tmp;
if (t_m <= 2.15e-55) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * t_3;
} else if (t_m <= 7.4e+136) {
tmp = t_3 * (pow((cbrt((l * (2.0 / sin(k)))) / t_m), 3.0) / (2.0 + t_2));
} else {
tmp = 2.0 / ((tan(k) * (sin(k) * pow((t_m / pow(cbrt(l), 2.0)), 3.0))) * (1.0 + (t_2 + 1.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((k / t_m), 2.0);
double t_3 = l / Math.tan(k);
double tmp;
if (t_m <= 2.15e-55) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * t_3;
} else if (t_m <= 7.4e+136) {
tmp = t_3 * (Math.pow((Math.cbrt((l * (2.0 / Math.sin(k)))) / t_m), 3.0) / (2.0 + t_2));
} else {
tmp = 2.0 / ((Math.tan(k) * (Math.sin(k) * Math.pow((t_m / Math.pow(Math.cbrt(l), 2.0)), 3.0))) * (1.0 + (t_2 + 1.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 t_3 = Float64(l / tan(k)) tmp = 0.0 if (t_m <= 2.15e-55) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * t_3); elseif (t_m <= 7.4e+136) tmp = Float64(t_3 * Float64((Float64(cbrt(Float64(l * Float64(2.0 / sin(k)))) / t_m) ^ 3.0) / Float64(2.0 + t_2))); else tmp = Float64(2.0 / Float64(Float64(tan(k) * Float64(sin(k) * (Float64(t_m / (cbrt(l) ^ 2.0)) ^ 3.0))) * Float64(1.0 + Float64(t_2 + 1.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[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$3 = N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 2.15e-55], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$3), $MachinePrecision], If[LessEqual[t$95$m, 7.4e+136], N[(t$95$3 * N[(N[Power[N[(N[Power[N[(l * N[(2.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] / t$95$m), $MachinePrecision], 3.0], $MachinePrecision] / N[(2.0 + t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[Power[N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(t$95$2 + 1.0), $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_3 := \frac{\ell}{\tan k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 2.15 \cdot 10^{-55}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot t\_3\\
\mathbf{elif}\;t\_m \leq 7.4 \cdot 10^{+136}:\\
\;\;\;\;t\_3 \cdot \frac{{\left(\frac{\sqrt[3]{\ell \cdot \frac{2}{\sin k}}}{t\_m}\right)}^{3}}{2 + t\_2}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot \left(\sin k \cdot {\left(\frac{t\_m}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}\right)\right) \cdot \left(1 + \left(t\_2 + 1\right)\right)}\\
\end{array}
\end{array}
\end{array}
if t < 2.15000000000000005e-55Initial program 46.4%
Simplified52.4%
associate-/l/52.4%
associate-/r/54.0%
times-frac53.9%
div-inv53.9%
clear-num53.9%
Applied egg-rr53.9%
Taylor expanded in k around inf 67.2%
if 2.15000000000000005e-55 < t < 7.4000000000000002e136Initial program 76.8%
Simplified80.3%
associate-/l/80.3%
associate-/r/82.5%
times-frac84.7%
div-inv84.7%
clear-num84.8%
Applied egg-rr84.8%
add-cube-cbrt84.3%
pow384.4%
associate-*r/86.4%
cbrt-div86.4%
unpow386.4%
add-cbrt-cube95.1%
Applied egg-rr95.1%
if 7.4000000000000002e136 < t Initial program 83.6%
associate-/r*88.2%
add-cube-cbrt88.2%
pow388.2%
associate-/r*83.6%
cbrt-div83.6%
rem-cbrt-cube83.9%
cbrt-prod92.1%
pow292.1%
Applied egg-rr92.1%
Final simplification74.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_3 (/ l (tan k)))
(t_4 (/ 2.0 (sin k))))
(*
t_s
(if (<= t_m 1.16e-55)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) t_3)
(if (<= t_m 6.6e+136)
(* t_3 (/ (pow (/ (cbrt (* l t_4)) t_m) 3.0) t_2))
(/ (/ t_4 (pow (* t_m (pow (cbrt l) -2.0)) 3.0)) (* (tan k) t_2)))))))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 = l / tan(k);
double t_4 = 2.0 / sin(k);
double tmp;
if (t_m <= 1.16e-55) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * t_3;
} else if (t_m <= 6.6e+136) {
tmp = t_3 * (pow((cbrt((l * t_4)) / t_m), 3.0) / t_2);
} else {
tmp = (t_4 / pow((t_m * pow(cbrt(l), -2.0)), 3.0)) / (tan(k) * t_2);
}
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 = l / Math.tan(k);
double t_4 = 2.0 / Math.sin(k);
double tmp;
if (t_m <= 1.16e-55) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * t_3;
} else if (t_m <= 6.6e+136) {
tmp = t_3 * (Math.pow((Math.cbrt((l * t_4)) / t_m), 3.0) / t_2);
} else {
tmp = (t_4 / Math.pow((t_m * Math.pow(Math.cbrt(l), -2.0)), 3.0)) / (Math.tan(k) * t_2);
}
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(l / tan(k)) t_4 = Float64(2.0 / sin(k)) tmp = 0.0 if (t_m <= 1.16e-55) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * t_3); elseif (t_m <= 6.6e+136) tmp = Float64(t_3 * Float64((Float64(cbrt(Float64(l * t_4)) / t_m) ^ 3.0) / t_2)); else tmp = Float64(Float64(t_4 / (Float64(t_m * (cbrt(l) ^ -2.0)) ^ 3.0)) / Float64(tan(k) * t_2)); 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[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(2.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 1.16e-55], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$3), $MachinePrecision], If[LessEqual[t$95$m, 6.6e+136], N[(t$95$3 * N[(N[Power[N[(N[Power[N[(l * t$95$4), $MachinePrecision], 1/3], $MachinePrecision] / t$95$m), $MachinePrecision], 3.0], $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision], N[(N[(t$95$4 / N[Power[N[(t$95$m * N[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision] / N[(N[Tan[k], $MachinePrecision] * t$95$2), $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{\ell}{\tan k}\\
t_4 := \frac{2}{\sin k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.16 \cdot 10^{-55}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot t\_3\\
\mathbf{elif}\;t\_m \leq 6.6 \cdot 10^{+136}:\\
\;\;\;\;t\_3 \cdot \frac{{\left(\frac{\sqrt[3]{\ell \cdot t\_4}}{t\_m}\right)}^{3}}{t\_2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{t\_4}{{\left(t\_m \cdot {\left(\sqrt[3]{\ell}\right)}^{-2}\right)}^{3}}}{\tan k \cdot t\_2}\\
\end{array}
\end{array}
\end{array}
if t < 1.15999999999999996e-55Initial program 46.4%
Simplified52.4%
associate-/l/52.4%
associate-/r/54.0%
times-frac53.9%
div-inv53.9%
clear-num53.9%
Applied egg-rr53.9%
Taylor expanded in k around inf 67.2%
if 1.15999999999999996e-55 < t < 6.59999999999999984e136Initial program 76.8%
Simplified80.3%
associate-/l/80.3%
associate-/r/82.5%
times-frac84.7%
div-inv84.7%
clear-num84.8%
Applied egg-rr84.8%
add-cube-cbrt84.3%
pow384.4%
associate-*r/86.4%
cbrt-div86.4%
unpow386.4%
add-cbrt-cube95.1%
Applied egg-rr95.1%
if 6.59999999999999984e136 < t Initial program 83.6%
Simplified88.2%
associate-/l/88.2%
add-cube-cbrt88.2%
times-frac88.2%
Applied egg-rr95.8%
frac-times95.8%
unpow295.8%
pow395.8%
div-inv95.7%
pow-flip95.7%
metadata-eval95.7%
Applied egg-rr95.7%
cube-div92.0%
rem-cube-cbrt92.1%
Simplified92.1%
Final simplification74.4%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (/ (pow t_m 1.5) l)))
(*
t_s
(if (<= t_m 2.8e-54)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) (/ l (tan k)))
(if (<= t_m 1.25e+207)
(/
(/ (* (/ 2.0 t_2) (/ (/ 1.0 (sin k)) t_2)) (tan k))
(+ 2.0 (pow (/ k t_m) 2.0)))
(pow
(/
(/ (cbrt (/ 2.0 (sin k))) (* t_m (pow (cbrt l) -2.0)))
(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 = pow(t_m, 1.5) / l;
double tmp;
if (t_m <= 2.8e-54) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * (l / tan(k));
} else if (t_m <= 1.25e+207) {
tmp = (((2.0 / t_2) * ((1.0 / sin(k)) / t_2)) / tan(k)) / (2.0 + pow((k / t_m), 2.0));
} else {
tmp = pow(((cbrt((2.0 / sin(k))) / (t_m * pow(cbrt(l), -2.0))) / 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.pow(t_m, 1.5) / l;
double tmp;
if (t_m <= 2.8e-54) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * (l / Math.tan(k));
} else if (t_m <= 1.25e+207) {
tmp = (((2.0 / t_2) * ((1.0 / Math.sin(k)) / t_2)) / Math.tan(k)) / (2.0 + Math.pow((k / t_m), 2.0));
} else {
tmp = Math.pow(((Math.cbrt((2.0 / Math.sin(k))) / (t_m * Math.pow(Math.cbrt(l), -2.0))) / 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((t_m ^ 1.5) / l) tmp = 0.0 if (t_m <= 2.8e-54) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * Float64(l / tan(k))); elseif (t_m <= 1.25e+207) tmp = Float64(Float64(Float64(Float64(2.0 / t_2) * Float64(Float64(1.0 / sin(k)) / t_2)) / tan(k)) / Float64(2.0 + (Float64(k / t_m) ^ 2.0))); else tmp = Float64(Float64(cbrt(Float64(2.0 / sin(k))) / Float64(t_m * (cbrt(l) ^ -2.0))) / 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[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 2.8e-54], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 1.25e+207], N[(N[(N[(N[(2.0 / t$95$2), $MachinePrecision] * N[(N[(1.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(N[Power[N[(2.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] / N[(t$95$m * N[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[(2.0 * k), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $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}^{1.5}}{\ell}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 2.8 \cdot 10^{-54}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot \frac{\ell}{\tan k}\\
\mathbf{elif}\;t\_m \leq 1.25 \cdot 10^{+207}:\\
\;\;\;\;\frac{\frac{\frac{2}{t\_2} \cdot \frac{\frac{1}{\sin k}}{t\_2}}{\tan k}}{2 + {\left(\frac{k}{t\_m}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{\frac{\sqrt[3]{\frac{2}{\sin k}}}{t\_m \cdot {\left(\sqrt[3]{\ell}\right)}^{-2}}}{\sqrt[3]{2 \cdot k}}\right)}^{3}\\
\end{array}
\end{array}
\end{array}
if t < 2.8000000000000002e-54Initial program 47.0%
Simplified52.9%
associate-/l/52.9%
associate-/r/54.5%
times-frac54.4%
div-inv54.3%
clear-num54.4%
Applied egg-rr54.4%
Taylor expanded in k around inf 67.6%
if 2.8000000000000002e-54 < t < 1.25e207Initial program 78.0%
Simplified81.2%
div-inv81.2%
add-sqr-sqrt81.1%
times-frac81.1%
associate-/r*78.1%
sqrt-div78.1%
sqrt-pow178.1%
metadata-eval78.1%
sqrt-prod47.6%
add-sqr-sqrt61.8%
associate-/r*60.6%
sqrt-div60.6%
sqrt-pow160.7%
metadata-eval60.7%
sqrt-prod51.4%
add-sqr-sqrt92.4%
Applied egg-rr92.4%
if 1.25e207 < t Initial program 80.5%
Simplified87.6%
associate-/l/87.6%
add-cube-cbrt87.6%
times-frac87.6%
Applied egg-rr93.6%
add-cube-cbrt93.6%
pow393.6%
Applied egg-rr99.4%
Taylor expanded in k around 0 99.4%
*-commutative99.4%
Simplified99.4%
Final simplification74.4%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (/ (pow t_m 1.5) l)))
(*
t_s
(if (<= t_m 4.2e-54)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) (/ l (tan k)))
(if (<= t_m 3.9e+189)
(/
(/ (* (/ 2.0 t_2) (/ (/ 1.0 (sin k)) t_2)) (tan k))
(+ 2.0 (pow (/ k t_m) 2.0)))
(pow (/ (/ (pow (cbrt l) 2.0) (pow (cbrt k) 2.0)) 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 t_2 = pow(t_m, 1.5) / l;
double tmp;
if (t_m <= 4.2e-54) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * (l / tan(k));
} else if (t_m <= 3.9e+189) {
tmp = (((2.0 / t_2) * ((1.0 / sin(k)) / t_2)) / tan(k)) / (2.0 + pow((k / t_m), 2.0));
} else {
tmp = pow(((pow(cbrt(l), 2.0) / pow(cbrt(k), 2.0)) / t_m), 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(t_m, 1.5) / l;
double tmp;
if (t_m <= 4.2e-54) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * (l / Math.tan(k));
} else if (t_m <= 3.9e+189) {
tmp = (((2.0 / t_2) * ((1.0 / Math.sin(k)) / t_2)) / Math.tan(k)) / (2.0 + Math.pow((k / t_m), 2.0));
} else {
tmp = Math.pow(((Math.pow(Math.cbrt(l), 2.0) / Math.pow(Math.cbrt(k), 2.0)) / 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) t_2 = Float64((t_m ^ 1.5) / l) tmp = 0.0 if (t_m <= 4.2e-54) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * Float64(l / tan(k))); elseif (t_m <= 3.9e+189) tmp = Float64(Float64(Float64(Float64(2.0 / t_2) * Float64(Float64(1.0 / sin(k)) / t_2)) / tan(k)) / Float64(2.0 + (Float64(k / t_m) ^ 2.0))); else tmp = Float64(Float64((cbrt(l) ^ 2.0) / (cbrt(k) ^ 2.0)) / t_m) ^ 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[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 4.2e-54], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 3.9e+189], N[(N[(N[(N[(2.0 / t$95$2), $MachinePrecision] * N[(N[(1.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision] / N[Power[N[Power[k, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision], 3.0], $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}^{1.5}}{\ell}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 4.2 \cdot 10^{-54}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot \frac{\ell}{\tan k}\\
\mathbf{elif}\;t\_m \leq 3.9 \cdot 10^{+189}:\\
\;\;\;\;\frac{\frac{\frac{2}{t\_2} \cdot \frac{\frac{1}{\sin k}}{t\_2}}{\tan k}}{2 + {\left(\frac{k}{t\_m}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{\frac{{\left(\sqrt[3]{\ell}\right)}^{2}}{{\left(\sqrt[3]{k}\right)}^{2}}}{t\_m}\right)}^{3}\\
\end{array}
\end{array}
\end{array}
if t < 4.2e-54Initial program 47.0%
Simplified52.9%
associate-/l/52.9%
associate-/r/54.5%
times-frac54.4%
div-inv54.3%
clear-num54.4%
Applied egg-rr54.4%
Taylor expanded in k around inf 67.6%
if 4.2e-54 < t < 3.9e189Initial program 76.7%
Simplified80.0%
div-inv80.0%
add-sqr-sqrt79.9%
times-frac80.0%
associate-/r*76.8%
sqrt-div76.8%
sqrt-pow176.8%
metadata-eval76.8%
sqrt-prod46.5%
add-sqr-sqrt59.4%
associate-/r*58.1%
sqrt-div58.1%
sqrt-pow158.2%
metadata-eval58.2%
sqrt-prod50.4%
add-sqr-sqrt91.9%
Applied egg-rr91.9%
if 3.9e189 < t Initial program 83.7%
Simplified89.6%
Taylor expanded in k around 0 55.6%
add-cube-cbrt55.6%
pow355.6%
associate-/r*55.6%
cbrt-div55.6%
unpow355.6%
add-cbrt-cube56.1%
Applied egg-rr56.1%
cbrt-div56.1%
unpow256.1%
cbrt-prod67.2%
pow267.2%
unpow267.2%
cbrt-prod99.4%
pow299.4%
Applied egg-rr99.4%
Final simplification74.4%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (/ l (tan k))))
(*
t_s
(if (<= t_m 1.35e-55)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) t_2)
(if (<= t_m 8e+136)
(*
t_2
(/
(pow (/ (cbrt (* l (/ 2.0 (sin k)))) t_m) 3.0)
(+ 2.0 (pow (/ k t_m) 2.0))))
(pow (/ (/ (pow (cbrt l) 2.0) (pow (cbrt k) 2.0)) 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 t_2 = l / tan(k);
double tmp;
if (t_m <= 1.35e-55) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * t_2;
} else if (t_m <= 8e+136) {
tmp = t_2 * (pow((cbrt((l * (2.0 / sin(k)))) / t_m), 3.0) / (2.0 + pow((k / t_m), 2.0)));
} else {
tmp = pow(((pow(cbrt(l), 2.0) / pow(cbrt(k), 2.0)) / t_m), 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 = l / Math.tan(k);
double tmp;
if (t_m <= 1.35e-55) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * t_2;
} else if (t_m <= 8e+136) {
tmp = t_2 * (Math.pow((Math.cbrt((l * (2.0 / Math.sin(k)))) / t_m), 3.0) / (2.0 + Math.pow((k / t_m), 2.0)));
} else {
tmp = Math.pow(((Math.pow(Math.cbrt(l), 2.0) / Math.pow(Math.cbrt(k), 2.0)) / 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) t_2 = Float64(l / tan(k)) tmp = 0.0 if (t_m <= 1.35e-55) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * t_2); elseif (t_m <= 8e+136) tmp = Float64(t_2 * Float64((Float64(cbrt(Float64(l * Float64(2.0 / sin(k)))) / t_m) ^ 3.0) / Float64(2.0 + (Float64(k / t_m) ^ 2.0)))); else tmp = Float64(Float64((cbrt(l) ^ 2.0) / (cbrt(k) ^ 2.0)) / t_m) ^ 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[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 1.35e-55], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision], If[LessEqual[t$95$m, 8e+136], N[(t$95$2 * N[(N[Power[N[(N[Power[N[(l * N[(2.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] / t$95$m), $MachinePrecision], 3.0], $MachinePrecision] / N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision] / N[Power[N[Power[k, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision], 3.0], $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
t_m = \left|t\right|
\\
t_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \frac{\ell}{\tan k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 1.35 \cdot 10^{-55}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot t\_2\\
\mathbf{elif}\;t\_m \leq 8 \cdot 10^{+136}:\\
\;\;\;\;t\_2 \cdot \frac{{\left(\frac{\sqrt[3]{\ell \cdot \frac{2}{\sin k}}}{t\_m}\right)}^{3}}{2 + {\left(\frac{k}{t\_m}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{\frac{{\left(\sqrt[3]{\ell}\right)}^{2}}{{\left(\sqrt[3]{k}\right)}^{2}}}{t\_m}\right)}^{3}\\
\end{array}
\end{array}
\end{array}
if t < 1.35000000000000002e-55Initial program 46.4%
Simplified52.4%
associate-/l/52.4%
associate-/r/54.0%
times-frac53.9%
div-inv53.9%
clear-num53.9%
Applied egg-rr53.9%
Taylor expanded in k around inf 67.2%
if 1.35000000000000002e-55 < t < 8.00000000000000047e136Initial program 76.8%
Simplified80.3%
associate-/l/80.3%
associate-/r/82.5%
times-frac84.7%
div-inv84.7%
clear-num84.8%
Applied egg-rr84.8%
add-cube-cbrt84.3%
pow384.4%
associate-*r/86.4%
cbrt-div86.4%
unpow386.4%
add-cbrt-cube95.1%
Applied egg-rr95.1%
if 8.00000000000000047e136 < t Initial program 83.6%
Simplified88.2%
Taylor expanded in k around 0 62.5%
add-cube-cbrt62.5%
pow362.5%
associate-/r*62.5%
cbrt-div62.5%
unpow362.5%
add-cbrt-cube63.0%
Applied egg-rr63.0%
cbrt-div63.0%
unpow263.0%
cbrt-prod75.2%
pow275.2%
unpow275.2%
cbrt-prod99.4%
pow299.4%
Applied egg-rr99.4%
Final simplification75.1%
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_3 (/ l (tan k))))
(*
t_s
(if (<= t_m 3.8e-54)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) t_3)
(if (<= t_m 2.9e+74)
(* t_3 (/ (* (/ 2.0 (sin k)) (/ l (pow t_m 3.0))) (+ 2.0 t_2)))
(if (<= t_m 5.8e+153)
(/
2.0
(*
(+ 1.0 (+ t_2 1.0))
(* (tan k) (* (sin k) (* (/ (pow t_m 2.0) l) (/ t_m l))))))
(pow (/ (/ (pow (cbrt l) 2.0) (pow (cbrt k) 2.0)) 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 t_2 = pow((k / t_m), 2.0);
double t_3 = l / tan(k);
double tmp;
if (t_m <= 3.8e-54) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * t_3;
} else if (t_m <= 2.9e+74) {
tmp = t_3 * (((2.0 / sin(k)) * (l / pow(t_m, 3.0))) / (2.0 + t_2));
} else if (t_m <= 5.8e+153) {
tmp = 2.0 / ((1.0 + (t_2 + 1.0)) * (tan(k) * (sin(k) * ((pow(t_m, 2.0) / l) * (t_m / l)))));
} else {
tmp = pow(((pow(cbrt(l), 2.0) / pow(cbrt(k), 2.0)) / t_m), 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((k / t_m), 2.0);
double t_3 = l / Math.tan(k);
double tmp;
if (t_m <= 3.8e-54) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * t_3;
} else if (t_m <= 2.9e+74) {
tmp = t_3 * (((2.0 / Math.sin(k)) * (l / Math.pow(t_m, 3.0))) / (2.0 + t_2));
} else if (t_m <= 5.8e+153) {
tmp = 2.0 / ((1.0 + (t_2 + 1.0)) * (Math.tan(k) * (Math.sin(k) * ((Math.pow(t_m, 2.0) / l) * (t_m / l)))));
} else {
tmp = Math.pow(((Math.pow(Math.cbrt(l), 2.0) / Math.pow(Math.cbrt(k), 2.0)) / 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) t_2 = Float64(k / t_m) ^ 2.0 t_3 = Float64(l / tan(k)) tmp = 0.0 if (t_m <= 3.8e-54) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * t_3); elseif (t_m <= 2.9e+74) tmp = Float64(t_3 * Float64(Float64(Float64(2.0 / sin(k)) * Float64(l / (t_m ^ 3.0))) / Float64(2.0 + t_2))); elseif (t_m <= 5.8e+153) tmp = Float64(2.0 / Float64(Float64(1.0 + Float64(t_2 + 1.0)) * Float64(tan(k) * Float64(sin(k) * Float64(Float64((t_m ^ 2.0) / l) * Float64(t_m / l)))))); else tmp = Float64(Float64((cbrt(l) ^ 2.0) / (cbrt(k) ^ 2.0)) / t_m) ^ 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[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$3 = N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 3.8e-54], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$3), $MachinePrecision], If[LessEqual[t$95$m, 2.9e+74], N[(t$95$3 * N[(N[(N[(2.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 5.8e+153], 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[(N[(N[Power[t$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision] / N[Power[N[Power[k, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision], 3.0], $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_3 := \frac{\ell}{\tan k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 3.8 \cdot 10^{-54}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot t\_3\\
\mathbf{elif}\;t\_m \leq 2.9 \cdot 10^{+74}:\\
\;\;\;\;t\_3 \cdot \frac{\frac{2}{\sin k} \cdot \frac{\ell}{{t\_m}^{3}}}{2 + t\_2}\\
\mathbf{elif}\;t\_m \leq 5.8 \cdot 10^{+153}:\\
\;\;\;\;\frac{2}{\left(1 + \left(t\_2 + 1\right)\right) \cdot \left(\tan k \cdot \left(\sin k \cdot \left(\frac{{t\_m}^{2}}{\ell} \cdot \frac{t\_m}{\ell}\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{\frac{{\left(\sqrt[3]{\ell}\right)}^{2}}{{\left(\sqrt[3]{k}\right)}^{2}}}{t\_m}\right)}^{3}\\
\end{array}
\end{array}
\end{array}
if t < 3.8000000000000002e-54Initial program 47.0%
Simplified52.9%
associate-/l/52.9%
associate-/r/54.5%
times-frac54.4%
div-inv54.3%
clear-num54.4%
Applied egg-rr54.4%
Taylor expanded in k around inf 67.6%
if 3.8000000000000002e-54 < t < 2.9000000000000002e74Initial program 80.8%
Simplified83.3%
associate-/l/83.2%
associate-/r/87.3%
times-frac91.6%
div-inv91.6%
clear-num91.7%
Applied egg-rr91.7%
if 2.9000000000000002e74 < t < 5.80000000000000004e153Initial program 74.9%
unpow374.9%
times-frac95.9%
pow295.9%
Applied egg-rr95.9%
if 5.80000000000000004e153 < t Initial program 80.4%
Simplified85.8%
Taylor expanded in k around 0 55.0%
add-cube-cbrt55.0%
pow355.0%
associate-/r*55.0%
cbrt-div55.0%
unpow355.0%
add-cbrt-cube55.7%
Applied egg-rr55.7%
cbrt-div55.7%
unpow255.7%
cbrt-prod70.2%
pow270.2%
unpow270.2%
cbrt-prod99.3%
pow299.3%
Applied egg-rr99.3%
Final simplification74.8%
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_3 (/ l (tan k))))
(*
t_s
(if (<= t_m 3.2e-55)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) t_3)
(if (<= t_m 1.12e+74)
(* t_3 (/ (* (/ 2.0 (sin k)) (/ l (pow t_m 3.0))) (+ 2.0 t_2)))
(/
2.0
(*
(+ 1.0 (+ t_2 1.0))
(* (tan k) (* (sin k) (* (/ (pow t_m 2.0) l) (/ t_m 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 = pow((k / t_m), 2.0);
double t_3 = l / tan(k);
double tmp;
if (t_m <= 3.2e-55) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * t_3;
} else if (t_m <= 1.12e+74) {
tmp = t_3 * (((2.0 / sin(k)) * (l / pow(t_m, 3.0))) / (2.0 + t_2));
} else {
tmp = 2.0 / ((1.0 + (t_2 + 1.0)) * (tan(k) * (sin(k) * ((pow(t_m, 2.0) / l) * (t_m / 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) :: t_3
real(8) :: tmp
t_2 = (k / t_m) ** 2.0d0
t_3 = l / tan(k)
if (t_m <= 3.2d-55) then
tmp = (2.0d0 * (l / ((k ** 2.0d0) * (t_m * sin(k))))) * t_3
else if (t_m <= 1.12d+74) then
tmp = t_3 * (((2.0d0 / sin(k)) * (l / (t_m ** 3.0d0))) / (2.0d0 + t_2))
else
tmp = 2.0d0 / ((1.0d0 + (t_2 + 1.0d0)) * (tan(k) * (sin(k) * (((t_m ** 2.0d0) / l) * (t_m / 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 = Math.pow((k / t_m), 2.0);
double t_3 = l / Math.tan(k);
double tmp;
if (t_m <= 3.2e-55) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * t_3;
} else if (t_m <= 1.12e+74) {
tmp = t_3 * (((2.0 / Math.sin(k)) * (l / Math.pow(t_m, 3.0))) / (2.0 + t_2));
} else {
tmp = 2.0 / ((1.0 + (t_2 + 1.0)) * (Math.tan(k) * (Math.sin(k) * ((Math.pow(t_m, 2.0) / l) * (t_m / 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 = math.pow((k / t_m), 2.0) t_3 = l / math.tan(k) tmp = 0 if t_m <= 3.2e-55: tmp = (2.0 * (l / (math.pow(k, 2.0) * (t_m * math.sin(k))))) * t_3 elif t_m <= 1.12e+74: tmp = t_3 * (((2.0 / math.sin(k)) * (l / math.pow(t_m, 3.0))) / (2.0 + t_2)) else: tmp = 2.0 / ((1.0 + (t_2 + 1.0)) * (math.tan(k) * (math.sin(k) * ((math.pow(t_m, 2.0) / l) * (t_m / 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(k / t_m) ^ 2.0 t_3 = Float64(l / tan(k)) tmp = 0.0 if (t_m <= 3.2e-55) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * t_3); elseif (t_m <= 1.12e+74) tmp = Float64(t_3 * Float64(Float64(Float64(2.0 / sin(k)) * Float64(l / (t_m ^ 3.0))) / 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(Float64((t_m ^ 2.0) / l) * Float64(t_m / 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 = (k / t_m) ^ 2.0; t_3 = l / tan(k); tmp = 0.0; if (t_m <= 3.2e-55) tmp = (2.0 * (l / ((k ^ 2.0) * (t_m * sin(k))))) * t_3; elseif (t_m <= 1.12e+74) tmp = t_3 * (((2.0 / sin(k)) * (l / (t_m ^ 3.0))) / (2.0 + t_2)); else tmp = 2.0 / ((1.0 + (t_2 + 1.0)) * (tan(k) * (sin(k) * (((t_m ^ 2.0) / l) * (t_m / 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[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$3 = N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 3.2e-55], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$3), $MachinePrecision], If[LessEqual[t$95$m, 1.12e+74], N[(t$95$3 * N[(N[(N[(2.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 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[(N[(N[Power[t$95$m, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t$95$m / l), $MachinePrecision]), $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_3 := \frac{\ell}{\tan k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 3.2 \cdot 10^{-55}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot t\_3\\
\mathbf{elif}\;t\_m \leq 1.12 \cdot 10^{+74}:\\
\;\;\;\;t\_3 \cdot \frac{\frac{2}{\sin k} \cdot \frac{\ell}{{t\_m}^{3}}}{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(\frac{{t\_m}^{2}}{\ell} \cdot \frac{t\_m}{\ell}\right)\right)\right)}\\
\end{array}
\end{array}
\end{array}
if t < 3.2000000000000001e-55Initial program 47.0%
Simplified52.9%
associate-/l/52.9%
associate-/r/54.5%
times-frac54.4%
div-inv54.3%
clear-num54.4%
Applied egg-rr54.4%
Taylor expanded in k around inf 67.6%
if 3.2000000000000001e-55 < t < 1.12000000000000003e74Initial program 80.8%
Simplified83.3%
associate-/l/83.2%
associate-/r/87.3%
times-frac91.6%
div-inv91.6%
clear-num91.7%
Applied egg-rr91.7%
if 1.12000000000000003e74 < t Initial program 77.4%
unpow377.4%
times-frac91.2%
pow291.2%
Applied egg-rr91.2%
Final simplification73.7%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (/ l (tan k))))
(*
t_s
(if (<= t_m 3.5e-54)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) t_2)
(if (<= t_m 2.85e+87)
(*
t_2
(/
(* (/ 2.0 (sin k)) (/ l (pow t_m 3.0)))
(+ 2.0 (pow (/ k t_m) 2.0))))
(/
2.0
(pow
(* (hypot 1.0 (hypot 1.0 (/ k t_m))) (* k (/ (pow t_m 1.5) l)))
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 t_2 = l / tan(k);
double tmp;
if (t_m <= 3.5e-54) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * t_2;
} else if (t_m <= 2.85e+87) {
tmp = t_2 * (((2.0 / sin(k)) * (l / pow(t_m, 3.0))) / (2.0 + pow((k / t_m), 2.0)));
} else {
tmp = 2.0 / pow((hypot(1.0, hypot(1.0, (k / t_m))) * (k * (pow(t_m, 1.5) / l))), 2.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 = l / Math.tan(k);
double tmp;
if (t_m <= 3.5e-54) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * t_2;
} else if (t_m <= 2.85e+87) {
tmp = t_2 * (((2.0 / Math.sin(k)) * (l / Math.pow(t_m, 3.0))) / (2.0 + Math.pow((k / t_m), 2.0)));
} else {
tmp = 2.0 / Math.pow((Math.hypot(1.0, Math.hypot(1.0, (k / t_m))) * (k * (Math.pow(t_m, 1.5) / l))), 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): t_2 = l / math.tan(k) tmp = 0 if t_m <= 3.5e-54: tmp = (2.0 * (l / (math.pow(k, 2.0) * (t_m * math.sin(k))))) * t_2 elif t_m <= 2.85e+87: tmp = t_2 * (((2.0 / math.sin(k)) * (l / math.pow(t_m, 3.0))) / (2.0 + math.pow((k / t_m), 2.0))) else: tmp = 2.0 / math.pow((math.hypot(1.0, math.hypot(1.0, (k / t_m))) * (k * (math.pow(t_m, 1.5) / l))), 2.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(l / tan(k)) tmp = 0.0 if (t_m <= 3.5e-54) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * t_2); elseif (t_m <= 2.85e+87) tmp = Float64(t_2 * Float64(Float64(Float64(2.0 / sin(k)) * Float64(l / (t_m ^ 3.0))) / Float64(2.0 + (Float64(k / t_m) ^ 2.0)))); else tmp = Float64(2.0 / (Float64(hypot(1.0, hypot(1.0, Float64(k / t_m))) * Float64(k * Float64((t_m ^ 1.5) / l))) ^ 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) t_2 = l / tan(k); tmp = 0.0; if (t_m <= 3.5e-54) tmp = (2.0 * (l / ((k ^ 2.0) * (t_m * sin(k))))) * t_2; elseif (t_m <= 2.85e+87) tmp = t_2 * (((2.0 / sin(k)) * (l / (t_m ^ 3.0))) / (2.0 + ((k / t_m) ^ 2.0))); else tmp = 2.0 / ((hypot(1.0, hypot(1.0, (k / t_m))) * (k * ((t_m ^ 1.5) / l))) ^ 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_] := Block[{t$95$2 = N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 3.5e-54], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision], If[LessEqual[t$95$m, 2.85e+87], N[(t$95$2 * N[(N[(N[(2.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[Sqrt[1.0 ^ 2 + N[Sqrt[1.0 ^ 2 + N[(k / t$95$m), $MachinePrecision] ^ 2], $MachinePrecision] ^ 2], $MachinePrecision] * N[(k * N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.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{\ell}{\tan k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 3.5 \cdot 10^{-54}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot t\_2\\
\mathbf{elif}\;t\_m \leq 2.85 \cdot 10^{+87}:\\
\;\;\;\;t\_2 \cdot \frac{\frac{2}{\sin k} \cdot \frac{\ell}{{t\_m}^{3}}}{2 + {\left(\frac{k}{t\_m}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\mathsf{hypot}\left(1, \mathsf{hypot}\left(1, \frac{k}{t\_m}\right)\right) \cdot \left(k \cdot \frac{{t\_m}^{1.5}}{\ell}\right)\right)}^{2}}\\
\end{array}
\end{array}
\end{array}
if t < 3.49999999999999982e-54Initial program 47.0%
Simplified52.9%
associate-/l/52.9%
associate-/r/54.5%
times-frac54.4%
div-inv54.3%
clear-num54.4%
Applied egg-rr54.4%
Taylor expanded in k around inf 67.6%
if 3.49999999999999982e-54 < t < 2.85000000000000019e87Initial program 80.5%
Simplified85.2%
associate-/l/85.2%
associate-/r/88.1%
times-frac91.2%
div-inv91.2%
clear-num91.3%
Applied egg-rr91.3%
if 2.85000000000000019e87 < t Initial program 76.8%
+-commutative76.8%
associate-+r+76.8%
metadata-eval76.8%
associate-*l*61.9%
associate-/r*65.3%
add-sqr-sqrt56.5%
pow256.5%
Applied egg-rr70.4%
Taylor expanded in k around 0 85.9%
Final simplification73.0%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 0.00012)
(/
2.0
(pow (* (hypot 1.0 (hypot 1.0 (/ k t_m))) (* k (/ (pow t_m 1.5) l))) 2.0))
(if (<= k 3.4e+151)
(* (/ l (tan k)) (/ (* 2.0 l) (* (sin k) (* t_m (pow k 2.0)))))
(/ 2.0 (pow (* (* k (/ (sin k) l)) (sqrt (/ t_m (cos 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 (k <= 0.00012) {
tmp = 2.0 / pow((hypot(1.0, hypot(1.0, (k / t_m))) * (k * (pow(t_m, 1.5) / l))), 2.0);
} else if (k <= 3.4e+151) {
tmp = (l / tan(k)) * ((2.0 * l) / (sin(k) * (t_m * pow(k, 2.0))));
} else {
tmp = 2.0 / pow(((k * (sin(k) / l)) * sqrt((t_m / cos(k)))), 2.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 (k <= 0.00012) {
tmp = 2.0 / Math.pow((Math.hypot(1.0, Math.hypot(1.0, (k / t_m))) * (k * (Math.pow(t_m, 1.5) / l))), 2.0);
} else if (k <= 3.4e+151) {
tmp = (l / Math.tan(k)) * ((2.0 * l) / (Math.sin(k) * (t_m * Math.pow(k, 2.0))));
} else {
tmp = 2.0 / Math.pow(((k * (Math.sin(k) / l)) * Math.sqrt((t_m / Math.cos(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 k <= 0.00012: tmp = 2.0 / math.pow((math.hypot(1.0, math.hypot(1.0, (k / t_m))) * (k * (math.pow(t_m, 1.5) / l))), 2.0) elif k <= 3.4e+151: tmp = (l / math.tan(k)) * ((2.0 * l) / (math.sin(k) * (t_m * math.pow(k, 2.0)))) else: tmp = 2.0 / math.pow(((k * (math.sin(k) / l)) * math.sqrt((t_m / math.cos(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 (k <= 0.00012) tmp = Float64(2.0 / (Float64(hypot(1.0, hypot(1.0, Float64(k / t_m))) * Float64(k * Float64((t_m ^ 1.5) / l))) ^ 2.0)); elseif (k <= 3.4e+151) tmp = Float64(Float64(l / tan(k)) * Float64(Float64(2.0 * l) / Float64(sin(k) * Float64(t_m * (k ^ 2.0))))); else tmp = Float64(2.0 / (Float64(Float64(k * Float64(sin(k) / l)) * sqrt(Float64(t_m / cos(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 (k <= 0.00012) tmp = 2.0 / ((hypot(1.0, hypot(1.0, (k / t_m))) * (k * ((t_m ^ 1.5) / l))) ^ 2.0); elseif (k <= 3.4e+151) tmp = (l / tan(k)) * ((2.0 * l) / (sin(k) * (t_m * (k ^ 2.0)))); else tmp = 2.0 / (((k * (sin(k) / l)) * sqrt((t_m / cos(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[k, 0.00012], N[(2.0 / N[Power[N[(N[Sqrt[1.0 ^ 2 + N[Sqrt[1.0 ^ 2 + N[(k / t$95$m), $MachinePrecision] ^ 2], $MachinePrecision] ^ 2], $MachinePrecision] * N[(k * N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 3.4e+151], N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 * l), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[(t$95$m * N[Power[k, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[(k * N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t$95$m / N[Cos[k], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.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}\;k \leq 0.00012:\\
\;\;\;\;\frac{2}{{\left(\mathsf{hypot}\left(1, \mathsf{hypot}\left(1, \frac{k}{t\_m}\right)\right) \cdot \left(k \cdot \frac{{t\_m}^{1.5}}{\ell}\right)\right)}^{2}}\\
\mathbf{elif}\;k \leq 3.4 \cdot 10^{+151}:\\
\;\;\;\;\frac{\ell}{\tan k} \cdot \frac{2 \cdot \ell}{\sin k \cdot \left(t\_m \cdot {k}^{2}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\left(k \cdot \frac{\sin k}{\ell}\right) \cdot \sqrt{\frac{t\_m}{\cos k}}\right)}^{2}}\\
\end{array}
\end{array}
if k < 1.20000000000000003e-4Initial program 59.6%
+-commutative59.6%
associate-+r+59.6%
metadata-eval59.6%
associate-*l*54.3%
associate-/r*59.4%
add-sqr-sqrt30.8%
pow230.8%
Applied egg-rr36.0%
Taylor expanded in k around 0 38.7%
if 1.20000000000000003e-4 < k < 3.4e151Initial program 39.4%
Simplified44.5%
associate-/l/44.6%
associate-/r/44.6%
times-frac46.8%
div-inv46.7%
clear-num46.7%
Applied egg-rr46.7%
Taylor expanded in k around inf 83.7%
associate-*r/83.7%
associate-*r*83.7%
Simplified83.7%
if 3.4e151 < k Initial program 46.5%
+-commutative46.5%
associate-+r+46.5%
metadata-eval46.5%
associate-*l*46.5%
associate-/r*55.2%
add-sqr-sqrt42.5%
pow242.5%
Applied egg-rr41.5%
Taylor expanded in k around inf 62.2%
associate-/l*62.2%
Simplified62.2%
Final simplification48.1%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (/ l (tan k))))
(*
t_s
(if (<= t_m 6e-54)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) t_2)
(if (<= t_m 5.5e+102)
(*
t_2
(/ (* 2.0 (/ l (* k (pow t_m 3.0)))) (+ 2.0 (pow (/ k t_m) 2.0))))
(/ (pow l 2.0) (pow (* t_m (pow (cbrt k) 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 = l / tan(k);
double tmp;
if (t_m <= 6e-54) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * t_2;
} else if (t_m <= 5.5e+102) {
tmp = t_2 * ((2.0 * (l / (k * pow(t_m, 3.0)))) / (2.0 + pow((k / t_m), 2.0)));
} else {
tmp = pow(l, 2.0) / pow((t_m * pow(cbrt(k), 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 = l / Math.tan(k);
double tmp;
if (t_m <= 6e-54) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * t_2;
} else if (t_m <= 5.5e+102) {
tmp = t_2 * ((2.0 * (l / (k * Math.pow(t_m, 3.0)))) / (2.0 + Math.pow((k / t_m), 2.0)));
} else {
tmp = Math.pow(l, 2.0) / Math.pow((t_m * Math.pow(Math.cbrt(k), 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 = Float64(l / tan(k)) tmp = 0.0 if (t_m <= 6e-54) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * t_2); elseif (t_m <= 5.5e+102) tmp = Float64(t_2 * Float64(Float64(2.0 * Float64(l / Float64(k * (t_m ^ 3.0)))) / Float64(2.0 + (Float64(k / t_m) ^ 2.0)))); else tmp = Float64((l ^ 2.0) / (Float64(t_m * (cbrt(k) ^ 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[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 6e-54], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision], If[LessEqual[t$95$m, 5.5e+102], N[(t$95$2 * N[(N[(2.0 * N[(l / N[(k * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[Power[N[(k / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[N[(t$95$m * N[Power[N[Power[k, 1/3], $MachinePrecision], 2.0], $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{\ell}{\tan k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 6 \cdot 10^{-54}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot t\_2\\
\mathbf{elif}\;t\_m \leq 5.5 \cdot 10^{+102}:\\
\;\;\;\;t\_2 \cdot \frac{2 \cdot \frac{\ell}{k \cdot {t\_m}^{3}}}{2 + {\left(\frac{k}{t\_m}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\ell}^{2}}{{\left(t\_m \cdot {\left(\sqrt[3]{k}\right)}^{2}\right)}^{3}}\\
\end{array}
\end{array}
\end{array}
if t < 6.00000000000000018e-54Initial program 47.0%
Simplified52.9%
associate-/l/52.9%
associate-/r/54.5%
times-frac54.4%
div-inv54.3%
clear-num54.4%
Applied egg-rr54.4%
Taylor expanded in k around inf 67.6%
if 6.00000000000000018e-54 < t < 5.49999999999999981e102Initial program 82.7%
Simplified86.8%
associate-/l/86.8%
associate-/r/89.4%
times-frac92.2%
div-inv92.2%
clear-num92.2%
Applied egg-rr92.2%
Taylor expanded in k around 0 84.5%
if 5.49999999999999981e102 < t Initial program 73.7%
Simplified77.6%
Taylor expanded in k around 0 56.8%
add-cube-cbrt56.8%
pow356.8%
*-commutative56.8%
cbrt-prod56.8%
unpow356.8%
add-cbrt-cube57.3%
unpow257.3%
cbrt-prod73.9%
pow273.9%
Applied egg-rr73.9%
Final simplification70.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 6.4e-54)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) (/ l (tan k)))
(/ 2.0 (pow (* (* k (/ (sqrt 2.0) l)) (fabs (pow t_m 1.5))) 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.4e-54) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * (l / tan(k));
} else {
tmp = 2.0 / pow(((k * (sqrt(2.0) / l)) * fabs(pow(t_m, 1.5))), 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.4d-54) then
tmp = (2.0d0 * (l / ((k ** 2.0d0) * (t_m * sin(k))))) * (l / tan(k))
else
tmp = 2.0d0 / (((k * (sqrt(2.0d0) / l)) * abs((t_m ** 1.5d0))) ** 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.4e-54) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * (l / Math.tan(k));
} else {
tmp = 2.0 / Math.pow(((k * (Math.sqrt(2.0) / l)) * Math.abs(Math.pow(t_m, 1.5))), 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.4e-54: tmp = (2.0 * (l / (math.pow(k, 2.0) * (t_m * math.sin(k))))) * (l / math.tan(k)) else: tmp = 2.0 / math.pow(((k * (math.sqrt(2.0) / l)) * math.fabs(math.pow(t_m, 1.5))), 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.4e-54) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * Float64(l / tan(k))); else tmp = Float64(2.0 / (Float64(Float64(k * Float64(sqrt(2.0) / l)) * abs((t_m ^ 1.5))) ^ 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.4e-54) tmp = (2.0 * (l / ((k ^ 2.0) * (t_m * sin(k))))) * (l / tan(k)); else tmp = 2.0 / (((k * (sqrt(2.0) / l)) * abs((t_m ^ 1.5))) ^ 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.4e-54], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[(k * N[(N[Sqrt[2.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[Abs[N[Power[t$95$m, 1.5], $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.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 6.4 \cdot 10^{-54}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot \frac{\ell}{\tan k}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\left(k \cdot \frac{\sqrt{2}}{\ell}\right) \cdot \left|{t\_m}^{1.5}\right|\right)}^{2}}\\
\end{array}
\end{array}
if t < 6.39999999999999997e-54Initial program 47.0%
Simplified52.9%
associate-/l/52.9%
associate-/r/54.5%
times-frac54.4%
div-inv54.3%
clear-num54.4%
Applied egg-rr54.4%
Taylor expanded in k around inf 67.6%
if 6.39999999999999997e-54 < t Initial program 78.6%
+-commutative78.6%
associate-+r+78.6%
metadata-eval78.6%
associate-*l*69.5%
associate-/r*73.5%
add-sqr-sqrt67.4%
pow267.4%
Applied egg-rr81.4%
Taylor expanded in k around 0 81.6%
associate-/l*81.6%
metadata-eval81.6%
pow-sqr81.6%
rem-sqrt-square84.5%
Simplified84.5%
Final simplification72.0%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (/ l (tan k))))
(*
t_s
(if (<= t_m 5.5e-54)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) t_2)
(*
t_2
(/ (* 2.0 (/ l (* k (pow t_m 3.0)))) (+ 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 t_2 = l / tan(k);
double tmp;
if (t_m <= 5.5e-54) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * t_2;
} else {
tmp = t_2 * ((2.0 * (l / (k * pow(t_m, 3.0)))) / (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) :: t_2
real(8) :: tmp
t_2 = l / tan(k)
if (t_m <= 5.5d-54) then
tmp = (2.0d0 * (l / ((k ** 2.0d0) * (t_m * sin(k))))) * t_2
else
tmp = t_2 * ((2.0d0 * (l / (k * (t_m ** 3.0d0)))) / (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 t_2 = l / Math.tan(k);
double tmp;
if (t_m <= 5.5e-54) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * t_2;
} else {
tmp = t_2 * ((2.0 * (l / (k * Math.pow(t_m, 3.0)))) / (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): t_2 = l / math.tan(k) tmp = 0 if t_m <= 5.5e-54: tmp = (2.0 * (l / (math.pow(k, 2.0) * (t_m * math.sin(k))))) * t_2 else: tmp = t_2 * ((2.0 * (l / (k * math.pow(t_m, 3.0)))) / (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) t_2 = Float64(l / tan(k)) tmp = 0.0 if (t_m <= 5.5e-54) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * t_2); else tmp = Float64(t_2 * Float64(Float64(2.0 * Float64(l / Float64(k * (t_m ^ 3.0)))) / 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) t_2 = l / tan(k); tmp = 0.0; if (t_m <= 5.5e-54) tmp = (2.0 * (l / ((k ^ 2.0) * (t_m * sin(k))))) * t_2; else tmp = t_2 * ((2.0 * (l / (k * (t_m ^ 3.0)))) / (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_] := Block[{t$95$2 = N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 5.5e-54], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision], N[(t$95$2 * N[(N[(2.0 * N[(l / N[(k * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 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)
\\
\begin{array}{l}
t_2 := \frac{\ell}{\tan k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 5.5 \cdot 10^{-54}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_2 \cdot \frac{2 \cdot \frac{\ell}{k \cdot {t\_m}^{3}}}{2 + {\left(\frac{k}{t\_m}\right)}^{2}}\\
\end{array}
\end{array}
\end{array}
if t < 5.50000000000000046e-54Initial program 47.0%
Simplified52.9%
associate-/l/52.9%
associate-/r/54.5%
times-frac54.4%
div-inv54.3%
clear-num54.4%
Applied egg-rr54.4%
Taylor expanded in k around inf 67.6%
if 5.50000000000000046e-54 < t Initial program 78.6%
Simplified82.6%
associate-/l/82.6%
associate-/r/84.0%
times-frac85.4%
div-inv85.4%
clear-num85.5%
Applied egg-rr85.5%
Taylor expanded in k around 0 81.3%
Final simplification71.1%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(let* ((t_2 (/ l (tan k))))
(*
t_s
(if (<= t_m 6.4e-54)
(* (* 2.0 (/ l (* (pow k 2.0) (* t_m (sin k))))) t_2)
(* t_2 (/ l (* 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 t_2 = l / tan(k);
double tmp;
if (t_m <= 6.4e-54) {
tmp = (2.0 * (l / (pow(k, 2.0) * (t_m * sin(k))))) * t_2;
} else {
tmp = t_2 * (l / (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) :: t_2
real(8) :: tmp
t_2 = l / tan(k)
if (t_m <= 6.4d-54) then
tmp = (2.0d0 * (l / ((k ** 2.0d0) * (t_m * sin(k))))) * t_2
else
tmp = t_2 * (l / (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 t_2 = l / Math.tan(k);
double tmp;
if (t_m <= 6.4e-54) {
tmp = (2.0 * (l / (Math.pow(k, 2.0) * (t_m * Math.sin(k))))) * t_2;
} else {
tmp = t_2 * (l / (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): t_2 = l / math.tan(k) tmp = 0 if t_m <= 6.4e-54: tmp = (2.0 * (l / (math.pow(k, 2.0) * (t_m * math.sin(k))))) * t_2 else: tmp = t_2 * (l / (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) t_2 = Float64(l / tan(k)) tmp = 0.0 if (t_m <= 6.4e-54) tmp = Float64(Float64(2.0 * Float64(l / Float64((k ^ 2.0) * Float64(t_m * sin(k))))) * t_2); else tmp = Float64(t_2 * Float64(l / Float64(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) t_2 = l / tan(k); tmp = 0.0; if (t_m <= 6.4e-54) tmp = (2.0 * (l / ((k ^ 2.0) * (t_m * sin(k))))) * t_2; else tmp = t_2 * (l / (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_] := Block[{t$95$2 = N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 6.4e-54], N[(N[(2.0 * N[(l / N[(N[Power[k, 2.0], $MachinePrecision] * N[(t$95$m * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision], N[(t$95$2 * N[(l / N[(k * N[Power[t$95$m, 3.0], $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 := \frac{\ell}{\tan k}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 6.4 \cdot 10^{-54}:\\
\;\;\;\;\left(2 \cdot \frac{\ell}{{k}^{2} \cdot \left(t\_m \cdot \sin k\right)}\right) \cdot t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_2 \cdot \frac{\ell}{k \cdot {t\_m}^{3}}\\
\end{array}
\end{array}
\end{array}
if t < 6.39999999999999997e-54Initial program 47.0%
Simplified52.9%
associate-/l/52.9%
associate-/r/54.5%
times-frac54.4%
div-inv54.3%
clear-num54.4%
Applied egg-rr54.4%
Taylor expanded in k around inf 67.6%
if 6.39999999999999997e-54 < t Initial program 78.6%
Simplified82.6%
associate-/l/82.6%
associate-/r/84.0%
times-frac85.4%
div-inv85.4%
clear-num85.5%
Applied egg-rr85.5%
Taylor expanded in k around 0 78.3%
Final simplification70.4%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 5.6e+139)
(* (/ l (tan k)) (/ l (* (sin k) (pow t_m 3.0))))
(/ (* 2.0 (* (cos k) (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) {
double tmp;
if (k <= 5.6e+139) {
tmp = (l / tan(k)) * (l / (sin(k) * pow(t_m, 3.0)));
} else {
tmp = (2.0 * (cos(k) * pow(l, 2.0))) / (t_m * pow(k, 4.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 (k <= 5.6d+139) then
tmp = (l / tan(k)) * (l / (sin(k) * (t_m ** 3.0d0)))
else
tmp = (2.0d0 * (cos(k) * (l ** 2.0d0))) / (t_m * (k ** 4.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 (k <= 5.6e+139) {
tmp = (l / Math.tan(k)) * (l / (Math.sin(k) * Math.pow(t_m, 3.0)));
} else {
tmp = (2.0 * (Math.cos(k) * Math.pow(l, 2.0))) / (t_m * Math.pow(k, 4.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 k <= 5.6e+139: tmp = (l / math.tan(k)) * (l / (math.sin(k) * math.pow(t_m, 3.0))) else: tmp = (2.0 * (math.cos(k) * math.pow(l, 2.0))) / (t_m * math.pow(k, 4.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 (k <= 5.6e+139) tmp = Float64(Float64(l / tan(k)) * Float64(l / Float64(sin(k) * (t_m ^ 3.0)))); else tmp = Float64(Float64(2.0 * Float64(cos(k) * (l ^ 2.0))) / Float64(t_m * (k ^ 4.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 (k <= 5.6e+139) tmp = (l / tan(k)) * (l / (sin(k) * (t_m ^ 3.0))); else tmp = (2.0 * (cos(k) * (l ^ 2.0))) / (t_m * (k ^ 4.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[k, 5.6e+139], N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(N[Cos[k], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[k, 4.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}\;k \leq 5.6 \cdot 10^{+139}:\\
\;\;\;\;\frac{\ell}{\tan k} \cdot \frac{\ell}{\sin k \cdot {t\_m}^{3}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \left(\cos k \cdot {\ell}^{2}\right)}{t\_m \cdot {k}^{4}}\\
\end{array}
\end{array}
if k < 5.5999999999999997e139Initial program 56.5%
Simplified61.6%
associate-/l/61.6%
associate-/r/63.3%
times-frac63.7%
div-inv63.6%
clear-num63.7%
Applied egg-rr63.7%
Taylor expanded in t around inf 66.2%
if 5.5999999999999997e139 < k Initial program 42.9%
Simplified51.0%
associate-/l/51.0%
associate-/r/50.9%
times-frac50.9%
div-inv51.0%
clear-num50.9%
Applied egg-rr50.9%
Taylor expanded in k around inf 55.0%
associate-*r/55.0%
*-commutative55.0%
associate-*r*55.0%
*-commutative55.0%
*-commutative55.0%
*-commutative55.0%
Simplified55.0%
Taylor expanded in k around 0 55.0%
Final simplification65.0%
t_m = (fabs.f64 t)
t_s = (copysign.f64 1 t)
(FPCore (t_s t_m l k)
:precision binary64
(*
t_s
(if (<= k 2.1e+145)
(* (/ l (tan k)) (/ l (* (sin k) (pow t_m 3.0))))
(/ 2.0 (/ (* t_m (pow k 4.0)) (pow l 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 (k <= 2.1e+145) {
tmp = (l / tan(k)) * (l / (sin(k) * pow(t_m, 3.0)));
} else {
tmp = 2.0 / ((t_m * pow(k, 4.0)) / pow(l, 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 (k <= 2.1d+145) then
tmp = (l / tan(k)) * (l / (sin(k) * (t_m ** 3.0d0)))
else
tmp = 2.0d0 / ((t_m * (k ** 4.0d0)) / (l ** 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 (k <= 2.1e+145) {
tmp = (l / Math.tan(k)) * (l / (Math.sin(k) * Math.pow(t_m, 3.0)));
} else {
tmp = 2.0 / ((t_m * Math.pow(k, 4.0)) / Math.pow(l, 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 k <= 2.1e+145: tmp = (l / math.tan(k)) * (l / (math.sin(k) * math.pow(t_m, 3.0))) else: tmp = 2.0 / ((t_m * math.pow(k, 4.0)) / math.pow(l, 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 (k <= 2.1e+145) tmp = Float64(Float64(l / tan(k)) * Float64(l / Float64(sin(k) * (t_m ^ 3.0)))); else tmp = Float64(2.0 / Float64(Float64(t_m * (k ^ 4.0)) / (l ^ 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 (k <= 2.1e+145) tmp = (l / tan(k)) * (l / (sin(k) * (t_m ^ 3.0))); else tmp = 2.0 / ((t_m * (k ^ 4.0)) / (l ^ 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[k, 2.1e+145], N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[(N[Sin[k], $MachinePrecision] * N[Power[t$95$m, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t$95$m * N[Power[k, 4.0], $MachinePrecision]), $MachinePrecision] / N[Power[l, 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}\;k \leq 2.1 \cdot 10^{+145}:\\
\;\;\;\;\frac{\ell}{\tan k} \cdot \frac{\ell}{\sin k \cdot {t\_m}^{3}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{t\_m \cdot {k}^{4}}{{\ell}^{2}}}\\
\end{array}
\end{array}
if k < 2.09999999999999989e145Initial program 56.5%
Simplified61.6%
associate-/l/61.6%
associate-/r/63.3%
times-frac63.7%
div-inv63.6%
clear-num63.7%
Applied egg-rr63.7%
Taylor expanded in t around inf 66.2%
if 2.09999999999999989e145 < k Initial program 42.9%
Taylor expanded in t around 0 55.0%
associate-*r*55.0%
times-frac55.0%
Simplified55.0%
Taylor expanded in k around 0 55.0%
Final simplification65.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 1.45e-58)
(/ 2.0 (* (pow k 4.0) (/ t_m (pow l 2.0))))
(* (/ l (tan k)) (/ l (* 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 <= 1.45e-58) {
tmp = 2.0 / (pow(k, 4.0) * (t_m / pow(l, 2.0)));
} else {
tmp = (l / tan(k)) * (l / (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 <= 1.45d-58) then
tmp = 2.0d0 / ((k ** 4.0d0) * (t_m / (l ** 2.0d0)))
else
tmp = (l / tan(k)) * (l / (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 <= 1.45e-58) {
tmp = 2.0 / (Math.pow(k, 4.0) * (t_m / Math.pow(l, 2.0)));
} else {
tmp = (l / Math.tan(k)) * (l / (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 <= 1.45e-58: tmp = 2.0 / (math.pow(k, 4.0) * (t_m / math.pow(l, 2.0))) else: tmp = (l / math.tan(k)) * (l / (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 <= 1.45e-58) tmp = Float64(2.0 / Float64((k ^ 4.0) * Float64(t_m / (l ^ 2.0)))); else tmp = Float64(Float64(l / tan(k)) * Float64(l / Float64(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 <= 1.45e-58) tmp = 2.0 / ((k ^ 4.0) * (t_m / (l ^ 2.0))); else tmp = (l / tan(k)) * (l / (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, 1.45e-58], N[(2.0 / N[(N[Power[k, 4.0], $MachinePrecision] * N[(t$95$m / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l / N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(l / N[(k * N[Power[t$95$m, 3.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.45 \cdot 10^{-58}:\\
\;\;\;\;\frac{2}{{k}^{4} \cdot \frac{t\_m}{{\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell}{\tan k} \cdot \frac{\ell}{k \cdot {t\_m}^{3}}\\
\end{array}
\end{array}
if t < 1.44999999999999995e-58Initial program 45.8%
Taylor expanded in t around 0 62.0%
associate-*r*62.0%
times-frac65.0%
Simplified65.0%
Taylor expanded in k around 0 45.6%
associate-/l*47.8%
Simplified47.8%
if 1.44999999999999995e-58 < t Initial program 79.8%
Simplified83.6%
associate-/l/83.6%
associate-/r/84.9%
times-frac86.2%
div-inv86.2%
clear-num86.3%
Applied egg-rr86.3%
Taylor expanded in k around 0 78.1%
Final simplification56.1%
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.12e-38)
(/ 2.0 (* (pow k 4.0) (/ t_m (pow l 2.0))))
(/ (pow l 2.0) (* (pow t_m 3.0) (* k k))))))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.12e-38) {
tmp = 2.0 / (pow(k, 4.0) * (t_m / pow(l, 2.0)));
} else {
tmp = pow(l, 2.0) / (pow(t_m, 3.0) * (k * k));
}
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.12d-38) then
tmp = 2.0d0 / ((k ** 4.0d0) * (t_m / (l ** 2.0d0)))
else
tmp = (l ** 2.0d0) / ((t_m ** 3.0d0) * (k * k))
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.12e-38) {
tmp = 2.0 / (Math.pow(k, 4.0) * (t_m / Math.pow(l, 2.0)));
} else {
tmp = Math.pow(l, 2.0) / (Math.pow(t_m, 3.0) * (k * k));
}
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.12e-38: tmp = 2.0 / (math.pow(k, 4.0) * (t_m / math.pow(l, 2.0))) else: tmp = math.pow(l, 2.0) / (math.pow(t_m, 3.0) * (k * k)) 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.12e-38) tmp = Float64(2.0 / Float64((k ^ 4.0) * Float64(t_m / (l ^ 2.0)))); else tmp = Float64((l ^ 2.0) / Float64((t_m ^ 3.0) * Float64(k * k))); 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.12e-38) tmp = 2.0 / ((k ^ 4.0) * (t_m / (l ^ 2.0))); else tmp = (l ^ 2.0) / ((t_m ^ 3.0) * (k * k)); 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.12e-38], N[(2.0 / N[(N[Power[k, 4.0], $MachinePrecision] * N[(t$95$m / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Power[l, 2.0], $MachinePrecision] / N[(N[Power[t$95$m, 3.0], $MachinePrecision] * N[(k * k), $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.12 \cdot 10^{-38}:\\
\;\;\;\;\frac{2}{{k}^{4} \cdot \frac{t\_m}{{\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\ell}^{2}}{{t\_m}^{3} \cdot \left(k \cdot k\right)}\\
\end{array}
\end{array}
if t < 1.1200000000000001e-38Initial program 47.0%
Taylor expanded in t around 0 62.7%
associate-*r*62.6%
times-frac65.6%
Simplified65.6%
Taylor expanded in k around 0 45.7%
associate-/l*47.9%
Simplified47.9%
if 1.1200000000000001e-38 < t Initial program 79.4%
Simplified83.6%
Taylor expanded in k around 0 67.4%
unpow267.4%
Applied egg-rr67.4%
Final simplification52.8%
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 55.1%
Taylor expanded in t around 0 59.9%
associate-*r*59.9%
times-frac62.8%
Simplified62.8%
Taylor expanded in k around 0 46.5%
associate-/l*48.4%
Simplified48.4%
Taylor expanded in k around 0 46.5%
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 k 4.0) (/ t_m (pow l 2.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(k, 4.0) * (t_m / pow(l, 2.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 / ((k ** 4.0d0) * (t_m / (l ** 2.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(k, 4.0) * (t_m / Math.pow(l, 2.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(k, 4.0) * (t_m / math.pow(l, 2.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((k ^ 4.0) * Float64(t_m / (l ^ 2.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 / ((k ^ 4.0) * (t_m / (l ^ 2.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[k, 4.0], $MachinePrecision] * N[(t$95$m / N[Power[l, 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 \frac{2}{{k}^{4} \cdot \frac{t\_m}{{\ell}^{2}}}
\end{array}
Initial program 55.1%
Taylor expanded in t around 0 59.9%
associate-*r*59.9%
times-frac62.8%
Simplified62.8%
Taylor expanded in k around 0 46.5%
associate-/l*48.4%
Simplified48.4%
Final simplification48.4%
herbie shell --seed 2024039
(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))))