
(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 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (- (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) - 1.0));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t) ** 2.0d0)) - 1.0d0))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t), 2.0)) - 1.0));
}
def code(t, l, k): return 2.0 / ((((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * math.tan(k)) * ((1.0 + math.pow((k / t), 2.0)) - 1.0))
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) - 1.0))) end
function tmp = code(t, l, k) tmp = 2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t) ^ 2.0)) - 1.0)); end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)}
\end{array}
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (pow (cbrt l) -2.0))
(t_3 (* (sin k_m) (tan k_m)))
(t_4 (cbrt t_3)))
(*
t_s
(if (<= k_m 70.0)
(/ 2.0 (pow (* (sqrt t_3) (* k_m (/ (sqrt t_m) l))) 2.0))
(*
(* (/ (sqrt 2.0) k_m) (/ t_m (pow (* t_4 (* t_m t_2)) 2.0)))
(/ (/ 1.0 (* k_m (/ t_2 (sqrt 2.0)))) t_4))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = pow(cbrt(l), -2.0);
double t_3 = sin(k_m) * tan(k_m);
double t_4 = cbrt(t_3);
double tmp;
if (k_m <= 70.0) {
tmp = 2.0 / pow((sqrt(t_3) * (k_m * (sqrt(t_m) / l))), 2.0);
} else {
tmp = ((sqrt(2.0) / k_m) * (t_m / pow((t_4 * (t_m * t_2)), 2.0))) * ((1.0 / (k_m * (t_2 / sqrt(2.0)))) / t_4);
}
return t_s * tmp;
}
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = Math.pow(Math.cbrt(l), -2.0);
double t_3 = Math.sin(k_m) * Math.tan(k_m);
double t_4 = Math.cbrt(t_3);
double tmp;
if (k_m <= 70.0) {
tmp = 2.0 / Math.pow((Math.sqrt(t_3) * (k_m * (Math.sqrt(t_m) / l))), 2.0);
} else {
tmp = ((Math.sqrt(2.0) / k_m) * (t_m / Math.pow((t_4 * (t_m * t_2)), 2.0))) * ((1.0 / (k_m * (t_2 / Math.sqrt(2.0)))) / t_4);
}
return t_s * tmp;
}
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = cbrt(l) ^ -2.0 t_3 = Float64(sin(k_m) * tan(k_m)) t_4 = cbrt(t_3) tmp = 0.0 if (k_m <= 70.0) tmp = Float64(2.0 / (Float64(sqrt(t_3) * Float64(k_m * Float64(sqrt(t_m) / l))) ^ 2.0)); else tmp = Float64(Float64(Float64(sqrt(2.0) / k_m) * Float64(t_m / (Float64(t_4 * Float64(t_m * t_2)) ^ 2.0))) * Float64(Float64(1.0 / Float64(k_m * Float64(t_2 / sqrt(2.0)))) / t_4)); end return Float64(t_s * tmp) end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]}, Block[{t$95$3 = N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[Power[t$95$3, 1/3], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 70.0], N[(2.0 / N[Power[N[(N[Sqrt[t$95$3], $MachinePrecision] * N[(k$95$m * N[(N[Sqrt[t$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] * N[(t$95$m / N[Power[N[(t$95$4 * N[(t$95$m * t$95$2), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 / N[(k$95$m * N[(t$95$2 / N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$4), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := {\left(\sqrt[3]{\ell}\right)}^{-2}\\
t_3 := \sin k\_m \cdot \tan k\_m\\
t_4 := \sqrt[3]{t\_3}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 70:\\
\;\;\;\;\frac{2}{{\left(\sqrt{t\_3} \cdot \left(k\_m \cdot \frac{\sqrt{t\_m}}{\ell}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\sqrt{2}}{k\_m} \cdot \frac{t\_m}{{\left(t\_4 \cdot \left(t\_m \cdot t\_2\right)\right)}^{2}}\right) \cdot \frac{\frac{1}{k\_m \cdot \frac{t\_2}{\sqrt{2}}}}{t\_4}\\
\end{array}
\end{array}
\end{array}
if k < 70Initial program 40.7%
*-commutative40.7%
associate-/r*40.7%
Simplified47.5%
add-sqr-sqrt21.7%
pow221.7%
sqrt-prod19.1%
sqrt-div20.6%
sqrt-pow124.6%
metadata-eval24.6%
sqrt-prod15.0%
add-sqr-sqrt29.1%
Applied egg-rr29.1%
*-un-lft-identity29.1%
associate-/l/29.1%
+-rgt-identity29.1%
pow-prod-down33.7%
*-commutative33.7%
Applied egg-rr33.7%
*-lft-identity33.7%
associate-*l*33.7%
associate-*l/34.1%
Simplified34.1%
Taylor expanded in t around 0 43.2%
associate-*l/42.8%
associate-/l*43.2%
Simplified43.2%
if 70 < k Initial program 29.1%
*-commutative29.1%
associate-/r*29.1%
Simplified39.8%
add-sqr-sqrt39.7%
add-cube-cbrt39.7%
times-frac39.7%
Applied egg-rr68.1%
associate-/r/68.1%
associate-/r*68.1%
associate-/r/70.0%
Simplified70.0%
clear-num69.9%
inv-pow69.9%
div-inv69.9%
pow-flip70.0%
metadata-eval70.0%
associate-*l/70.0%
Applied egg-rr70.0%
unpow-170.0%
associate-/r/70.0%
*-commutative70.0%
times-frac70.0%
*-inverses70.0%
Simplified70.0%
associate-/l*74.2%
*-commutative74.2%
div-inv74.2%
pow-flip74.2%
metadata-eval74.2%
Applied egg-rr74.2%
Final simplification50.0%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (* (sin k_m) (tan k_m))) (t_3 (cbrt t_2)))
(*
t_s
(if (<= k_m 0.0152)
(/ 2.0 (pow (* (sqrt t_2) (* k_m (/ (sqrt t_m) l))) 2.0))
(if (<= k_m 8.2e+150)
(*
2.0
(/
1.0
(*
(pow k_m 2.0)
(* (/ (pow (sin k_m) 2.0) (pow l 2.0)) (/ t_m (cos k_m))))))
(*
(*
(* (sqrt 2.0) (/ t_m k_m))
(pow (* t_m (* t_3 (pow (cbrt l) -2.0))) -2.0))
(/
(/ (* t_m (/ (sqrt 2.0) k_m)) (/ t_m (pow (cbrt l) 2.0)))
t_3)))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = sin(k_m) * tan(k_m);
double t_3 = cbrt(t_2);
double tmp;
if (k_m <= 0.0152) {
tmp = 2.0 / pow((sqrt(t_2) * (k_m * (sqrt(t_m) / l))), 2.0);
} else if (k_m <= 8.2e+150) {
tmp = 2.0 * (1.0 / (pow(k_m, 2.0) * ((pow(sin(k_m), 2.0) / pow(l, 2.0)) * (t_m / cos(k_m)))));
} else {
tmp = ((sqrt(2.0) * (t_m / k_m)) * pow((t_m * (t_3 * pow(cbrt(l), -2.0))), -2.0)) * (((t_m * (sqrt(2.0) / k_m)) / (t_m / pow(cbrt(l), 2.0))) / t_3);
}
return t_s * tmp;
}
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = Math.sin(k_m) * Math.tan(k_m);
double t_3 = Math.cbrt(t_2);
double tmp;
if (k_m <= 0.0152) {
tmp = 2.0 / Math.pow((Math.sqrt(t_2) * (k_m * (Math.sqrt(t_m) / l))), 2.0);
} else if (k_m <= 8.2e+150) {
tmp = 2.0 * (1.0 / (Math.pow(k_m, 2.0) * ((Math.pow(Math.sin(k_m), 2.0) / Math.pow(l, 2.0)) * (t_m / Math.cos(k_m)))));
} else {
tmp = ((Math.sqrt(2.0) * (t_m / k_m)) * Math.pow((t_m * (t_3 * Math.pow(Math.cbrt(l), -2.0))), -2.0)) * (((t_m * (Math.sqrt(2.0) / k_m)) / (t_m / Math.pow(Math.cbrt(l), 2.0))) / t_3);
}
return t_s * tmp;
}
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = Float64(sin(k_m) * tan(k_m)) t_3 = cbrt(t_2) tmp = 0.0 if (k_m <= 0.0152) tmp = Float64(2.0 / (Float64(sqrt(t_2) * Float64(k_m * Float64(sqrt(t_m) / l))) ^ 2.0)); elseif (k_m <= 8.2e+150) tmp = Float64(2.0 * Float64(1.0 / Float64((k_m ^ 2.0) * Float64(Float64((sin(k_m) ^ 2.0) / (l ^ 2.0)) * Float64(t_m / cos(k_m)))))); else tmp = Float64(Float64(Float64(sqrt(2.0) * Float64(t_m / k_m)) * (Float64(t_m * Float64(t_3 * (cbrt(l) ^ -2.0))) ^ -2.0)) * Float64(Float64(Float64(t_m * Float64(sqrt(2.0) / k_m)) / Float64(t_m / (cbrt(l) ^ 2.0))) / t_3)); end return Float64(t_s * tmp) end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[Power[t$95$2, 1/3], $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 0.0152], N[(2.0 / N[Power[N[(N[Sqrt[t$95$2], $MachinePrecision] * N[(k$95$m * N[(N[Sqrt[t$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 8.2e+150], N[(2.0 * N[(1.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] * N[(t$95$m / k$95$m), $MachinePrecision]), $MachinePrecision] * N[Power[N[(t$95$m * N[(t$95$3 * N[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(t$95$m * N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision]), $MachinePrecision] / N[(t$95$m / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$3), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sin k\_m \cdot \tan k\_m\\
t_3 := \sqrt[3]{t\_2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 0.0152:\\
\;\;\;\;\frac{2}{{\left(\sqrt{t\_2} \cdot \left(k\_m \cdot \frac{\sqrt{t\_m}}{\ell}\right)\right)}^{2}}\\
\mathbf{elif}\;k\_m \leq 8.2 \cdot 10^{+150}:\\
\;\;\;\;2 \cdot \frac{1}{{k\_m}^{2} \cdot \left(\frac{{\sin k\_m}^{2}}{{\ell}^{2}} \cdot \frac{t\_m}{\cos k\_m}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\sqrt{2} \cdot \frac{t\_m}{k\_m}\right) \cdot {\left(t\_m \cdot \left(t\_3 \cdot {\left(\sqrt[3]{\ell}\right)}^{-2}\right)\right)}^{-2}\right) \cdot \frac{\frac{t\_m \cdot \frac{\sqrt{2}}{k\_m}}{\frac{t\_m}{{\left(\sqrt[3]{\ell}\right)}^{2}}}}{t\_3}\\
\end{array}
\end{array}
\end{array}
if k < 0.0152Initial program 40.7%
*-commutative40.7%
associate-/r*40.7%
Simplified47.5%
add-sqr-sqrt21.7%
pow221.7%
sqrt-prod19.1%
sqrt-div20.6%
sqrt-pow124.6%
metadata-eval24.6%
sqrt-prod15.0%
add-sqr-sqrt29.1%
Applied egg-rr29.1%
*-un-lft-identity29.1%
associate-/l/29.1%
+-rgt-identity29.1%
pow-prod-down33.7%
*-commutative33.7%
Applied egg-rr33.7%
*-lft-identity33.7%
associate-*l*33.7%
associate-*l/34.1%
Simplified34.1%
Taylor expanded in t around 0 43.2%
associate-*l/42.8%
associate-/l*43.2%
Simplified43.2%
if 0.0152 < k < 8.19999999999999988e150Initial program 17.6%
Simplified28.0%
Taylor expanded in t around 0 79.9%
clear-num80.0%
inv-pow80.0%
Applied egg-rr80.0%
unpow-180.0%
associate-/l*83.2%
*-commutative83.2%
times-frac83.4%
Simplified83.4%
if 8.19999999999999988e150 < k Initial program 41.4%
*-commutative41.4%
associate-/r*41.4%
Simplified52.5%
add-sqr-sqrt52.5%
add-cube-cbrt52.5%
times-frac52.5%
Applied egg-rr70.0%
associate-/r/70.0%
associate-/r*70.0%
associate-/r/70.0%
Simplified70.0%
div-inv70.0%
associate-*l/70.0%
pow-flip70.0%
*-commutative70.0%
div-inv70.1%
pow-flip70.1%
metadata-eval70.1%
metadata-eval70.1%
Applied egg-rr70.1%
associate-/l*70.1%
associate-*r*70.1%
*-commutative70.1%
associate-*l*70.2%
Simplified70.2%
Final simplification50.6%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 0.024)
(/
2.0
(pow (* (sqrt (* (sin k_m) (tan k_m))) (* k_m (/ (sqrt t_m) l))) 2.0))
(if (<= k_m 2.5e+165)
(*
2.0
(/
1.0
(*
(pow k_m 2.0)
(* (/ (pow (sin k_m) 2.0) (pow l 2.0)) (/ t_m (cos k_m))))))
(/
2.0
(*
(sin k_m)
(* (tan k_m) (pow (* (/ (pow t_m 1.5) l) (/ k_m t_m)) 2.0))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 0.024) {
tmp = 2.0 / pow((sqrt((sin(k_m) * tan(k_m))) * (k_m * (sqrt(t_m) / l))), 2.0);
} else if (k_m <= 2.5e+165) {
tmp = 2.0 * (1.0 / (pow(k_m, 2.0) * ((pow(sin(k_m), 2.0) / pow(l, 2.0)) * (t_m / cos(k_m)))));
} else {
tmp = 2.0 / (sin(k_m) * (tan(k_m) * pow(((pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0)));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 0.024d0) then
tmp = 2.0d0 / ((sqrt((sin(k_m) * tan(k_m))) * (k_m * (sqrt(t_m) / l))) ** 2.0d0)
else if (k_m <= 2.5d+165) then
tmp = 2.0d0 * (1.0d0 / ((k_m ** 2.0d0) * (((sin(k_m) ** 2.0d0) / (l ** 2.0d0)) * (t_m / cos(k_m)))))
else
tmp = 2.0d0 / (sin(k_m) * (tan(k_m) * ((((t_m ** 1.5d0) / l) * (k_m / t_m)) ** 2.0d0)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 0.024) {
tmp = 2.0 / Math.pow((Math.sqrt((Math.sin(k_m) * Math.tan(k_m))) * (k_m * (Math.sqrt(t_m) / l))), 2.0);
} else if (k_m <= 2.5e+165) {
tmp = 2.0 * (1.0 / (Math.pow(k_m, 2.0) * ((Math.pow(Math.sin(k_m), 2.0) / Math.pow(l, 2.0)) * (t_m / Math.cos(k_m)))));
} else {
tmp = 2.0 / (Math.sin(k_m) * (Math.tan(k_m) * Math.pow(((Math.pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 0.024: tmp = 2.0 / math.pow((math.sqrt((math.sin(k_m) * math.tan(k_m))) * (k_m * (math.sqrt(t_m) / l))), 2.0) elif k_m <= 2.5e+165: tmp = 2.0 * (1.0 / (math.pow(k_m, 2.0) * ((math.pow(math.sin(k_m), 2.0) / math.pow(l, 2.0)) * (t_m / math.cos(k_m))))) else: tmp = 2.0 / (math.sin(k_m) * (math.tan(k_m) * math.pow(((math.pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 0.024) tmp = Float64(2.0 / (Float64(sqrt(Float64(sin(k_m) * tan(k_m))) * Float64(k_m * Float64(sqrt(t_m) / l))) ^ 2.0)); elseif (k_m <= 2.5e+165) tmp = Float64(2.0 * Float64(1.0 / Float64((k_m ^ 2.0) * Float64(Float64((sin(k_m) ^ 2.0) / (l ^ 2.0)) * Float64(t_m / cos(k_m)))))); else tmp = Float64(2.0 / Float64(sin(k_m) * Float64(tan(k_m) * (Float64(Float64((t_m ^ 1.5) / l) * Float64(k_m / t_m)) ^ 2.0)))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 0.024) tmp = 2.0 / ((sqrt((sin(k_m) * tan(k_m))) * (k_m * (sqrt(t_m) / l))) ^ 2.0); elseif (k_m <= 2.5e+165) tmp = 2.0 * (1.0 / ((k_m ^ 2.0) * (((sin(k_m) ^ 2.0) / (l ^ 2.0)) * (t_m / cos(k_m))))); else tmp = 2.0 / (sin(k_m) * (tan(k_m) * ((((t_m ^ 1.5) / l) * (k_m / t_m)) ^ 2.0))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 0.024], N[(2.0 / N[Power[N[(N[Sqrt[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(k$95$m * N[(N[Sqrt[t$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2.5e+165], N[(2.0 * N[(1.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(t$95$m / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Tan[k$95$m], $MachinePrecision] * N[Power[N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 0.024:\\
\;\;\;\;\frac{2}{{\left(\sqrt{\sin k\_m \cdot \tan k\_m} \cdot \left(k\_m \cdot \frac{\sqrt{t\_m}}{\ell}\right)\right)}^{2}}\\
\mathbf{elif}\;k\_m \leq 2.5 \cdot 10^{+165}:\\
\;\;\;\;2 \cdot \frac{1}{{k\_m}^{2} \cdot \left(\frac{{\sin k\_m}^{2}}{{\ell}^{2}} \cdot \frac{t\_m}{\cos k\_m}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\sin k\_m \cdot \left(\tan k\_m \cdot {\left(\frac{{t\_m}^{1.5}}{\ell} \cdot \frac{k\_m}{t\_m}\right)}^{2}\right)}\\
\end{array}
\end{array}
if k < 0.024Initial program 40.7%
*-commutative40.7%
associate-/r*40.7%
Simplified47.5%
add-sqr-sqrt21.7%
pow221.7%
sqrt-prod19.1%
sqrt-div20.6%
sqrt-pow124.6%
metadata-eval24.6%
sqrt-prod15.0%
add-sqr-sqrt29.1%
Applied egg-rr29.1%
*-un-lft-identity29.1%
associate-/l/29.1%
+-rgt-identity29.1%
pow-prod-down33.7%
*-commutative33.7%
Applied egg-rr33.7%
*-lft-identity33.7%
associate-*l*33.7%
associate-*l/34.1%
Simplified34.1%
Taylor expanded in t around 0 43.2%
associate-*l/42.8%
associate-/l*43.2%
Simplified43.2%
if 0.024 < k < 2.49999999999999985e165Initial program 19.1%
Simplified34.7%
Taylor expanded in t around 0 81.8%
clear-num81.8%
inv-pow81.8%
Applied egg-rr81.8%
unpow-181.8%
associate-/l*84.8%
*-commutative84.8%
times-frac85.0%
Simplified85.0%
if 2.49999999999999985e165 < k Initial program 42.4%
*-commutative42.4%
associate-/r*42.4%
Simplified46.5%
add-sqr-sqrt12.8%
pow212.8%
sqrt-prod12.8%
sqrt-div12.8%
sqrt-pow116.7%
metadata-eval16.7%
sqrt-prod8.3%
add-sqr-sqrt21.0%
Applied egg-rr21.0%
*-un-lft-identity21.0%
associate-/l/21.0%
+-rgt-identity21.0%
pow-prod-down21.0%
*-commutative21.0%
Applied egg-rr21.0%
*-lft-identity21.0%
associate-*l*21.0%
associate-*l/21.0%
Simplified21.0%
*-un-lft-identity21.0%
unpow-prod-down21.0%
pow221.0%
add-sqr-sqrt47.8%
associate-/l*47.9%
unpow-prod-down33.7%
pow233.7%
pow-prod-up57.8%
metadata-eval57.8%
Applied egg-rr57.8%
*-lft-identity57.8%
associate-*l*57.8%
unpow257.8%
metadata-eval57.8%
pow-sqr33.7%
swap-sqr47.8%
unpow247.8%
associate-*r/47.7%
associate-*l/47.8%
Simplified47.8%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 0.0195)
(/
2.0
(pow (* (sqrt (* (sin k_m) (tan k_m))) (* k_m (/ (sqrt t_m) l))) 2.0))
(if (<= k_m 2.5e+165)
(*
(/ 2.0 (pow k_m 2.0))
(/ (pow l 2.0) (* t_m (/ (pow (sin k_m) 2.0) (cos k_m)))))
(/
2.0
(*
(sin k_m)
(* (tan k_m) (pow (* (/ (pow t_m 1.5) l) (/ k_m t_m)) 2.0))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 0.0195) {
tmp = 2.0 / pow((sqrt((sin(k_m) * tan(k_m))) * (k_m * (sqrt(t_m) / l))), 2.0);
} else if (k_m <= 2.5e+165) {
tmp = (2.0 / pow(k_m, 2.0)) * (pow(l, 2.0) / (t_m * (pow(sin(k_m), 2.0) / cos(k_m))));
} else {
tmp = 2.0 / (sin(k_m) * (tan(k_m) * pow(((pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0)));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 0.0195d0) then
tmp = 2.0d0 / ((sqrt((sin(k_m) * tan(k_m))) * (k_m * (sqrt(t_m) / l))) ** 2.0d0)
else if (k_m <= 2.5d+165) then
tmp = (2.0d0 / (k_m ** 2.0d0)) * ((l ** 2.0d0) / (t_m * ((sin(k_m) ** 2.0d0) / cos(k_m))))
else
tmp = 2.0d0 / (sin(k_m) * (tan(k_m) * ((((t_m ** 1.5d0) / l) * (k_m / t_m)) ** 2.0d0)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 0.0195) {
tmp = 2.0 / Math.pow((Math.sqrt((Math.sin(k_m) * Math.tan(k_m))) * (k_m * (Math.sqrt(t_m) / l))), 2.0);
} else if (k_m <= 2.5e+165) {
tmp = (2.0 / Math.pow(k_m, 2.0)) * (Math.pow(l, 2.0) / (t_m * (Math.pow(Math.sin(k_m), 2.0) / Math.cos(k_m))));
} else {
tmp = 2.0 / (Math.sin(k_m) * (Math.tan(k_m) * Math.pow(((Math.pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 0.0195: tmp = 2.0 / math.pow((math.sqrt((math.sin(k_m) * math.tan(k_m))) * (k_m * (math.sqrt(t_m) / l))), 2.0) elif k_m <= 2.5e+165: tmp = (2.0 / math.pow(k_m, 2.0)) * (math.pow(l, 2.0) / (t_m * (math.pow(math.sin(k_m), 2.0) / math.cos(k_m)))) else: tmp = 2.0 / (math.sin(k_m) * (math.tan(k_m) * math.pow(((math.pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 0.0195) tmp = Float64(2.0 / (Float64(sqrt(Float64(sin(k_m) * tan(k_m))) * Float64(k_m * Float64(sqrt(t_m) / l))) ^ 2.0)); elseif (k_m <= 2.5e+165) tmp = Float64(Float64(2.0 / (k_m ^ 2.0)) * Float64((l ^ 2.0) / Float64(t_m * Float64((sin(k_m) ^ 2.0) / cos(k_m))))); else tmp = Float64(2.0 / Float64(sin(k_m) * Float64(tan(k_m) * (Float64(Float64((t_m ^ 1.5) / l) * Float64(k_m / t_m)) ^ 2.0)))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 0.0195) tmp = 2.0 / ((sqrt((sin(k_m) * tan(k_m))) * (k_m * (sqrt(t_m) / l))) ^ 2.0); elseif (k_m <= 2.5e+165) tmp = (2.0 / (k_m ^ 2.0)) * ((l ^ 2.0) / (t_m * ((sin(k_m) ^ 2.0) / cos(k_m)))); else tmp = 2.0 / (sin(k_m) * (tan(k_m) * ((((t_m ^ 1.5) / l) * (k_m / t_m)) ^ 2.0))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 0.0195], N[(2.0 / N[Power[N[(N[Sqrt[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(k$95$m * N[(N[Sqrt[t$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2.5e+165], N[(N[(2.0 / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Power[l, 2.0], $MachinePrecision] / N[(t$95$m * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Tan[k$95$m], $MachinePrecision] * N[Power[N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 0.0195:\\
\;\;\;\;\frac{2}{{\left(\sqrt{\sin k\_m \cdot \tan k\_m} \cdot \left(k\_m \cdot \frac{\sqrt{t\_m}}{\ell}\right)\right)}^{2}}\\
\mathbf{elif}\;k\_m \leq 2.5 \cdot 10^{+165}:\\
\;\;\;\;\frac{2}{{k\_m}^{2}} \cdot \frac{{\ell}^{2}}{t\_m \cdot \frac{{\sin k\_m}^{2}}{\cos k\_m}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\sin k\_m \cdot \left(\tan k\_m \cdot {\left(\frac{{t\_m}^{1.5}}{\ell} \cdot \frac{k\_m}{t\_m}\right)}^{2}\right)}\\
\end{array}
\end{array}
if k < 0.0195Initial program 40.7%
*-commutative40.7%
associate-/r*40.7%
Simplified47.5%
add-sqr-sqrt21.7%
pow221.7%
sqrt-prod19.1%
sqrt-div20.6%
sqrt-pow124.6%
metadata-eval24.6%
sqrt-prod15.0%
add-sqr-sqrt29.1%
Applied egg-rr29.1%
*-un-lft-identity29.1%
associate-/l/29.1%
+-rgt-identity29.1%
pow-prod-down33.7%
*-commutative33.7%
Applied egg-rr33.7%
*-lft-identity33.7%
associate-*l*33.7%
associate-*l/34.1%
Simplified34.1%
Taylor expanded in t around 0 43.2%
associate-*l/42.8%
associate-/l*43.2%
Simplified43.2%
if 0.0195 < k < 2.49999999999999985e165Initial program 19.1%
Simplified34.7%
Taylor expanded in t around 0 81.9%
associate-/l*81.8%
Simplified81.8%
associate-*l/81.8%
pow281.8%
associate-/l*81.9%
Applied egg-rr81.9%
pow281.9%
times-frac84.7%
pow284.7%
Applied egg-rr84.7%
if 2.49999999999999985e165 < k Initial program 42.4%
*-commutative42.4%
associate-/r*42.4%
Simplified46.5%
add-sqr-sqrt12.8%
pow212.8%
sqrt-prod12.8%
sqrt-div12.8%
sqrt-pow116.7%
metadata-eval16.7%
sqrt-prod8.3%
add-sqr-sqrt21.0%
Applied egg-rr21.0%
*-un-lft-identity21.0%
associate-/l/21.0%
+-rgt-identity21.0%
pow-prod-down21.0%
*-commutative21.0%
Applied egg-rr21.0%
*-lft-identity21.0%
associate-*l*21.0%
associate-*l/21.0%
Simplified21.0%
*-un-lft-identity21.0%
unpow-prod-down21.0%
pow221.0%
add-sqr-sqrt47.8%
associate-/l*47.9%
unpow-prod-down33.7%
pow233.7%
pow-prod-up57.8%
metadata-eval57.8%
Applied egg-rr57.8%
*-lft-identity57.8%
associate-*l*57.8%
unpow257.8%
metadata-eval57.8%
pow-sqr33.7%
swap-sqr47.8%
unpow247.8%
associate-*r/47.7%
associate-*l/47.8%
Simplified47.8%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 70.0)
(/
2.0
(pow (* (sqrt (* (sin k_m) (tan k_m))) (* k_m (/ (sqrt t_m) l))) 2.0))
(if (<= k_m 2.5e+165)
(*
(* 2.0 (/ (/ (/ (cos k_m) (pow k_m 2.0)) t_m) (pow (sin k_m) 2.0)))
(* l l))
(/
2.0
(*
(sin k_m)
(* (tan k_m) (pow (* (/ (pow t_m 1.5) l) (/ k_m t_m)) 2.0))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 70.0) {
tmp = 2.0 / pow((sqrt((sin(k_m) * tan(k_m))) * (k_m * (sqrt(t_m) / l))), 2.0);
} else if (k_m <= 2.5e+165) {
tmp = (2.0 * (((cos(k_m) / pow(k_m, 2.0)) / t_m) / pow(sin(k_m), 2.0))) * (l * l);
} else {
tmp = 2.0 / (sin(k_m) * (tan(k_m) * pow(((pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0)));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 70.0d0) then
tmp = 2.0d0 / ((sqrt((sin(k_m) * tan(k_m))) * (k_m * (sqrt(t_m) / l))) ** 2.0d0)
else if (k_m <= 2.5d+165) then
tmp = (2.0d0 * (((cos(k_m) / (k_m ** 2.0d0)) / t_m) / (sin(k_m) ** 2.0d0))) * (l * l)
else
tmp = 2.0d0 / (sin(k_m) * (tan(k_m) * ((((t_m ** 1.5d0) / l) * (k_m / t_m)) ** 2.0d0)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 70.0) {
tmp = 2.0 / Math.pow((Math.sqrt((Math.sin(k_m) * Math.tan(k_m))) * (k_m * (Math.sqrt(t_m) / l))), 2.0);
} else if (k_m <= 2.5e+165) {
tmp = (2.0 * (((Math.cos(k_m) / Math.pow(k_m, 2.0)) / t_m) / Math.pow(Math.sin(k_m), 2.0))) * (l * l);
} else {
tmp = 2.0 / (Math.sin(k_m) * (Math.tan(k_m) * Math.pow(((Math.pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if k_m <= 70.0: tmp = 2.0 / math.pow((math.sqrt((math.sin(k_m) * math.tan(k_m))) * (k_m * (math.sqrt(t_m) / l))), 2.0) elif k_m <= 2.5e+165: tmp = (2.0 * (((math.cos(k_m) / math.pow(k_m, 2.0)) / t_m) / math.pow(math.sin(k_m), 2.0))) * (l * l) else: tmp = 2.0 / (math.sin(k_m) * (math.tan(k_m) * math.pow(((math.pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 70.0) tmp = Float64(2.0 / (Float64(sqrt(Float64(sin(k_m) * tan(k_m))) * Float64(k_m * Float64(sqrt(t_m) / l))) ^ 2.0)); elseif (k_m <= 2.5e+165) tmp = Float64(Float64(2.0 * Float64(Float64(Float64(cos(k_m) / (k_m ^ 2.0)) / t_m) / (sin(k_m) ^ 2.0))) * Float64(l * l)); else tmp = Float64(2.0 / Float64(sin(k_m) * Float64(tan(k_m) * (Float64(Float64((t_m ^ 1.5) / l) * Float64(k_m / t_m)) ^ 2.0)))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (k_m <= 70.0) tmp = 2.0 / ((sqrt((sin(k_m) * tan(k_m))) * (k_m * (sqrt(t_m) / l))) ^ 2.0); elseif (k_m <= 2.5e+165) tmp = (2.0 * (((cos(k_m) / (k_m ^ 2.0)) / t_m) / (sin(k_m) ^ 2.0))) * (l * l); else tmp = 2.0 / (sin(k_m) * (tan(k_m) * ((((t_m ^ 1.5) / l) * (k_m / t_m)) ^ 2.0))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 70.0], N[(2.0 / N[Power[N[(N[Sqrt[N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(k$95$m * N[(N[Sqrt[t$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2.5e+165], N[(N[(2.0 * N[(N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l * l), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Tan[k$95$m], $MachinePrecision] * N[Power[N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 70:\\
\;\;\;\;\frac{2}{{\left(\sqrt{\sin k\_m \cdot \tan k\_m} \cdot \left(k\_m \cdot \frac{\sqrt{t\_m}}{\ell}\right)\right)}^{2}}\\
\mathbf{elif}\;k\_m \leq 2.5 \cdot 10^{+165}:\\
\;\;\;\;\left(2 \cdot \frac{\frac{\frac{\cos k\_m}{{k\_m}^{2}}}{t\_m}}{{\sin k\_m}^{2}}\right) \cdot \left(\ell \cdot \ell\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\sin k\_m \cdot \left(\tan k\_m \cdot {\left(\frac{{t\_m}^{1.5}}{\ell} \cdot \frac{k\_m}{t\_m}\right)}^{2}\right)}\\
\end{array}
\end{array}
if k < 70Initial program 40.7%
*-commutative40.7%
associate-/r*40.7%
Simplified47.5%
add-sqr-sqrt21.7%
pow221.7%
sqrt-prod19.1%
sqrt-div20.6%
sqrt-pow124.6%
metadata-eval24.6%
sqrt-prod15.0%
add-sqr-sqrt29.1%
Applied egg-rr29.1%
*-un-lft-identity29.1%
associate-/l/29.1%
+-rgt-identity29.1%
pow-prod-down33.7%
*-commutative33.7%
Applied egg-rr33.7%
*-lft-identity33.7%
associate-*l*33.7%
associate-*l/34.1%
Simplified34.1%
Taylor expanded in t around 0 43.2%
associate-*l/42.8%
associate-/l*43.2%
Simplified43.2%
if 70 < k < 2.49999999999999985e165Initial program 19.1%
Simplified34.7%
Taylor expanded in t around 0 81.9%
associate-/r*81.9%
associate-/r*81.9%
Simplified81.9%
if 2.49999999999999985e165 < k Initial program 42.4%
*-commutative42.4%
associate-/r*42.4%
Simplified46.5%
add-sqr-sqrt12.8%
pow212.8%
sqrt-prod12.8%
sqrt-div12.8%
sqrt-pow116.7%
metadata-eval16.7%
sqrt-prod8.3%
add-sqr-sqrt21.0%
Applied egg-rr21.0%
*-un-lft-identity21.0%
associate-/l/21.0%
+-rgt-identity21.0%
pow-prod-down21.0%
*-commutative21.0%
Applied egg-rr21.0%
*-lft-identity21.0%
associate-*l*21.0%
associate-*l/21.0%
Simplified21.0%
*-un-lft-identity21.0%
unpow-prod-down21.0%
pow221.0%
add-sqr-sqrt47.8%
associate-/l*47.9%
unpow-prod-down33.7%
pow233.7%
pow-prod-up57.8%
metadata-eval57.8%
Applied egg-rr57.8%
*-lft-identity57.8%
associate-*l*57.8%
unpow257.8%
metadata-eval57.8%
pow-sqr33.7%
swap-sqr47.8%
unpow247.8%
associate-*r/47.7%
associate-*l/47.8%
Simplified47.8%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (pow (sin k_m) 2.0)))
(*
t_s
(if (<= t_m 2.1e-209)
(* (* 2.0 (/ (/ (/ (cos k_m) (pow k_m 2.0)) t_m) t_2)) (* l l))
(if (<= t_m 7.1e+171)
(/
2.0
(*
(sin k_m)
(* (tan k_m) (pow (* (/ (pow t_m 1.5) l) (/ k_m t_m)) 2.0))))
(* (* l l) (/ 2.0 (* (pow k_m 2.0) (/ (* t_m t_2) (cos k_m))))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = pow(sin(k_m), 2.0);
double tmp;
if (t_m <= 2.1e-209) {
tmp = (2.0 * (((cos(k_m) / pow(k_m, 2.0)) / t_m) / t_2)) * (l * l);
} else if (t_m <= 7.1e+171) {
tmp = 2.0 / (sin(k_m) * (tan(k_m) * pow(((pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0)));
} else {
tmp = (l * l) * (2.0 / (pow(k_m, 2.0) * ((t_m * t_2) / cos(k_m))));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: t_2
real(8) :: tmp
t_2 = sin(k_m) ** 2.0d0
if (t_m <= 2.1d-209) then
tmp = (2.0d0 * (((cos(k_m) / (k_m ** 2.0d0)) / t_m) / t_2)) * (l * l)
else if (t_m <= 7.1d+171) then
tmp = 2.0d0 / (sin(k_m) * (tan(k_m) * ((((t_m ** 1.5d0) / l) * (k_m / t_m)) ** 2.0d0)))
else
tmp = (l * l) * (2.0d0 / ((k_m ** 2.0d0) * ((t_m * t_2) / cos(k_m))))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = Math.pow(Math.sin(k_m), 2.0);
double tmp;
if (t_m <= 2.1e-209) {
tmp = (2.0 * (((Math.cos(k_m) / Math.pow(k_m, 2.0)) / t_m) / t_2)) * (l * l);
} else if (t_m <= 7.1e+171) {
tmp = 2.0 / (Math.sin(k_m) * (Math.tan(k_m) * Math.pow(((Math.pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0)));
} else {
tmp = (l * l) * (2.0 / (Math.pow(k_m, 2.0) * ((t_m * t_2) / Math.cos(k_m))));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): t_2 = math.pow(math.sin(k_m), 2.0) tmp = 0 if t_m <= 2.1e-209: tmp = (2.0 * (((math.cos(k_m) / math.pow(k_m, 2.0)) / t_m) / t_2)) * (l * l) elif t_m <= 7.1e+171: tmp = 2.0 / (math.sin(k_m) * (math.tan(k_m) * math.pow(((math.pow(t_m, 1.5) / l) * (k_m / t_m)), 2.0))) else: tmp = (l * l) * (2.0 / (math.pow(k_m, 2.0) * ((t_m * t_2) / math.cos(k_m)))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = sin(k_m) ^ 2.0 tmp = 0.0 if (t_m <= 2.1e-209) tmp = Float64(Float64(2.0 * Float64(Float64(Float64(cos(k_m) / (k_m ^ 2.0)) / t_m) / t_2)) * Float64(l * l)); elseif (t_m <= 7.1e+171) tmp = Float64(2.0 / Float64(sin(k_m) * Float64(tan(k_m) * (Float64(Float64((t_m ^ 1.5) / l) * Float64(k_m / t_m)) ^ 2.0)))); else tmp = Float64(Float64(l * l) * Float64(2.0 / Float64((k_m ^ 2.0) * Float64(Float64(t_m * t_2) / cos(k_m))))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) t_2 = sin(k_m) ^ 2.0; tmp = 0.0; if (t_m <= 2.1e-209) tmp = (2.0 * (((cos(k_m) / (k_m ^ 2.0)) / t_m) / t_2)) * (l * l); elseif (t_m <= 7.1e+171) tmp = 2.0 / (sin(k_m) * (tan(k_m) * ((((t_m ^ 1.5) / l) * (k_m / t_m)) ^ 2.0))); else tmp = (l * l) * (2.0 / ((k_m ^ 2.0) * ((t_m * t_2) / cos(k_m)))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]}, N[(t$95$s * If[LessEqual[t$95$m, 2.1e-209], N[(N[(2.0 * N[(N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision] * N[(l * l), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 7.1e+171], N[(2.0 / N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Tan[k$95$m], $MachinePrecision] * N[Power[N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(N[(t$95$m * t$95$2), $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := {\sin k\_m}^{2}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 2.1 \cdot 10^{-209}:\\
\;\;\;\;\left(2 \cdot \frac{\frac{\frac{\cos k\_m}{{k\_m}^{2}}}{t\_m}}{t\_2}\right) \cdot \left(\ell \cdot \ell\right)\\
\mathbf{elif}\;t\_m \leq 7.1 \cdot 10^{+171}:\\
\;\;\;\;\frac{2}{\sin k\_m \cdot \left(\tan k\_m \cdot {\left(\frac{{t\_m}^{1.5}}{\ell} \cdot \frac{k\_m}{t\_m}\right)}^{2}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{{k\_m}^{2} \cdot \frac{t\_m \cdot t\_2}{\cos k\_m}}\\
\end{array}
\end{array}
\end{array}
if t < 2.09999999999999996e-209Initial program 35.4%
Simplified43.1%
Taylor expanded in t around 0 75.5%
associate-/r*75.5%
associate-/r*75.5%
Simplified75.5%
if 2.09999999999999996e-209 < t < 7.09999999999999969e171Initial program 53.7%
*-commutative53.7%
associate-/r*53.7%
Simplified58.8%
add-sqr-sqrt43.5%
pow243.5%
sqrt-prod43.5%
sqrt-div47.2%
sqrt-pow158.4%
metadata-eval58.4%
sqrt-prod35.4%
add-sqr-sqrt69.7%
Applied egg-rr69.7%
*-un-lft-identity69.7%
associate-/l/69.7%
+-rgt-identity69.7%
pow-prod-down74.8%
*-commutative74.8%
Applied egg-rr74.8%
*-lft-identity74.8%
associate-*l*74.7%
associate-*l/75.8%
Simplified75.8%
*-un-lft-identity75.8%
unpow-prod-down74.6%
pow274.6%
add-sqr-sqrt94.8%
associate-/l*94.8%
unpow-prod-down71.2%
pow271.2%
pow-prod-up71.2%
metadata-eval71.2%
Applied egg-rr71.2%
*-lft-identity71.2%
associate-*l*72.4%
unpow272.4%
metadata-eval72.4%
pow-sqr72.4%
swap-sqr96.0%
unpow296.0%
associate-*r/96.0%
associate-*l/96.0%
Simplified96.0%
if 7.09999999999999969e171 < t Initial program 4.2%
Simplified29.2%
Taylor expanded in t around 0 75.5%
associate-/l*75.5%
Simplified75.5%
Final simplification81.8%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 0.0)
(/ 2.0 (pow (* k_m (/ (* (pow t_m 1.5) (/ k_m t_m)) l)) 2.0))
(*
(* 2.0 (/ (/ (/ (cos k_m) (pow k_m 2.0)) t_m) (pow (sin k_m) 2.0)))
(* l l)))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / pow((k_m * ((pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else {
tmp = (2.0 * (((cos(k_m) / pow(k_m, 2.0)) / t_m) / pow(sin(k_m), 2.0))) * (l * l);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 0.0d0) then
tmp = 2.0d0 / ((k_m * (((t_m ** 1.5d0) * (k_m / t_m)) / l)) ** 2.0d0)
else
tmp = (2.0d0 * (((cos(k_m) / (k_m ** 2.0d0)) / t_m) / (sin(k_m) ** 2.0d0))) * (l * l)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / Math.pow((k_m * ((Math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else {
tmp = (2.0 * (((Math.cos(k_m) / Math.pow(k_m, 2.0)) / t_m) / Math.pow(Math.sin(k_m), 2.0))) * (l * l);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 0.0: tmp = 2.0 / math.pow((k_m * ((math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0) else: tmp = (2.0 * (((math.cos(k_m) / math.pow(k_m, 2.0)) / t_m) / math.pow(math.sin(k_m), 2.0))) * (l * l) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 0.0) tmp = Float64(2.0 / (Float64(k_m * Float64(Float64((t_m ^ 1.5) * Float64(k_m / t_m)) / l)) ^ 2.0)); else tmp = Float64(Float64(2.0 * Float64(Float64(Float64(cos(k_m) / (k_m ^ 2.0)) / t_m) / (sin(k_m) ^ 2.0))) * Float64(l * l)); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 0.0) tmp = 2.0 / ((k_m * (((t_m ^ 1.5) * (k_m / t_m)) / l)) ^ 2.0); else tmp = (2.0 * (((cos(k_m) / (k_m ^ 2.0)) / t_m) / (sin(k_m) ^ 2.0))) * (l * l); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 0.0], N[(2.0 / N[Power[N[(k$95$m * N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(N[(N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l * l), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 0:\\
\;\;\;\;\frac{2}{{\left(k\_m \cdot \frac{{t\_m}^{1.5} \cdot \frac{k\_m}{t\_m}}{\ell}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \frac{\frac{\frac{\cos k\_m}{{k\_m}^{2}}}{t\_m}}{{\sin k\_m}^{2}}\right) \cdot \left(\ell \cdot \ell\right)\\
\end{array}
\end{array}
if (*.f64 l l) < 0.0Initial program 18.6%
*-commutative18.6%
associate-/r*18.6%
Simplified36.3%
add-sqr-sqrt13.0%
pow213.0%
sqrt-prod9.3%
sqrt-div9.3%
sqrt-pow111.2%
metadata-eval11.2%
sqrt-prod16.6%
add-sqr-sqrt27.8%
Applied egg-rr27.8%
*-un-lft-identity27.8%
associate-/l/27.9%
+-rgt-identity27.9%
pow-prod-down28.1%
*-commutative28.1%
Applied egg-rr28.1%
*-lft-identity28.1%
associate-*l*28.1%
associate-*l/29.8%
Simplified29.8%
Taylor expanded in k around 0 42.6%
if 0.0 < (*.f64 l l) Initial program 43.4%
Simplified50.3%
Taylor expanded in t around 0 81.9%
associate-/r*81.9%
associate-/r*81.9%
Simplified81.9%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 2e-92)
(/ 2.0 (pow (* k_m (/ (* (pow t_m 1.5) (/ k_m t_m)) l)) 2.0))
(*
(* l l)
(/
2.0
(*
(pow k_m 2.0)
(/ (* t_m (- 0.5 (/ (cos (* k_m 2.0)) 2.0))) (cos k_m))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 2e-92) {
tmp = 2.0 / pow((k_m * ((pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else {
tmp = (l * l) * (2.0 / (pow(k_m, 2.0) * ((t_m * (0.5 - (cos((k_m * 2.0)) / 2.0))) / cos(k_m))));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 2d-92) then
tmp = 2.0d0 / ((k_m * (((t_m ** 1.5d0) * (k_m / t_m)) / l)) ** 2.0d0)
else
tmp = (l * l) * (2.0d0 / ((k_m ** 2.0d0) * ((t_m * (0.5d0 - (cos((k_m * 2.0d0)) / 2.0d0))) / cos(k_m))))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 2e-92) {
tmp = 2.0 / Math.pow((k_m * ((Math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else {
tmp = (l * l) * (2.0 / (Math.pow(k_m, 2.0) * ((t_m * (0.5 - (Math.cos((k_m * 2.0)) / 2.0))) / Math.cos(k_m))));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 2e-92: tmp = 2.0 / math.pow((k_m * ((math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0) else: tmp = (l * l) * (2.0 / (math.pow(k_m, 2.0) * ((t_m * (0.5 - (math.cos((k_m * 2.0)) / 2.0))) / math.cos(k_m)))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 2e-92) tmp = Float64(2.0 / (Float64(k_m * Float64(Float64((t_m ^ 1.5) * Float64(k_m / t_m)) / l)) ^ 2.0)); else tmp = Float64(Float64(l * l) * Float64(2.0 / Float64((k_m ^ 2.0) * Float64(Float64(t_m * Float64(0.5 - Float64(cos(Float64(k_m * 2.0)) / 2.0))) / cos(k_m))))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 2e-92) tmp = 2.0 / ((k_m * (((t_m ^ 1.5) * (k_m / t_m)) / l)) ^ 2.0); else tmp = (l * l) * (2.0 / ((k_m ^ 2.0) * ((t_m * (0.5 - (cos((k_m * 2.0)) / 2.0))) / cos(k_m)))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 2e-92], N[(2.0 / N[Power[N[(k$95$m * N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(N[(t$95$m * N[(0.5 - N[(N[Cos[N[(k$95$m * 2.0), $MachinePrecision]], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 2 \cdot 10^{-92}:\\
\;\;\;\;\frac{2}{{\left(k\_m \cdot \frac{{t\_m}^{1.5} \cdot \frac{k\_m}{t\_m}}{\ell}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{{k\_m}^{2} \cdot \frac{t\_m \cdot \left(0.5 - \frac{\cos \left(k\_m \cdot 2\right)}{2}\right)}{\cos k\_m}}\\
\end{array}
\end{array}
if (*.f64 l l) < 1.99999999999999998e-92Initial program 33.5%
*-commutative33.5%
associate-/r*33.5%
Simplified45.8%
add-sqr-sqrt18.7%
pow218.7%
sqrt-prod13.8%
sqrt-div15.7%
sqrt-pow117.7%
metadata-eval17.7%
sqrt-prod16.7%
add-sqr-sqrt26.5%
Applied egg-rr26.5%
*-un-lft-identity26.5%
associate-/l/26.5%
+-rgt-identity26.5%
pow-prod-down29.7%
*-commutative29.7%
Applied egg-rr29.7%
*-lft-identity29.7%
associate-*l*29.7%
associate-*l/30.6%
Simplified30.6%
Taylor expanded in k around 0 40.2%
if 1.99999999999999998e-92 < (*.f64 l l) Initial program 41.2%
Simplified47.1%
Taylor expanded in t around 0 79.6%
associate-/l*79.6%
Simplified79.6%
unpow279.6%
sin-mult76.1%
Applied egg-rr76.1%
div-sub76.1%
+-inverses76.1%
cos-076.1%
metadata-eval76.1%
count-276.1%
*-commutative76.1%
Simplified76.1%
Final simplification61.8%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= t_m 4.4e-220)
(* -2.0 (/ (/ (pow l 2.0) (pow k_m 4.0)) (- t_m)))
(if (<= t_m 1.25e+172)
(/ 2.0 (pow (* k_m (/ (* (pow t_m 1.5) (/ k_m t_m)) l)) 2.0))
(/ (* 2.0 (pow l 2.0)) (* (pow k_m 2.0) (* t_m (pow k_m 2.0))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (t_m <= 4.4e-220) {
tmp = -2.0 * ((pow(l, 2.0) / pow(k_m, 4.0)) / -t_m);
} else if (t_m <= 1.25e+172) {
tmp = 2.0 / pow((k_m * ((pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else {
tmp = (2.0 * pow(l, 2.0)) / (pow(k_m, 2.0) * (t_m * pow(k_m, 2.0)));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (t_m <= 4.4d-220) then
tmp = (-2.0d0) * (((l ** 2.0d0) / (k_m ** 4.0d0)) / -t_m)
else if (t_m <= 1.25d+172) then
tmp = 2.0d0 / ((k_m * (((t_m ** 1.5d0) * (k_m / t_m)) / l)) ** 2.0d0)
else
tmp = (2.0d0 * (l ** 2.0d0)) / ((k_m ** 2.0d0) * (t_m * (k_m ** 2.0d0)))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (t_m <= 4.4e-220) {
tmp = -2.0 * ((Math.pow(l, 2.0) / Math.pow(k_m, 4.0)) / -t_m);
} else if (t_m <= 1.25e+172) {
tmp = 2.0 / Math.pow((k_m * ((Math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else {
tmp = (2.0 * Math.pow(l, 2.0)) / (Math.pow(k_m, 2.0) * (t_m * Math.pow(k_m, 2.0)));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if t_m <= 4.4e-220: tmp = -2.0 * ((math.pow(l, 2.0) / math.pow(k_m, 4.0)) / -t_m) elif t_m <= 1.25e+172: tmp = 2.0 / math.pow((k_m * ((math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0) else: tmp = (2.0 * math.pow(l, 2.0)) / (math.pow(k_m, 2.0) * (t_m * math.pow(k_m, 2.0))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (t_m <= 4.4e-220) tmp = Float64(-2.0 * Float64(Float64((l ^ 2.0) / (k_m ^ 4.0)) / Float64(-t_m))); elseif (t_m <= 1.25e+172) tmp = Float64(2.0 / (Float64(k_m * Float64(Float64((t_m ^ 1.5) * Float64(k_m / t_m)) / l)) ^ 2.0)); else tmp = Float64(Float64(2.0 * (l ^ 2.0)) / Float64((k_m ^ 2.0) * Float64(t_m * (k_m ^ 2.0)))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (t_m <= 4.4e-220) tmp = -2.0 * (((l ^ 2.0) / (k_m ^ 4.0)) / -t_m); elseif (t_m <= 1.25e+172) tmp = 2.0 / ((k_m * (((t_m ^ 1.5) * (k_m / t_m)) / l)) ^ 2.0); else tmp = (2.0 * (l ^ 2.0)) / ((k_m ^ 2.0) * (t_m * (k_m ^ 2.0))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[t$95$m, 4.4e-220], N[(-2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision] / (-t$95$m)), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$m, 1.25e+172], N[(2.0 / N[Power[N[(k$95$m * N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;t\_m \leq 4.4 \cdot 10^{-220}:\\
\;\;\;\;-2 \cdot \frac{\frac{{\ell}^{2}}{{k\_m}^{4}}}{-t\_m}\\
\mathbf{elif}\;t\_m \leq 1.25 \cdot 10^{+172}:\\
\;\;\;\;\frac{2}{{\left(k\_m \cdot \frac{{t\_m}^{1.5} \cdot \frac{k\_m}{t\_m}}{\ell}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot {\ell}^{2}}{{k\_m}^{2} \cdot \left(t\_m \cdot {k\_m}^{2}\right)}\\
\end{array}
\end{array}
if t < 4.39999999999999973e-220Initial program 34.8%
*-commutative34.8%
associate-/r*34.8%
Simplified41.2%
add-sqr-sqrt8.1%
pow28.1%
sqrt-prod3.3%
sqrt-div3.3%
sqrt-pow13.3%
metadata-eval3.3%
sqrt-prod1.3%
add-sqr-sqrt4.0%
Applied egg-rr4.0%
Taylor expanded in k around 0 4.0%
Taylor expanded in t around -inf 0.0%
associate-/r*0.0%
*-commutative0.0%
unpow20.0%
rem-square-sqrt64.3%
mul-1-neg64.3%
Simplified64.3%
if 4.39999999999999973e-220 < t < 1.25e172Initial program 54.1%
*-commutative54.1%
associate-/r*54.2%
Simplified59.0%
add-sqr-sqrt44.3%
pow244.3%
sqrt-prod44.3%
sqrt-div47.9%
sqrt-pow158.7%
metadata-eval58.7%
sqrt-prod36.5%
add-sqr-sqrt69.6%
Applied egg-rr69.6%
*-un-lft-identity69.6%
associate-/l/69.6%
+-rgt-identity69.6%
pow-prod-down74.5%
*-commutative74.5%
Applied egg-rr74.5%
*-lft-identity74.5%
associate-*l*74.4%
associate-*l/75.5%
Simplified75.5%
Taylor expanded in k around 0 77.9%
if 1.25e172 < t Initial program 4.2%
Simplified29.2%
Taylor expanded in t around 0 75.5%
associate-/l*75.5%
Simplified75.5%
associate-*l/75.5%
pow275.5%
associate-/l*75.5%
Applied egg-rr75.5%
Taylor expanded in k around 0 71.3%
Final simplification69.3%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 0.0)
(/ 2.0 (pow (* k_m (/ (* (pow t_m 1.5) (/ k_m t_m)) l)) 2.0))
(*
(* l l)
(/ 2.0 (* (pow k_m 2.0) (/ (* t_m (pow k_m 2.0)) (cos k_m))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / pow((k_m * ((pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else {
tmp = (l * l) * (2.0 / (pow(k_m, 2.0) * ((t_m * pow(k_m, 2.0)) / cos(k_m))));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 0.0d0) then
tmp = 2.0d0 / ((k_m * (((t_m ** 1.5d0) * (k_m / t_m)) / l)) ** 2.0d0)
else
tmp = (l * l) * (2.0d0 / ((k_m ** 2.0d0) * ((t_m * (k_m ** 2.0d0)) / cos(k_m))))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / Math.pow((k_m * ((Math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else {
tmp = (l * l) * (2.0 / (Math.pow(k_m, 2.0) * ((t_m * Math.pow(k_m, 2.0)) / Math.cos(k_m))));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 0.0: tmp = 2.0 / math.pow((k_m * ((math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0) else: tmp = (l * l) * (2.0 / (math.pow(k_m, 2.0) * ((t_m * math.pow(k_m, 2.0)) / math.cos(k_m)))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 0.0) tmp = Float64(2.0 / (Float64(k_m * Float64(Float64((t_m ^ 1.5) * Float64(k_m / t_m)) / l)) ^ 2.0)); else tmp = Float64(Float64(l * l) * Float64(2.0 / Float64((k_m ^ 2.0) * Float64(Float64(t_m * (k_m ^ 2.0)) / cos(k_m))))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 0.0) tmp = 2.0 / ((k_m * (((t_m ^ 1.5) * (k_m / t_m)) / l)) ^ 2.0); else tmp = (l * l) * (2.0 / ((k_m ^ 2.0) * ((t_m * (k_m ^ 2.0)) / cos(k_m)))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 0.0], N[(2.0 / N[Power[N[(k$95$m * N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(N[Power[k$95$m, 2.0], $MachinePrecision] * N[(N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 0:\\
\;\;\;\;\frac{2}{{\left(k\_m \cdot \frac{{t\_m}^{1.5} \cdot \frac{k\_m}{t\_m}}{\ell}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{2}{{k\_m}^{2} \cdot \frac{t\_m \cdot {k\_m}^{2}}{\cos k\_m}}\\
\end{array}
\end{array}
if (*.f64 l l) < 0.0Initial program 18.6%
*-commutative18.6%
associate-/r*18.6%
Simplified36.3%
add-sqr-sqrt13.0%
pow213.0%
sqrt-prod9.3%
sqrt-div9.3%
sqrt-pow111.2%
metadata-eval11.2%
sqrt-prod16.6%
add-sqr-sqrt27.8%
Applied egg-rr27.8%
*-un-lft-identity27.8%
associate-/l/27.9%
+-rgt-identity27.9%
pow-prod-down28.1%
*-commutative28.1%
Applied egg-rr28.1%
*-lft-identity28.1%
associate-*l*28.1%
associate-*l/29.8%
Simplified29.8%
Taylor expanded in k around 0 42.6%
if 0.0 < (*.f64 l l) Initial program 43.4%
Simplified50.3%
Taylor expanded in t around 0 81.9%
associate-/l*81.8%
Simplified81.8%
Taylor expanded in k around 0 73.2%
Final simplification66.8%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 0.0)
(/ 2.0 (pow (* k_m (/ (* (pow t_m 1.5) (/ k_m t_m)) l)) 2.0))
(if (<= (* l l) 2e+304)
(* -2.0 (/ (/ (pow l 2.0) (pow k_m 4.0)) (- t_m)))
(*
(* l l)
(/
(+ (* -0.3333333333333333 (/ (pow k_m 2.0) t_m)) (* 2.0 (/ 1.0 t_m)))
(pow k_m 4.0)))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / pow((k_m * ((pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else if ((l * l) <= 2e+304) {
tmp = -2.0 * ((pow(l, 2.0) / pow(k_m, 4.0)) / -t_m);
} else {
tmp = (l * l) * (((-0.3333333333333333 * (pow(k_m, 2.0) / t_m)) + (2.0 * (1.0 / t_m))) / pow(k_m, 4.0));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 0.0d0) then
tmp = 2.0d0 / ((k_m * (((t_m ** 1.5d0) * (k_m / t_m)) / l)) ** 2.0d0)
else if ((l * l) <= 2d+304) then
tmp = (-2.0d0) * (((l ** 2.0d0) / (k_m ** 4.0d0)) / -t_m)
else
tmp = (l * l) * ((((-0.3333333333333333d0) * ((k_m ** 2.0d0) / t_m)) + (2.0d0 * (1.0d0 / t_m))) / (k_m ** 4.0d0))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / Math.pow((k_m * ((Math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else if ((l * l) <= 2e+304) {
tmp = -2.0 * ((Math.pow(l, 2.0) / Math.pow(k_m, 4.0)) / -t_m);
} else {
tmp = (l * l) * (((-0.3333333333333333 * (Math.pow(k_m, 2.0) / t_m)) + (2.0 * (1.0 / t_m))) / Math.pow(k_m, 4.0));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 0.0: tmp = 2.0 / math.pow((k_m * ((math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0) elif (l * l) <= 2e+304: tmp = -2.0 * ((math.pow(l, 2.0) / math.pow(k_m, 4.0)) / -t_m) else: tmp = (l * l) * (((-0.3333333333333333 * (math.pow(k_m, 2.0) / t_m)) + (2.0 * (1.0 / t_m))) / math.pow(k_m, 4.0)) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 0.0) tmp = Float64(2.0 / (Float64(k_m * Float64(Float64((t_m ^ 1.5) * Float64(k_m / t_m)) / l)) ^ 2.0)); elseif (Float64(l * l) <= 2e+304) tmp = Float64(-2.0 * Float64(Float64((l ^ 2.0) / (k_m ^ 4.0)) / Float64(-t_m))); else tmp = Float64(Float64(l * l) * Float64(Float64(Float64(-0.3333333333333333 * Float64((k_m ^ 2.0) / t_m)) + Float64(2.0 * Float64(1.0 / t_m))) / (k_m ^ 4.0))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 0.0) tmp = 2.0 / ((k_m * (((t_m ^ 1.5) * (k_m / t_m)) / l)) ^ 2.0); elseif ((l * l) <= 2e+304) tmp = -2.0 * (((l ^ 2.0) / (k_m ^ 4.0)) / -t_m); else tmp = (l * l) * (((-0.3333333333333333 * ((k_m ^ 2.0) / t_m)) + (2.0 * (1.0 / t_m))) / (k_m ^ 4.0)); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 0.0], N[(2.0 / N[Power[N[(k$95$m * N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(l * l), $MachinePrecision], 2e+304], N[(-2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision] / (-t$95$m)), $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] * N[(N[(N[(-0.3333333333333333 * N[(N[Power[k$95$m, 2.0], $MachinePrecision] / t$95$m), $MachinePrecision]), $MachinePrecision] + N[(2.0 * N[(1.0 / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 0:\\
\;\;\;\;\frac{2}{{\left(k\_m \cdot \frac{{t\_m}^{1.5} \cdot \frac{k\_m}{t\_m}}{\ell}\right)}^{2}}\\
\mathbf{elif}\;\ell \cdot \ell \leq 2 \cdot 10^{+304}:\\
\;\;\;\;-2 \cdot \frac{\frac{{\ell}^{2}}{{k\_m}^{4}}}{-t\_m}\\
\mathbf{else}:\\
\;\;\;\;\left(\ell \cdot \ell\right) \cdot \frac{-0.3333333333333333 \cdot \frac{{k\_m}^{2}}{t\_m} + 2 \cdot \frac{1}{t\_m}}{{k\_m}^{4}}\\
\end{array}
\end{array}
if (*.f64 l l) < 0.0Initial program 18.6%
*-commutative18.6%
associate-/r*18.6%
Simplified36.3%
add-sqr-sqrt13.0%
pow213.0%
sqrt-prod9.3%
sqrt-div9.3%
sqrt-pow111.2%
metadata-eval11.2%
sqrt-prod16.6%
add-sqr-sqrt27.8%
Applied egg-rr27.8%
*-un-lft-identity27.8%
associate-/l/27.9%
+-rgt-identity27.9%
pow-prod-down28.1%
*-commutative28.1%
Applied egg-rr28.1%
*-lft-identity28.1%
associate-*l*28.1%
associate-*l/29.8%
Simplified29.8%
Taylor expanded in k around 0 42.6%
if 0.0 < (*.f64 l l) < 1.9999999999999999e304Initial program 46.8%
*-commutative46.8%
associate-/r*46.8%
Simplified54.1%
add-sqr-sqrt24.2%
pow224.2%
sqrt-prod21.2%
sqrt-div23.4%
sqrt-pow127.6%
metadata-eval27.6%
sqrt-prod13.8%
add-sqr-sqrt27.6%
Applied egg-rr27.6%
Taylor expanded in k around 0 31.2%
Taylor expanded in t around -inf 0.0%
associate-/r*0.0%
*-commutative0.0%
unpow20.0%
rem-square-sqrt74.4%
mul-1-neg74.4%
Simplified74.4%
if 1.9999999999999999e304 < (*.f64 l l) Initial program 36.0%
Simplified36.0%
Taylor expanded in k around 0 61.1%
Final simplification64.4%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= (* l l) 0.0)
(/ 2.0 (pow (* k_m (/ (* (pow t_m 1.5) (/ k_m t_m)) l)) 2.0))
(* -2.0 (/ (/ (pow l 2.0) (pow k_m 4.0)) (- t_m))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / pow((k_m * ((pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else {
tmp = -2.0 * ((pow(l, 2.0) / pow(k_m, 4.0)) / -t_m);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l * l) <= 0.0d0) then
tmp = 2.0d0 / ((k_m * (((t_m ** 1.5d0) * (k_m / t_m)) / l)) ** 2.0d0)
else
tmp = (-2.0d0) * (((l ** 2.0d0) / (k_m ** 4.0d0)) / -t_m)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l * l) <= 0.0) {
tmp = 2.0 / Math.pow((k_m * ((Math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0);
} else {
tmp = -2.0 * ((Math.pow(l, 2.0) / Math.pow(k_m, 4.0)) / -t_m);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l * l) <= 0.0: tmp = 2.0 / math.pow((k_m * ((math.pow(t_m, 1.5) * (k_m / t_m)) / l)), 2.0) else: tmp = -2.0 * ((math.pow(l, 2.0) / math.pow(k_m, 4.0)) / -t_m) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (Float64(l * l) <= 0.0) tmp = Float64(2.0 / (Float64(k_m * Float64(Float64((t_m ^ 1.5) * Float64(k_m / t_m)) / l)) ^ 2.0)); else tmp = Float64(-2.0 * Float64(Float64((l ^ 2.0) / (k_m ^ 4.0)) / Float64(-t_m))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l * l) <= 0.0) tmp = 2.0 / ((k_m * (((t_m ^ 1.5) * (k_m / t_m)) / l)) ^ 2.0); else tmp = -2.0 * (((l ^ 2.0) / (k_m ^ 4.0)) / -t_m); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[N[(l * l), $MachinePrecision], 0.0], N[(2.0 / N[Power[N[(k$95$m * N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(-2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision] / (-t$95$m)), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \cdot \ell \leq 0:\\
\;\;\;\;\frac{2}{{\left(k\_m \cdot \frac{{t\_m}^{1.5} \cdot \frac{k\_m}{t\_m}}{\ell}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;-2 \cdot \frac{\frac{{\ell}^{2}}{{k\_m}^{4}}}{-t\_m}\\
\end{array}
\end{array}
if (*.f64 l l) < 0.0Initial program 18.6%
*-commutative18.6%
associate-/r*18.6%
Simplified36.3%
add-sqr-sqrt13.0%
pow213.0%
sqrt-prod9.3%
sqrt-div9.3%
sqrt-pow111.2%
metadata-eval11.2%
sqrt-prod16.6%
add-sqr-sqrt27.8%
Applied egg-rr27.8%
*-un-lft-identity27.8%
associate-/l/27.9%
+-rgt-identity27.9%
pow-prod-down28.1%
*-commutative28.1%
Applied egg-rr28.1%
*-lft-identity28.1%
associate-*l*28.1%
associate-*l/29.8%
Simplified29.8%
Taylor expanded in k around 0 42.6%
if 0.0 < (*.f64 l l) Initial program 43.4%
*-commutative43.4%
associate-/r*43.4%
Simplified48.4%
add-sqr-sqrt22.0%
pow222.0%
sqrt-prod19.5%
sqrt-div20.9%
sqrt-pow125.8%
metadata-eval25.8%
sqrt-prod12.4%
add-sqr-sqrt26.3%
Applied egg-rr26.3%
Taylor expanded in k around 0 28.3%
Taylor expanded in t around -inf 0.0%
associate-/r*0.0%
*-commutative0.0%
unpow20.0%
rem-square-sqrt68.2%
mul-1-neg68.2%
Simplified68.2%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= l 1.95e-162)
(/ 2.0 (pow (* k_m (* (/ (pow t_m 1.5) l) (/ k_m t_m))) 2.0))
(* -2.0 (/ (/ (pow l 2.0) (pow k_m 4.0)) (- t_m))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (l <= 1.95e-162) {
tmp = 2.0 / pow((k_m * ((pow(t_m, 1.5) / l) * (k_m / t_m))), 2.0);
} else {
tmp = -2.0 * ((pow(l, 2.0) / pow(k_m, 4.0)) / -t_m);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (l <= 1.95d-162) then
tmp = 2.0d0 / ((k_m * (((t_m ** 1.5d0) / l) * (k_m / t_m))) ** 2.0d0)
else
tmp = (-2.0d0) * (((l ** 2.0d0) / (k_m ** 4.0d0)) / -t_m)
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (l <= 1.95e-162) {
tmp = 2.0 / Math.pow((k_m * ((Math.pow(t_m, 1.5) / l) * (k_m / t_m))), 2.0);
} else {
tmp = -2.0 * ((Math.pow(l, 2.0) / Math.pow(k_m, 4.0)) / -t_m);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if l <= 1.95e-162: tmp = 2.0 / math.pow((k_m * ((math.pow(t_m, 1.5) / l) * (k_m / t_m))), 2.0) else: tmp = -2.0 * ((math.pow(l, 2.0) / math.pow(k_m, 4.0)) / -t_m) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (l <= 1.95e-162) tmp = Float64(2.0 / (Float64(k_m * Float64(Float64((t_m ^ 1.5) / l) * Float64(k_m / t_m))) ^ 2.0)); else tmp = Float64(-2.0 * Float64(Float64((l ^ 2.0) / (k_m ^ 4.0)) / Float64(-t_m))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if (l <= 1.95e-162) tmp = 2.0 / ((k_m * (((t_m ^ 1.5) / l) * (k_m / t_m))) ^ 2.0); else tmp = -2.0 * (((l ^ 2.0) / (k_m ^ 4.0)) / -t_m); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[l, 1.95e-162], N[(2.0 / N[Power[N[(k$95$m * N[(N[(N[Power[t$95$m, 1.5], $MachinePrecision] / l), $MachinePrecision] * N[(k$95$m / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], N[(-2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision] / (-t$95$m)), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \leq 1.95 \cdot 10^{-162}:\\
\;\;\;\;\frac{2}{{\left(k\_m \cdot \left(\frac{{t\_m}^{1.5}}{\ell} \cdot \frac{k\_m}{t\_m}\right)\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;-2 \cdot \frac{\frac{{\ell}^{2}}{{k\_m}^{4}}}{-t\_m}\\
\end{array}
\end{array}
if l < 1.95e-162Initial program 33.8%
*-commutative33.8%
associate-/r*33.8%
Simplified43.8%
add-sqr-sqrt18.2%
pow218.2%
sqrt-prod15.6%
sqrt-div16.9%
sqrt-pow121.4%
metadata-eval21.4%
sqrt-prod5.8%
add-sqr-sqrt27.8%
Applied egg-rr27.8%
Taylor expanded in k around 0 31.6%
*-un-lft-identity31.6%
associate-/l/31.6%
+-rgt-identity31.6%
pow-prod-down34.3%
*-commutative34.3%
Applied egg-rr34.3%
*-lft-identity34.3%
associate-*l*34.3%
Simplified34.3%
if 1.95e-162 < l Initial program 44.8%
*-commutative44.8%
associate-/r*44.8%
Simplified48.9%
add-sqr-sqrt23.1%
pow223.1%
sqrt-prod19.9%
sqrt-div20.9%
sqrt-pow124.8%
metadata-eval24.8%
sqrt-prod24.8%
add-sqr-sqrt24.8%
Applied egg-rr24.8%
Taylor expanded in k around 0 28.0%
Taylor expanded in t around -inf 0.0%
associate-/r*0.0%
*-commutative0.0%
unpow20.0%
rem-square-sqrt72.7%
mul-1-neg72.7%
Simplified72.7%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* -2.0 (/ (/ (pow l 2.0) (pow k_m 4.0)) (- t_m)))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (-2.0 * ((pow(l, 2.0) / pow(k_m, 4.0)) / -t_m));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((-2.0d0) * (((l ** 2.0d0) / (k_m ** 4.0d0)) / -t_m))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (-2.0 * ((Math.pow(l, 2.0) / Math.pow(k_m, 4.0)) / -t_m));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (-2.0 * ((math.pow(l, 2.0) / math.pow(k_m, 4.0)) / -t_m))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(-2.0 * Float64(Float64((l ^ 2.0) / (k_m ^ 4.0)) / Float64(-t_m)))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (-2.0 * (((l ^ 2.0) / (k_m ^ 4.0)) / -t_m)); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(-2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision] / (-t$95$m)), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(-2 \cdot \frac{\frac{{\ell}^{2}}{{k\_m}^{4}}}{-t\_m}\right)
\end{array}
Initial program 38.1%
*-commutative38.1%
associate-/r*38.1%
Simplified45.8%
add-sqr-sqrt20.1%
pow220.1%
sqrt-prod17.3%
sqrt-div18.5%
sqrt-pow122.7%
metadata-eval22.7%
sqrt-prod13.3%
add-sqr-sqrt26.6%
Applied egg-rr26.6%
Taylor expanded in k around 0 30.2%
Taylor expanded in t around -inf 0.0%
associate-/r*0.0%
*-commutative0.0%
unpow20.0%
rem-square-sqrt64.6%
mul-1-neg64.6%
Simplified64.6%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* (/ (pow l 2.0) (pow k_m 4.0)) (/ 2.0 t_m))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((pow(l, 2.0) / pow(k_m, 4.0)) * (2.0 / t_m));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (((l ** 2.0d0) / (k_m ** 4.0d0)) * (2.0d0 / t_m))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((Math.pow(l, 2.0) / Math.pow(k_m, 4.0)) * (2.0 / t_m));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * ((math.pow(l, 2.0) / math.pow(k_m, 4.0)) * (2.0 / t_m))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(Float64((l ^ 2.0) / (k_m ^ 4.0)) * Float64(2.0 / t_m))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (((l ^ 2.0) / (k_m ^ 4.0)) * (2.0 / t_m)); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision] * N[(2.0 / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\frac{{\ell}^{2}}{{k\_m}^{4}} \cdot \frac{2}{t\_m}\right)
\end{array}
Initial program 38.1%
*-commutative38.1%
associate-/r*38.1%
Simplified45.8%
add-sqr-sqrt20.1%
pow220.1%
sqrt-prod17.3%
sqrt-div18.5%
sqrt-pow122.7%
metadata-eval22.7%
sqrt-prod13.3%
add-sqr-sqrt26.6%
Applied egg-rr26.6%
Taylor expanded in k around 0 30.2%
Taylor expanded in k around 0 63.5%
associate-*r/63.5%
*-commutative63.5%
times-frac64.2%
Simplified64.2%
Final simplification64.2%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* (* l l) (/ 1.0 (* t_m (/ (pow k_m 4.0) 2.0))))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (1.0 / (t_m * (pow(k_m, 4.0) / 2.0))));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((l * l) * (1.0d0 / (t_m * ((k_m ** 4.0d0) / 2.0d0))))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (1.0 / (t_m * (Math.pow(k_m, 4.0) / 2.0))));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * ((l * l) * (1.0 / (t_m * (math.pow(k_m, 4.0) / 2.0))))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(Float64(l * l) * Float64(1.0 / Float64(t_m * Float64((k_m ^ 4.0) / 2.0))))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((l * l) * (1.0 / (t_m * ((k_m ^ 4.0) / 2.0)))); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(N[(l * l), $MachinePrecision] * N[(1.0 / N[(t$95$m * N[(N[Power[k$95$m, 4.0], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(\ell \cdot \ell\right) \cdot \frac{1}{t\_m \cdot \frac{{k\_m}^{4}}{2}}\right)
\end{array}
Initial program 38.1%
Simplified47.8%
Taylor expanded in k around 0 63.5%
clear-num63.5%
inv-pow63.5%
*-commutative63.5%
Applied egg-rr63.5%
unpow-163.5%
associate-/l*63.5%
Simplified63.5%
Final simplification63.5%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* (* l l) (/ 2.0 (* t_m (pow k_m 4.0))))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (2.0 / (t_m * pow(k_m, 4.0))));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((l * l) * (2.0d0 / (t_m * (k_m ** 4.0d0))))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (2.0 / (t_m * Math.pow(k_m, 4.0))));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * ((l * l) * (2.0 / (t_m * math.pow(k_m, 4.0))))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(Float64(l * l) * Float64(2.0 / Float64(t_m * (k_m ^ 4.0))))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((l * l) * (2.0 / (t_m * (k_m ^ 4.0)))); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(N[(l * l), $MachinePrecision] * N[(2.0 / N[(t$95$m * N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(\ell \cdot \ell\right) \cdot \frac{2}{t\_m \cdot {k\_m}^{4}}\right)
\end{array}
Initial program 38.1%
Simplified47.8%
Taylor expanded in k around 0 63.5%
Final simplification63.5%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* (* l l) (/ -0.11666666666666667 t_m))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (-0.11666666666666667 / t_m));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((l * l) * ((-0.11666666666666667d0) / t_m))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * ((l * l) * (-0.11666666666666667 / t_m));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * ((l * l) * (-0.11666666666666667 / t_m))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(Float64(l * l) * Float64(-0.11666666666666667 / t_m))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((l * l) * (-0.11666666666666667 / t_m)); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(N[(l * l), $MachinePrecision] * N[(-0.11666666666666667 / t$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\left(\ell \cdot \ell\right) \cdot \frac{-0.11666666666666667}{t\_m}\right)
\end{array}
Initial program 38.1%
Simplified47.8%
Taylor expanded in k around 0 48.4%
Taylor expanded in k around inf 20.4%
Final simplification20.4%
herbie shell --seed 2024088
(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))))