
(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 13 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 (* (sin k_m) (tan k_m))) (t_3 (/ (sqrt 2.0) k_m)))
(*
t_s
(if (<= k_m 3.2e-43)
(pow (* (* t_3 (/ l (sin k_m))) (sqrt (/ (cos k_m) t_m))) 2.0)
(if (<= k_m 5.4e+142)
(*
l
(*
2.0
(* (/ l (pow k_m 2.0)) (/ (cos k_m) (* t_m (pow (sin k_m) 2.0))))))
(*
(/ (* t_3 t_m) (pow (/ (* t_m (cbrt (/ t_2 l))) (cbrt l)) 2.0))
(/
(/ 1.0 (* k_m (/ (pow (cbrt l) -2.0) (sqrt 2.0))))
(cbrt t_2))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = sin(k_m) * tan(k_m);
double t_3 = sqrt(2.0) / k_m;
double tmp;
if (k_m <= 3.2e-43) {
tmp = pow(((t_3 * (l / sin(k_m))) * sqrt((cos(k_m) / t_m))), 2.0);
} else if (k_m <= 5.4e+142) {
tmp = l * (2.0 * ((l / pow(k_m, 2.0)) * (cos(k_m) / (t_m * pow(sin(k_m), 2.0)))));
} else {
tmp = ((t_3 * t_m) / pow(((t_m * cbrt((t_2 / l))) / cbrt(l)), 2.0)) * ((1.0 / (k_m * (pow(cbrt(l), -2.0) / sqrt(2.0)))) / cbrt(t_2));
}
return t_s * tmp;
}
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = Math.sin(k_m) * Math.tan(k_m);
double t_3 = Math.sqrt(2.0) / k_m;
double tmp;
if (k_m <= 3.2e-43) {
tmp = Math.pow(((t_3 * (l / Math.sin(k_m))) * Math.sqrt((Math.cos(k_m) / t_m))), 2.0);
} else if (k_m <= 5.4e+142) {
tmp = l * (2.0 * ((l / Math.pow(k_m, 2.0)) * (Math.cos(k_m) / (t_m * Math.pow(Math.sin(k_m), 2.0)))));
} else {
tmp = ((t_3 * t_m) / Math.pow(((t_m * Math.cbrt((t_2 / l))) / Math.cbrt(l)), 2.0)) * ((1.0 / (k_m * (Math.pow(Math.cbrt(l), -2.0) / Math.sqrt(2.0)))) / Math.cbrt(t_2));
}
return t_s * tmp;
}
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = Float64(sin(k_m) * tan(k_m)) t_3 = Float64(sqrt(2.0) / k_m) tmp = 0.0 if (k_m <= 3.2e-43) tmp = Float64(Float64(t_3 * Float64(l / sin(k_m))) * sqrt(Float64(cos(k_m) / t_m))) ^ 2.0; elseif (k_m <= 5.4e+142) tmp = Float64(l * Float64(2.0 * Float64(Float64(l / (k_m ^ 2.0)) * Float64(cos(k_m) / Float64(t_m * (sin(k_m) ^ 2.0)))))); else tmp = Float64(Float64(Float64(t_3 * t_m) / (Float64(Float64(t_m * cbrt(Float64(t_2 / l))) / cbrt(l)) ^ 2.0)) * Float64(Float64(1.0 / Float64(k_m * Float64((cbrt(l) ^ -2.0) / sqrt(2.0)))) / cbrt(t_2))); end return Float64(t_s * tmp) end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 3.2e-43], N[Power[N[(N[(t$95$3 * N[(l / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[k$95$m, 5.4e+142], N[(l * N[(2.0 * N[(N[(l / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k$95$m], $MachinePrecision] / N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(t$95$3 * t$95$m), $MachinePrecision] / N[Power[N[(N[(t$95$m * N[Power[N[(t$95$2 / l), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] / N[Power[l, 1/3], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 / N[(k$95$m * N[(N[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision] / N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[t$95$2, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sin k\_m \cdot \tan k\_m\\
t_3 := \frac{\sqrt{2}}{k\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 3.2 \cdot 10^{-43}:\\
\;\;\;\;{\left(\left(t\_3 \cdot \frac{\ell}{\sin k\_m}\right) \cdot \sqrt{\frac{\cos k\_m}{t\_m}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 5.4 \cdot 10^{+142}:\\
\;\;\;\;\ell \cdot \left(2 \cdot \left(\frac{\ell}{{k\_m}^{2}} \cdot \frac{\cos k\_m}{t\_m \cdot {\sin k\_m}^{2}}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_3 \cdot t\_m}{{\left(\frac{t\_m \cdot \sqrt[3]{\frac{t\_2}{\ell}}}{\sqrt[3]{\ell}}\right)}^{2}} \cdot \frac{\frac{1}{k\_m \cdot \frac{{\left(\sqrt[3]{\ell}\right)}^{-2}}{\sqrt{2}}}}{\sqrt[3]{t\_2}}\\
\end{array}
\end{array}
\end{array}
if k < 3.19999999999999985e-43Initial program 34.7%
*-commutative34.7%
associate-/r*34.7%
Simplified40.4%
add-sqr-sqrt28.7%
Applied egg-rr32.6%
unpow232.6%
associate-/r/32.6%
*-commutative32.6%
times-frac33.2%
Simplified33.2%
Taylor expanded in k around inf 47.4%
*-commutative47.4%
times-frac48.5%
Simplified48.5%
if 3.19999999999999985e-43 < k < 5.39999999999999965e142Initial program 26.5%
*-commutative26.5%
associate-/r*26.5%
Simplified49.1%
associate-*l/52.1%
associate-/r*54.9%
Applied egg-rr54.9%
associate-/r/57.9%
div-inv57.9%
+-rgt-identity57.9%
pow-flip57.9%
metadata-eval57.9%
associate-/l*57.9%
Applied egg-rr57.9%
Taylor expanded in k around inf 89.6%
times-frac99.0%
Applied egg-rr99.0%
if 5.39999999999999965e142 < k Initial program 46.4%
*-commutative46.4%
associate-/r*46.4%
Simplified53.6%
add-sqr-sqrt53.6%
add-cube-cbrt53.6%
times-frac53.6%
Applied egg-rr74.5%
associate-/r/74.5%
associate-/r*74.6%
associate-/r/74.6%
Simplified74.6%
add-cbrt-cube61.1%
unpow361.1%
unpow261.1%
cbrt-prod57.2%
cbrt-div57.2%
cbrt-prod57.2%
associate-*l/57.2%
associate-/l/61.1%
cbrt-div61.1%
associate-/l*61.1%
cbrt-prod61.1%
unpow361.1%
add-cbrt-cube74.7%
Applied egg-rr74.7%
clear-num74.7%
inv-pow74.7%
div-inv74.7%
pow-flip74.7%
metadata-eval74.7%
associate-*l/74.7%
Applied egg-rr74.7%
unpow-174.7%
associate-/r/74.8%
*-commutative74.8%
times-frac74.8%
*-inverses74.8%
Simplified74.8%
Final simplification57.5%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(let* ((t_2 (* (sin k_m) (tan k_m))) (t_3 (/ (sqrt 2.0) k_m)))
(*
t_s
(if (<= k_m 3.1e-43)
(pow (* (* t_3 (/ l (sin k_m))) (sqrt (/ (cos k_m) t_m))) 2.0)
(if (<= k_m 3.3e+142)
(*
l
(*
2.0
(* (/ l (pow k_m 2.0)) (/ (cos k_m) (* t_m (pow (sin k_m) 2.0))))))
(*
(/ (* t_3 t_m) (pow (/ (* t_m (cbrt (/ t_2 l))) (cbrt l)) 2.0))
(/ (/ t_3 (pow (cbrt l) -2.0)) (cbrt t_2))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = sin(k_m) * tan(k_m);
double t_3 = sqrt(2.0) / k_m;
double tmp;
if (k_m <= 3.1e-43) {
tmp = pow(((t_3 * (l / sin(k_m))) * sqrt((cos(k_m) / t_m))), 2.0);
} else if (k_m <= 3.3e+142) {
tmp = l * (2.0 * ((l / pow(k_m, 2.0)) * (cos(k_m) / (t_m * pow(sin(k_m), 2.0)))));
} else {
tmp = ((t_3 * t_m) / pow(((t_m * cbrt((t_2 / l))) / cbrt(l)), 2.0)) * ((t_3 / pow(cbrt(l), -2.0)) / cbrt(t_2));
}
return t_s * tmp;
}
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = Math.sin(k_m) * Math.tan(k_m);
double t_3 = Math.sqrt(2.0) / k_m;
double tmp;
if (k_m <= 3.1e-43) {
tmp = Math.pow(((t_3 * (l / Math.sin(k_m))) * Math.sqrt((Math.cos(k_m) / t_m))), 2.0);
} else if (k_m <= 3.3e+142) {
tmp = l * (2.0 * ((l / Math.pow(k_m, 2.0)) * (Math.cos(k_m) / (t_m * Math.pow(Math.sin(k_m), 2.0)))));
} else {
tmp = ((t_3 * t_m) / Math.pow(((t_m * Math.cbrt((t_2 / l))) / Math.cbrt(l)), 2.0)) * ((t_3 / Math.pow(Math.cbrt(l), -2.0)) / Math.cbrt(t_2));
}
return t_s * tmp;
}
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = Float64(sin(k_m) * tan(k_m)) t_3 = Float64(sqrt(2.0) / k_m) tmp = 0.0 if (k_m <= 3.1e-43) tmp = Float64(Float64(t_3 * Float64(l / sin(k_m))) * sqrt(Float64(cos(k_m) / t_m))) ^ 2.0; elseif (k_m <= 3.3e+142) tmp = Float64(l * Float64(2.0 * Float64(Float64(l / (k_m ^ 2.0)) * Float64(cos(k_m) / Float64(t_m * (sin(k_m) ^ 2.0)))))); else tmp = Float64(Float64(Float64(t_3 * t_m) / (Float64(Float64(t_m * cbrt(Float64(t_2 / l))) / cbrt(l)) ^ 2.0)) * Float64(Float64(t_3 / (cbrt(l) ^ -2.0)) / cbrt(t_2))); end return Float64(t_s * tmp) end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision]}, N[(t$95$s * If[LessEqual[k$95$m, 3.1e-43], N[Power[N[(N[(t$95$3 * N[(l / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[k$95$m, 3.3e+142], N[(l * N[(2.0 * N[(N[(l / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k$95$m], $MachinePrecision] / N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(t$95$3 * t$95$m), $MachinePrecision] / N[Power[N[(N[(t$95$m * N[Power[N[(t$95$2 / l), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] / N[Power[l, 1/3], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[(t$95$3 / N[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision] / N[Power[t$95$2, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sin k\_m \cdot \tan k\_m\\
t_3 := \frac{\sqrt{2}}{k\_m}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 3.1 \cdot 10^{-43}:\\
\;\;\;\;{\left(\left(t\_3 \cdot \frac{\ell}{\sin k\_m}\right) \cdot \sqrt{\frac{\cos k\_m}{t\_m}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 3.3 \cdot 10^{+142}:\\
\;\;\;\;\ell \cdot \left(2 \cdot \left(\frac{\ell}{{k\_m}^{2}} \cdot \frac{\cos k\_m}{t\_m \cdot {\sin k\_m}^{2}}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_3 \cdot t\_m}{{\left(\frac{t\_m \cdot \sqrt[3]{\frac{t\_2}{\ell}}}{\sqrt[3]{\ell}}\right)}^{2}} \cdot \frac{\frac{t\_3}{{\left(\sqrt[3]{\ell}\right)}^{-2}}}{\sqrt[3]{t\_2}}\\
\end{array}
\end{array}
\end{array}
if k < 3.0999999999999999e-43Initial program 34.7%
*-commutative34.7%
associate-/r*34.7%
Simplified40.4%
add-sqr-sqrt28.7%
Applied egg-rr32.6%
unpow232.6%
associate-/r/32.6%
*-commutative32.6%
times-frac33.2%
Simplified33.2%
Taylor expanded in k around inf 47.4%
*-commutative47.4%
times-frac48.5%
Simplified48.5%
if 3.0999999999999999e-43 < k < 3.3000000000000002e142Initial program 26.5%
*-commutative26.5%
associate-/r*26.5%
Simplified49.1%
associate-*l/52.1%
associate-/r*54.9%
Applied egg-rr54.9%
associate-/r/57.9%
div-inv57.9%
+-rgt-identity57.9%
pow-flip57.9%
metadata-eval57.9%
associate-/l*57.9%
Applied egg-rr57.9%
Taylor expanded in k around inf 89.6%
times-frac99.0%
Applied egg-rr99.0%
if 3.3000000000000002e142 < k Initial program 46.4%
*-commutative46.4%
associate-/r*46.4%
Simplified53.6%
add-sqr-sqrt53.6%
add-cube-cbrt53.6%
times-frac53.6%
Applied egg-rr74.5%
associate-/r/74.5%
associate-/r*74.6%
associate-/r/74.6%
Simplified74.6%
add-cbrt-cube61.1%
unpow361.1%
unpow261.1%
cbrt-prod57.2%
cbrt-div57.2%
cbrt-prod57.2%
associate-*l/57.2%
associate-/l/61.1%
cbrt-div61.1%
associate-/l*61.1%
cbrt-prod61.1%
unpow361.1%
add-cbrt-cube74.7%
Applied egg-rr74.7%
associate-/l*74.7%
div-inv74.7%
pow-flip74.7%
metadata-eval74.7%
Applied egg-rr74.7%
*-commutative74.7%
associate-/r*74.7%
*-inverses74.7%
associate-*l/74.7%
*-lft-identity74.7%
Simplified74.7%
Final simplification57.5%
k_m = (fabs.f64 k)
t\_m = (fabs.f64 t)
t\_s = (copysign.f64 #s(literal 1 binary64) t)
(FPCore (t_s t_m l k_m)
:precision binary64
(*
t_s
(if (<= k_m 3.2e-43)
(pow
(* (* (/ (sqrt 2.0) k_m) (/ l (sin k_m))) (sqrt (/ (cos k_m) t_m)))
2.0)
(if (<= k_m 5e+151)
(*
l
(*
2.0
(* (/ l (pow k_m 2.0)) (/ (cos k_m) (* t_m (pow (sin k_m) 2.0))))))
(pow
(/
(* (pow (cbrt l) 2.0) (cbrt 2.0))
(* t_m (cbrt (* (* (sin k_m) (tan k_m)) (pow (/ k_m t_m) 2.0)))))
3.0)))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 3.2e-43) {
tmp = pow((((sqrt(2.0) / k_m) * (l / sin(k_m))) * sqrt((cos(k_m) / t_m))), 2.0);
} else if (k_m <= 5e+151) {
tmp = l * (2.0 * ((l / pow(k_m, 2.0)) * (cos(k_m) / (t_m * pow(sin(k_m), 2.0)))));
} else {
tmp = pow(((pow(cbrt(l), 2.0) * cbrt(2.0)) / (t_m * cbrt(((sin(k_m) * tan(k_m)) * pow((k_m / t_m), 2.0))))), 3.0);
}
return t_s * tmp;
}
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 3.2e-43) {
tmp = Math.pow((((Math.sqrt(2.0) / k_m) * (l / Math.sin(k_m))) * Math.sqrt((Math.cos(k_m) / t_m))), 2.0);
} else if (k_m <= 5e+151) {
tmp = l * (2.0 * ((l / Math.pow(k_m, 2.0)) * (Math.cos(k_m) / (t_m * Math.pow(Math.sin(k_m), 2.0)))));
} else {
tmp = Math.pow(((Math.pow(Math.cbrt(l), 2.0) * Math.cbrt(2.0)) / (t_m * Math.cbrt(((Math.sin(k_m) * Math.tan(k_m)) * Math.pow((k_m / t_m), 2.0))))), 3.0);
}
return t_s * tmp;
}
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 3.2e-43) tmp = Float64(Float64(Float64(sqrt(2.0) / k_m) * Float64(l / sin(k_m))) * sqrt(Float64(cos(k_m) / t_m))) ^ 2.0; elseif (k_m <= 5e+151) tmp = Float64(l * Float64(2.0 * Float64(Float64(l / (k_m ^ 2.0)) * Float64(cos(k_m) / Float64(t_m * (sin(k_m) ^ 2.0)))))); else tmp = Float64(Float64((cbrt(l) ^ 2.0) * cbrt(2.0)) / Float64(t_m * cbrt(Float64(Float64(sin(k_m) * tan(k_m)) * (Float64(k_m / t_m) ^ 2.0))))) ^ 3.0; end return Float64(t_s * tmp) end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 3.2e-43], N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] * N[(l / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[k$95$m, 5e+151], N[(l * N[(2.0 * N[(N[(l / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k$95$m], $MachinePrecision] / N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[(N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision] * N[Power[2.0, 1/3], $MachinePrecision]), $MachinePrecision] / N[(t$95$m * N[Power[N[(N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision] * N[Power[N[(k$95$m / t$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 3.2 \cdot 10^{-43}:\\
\;\;\;\;{\left(\left(\frac{\sqrt{2}}{k\_m} \cdot \frac{\ell}{\sin k\_m}\right) \cdot \sqrt{\frac{\cos k\_m}{t\_m}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 5 \cdot 10^{+151}:\\
\;\;\;\;\ell \cdot \left(2 \cdot \left(\frac{\ell}{{k\_m}^{2}} \cdot \frac{\cos k\_m}{t\_m \cdot {\sin k\_m}^{2}}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{{\left(\sqrt[3]{\ell}\right)}^{2} \cdot \sqrt[3]{2}}{t\_m \cdot \sqrt[3]{\left(\sin k\_m \cdot \tan k\_m\right) \cdot {\left(\frac{k\_m}{t\_m}\right)}^{2}}}\right)}^{3}\\
\end{array}
\end{array}
if k < 3.19999999999999985e-43Initial program 34.7%
*-commutative34.7%
associate-/r*34.7%
Simplified40.4%
add-sqr-sqrt28.7%
Applied egg-rr32.6%
unpow232.6%
associate-/r/32.6%
*-commutative32.6%
times-frac33.2%
Simplified33.2%
Taylor expanded in k around inf 47.4%
*-commutative47.4%
times-frac48.5%
Simplified48.5%
if 3.19999999999999985e-43 < k < 5.0000000000000002e151Initial program 27.9%
*-commutative27.9%
associate-/r*28.0%
Simplified49.1%
associate-*l/52.0%
associate-/r*54.6%
Applied egg-rr54.6%
associate-/r/57.4%
div-inv57.4%
+-rgt-identity57.4%
pow-flip57.4%
metadata-eval57.4%
associate-/l*57.4%
Applied egg-rr57.4%
Taylor expanded in k around inf 90.2%
times-frac99.0%
Applied egg-rr99.0%
if 5.0000000000000002e151 < k Initial program 46.2%
Simplified53.8%
add-cube-cbrt53.8%
pow353.8%
Applied egg-rr58.6%
*-commutative58.6%
cbrt-prod58.6%
pow258.6%
cbrt-prod69.9%
pow269.9%
Applied egg-rr69.9%
Final simplification57.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 (<= k_m 3.2e-43)
(pow
(* (* (/ (sqrt 2.0) k_m) (/ l (sin k_m))) (sqrt (/ (cos k_m) t_m)))
2.0)
(if (<= k_m 1.8e+152)
(*
l
(*
2.0
(* (/ l (pow k_m 2.0)) (/ (cos k_m) (* t_m (pow (sin k_m) 2.0))))))
(pow
(/
(cbrt (* l (* 2.0 (pow (/ t_m k_m) 2.0))))
(* t_m (cbrt (* (sin k_m) (/ (tan k_m) l)))))
3.0)))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 3.2e-43) {
tmp = pow((((sqrt(2.0) / k_m) * (l / sin(k_m))) * sqrt((cos(k_m) / t_m))), 2.0);
} else if (k_m <= 1.8e+152) {
tmp = l * (2.0 * ((l / pow(k_m, 2.0)) * (cos(k_m) / (t_m * pow(sin(k_m), 2.0)))));
} else {
tmp = pow((cbrt((l * (2.0 * pow((t_m / k_m), 2.0)))) / (t_m * cbrt((sin(k_m) * (tan(k_m) / l))))), 3.0);
}
return t_s * tmp;
}
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 3.2e-43) {
tmp = Math.pow((((Math.sqrt(2.0) / k_m) * (l / Math.sin(k_m))) * Math.sqrt((Math.cos(k_m) / t_m))), 2.0);
} else if (k_m <= 1.8e+152) {
tmp = l * (2.0 * ((l / Math.pow(k_m, 2.0)) * (Math.cos(k_m) / (t_m * Math.pow(Math.sin(k_m), 2.0)))));
} else {
tmp = Math.pow((Math.cbrt((l * (2.0 * Math.pow((t_m / k_m), 2.0)))) / (t_m * Math.cbrt((Math.sin(k_m) * (Math.tan(k_m) / l))))), 3.0);
}
return t_s * tmp;
}
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if (k_m <= 3.2e-43) tmp = Float64(Float64(Float64(sqrt(2.0) / k_m) * Float64(l / sin(k_m))) * sqrt(Float64(cos(k_m) / t_m))) ^ 2.0; elseif (k_m <= 1.8e+152) tmp = Float64(l * Float64(2.0 * Float64(Float64(l / (k_m ^ 2.0)) * Float64(cos(k_m) / Float64(t_m * (sin(k_m) ^ 2.0)))))); else tmp = Float64(cbrt(Float64(l * Float64(2.0 * (Float64(t_m / k_m) ^ 2.0)))) / Float64(t_m * cbrt(Float64(sin(k_m) * Float64(tan(k_m) / l))))) ^ 3.0; end return Float64(t_s * tmp) end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[LessEqual[k$95$m, 3.2e-43], N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] * N[(l / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[k$95$m, 1.8e+152], N[(l * N[(2.0 * N[(N[(l / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k$95$m], $MachinePrecision] / N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(N[Power[N[(l * N[(2.0 * N[Power[N[(t$95$m / k$95$m), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] / N[(t$95$m * N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Tan[k$95$m], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;k\_m \leq 3.2 \cdot 10^{-43}:\\
\;\;\;\;{\left(\left(\frac{\sqrt{2}}{k\_m} \cdot \frac{\ell}{\sin k\_m}\right) \cdot \sqrt{\frac{\cos k\_m}{t\_m}}\right)}^{2}\\
\mathbf{elif}\;k\_m \leq 1.8 \cdot 10^{+152}:\\
\;\;\;\;\ell \cdot \left(2 \cdot \left(\frac{\ell}{{k\_m}^{2}} \cdot \frac{\cos k\_m}{t\_m \cdot {\sin k\_m}^{2}}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;{\left(\frac{\sqrt[3]{\ell \cdot \left(2 \cdot {\left(\frac{t\_m}{k\_m}\right)}^{2}\right)}}{t\_m \cdot \sqrt[3]{\sin k\_m \cdot \frac{\tan k\_m}{\ell}}}\right)}^{3}\\
\end{array}
\end{array}
if k < 3.19999999999999985e-43Initial program 34.7%
*-commutative34.7%
associate-/r*34.7%
Simplified40.4%
add-sqr-sqrt28.7%
Applied egg-rr32.6%
unpow232.6%
associate-/r/32.6%
*-commutative32.6%
times-frac33.2%
Simplified33.2%
Taylor expanded in k around inf 47.4%
*-commutative47.4%
times-frac48.5%
Simplified48.5%
if 3.19999999999999985e-43 < k < 1.7999999999999999e152Initial program 27.9%
*-commutative27.9%
associate-/r*28.0%
Simplified49.1%
associate-*l/52.0%
associate-/r*54.6%
Applied egg-rr54.6%
associate-/r/57.4%
div-inv57.4%
+-rgt-identity57.4%
pow-flip57.4%
metadata-eval57.4%
associate-/l*57.4%
Applied egg-rr57.4%
Taylor expanded in k around inf 90.2%
times-frac99.0%
Applied egg-rr99.0%
if 1.7999999999999999e152 < k Initial program 46.2%
*-commutative46.2%
associate-/r*46.2%
Simplified53.8%
associate-*l/53.8%
associate-/r*58.1%
Applied egg-rr58.1%
associate-/r/58.1%
div-inv58.1%
+-rgt-identity58.1%
pow-flip58.1%
metadata-eval58.1%
associate-/l*58.1%
Applied egg-rr58.1%
add-cube-cbrt58.1%
pow358.1%
Applied egg-rr69.9%
Final simplification57.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 (or (<= l 1.45e-259) (not (<= l 1.28e+154)))
(pow
(* (* (/ (sqrt 2.0) k_m) (/ l (sin k_m))) (sqrt (/ (cos k_m) t_m)))
2.0)
(*
l
(*
2.0
(* (/ l (pow k_m 2.0)) (/ (cos k_m) (* t_m (pow (sin k_m) 2.0)))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l <= 1.45e-259) || !(l <= 1.28e+154)) {
tmp = pow((((sqrt(2.0) / k_m) * (l / sin(k_m))) * sqrt((cos(k_m) / t_m))), 2.0);
} else {
tmp = l * (2.0 * ((l / pow(k_m, 2.0)) * (cos(k_m) / (t_m * pow(sin(k_m), 2.0)))));
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if ((l <= 1.45d-259) .or. (.not. (l <= 1.28d+154))) then
tmp = (((sqrt(2.0d0) / k_m) * (l / sin(k_m))) * sqrt((cos(k_m) / t_m))) ** 2.0d0
else
tmp = l * (2.0d0 * ((l / (k_m ** 2.0d0)) * (cos(k_m) / (t_m * (sin(k_m) ** 2.0d0)))))
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if ((l <= 1.45e-259) || !(l <= 1.28e+154)) {
tmp = Math.pow((((Math.sqrt(2.0) / k_m) * (l / Math.sin(k_m))) * Math.sqrt((Math.cos(k_m) / t_m))), 2.0);
} else {
tmp = l * (2.0 * ((l / Math.pow(k_m, 2.0)) * (Math.cos(k_m) / (t_m * Math.pow(Math.sin(k_m), 2.0)))));
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): tmp = 0 if (l <= 1.45e-259) or not (l <= 1.28e+154): tmp = math.pow((((math.sqrt(2.0) / k_m) * (l / math.sin(k_m))) * math.sqrt((math.cos(k_m) / t_m))), 2.0) else: tmp = l * (2.0 * ((l / math.pow(k_m, 2.0)) * (math.cos(k_m) / (t_m * math.pow(math.sin(k_m), 2.0))))) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) tmp = 0.0 if ((l <= 1.45e-259) || !(l <= 1.28e+154)) tmp = Float64(Float64(Float64(sqrt(2.0) / k_m) * Float64(l / sin(k_m))) * sqrt(Float64(cos(k_m) / t_m))) ^ 2.0; else tmp = Float64(l * Float64(2.0 * Float64(Float64(l / (k_m ^ 2.0)) * Float64(cos(k_m) / Float64(t_m * (sin(k_m) ^ 2.0)))))); end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) tmp = 0.0; if ((l <= 1.45e-259) || ~((l <= 1.28e+154))) tmp = (((sqrt(2.0) / k_m) * (l / sin(k_m))) * sqrt((cos(k_m) / t_m))) ^ 2.0; else tmp = l * (2.0 * ((l / (k_m ^ 2.0)) * (cos(k_m) / (t_m * (sin(k_m) ^ 2.0))))); end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * If[Or[LessEqual[l, 1.45e-259], N[Not[LessEqual[l, 1.28e+154]], $MachinePrecision]], N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] * N[(l / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(l * N[(2.0 * N[(N[(l / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k$95$m], $MachinePrecision] / N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \leq 1.45 \cdot 10^{-259} \lor \neg \left(\ell \leq 1.28 \cdot 10^{+154}\right):\\
\;\;\;\;{\left(\left(\frac{\sqrt{2}}{k\_m} \cdot \frac{\ell}{\sin k\_m}\right) \cdot \sqrt{\frac{\cos k\_m}{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \left(2 \cdot \left(\frac{\ell}{{k\_m}^{2}} \cdot \frac{\cos k\_m}{t\_m \cdot {\sin k\_m}^{2}}\right)\right)\\
\end{array}
\end{array}
if l < 1.45000000000000004e-259 or 1.2800000000000001e154 < l Initial program 35.7%
*-commutative35.7%
associate-/r*35.7%
Simplified42.0%
add-sqr-sqrt31.6%
Applied egg-rr27.4%
unpow227.4%
associate-/r/27.4%
*-commutative27.4%
times-frac27.6%
Simplified27.6%
Taylor expanded in k around inf 48.2%
*-commutative48.2%
times-frac49.4%
Simplified49.4%
if 1.45000000000000004e-259 < l < 1.2800000000000001e154Initial program 33.3%
*-commutative33.3%
associate-/r*33.3%
Simplified44.8%
associate-*l/46.1%
associate-/r*51.0%
Applied egg-rr51.0%
associate-/r/52.2%
div-inv52.2%
+-rgt-identity52.2%
pow-flip52.2%
metadata-eval52.2%
associate-/l*53.4%
Applied egg-rr53.4%
Taylor expanded in k around inf 92.4%
times-frac94.1%
Applied egg-rr94.1%
Final simplification63.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
(let* ((t_2 (sqrt (/ (cos k_m) t_m))))
(*
t_s
(if (<= l 9.5e-260)
(pow (* (* (/ (sqrt 2.0) k_m) (/ l (sin k_m))) t_2) 2.0)
(if (<= l 2.45e+153)
(*
l
(*
2.0
(* (/ l (pow k_m 2.0)) (/ (cos k_m) (* t_m (pow (sin k_m) 2.0))))))
(pow (* t_2 (/ (* (sqrt 2.0) l) (* k_m (sin k_m)))) 2.0))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double t_2 = sqrt((cos(k_m) / t_m));
double tmp;
if (l <= 9.5e-260) {
tmp = pow((((sqrt(2.0) / k_m) * (l / sin(k_m))) * t_2), 2.0);
} else if (l <= 2.45e+153) {
tmp = l * (2.0 * ((l / pow(k_m, 2.0)) * (cos(k_m) / (t_m * pow(sin(k_m), 2.0)))));
} else {
tmp = pow((t_2 * ((sqrt(2.0) * l) / (k_m * sin(k_m)))), 2.0);
}
return t_s * tmp;
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: t_2
real(8) :: tmp
t_2 = sqrt((cos(k_m) / t_m))
if (l <= 9.5d-260) then
tmp = (((sqrt(2.0d0) / k_m) * (l / sin(k_m))) * t_2) ** 2.0d0
else if (l <= 2.45d+153) then
tmp = l * (2.0d0 * ((l / (k_m ** 2.0d0)) * (cos(k_m) / (t_m * (sin(k_m) ** 2.0d0)))))
else
tmp = (t_2 * ((sqrt(2.0d0) * l) / (k_m * sin(k_m)))) ** 2.0d0
end if
code = t_s * tmp
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
double t_2 = Math.sqrt((Math.cos(k_m) / t_m));
double tmp;
if (l <= 9.5e-260) {
tmp = Math.pow((((Math.sqrt(2.0) / k_m) * (l / Math.sin(k_m))) * t_2), 2.0);
} else if (l <= 2.45e+153) {
tmp = l * (2.0 * ((l / Math.pow(k_m, 2.0)) * (Math.cos(k_m) / (t_m * Math.pow(Math.sin(k_m), 2.0)))));
} else {
tmp = Math.pow((t_2 * ((Math.sqrt(2.0) * l) / (k_m * Math.sin(k_m)))), 2.0);
}
return t_s * tmp;
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): t_2 = math.sqrt((math.cos(k_m) / t_m)) tmp = 0 if l <= 9.5e-260: tmp = math.pow((((math.sqrt(2.0) / k_m) * (l / math.sin(k_m))) * t_2), 2.0) elif l <= 2.45e+153: tmp = l * (2.0 * ((l / math.pow(k_m, 2.0)) * (math.cos(k_m) / (t_m * math.pow(math.sin(k_m), 2.0))))) else: tmp = math.pow((t_2 * ((math.sqrt(2.0) * l) / (k_m * math.sin(k_m)))), 2.0) return t_s * tmp
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) t_2 = sqrt(Float64(cos(k_m) / t_m)) tmp = 0.0 if (l <= 9.5e-260) tmp = Float64(Float64(Float64(sqrt(2.0) / k_m) * Float64(l / sin(k_m))) * t_2) ^ 2.0; elseif (l <= 2.45e+153) tmp = Float64(l * Float64(2.0 * Float64(Float64(l / (k_m ^ 2.0)) * Float64(cos(k_m) / Float64(t_m * (sin(k_m) ^ 2.0)))))); else tmp = Float64(t_2 * Float64(Float64(sqrt(2.0) * l) / Float64(k_m * sin(k_m)))) ^ 2.0; end return Float64(t_s * tmp) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp_2 = code(t_s, t_m, l, k_m) t_2 = sqrt((cos(k_m) / t_m)); tmp = 0.0; if (l <= 9.5e-260) tmp = (((sqrt(2.0) / k_m) * (l / sin(k_m))) * t_2) ^ 2.0; elseif (l <= 2.45e+153) tmp = l * (2.0 * ((l / (k_m ^ 2.0)) * (cos(k_m) / (t_m * (sin(k_m) ^ 2.0))))); else tmp = (t_2 * ((sqrt(2.0) * l) / (k_m * sin(k_m)))) ^ 2.0; end tmp_2 = t_s * tmp; end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := Block[{t$95$2 = N[Sqrt[N[(N[Cos[k$95$m], $MachinePrecision] / t$95$m), $MachinePrecision]], $MachinePrecision]}, N[(t$95$s * If[LessEqual[l, 9.5e-260], N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] / k$95$m), $MachinePrecision] * N[(l / N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision], 2.0], $MachinePrecision], If[LessEqual[l, 2.45e+153], N[(l * N[(2.0 * N[(N[(l / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k$95$m], $MachinePrecision] / N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Power[N[(t$95$2 * N[(N[(N[Sqrt[2.0], $MachinePrecision] * l), $MachinePrecision] / N[(k$95$m * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision]]]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
\begin{array}{l}
t_2 := \sqrt{\frac{\cos k\_m}{t\_m}}\\
t\_s \cdot \begin{array}{l}
\mathbf{if}\;\ell \leq 9.5 \cdot 10^{-260}:\\
\;\;\;\;{\left(\left(\frac{\sqrt{2}}{k\_m} \cdot \frac{\ell}{\sin k\_m}\right) \cdot t\_2\right)}^{2}\\
\mathbf{elif}\;\ell \leq 2.45 \cdot 10^{+153}:\\
\;\;\;\;\ell \cdot \left(2 \cdot \left(\frac{\ell}{{k\_m}^{2}} \cdot \frac{\cos k\_m}{t\_m \cdot {\sin k\_m}^{2}}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;{\left(t\_2 \cdot \frac{\sqrt{2} \cdot \ell}{k\_m \cdot \sin k\_m}\right)}^{2}\\
\end{array}
\end{array}
\end{array}
if l < 9.5000000000000001e-260Initial program 35.3%
*-commutative35.3%
associate-/r*35.3%
Simplified43.1%
add-sqr-sqrt33.1%
Applied egg-rr25.7%
unpow225.7%
associate-/r/25.7%
*-commutative25.7%
times-frac25.8%
Simplified25.8%
Taylor expanded in k around inf 46.0%
*-commutative46.0%
times-frac47.6%
Simplified47.6%
if 9.5000000000000001e-260 < l < 2.45000000000000001e153Initial program 33.3%
*-commutative33.3%
associate-/r*33.3%
Simplified44.8%
associate-*l/46.1%
associate-/r*51.0%
Applied egg-rr51.0%
associate-/r/52.2%
div-inv52.2%
+-rgt-identity52.2%
pow-flip52.2%
metadata-eval52.2%
associate-/l*53.4%
Applied egg-rr53.4%
Taylor expanded in k around inf 92.4%
times-frac94.1%
Applied egg-rr94.1%
if 2.45000000000000001e153 < l Initial program 37.2%
*-commutative37.2%
associate-/r*37.2%
Simplified37.2%
add-sqr-sqrt25.7%
Applied egg-rr34.6%
unpow234.6%
associate-/r/34.6%
*-commutative34.6%
times-frac34.6%
Simplified34.6%
Taylor expanded in k around inf 56.8%
Final simplification63.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 1.6e-40)
(pow (* (/ (* (sqrt 2.0) l) (pow k_m 2.0)) (sqrt (/ 1.0 t_m))) 2.0)
(*
l
(*
2.0
(* l (/ (cos k_m) (* (pow (sin k_m) 2.0) (* t_m (pow k_m 2.0))))))))))k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
double tmp;
if (k_m <= 1.6e-40) {
tmp = pow((((sqrt(2.0) * l) / pow(k_m, 2.0)) * sqrt((1.0 / t_m))), 2.0);
} else {
tmp = l * (2.0 * (l * (cos(k_m) / (pow(sin(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 (k_m <= 1.6d-40) then
tmp = (((sqrt(2.0d0) * l) / (k_m ** 2.0d0)) * sqrt((1.0d0 / t_m))) ** 2.0d0
else
tmp = l * (2.0d0 * (l * (cos(k_m) / ((sin(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 (k_m <= 1.6e-40) {
tmp = Math.pow((((Math.sqrt(2.0) * l) / Math.pow(k_m, 2.0)) * Math.sqrt((1.0 / t_m))), 2.0);
} else {
tmp = l * (2.0 * (l * (Math.cos(k_m) / (Math.pow(Math.sin(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 k_m <= 1.6e-40: tmp = math.pow((((math.sqrt(2.0) * l) / math.pow(k_m, 2.0)) * math.sqrt((1.0 / t_m))), 2.0) else: tmp = l * (2.0 * (l * (math.cos(k_m) / (math.pow(math.sin(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 (k_m <= 1.6e-40) tmp = Float64(Float64(Float64(sqrt(2.0) * l) / (k_m ^ 2.0)) * sqrt(Float64(1.0 / t_m))) ^ 2.0; else tmp = Float64(l * Float64(2.0 * Float64(l * Float64(cos(k_m) / Float64((sin(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 (k_m <= 1.6e-40) tmp = (((sqrt(2.0) * l) / (k_m ^ 2.0)) * sqrt((1.0 / t_m))) ^ 2.0; else tmp = l * (2.0 * (l * (cos(k_m) / ((sin(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[k$95$m, 1.6e-40], N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] * l), $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision], N[(l * N[(2.0 * N[(l * N[(N[Cos[k$95$m], $MachinePrecision] / N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] * N[(t$95$m * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $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 1.6 \cdot 10^{-40}:\\
\;\;\;\;{\left(\frac{\sqrt{2} \cdot \ell}{{k\_m}^{2}} \cdot \sqrt{\frac{1}{t\_m}}\right)}^{2}\\
\mathbf{else}:\\
\;\;\;\;\ell \cdot \left(2 \cdot \left(\ell \cdot \frac{\cos k\_m}{{\sin k\_m}^{2} \cdot \left(t\_m \cdot {k\_m}^{2}\right)}\right)\right)\\
\end{array}
\end{array}
if k < 1.60000000000000001e-40Initial program 34.5%
*-commutative34.5%
associate-/r*34.5%
Simplified40.2%
add-sqr-sqrt28.5%
Applied egg-rr32.4%
unpow232.4%
associate-/r/32.4%
*-commutative32.4%
times-frac33.0%
Simplified33.0%
Taylor expanded in k around 0 40.9%
if 1.60000000000000001e-40 < k Initial program 36.6%
*-commutative36.6%
associate-/r*36.6%
Simplified52.1%
associate-*l/53.7%
associate-/r*57.1%
Applied egg-rr57.1%
associate-/r/58.7%
div-inv58.7%
+-rgt-identity58.7%
pow-flip58.7%
metadata-eval58.7%
associate-/l*58.7%
Applied egg-rr58.7%
Taylor expanded in k around inf 77.8%
associate-/l*77.8%
*-commutative77.8%
*-commutative77.8%
associate-*l*77.8%
*-commutative77.8%
Simplified77.8%
Final simplification49.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 (* l (* 2.0 (* (/ l (pow k_m 2.0)) (/ (cos k_m) (* t_m (pow (sin k_m) 2.0))))))))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * (l * (2.0 * ((l / pow(k_m, 2.0)) * (cos(k_m) / (t_m * pow(sin(k_m), 2.0))))));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (l * (2.0d0 * ((l / (k_m ** 2.0d0)) * (cos(k_m) / (t_m * (sin(k_m) ** 2.0d0))))))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (l * (2.0 * ((l / Math.pow(k_m, 2.0)) * (Math.cos(k_m) / (t_m * Math.pow(Math.sin(k_m), 2.0))))));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (l * (2.0 * ((l / math.pow(k_m, 2.0)) * (math.cos(k_m) / (t_m * math.pow(math.sin(k_m), 2.0))))))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(l * Float64(2.0 * Float64(Float64(l / (k_m ^ 2.0)) * Float64(cos(k_m) / Float64(t_m * (sin(k_m) ^ 2.0))))))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (l * (2.0 * ((l / (k_m ^ 2.0)) * (cos(k_m) / (t_m * (sin(k_m) ^ 2.0)))))); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(l * N[(2.0 * N[(N[(l / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k$95$m], $MachinePrecision] / N[(t$95$m * N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\ell \cdot \left(2 \cdot \left(\frac{\ell}{{k\_m}^{2}} \cdot \frac{\cos k\_m}{t\_m \cdot {\sin k\_m}^{2}}\right)\right)\right)
\end{array}
Initial program 35.0%
*-commutative35.0%
associate-/r*35.0%
Simplified42.9%
associate-*l/44.1%
associate-/r*49.7%
Applied egg-rr49.7%
associate-/r/50.4%
div-inv50.4%
+-rgt-identity50.4%
pow-flip50.4%
metadata-eval50.4%
associate-/l*50.8%
Applied egg-rr50.8%
Taylor expanded in k around inf 84.3%
times-frac87.4%
Applied egg-rr87.4%
Final simplification87.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 (pow (* (/ (* (sqrt 2.0) l) (pow k_m 2.0)) (sqrt (/ 1.0 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) {
return t_s * pow((((sqrt(2.0) * l) / pow(k_m, 2.0)) * sqrt((1.0 / t_m))), 2.0);
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * ((((sqrt(2.0d0) * l) / (k_m ** 2.0d0)) * sqrt((1.0d0 / t_m))) ** 2.0d0)
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * Math.pow((((Math.sqrt(2.0) * l) / Math.pow(k_m, 2.0)) * Math.sqrt((1.0 / t_m))), 2.0);
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * math.pow((((math.sqrt(2.0) * l) / math.pow(k_m, 2.0)) * math.sqrt((1.0 / t_m))), 2.0)
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * (Float64(Float64(Float64(sqrt(2.0) * l) / (k_m ^ 2.0)) * sqrt(Float64(1.0 / t_m))) ^ 2.0)) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * ((((sqrt(2.0) * l) / (k_m ^ 2.0)) * sqrt((1.0 / t_m))) ^ 2.0); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[Power[N[(N[(N[(N[Sqrt[2.0], $MachinePrecision] * l), $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(1.0 / t$95$m), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 2.0], $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{\sqrt{2} \cdot \ell}{{k\_m}^{2}} \cdot \sqrt{\frac{1}{t\_m}}\right)}^{2}
\end{array}
Initial program 35.0%
*-commutative35.0%
associate-/r*35.0%
Simplified42.9%
add-sqr-sqrt32.7%
Applied egg-rr29.0%
unpow229.0%
associate-/r/29.0%
*-commutative29.0%
times-frac29.5%
Simplified29.5%
Taylor expanded in k around 0 39.2%
Final simplification39.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 (pow (* l (* (sqrt (/ 1.0 t_m)) (/ (sqrt 2.0) (pow k_m 2.0)))) 2.0)))
k_m = fabs(k);
t\_m = fabs(t);
t\_s = copysign(1.0, t);
double code(double t_s, double t_m, double l, double k_m) {
return t_s * pow((l * (sqrt((1.0 / t_m)) * (sqrt(2.0) / pow(k_m, 2.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 * (sqrt((1.0d0 / t_m)) * (sqrt(2.0d0) / (k_m ** 2.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 * Math.pow((l * (Math.sqrt((1.0 / t_m)) * (Math.sqrt(2.0) / Math.pow(k_m, 2.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 * math.pow((l * (math.sqrt((1.0 / t_m)) * (math.sqrt(2.0) / math.pow(k_m, 2.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(l * Float64(sqrt(Float64(1.0 / t_m)) * Float64(sqrt(2.0) / (k_m ^ 2.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 * (sqrt((1.0 / t_m)) * (sqrt(2.0) / (k_m ^ 2.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[Power[N[(l * N[(N[Sqrt[N[(1.0 / t$95$m), $MachinePrecision]], $MachinePrecision] * N[(N[Sqrt[2.0], $MachinePrecision] / N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 2.0], $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(\ell \cdot \left(\sqrt{\frac{1}{t\_m}} \cdot \frac{\sqrt{2}}{{k\_m}^{2}}\right)\right)}^{2}
\end{array}
Initial program 35.0%
*-commutative35.0%
associate-/r*35.0%
Simplified42.9%
add-sqr-sqrt32.7%
Applied egg-rr29.0%
unpow229.0%
associate-/r/29.0%
*-commutative29.0%
times-frac29.5%
Simplified29.5%
Taylor expanded in k around 0 39.2%
associate-/l*38.9%
associate-*l*38.9%
Simplified38.9%
Final simplification38.9%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* l (* 2.0 (/ (* l (cos k_m)) (* (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) {
return t_s * (l * (2.0 * ((l * cos(k_m)) / (pow(k_m, 2.0) * (t_m * pow(k_m, 2.0))))));
}
k_m = abs(k)
t\_m = abs(t)
t\_s = copysign(1.0d0, t)
real(8) function code(t_s, t_m, l, k_m)
real(8), intent (in) :: t_s
real(8), intent (in) :: t_m
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = t_s * (l * (2.0d0 * ((l * cos(k_m)) / ((k_m ** 2.0d0) * (t_m * (k_m ** 2.0d0))))))
end function
k_m = Math.abs(k);
t\_m = Math.abs(t);
t\_s = Math.copySign(1.0, t);
public static double code(double t_s, double t_m, double l, double k_m) {
return t_s * (l * (2.0 * ((l * Math.cos(k_m)) / (Math.pow(k_m, 2.0) * (t_m * Math.pow(k_m, 2.0))))));
}
k_m = math.fabs(k) t\_m = math.fabs(t) t\_s = math.copysign(1.0, t) def code(t_s, t_m, l, k_m): return t_s * (l * (2.0 * ((l * math.cos(k_m)) / (math.pow(k_m, 2.0) * (t_m * math.pow(k_m, 2.0))))))
k_m = abs(k) t\_m = abs(t) t\_s = copysign(1.0, t) function code(t_s, t_m, l, k_m) return Float64(t_s * Float64(l * Float64(2.0 * Float64(Float64(l * cos(k_m)) / Float64((k_m ^ 2.0) * Float64(t_m * (k_m ^ 2.0))))))) end
k_m = abs(k); t\_m = abs(t); t\_s = sign(t) * abs(1.0); function tmp = code(t_s, t_m, l, k_m) tmp = t_s * (l * (2.0 * ((l * cos(k_m)) / ((k_m ^ 2.0) * (t_m * (k_m ^ 2.0)))))); end
k_m = N[Abs[k], $MachinePrecision]
t\_m = N[Abs[t], $MachinePrecision]
t\_s = N[With[{TMP1 = Abs[1.0], TMP2 = Sign[t]}, TMP1 * If[TMP2 == 0, 1, TMP2]], $MachinePrecision]
code[t$95$s_, t$95$m_, l_, k$95$m_] := N[(t$95$s * N[(l * N[(2.0 * N[(N[(l * N[Cos[k$95$m], $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]), $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(\ell \cdot \left(2 \cdot \frac{\ell \cdot \cos k\_m}{{k\_m}^{2} \cdot \left(t\_m \cdot {k\_m}^{2}\right)}\right)\right)
\end{array}
Initial program 35.0%
*-commutative35.0%
associate-/r*35.0%
Simplified42.9%
associate-*l/44.1%
associate-/r*49.7%
Applied egg-rr49.7%
associate-/r/50.4%
div-inv50.4%
+-rgt-identity50.4%
pow-flip50.4%
metadata-eval50.4%
associate-/l*50.8%
Applied egg-rr50.8%
Taylor expanded in k around inf 84.3%
Taylor expanded in k around 0 75.1%
Final simplification75.1%
k_m = (fabs.f64 k) t\_m = (fabs.f64 t) t\_s = (copysign.f64 #s(literal 1 binary64) t) (FPCore (t_s t_m l k_m) :precision binary64 (* t_s (* l (* 2.0 (/ (/ l (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 * (l * (2.0 * ((l / 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 * (l * (2.0d0 * ((l / (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 * (l * (2.0 * ((l / 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 * (l * (2.0 * ((l / 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(l * Float64(2.0 * Float64(Float64(l / (k_m ^ 4.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 * ((l / (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[(l * N[(2.0 * N[(N[(l / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision] / t$95$m), $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(\ell \cdot \left(2 \cdot \frac{\frac{\ell}{{k\_m}^{4}}}{t\_m}\right)\right)
\end{array}
Initial program 35.0%
*-commutative35.0%
associate-/r*35.0%
Simplified42.9%
associate-*l/44.1%
associate-/r*49.7%
Applied egg-rr49.7%
associate-/r/50.4%
div-inv50.4%
+-rgt-identity50.4%
pow-flip50.4%
metadata-eval50.4%
associate-/l*50.8%
Applied egg-rr50.8%
Taylor expanded in k around 0 72.3%
associate-/r*72.8%
Simplified72.8%
Final simplification72.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 (* l (* 2.0 (/ l (* 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 * (2.0 * (l / (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 * (2.0d0 * (l / (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 * (2.0 * (l / (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 * (2.0 * (l / (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(l * Float64(2.0 * Float64(l / 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 * (2.0 * (l / (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[(l * N[(2.0 * N[(l / N[(t$95$m * N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
t\_m = \left|t\right|
\\
t\_s = \mathsf{copysign}\left(1, t\right)
\\
t\_s \cdot \left(\ell \cdot \left(2 \cdot \frac{\ell}{t\_m \cdot {k\_m}^{4}}\right)\right)
\end{array}
Initial program 35.0%
*-commutative35.0%
associate-/r*35.0%
Simplified42.9%
associate-*l/44.1%
associate-/r*49.7%
Applied egg-rr49.7%
associate-/r/50.4%
div-inv50.4%
+-rgt-identity50.4%
pow-flip50.4%
metadata-eval50.4%
associate-/l*50.8%
Applied egg-rr50.8%
Taylor expanded in k around 0 72.3%
Final simplification72.3%
herbie shell --seed 2024101
(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))))