
(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 20 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)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (pow (cbrt l) 2.0)))
(if (<= k_m 1.9e-85)
(/ 2.0 (* (pow (/ (* t (cbrt (sin k_m))) t_1) 3.0) (* 2.0 k_m)))
(if (or (<= k_m 0.02) (not (<= k_m 1.12e+154)))
(/
2.0
(pow
(*
(/ t t_1)
(cbrt (* (sin k_m) (* (tan k_m) (+ 2.0 (pow (/ k_m t) 2.0))))))
3.0))
(*
(/ 2.0 (* t (pow k_m 2.0)))
(/ (* (pow l 2.0) (cos k_m)) (pow (sin k_m) 2.0)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = pow(cbrt(l), 2.0);
double tmp;
if (k_m <= 1.9e-85) {
tmp = 2.0 / (pow(((t * cbrt(sin(k_m))) / t_1), 3.0) * (2.0 * k_m));
} else if ((k_m <= 0.02) || !(k_m <= 1.12e+154)) {
tmp = 2.0 / pow(((t / t_1) * cbrt((sin(k_m) * (tan(k_m) * (2.0 + pow((k_m / t), 2.0)))))), 3.0);
} else {
tmp = (2.0 / (t * pow(k_m, 2.0))) * ((pow(l, 2.0) * cos(k_m)) / pow(sin(k_m), 2.0));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double t_1 = Math.pow(Math.cbrt(l), 2.0);
double tmp;
if (k_m <= 1.9e-85) {
tmp = 2.0 / (Math.pow(((t * Math.cbrt(Math.sin(k_m))) / t_1), 3.0) * (2.0 * k_m));
} else if ((k_m <= 0.02) || !(k_m <= 1.12e+154)) {
tmp = 2.0 / Math.pow(((t / t_1) * Math.cbrt((Math.sin(k_m) * (Math.tan(k_m) * (2.0 + Math.pow((k_m / t), 2.0)))))), 3.0);
} else {
tmp = (2.0 / (t * Math.pow(k_m, 2.0))) * ((Math.pow(l, 2.0) * Math.cos(k_m)) / Math.pow(Math.sin(k_m), 2.0));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = cbrt(l) ^ 2.0 tmp = 0.0 if (k_m <= 1.9e-85) tmp = Float64(2.0 / Float64((Float64(Float64(t * cbrt(sin(k_m))) / t_1) ^ 3.0) * Float64(2.0 * k_m))); elseif ((k_m <= 0.02) || !(k_m <= 1.12e+154)) tmp = Float64(2.0 / (Float64(Float64(t / t_1) * cbrt(Float64(sin(k_m) * Float64(tan(k_m) * Float64(2.0 + (Float64(k_m / t) ^ 2.0)))))) ^ 3.0)); else tmp = Float64(Float64(2.0 / Float64(t * (k_m ^ 2.0))) * Float64(Float64((l ^ 2.0) * cos(k_m)) / (sin(k_m) ^ 2.0))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[k$95$m, 1.9e-85], N[(2.0 / N[(N[Power[N[(N[(t * N[Power[N[Sin[k$95$m], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision], 3.0], $MachinePrecision] * N[(2.0 * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[k$95$m, 0.02], N[Not[LessEqual[k$95$m, 1.12e+154]], $MachinePrecision]], N[(2.0 / N[Power[N[(N[(t / t$95$1), $MachinePrecision] * N[Power[N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Tan[k$95$m], $MachinePrecision] * N[(2.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(t * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := {\left(\sqrt[3]{\ell}\right)}^{2}\\
\mathbf{if}\;k\_m \leq 1.9 \cdot 10^{-85}:\\
\;\;\;\;\frac{2}{{\left(\frac{t \cdot \sqrt[3]{\sin k\_m}}{t\_1}\right)}^{3} \cdot \left(2 \cdot k\_m\right)}\\
\mathbf{elif}\;k\_m \leq 0.02 \lor \neg \left(k\_m \leq 1.12 \cdot 10^{+154}\right):\\
\;\;\;\;\frac{2}{{\left(\frac{t}{t\_1} \cdot \sqrt[3]{\sin k\_m \cdot \left(\tan k\_m \cdot \left(2 + {\left(\frac{k\_m}{t}\right)}^{2}\right)\right)}\right)}^{3}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{t \cdot {k\_m}^{2}} \cdot \frac{{\ell}^{2} \cdot \cos k\_m}{{\sin k\_m}^{2}}\\
\end{array}
\end{array}
if k < 1.8999999999999999e-85Initial program 52.6%
Simplified52.6%
associate-/r*60.0%
add-cube-cbrt59.9%
*-un-lft-identity59.9%
times-frac59.9%
pow259.9%
cbrt-div59.9%
rem-cbrt-cube59.9%
cbrt-div59.9%
rem-cbrt-cube71.8%
Applied egg-rr71.8%
add-cube-cbrt71.7%
pow371.7%
Applied egg-rr78.7%
*-commutative78.7%
Simplified78.7%
associate-*l/78.7%
Applied egg-rr78.7%
Taylor expanded in k around 0 74.3%
if 1.8999999999999999e-85 < k < 0.0200000000000000004 or 1.11999999999999994e154 < k Initial program 53.9%
Simplified53.9%
associate-/r*57.2%
add-cube-cbrt57.1%
*-un-lft-identity57.1%
times-frac57.1%
pow257.1%
cbrt-div57.2%
rem-cbrt-cube57.1%
cbrt-div57.1%
rem-cbrt-cube72.8%
Applied egg-rr72.8%
add-cube-cbrt72.8%
pow372.8%
Applied egg-rr80.1%
if 0.0200000000000000004 < k < 1.11999999999999994e154Initial program 63.6%
associate-/r*63.6%
+-commutative63.6%
unpow263.6%
sqr-neg63.6%
distribute-frac-neg263.6%
distribute-frac-neg263.6%
unpow263.6%
+-commutative63.6%
associate-*l*63.7%
associate-*l/63.7%
associate-/r/63.7%
+-commutative63.7%
associate-+r+63.7%
Simplified63.7%
add-cube-cbrt63.6%
pow363.5%
cbrt-prod63.6%
rem-cbrt-cube63.5%
Applied egg-rr63.5%
*-un-lft-identity63.5%
unpow-prod-down63.5%
pow363.5%
add-cube-cbrt63.7%
Applied egg-rr63.7%
*-lft-identity63.7%
*-commutative63.7%
associate-*l*63.7%
Simplified63.7%
Taylor expanded in k around inf 77.4%
associate-*r/77.4%
associate-*r*77.5%
times-frac77.5%
Simplified77.5%
Final simplification75.8%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (pow (/ k_m t) 2.0)))
(if (<=
(*
(* (tan k_m) (* (sin k_m) (/ (pow t 3.0) (* l l))))
(+ 1.0 (+ 1.0 t_1)))
INFINITY)
(/
2.0
(*
(* (tan k_m) (+ 2.0 t_1))
(* (sin k_m) (* (/ (pow t 2.0) l) (/ t l)))))
(/
2.0
(pow (* (/ t (pow (cbrt l) 2.0)) (cbrt (* 2.0 (pow k_m 2.0)))) 3.0)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = pow((k_m / t), 2.0);
double tmp;
if (((tan(k_m) * (sin(k_m) * (pow(t, 3.0) / (l * l)))) * (1.0 + (1.0 + t_1))) <= ((double) INFINITY)) {
tmp = 2.0 / ((tan(k_m) * (2.0 + t_1)) * (sin(k_m) * ((pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 / pow(((t / pow(cbrt(l), 2.0)) * cbrt((2.0 * pow(k_m, 2.0)))), 3.0);
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double t_1 = Math.pow((k_m / t), 2.0);
double tmp;
if (((Math.tan(k_m) * (Math.sin(k_m) * (Math.pow(t, 3.0) / (l * l)))) * (1.0 + (1.0 + t_1))) <= Double.POSITIVE_INFINITY) {
tmp = 2.0 / ((Math.tan(k_m) * (2.0 + t_1)) * (Math.sin(k_m) * ((Math.pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 / Math.pow(((t / Math.pow(Math.cbrt(l), 2.0)) * Math.cbrt((2.0 * Math.pow(k_m, 2.0)))), 3.0);
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = Float64(k_m / t) ^ 2.0 tmp = 0.0 if (Float64(Float64(tan(k_m) * Float64(sin(k_m) * Float64((t ^ 3.0) / Float64(l * l)))) * Float64(1.0 + Float64(1.0 + t_1))) <= Inf) tmp = Float64(2.0 / Float64(Float64(tan(k_m) * Float64(2.0 + t_1)) * Float64(sin(k_m) * Float64(Float64((t ^ 2.0) / l) * Float64(t / l))))); else tmp = Float64(2.0 / (Float64(Float64(t / (cbrt(l) ^ 2.0)) * cbrt(Float64(2.0 * (k_m ^ 2.0)))) ^ 3.0)); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(1.0 + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(2.0 / N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(2.0 + t$95$1), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[Power[N[(N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[Power[N[(2.0 * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := {\left(\frac{k\_m}{t}\right)}^{2}\\
\mathbf{if}\;\left(\tan k\_m \cdot \left(\sin k\_m \cdot \frac{{t}^{3}}{\ell \cdot \ell}\right)\right) \cdot \left(1 + \left(1 + t\_1\right)\right) \leq \infty:\\
\;\;\;\;\frac{2}{\left(\tan k\_m \cdot \left(2 + t\_1\right)\right) \cdot \left(\sin k\_m \cdot \left(\frac{{t}^{2}}{\ell} \cdot \frac{t}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}} \cdot \sqrt[3]{2 \cdot {k\_m}^{2}}\right)}^{3}}\\
\end{array}
\end{array}
if (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t #s(literal 3 binary64)) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 k t) #s(literal 2 binary64))) #s(literal 1 binary64))) < +inf.0Initial program 82.5%
Simplified82.5%
unpow382.4%
times-frac90.8%
pow290.8%
Applied egg-rr90.8%
associate-+r+90.8%
metadata-eval90.8%
+-commutative90.8%
distribute-lft-in90.8%
Applied egg-rr90.8%
+-commutative90.8%
distribute-lft-in90.8%
Simplified90.8%
if +inf.0 < (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t #s(literal 3 binary64)) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 #s(literal 1 binary64) (pow.f64 (/.f64 k t) #s(literal 2 binary64))) #s(literal 1 binary64))) Initial program 0.0%
Simplified14.1%
Taylor expanded in k around 0 25.9%
add-cube-cbrt25.9%
pow325.9%
*-commutative25.9%
cbrt-prod25.9%
associate-/l/11.7%
cbrt-div11.8%
unpow311.8%
add-cbrt-cube27.1%
cbrt-prod48.5%
pow248.5%
Applied egg-rr48.5%
Final simplification76.4%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= t 1.2e-145)
(*
(/ 2.0 (* t (pow k_m 2.0)))
(/ (* (pow l 2.0) (cos k_m)) (pow (sin k_m) 2.0)))
(/
2.0
(*
(pow (/ (* t (cbrt (sin k_m))) (pow (cbrt l) 2.0)) 3.0)
(* (tan k_m) (+ 1.0 (+ 1.0 (pow (/ k_m t) 2.0))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 1.2e-145) {
tmp = (2.0 / (t * pow(k_m, 2.0))) * ((pow(l, 2.0) * cos(k_m)) / pow(sin(k_m), 2.0));
} else {
tmp = 2.0 / (pow(((t * cbrt(sin(k_m))) / pow(cbrt(l), 2.0)), 3.0) * (tan(k_m) * (1.0 + (1.0 + pow((k_m / t), 2.0)))));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (t <= 1.2e-145) {
tmp = (2.0 / (t * Math.pow(k_m, 2.0))) * ((Math.pow(l, 2.0) * Math.cos(k_m)) / Math.pow(Math.sin(k_m), 2.0));
} else {
tmp = 2.0 / (Math.pow(((t * Math.cbrt(Math.sin(k_m))) / Math.pow(Math.cbrt(l), 2.0)), 3.0) * (Math.tan(k_m) * (1.0 + (1.0 + Math.pow((k_m / t), 2.0)))));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 1.2e-145) tmp = Float64(Float64(2.0 / Float64(t * (k_m ^ 2.0))) * Float64(Float64((l ^ 2.0) * cos(k_m)) / (sin(k_m) ^ 2.0))); else tmp = Float64(2.0 / Float64((Float64(Float64(t * cbrt(sin(k_m))) / (cbrt(l) ^ 2.0)) ^ 3.0) * Float64(tan(k_m) * Float64(1.0 + Float64(1.0 + (Float64(k_m / t) ^ 2.0)))))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 1.2e-145], N[(N[(2.0 / N[(t * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Power[N[(N[(t * N[Power[N[Sin[k$95$m], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(N[Tan[k$95$m], $MachinePrecision] * N[(1.0 + N[(1.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 1.2 \cdot 10^{-145}:\\
\;\;\;\;\frac{2}{t \cdot {k\_m}^{2}} \cdot \frac{{\ell}^{2} \cdot \cos k\_m}{{\sin k\_m}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{t \cdot \sqrt[3]{\sin k\_m}}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3} \cdot \left(\tan k\_m \cdot \left(1 + \left(1 + {\left(\frac{k\_m}{t}\right)}^{2}\right)\right)\right)}\\
\end{array}
\end{array}
if t < 1.20000000000000008e-145Initial program 51.1%
associate-/r*50.5%
+-commutative50.5%
unpow250.5%
sqr-neg50.5%
distribute-frac-neg250.5%
distribute-frac-neg250.5%
unpow250.5%
+-commutative50.5%
associate-*l*49.2%
associate-*l/49.2%
associate-/r/49.2%
+-commutative49.2%
associate-+r+49.2%
Simplified49.2%
add-cube-cbrt49.2%
pow349.2%
cbrt-prod49.2%
rem-cbrt-cube51.6%
Applied egg-rr51.6%
*-un-lft-identity51.6%
unpow-prod-down49.2%
pow349.2%
add-cube-cbrt49.2%
Applied egg-rr49.2%
*-lft-identity49.2%
*-commutative49.2%
associate-*l*50.4%
Simplified50.4%
Taylor expanded in k around inf 63.3%
associate-*r/63.3%
associate-*r*63.3%
times-frac63.8%
Simplified63.8%
if 1.20000000000000008e-145 < t Initial program 60.6%
Simplified60.6%
associate-/r*67.0%
add-cube-cbrt66.8%
*-un-lft-identity66.8%
times-frac66.8%
pow266.8%
cbrt-div66.9%
rem-cbrt-cube66.9%
cbrt-div66.8%
rem-cbrt-cube83.3%
Applied egg-rr83.3%
add-cube-cbrt83.3%
pow383.2%
Applied egg-rr87.2%
*-commutative87.2%
Simplified87.2%
associate-*l/87.2%
Applied egg-rr87.2%
Final simplification72.0%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= t 6.2e-141)
(*
(/ 2.0 (* t (pow k_m 2.0)))
(/ (* (pow l 2.0) (cos k_m)) (pow (sin k_m) 2.0)))
(/
2.0
(*
(* (tan k_m) (+ 1.0 (+ 1.0 (pow (/ k_m t) 2.0))))
(pow (* (cbrt (sin k_m)) (* t (pow (cbrt l) -2.0))) 3.0)))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 6.2e-141) {
tmp = (2.0 / (t * pow(k_m, 2.0))) * ((pow(l, 2.0) * cos(k_m)) / pow(sin(k_m), 2.0));
} else {
tmp = 2.0 / ((tan(k_m) * (1.0 + (1.0 + pow((k_m / t), 2.0)))) * pow((cbrt(sin(k_m)) * (t * pow(cbrt(l), -2.0))), 3.0));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (t <= 6.2e-141) {
tmp = (2.0 / (t * Math.pow(k_m, 2.0))) * ((Math.pow(l, 2.0) * Math.cos(k_m)) / Math.pow(Math.sin(k_m), 2.0));
} else {
tmp = 2.0 / ((Math.tan(k_m) * (1.0 + (1.0 + Math.pow((k_m / t), 2.0)))) * Math.pow((Math.cbrt(Math.sin(k_m)) * (t * Math.pow(Math.cbrt(l), -2.0))), 3.0));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 6.2e-141) tmp = Float64(Float64(2.0 / Float64(t * (k_m ^ 2.0))) * Float64(Float64((l ^ 2.0) * cos(k_m)) / (sin(k_m) ^ 2.0))); else tmp = Float64(2.0 / Float64(Float64(tan(k_m) * Float64(1.0 + Float64(1.0 + (Float64(k_m / t) ^ 2.0)))) * (Float64(cbrt(sin(k_m)) * Float64(t * (cbrt(l) ^ -2.0))) ^ 3.0))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 6.2e-141], N[(N[(2.0 / N[(t * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(1.0 + N[(1.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Power[N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 1/3], $MachinePrecision] * N[(t * N[Power[N[Power[l, 1/3], $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 6.2 \cdot 10^{-141}:\\
\;\;\;\;\frac{2}{t \cdot {k\_m}^{2}} \cdot \frac{{\ell}^{2} \cdot \cos k\_m}{{\sin k\_m}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\tan k\_m \cdot \left(1 + \left(1 + {\left(\frac{k\_m}{t}\right)}^{2}\right)\right)\right) \cdot {\left(\sqrt[3]{\sin k\_m} \cdot \left(t \cdot {\left(\sqrt[3]{\ell}\right)}^{-2}\right)\right)}^{3}}\\
\end{array}
\end{array}
if t < 6.20000000000000055e-141Initial program 51.4%
associate-/r*50.8%
+-commutative50.8%
unpow250.8%
sqr-neg50.8%
distribute-frac-neg250.8%
distribute-frac-neg250.8%
unpow250.8%
+-commutative50.8%
associate-*l*49.5%
associate-*l/49.5%
associate-/r/49.5%
+-commutative49.5%
associate-+r+49.5%
Simplified49.5%
add-cube-cbrt49.5%
pow349.5%
cbrt-prod49.5%
rem-cbrt-cube51.9%
Applied egg-rr51.9%
*-un-lft-identity51.9%
unpow-prod-down49.5%
pow349.5%
add-cube-cbrt49.5%
Applied egg-rr49.5%
*-lft-identity49.5%
*-commutative49.5%
associate-*l*50.7%
Simplified50.7%
Taylor expanded in k around inf 63.6%
associate-*r/63.6%
associate-*r*63.6%
times-frac64.0%
Simplified64.0%
if 6.20000000000000055e-141 < t Initial program 60.2%
Simplified60.2%
unpow360.1%
times-frac78.2%
pow278.2%
Applied egg-rr78.2%
add-cube-cbrt77.9%
pow377.9%
frac-times60.0%
unpow260.0%
unpow360.0%
cbrt-div60.1%
unpow360.0%
add-cbrt-cube70.1%
cbrt-unprod84.0%
unpow284.0%
add-cube-cbrt83.9%
pow383.8%
Applied egg-rr87.0%
Final simplification72.0%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= t 9.2e-145)
(*
(/ 2.0 (* t (pow k_m 2.0)))
(/ (* (pow l 2.0) (cos k_m)) (pow (sin k_m) 2.0)))
(if (<= t 1.18e+154)
(/
2.0
(*
(* (tan k_m) (+ 2.0 (pow (/ k_m t) 2.0)))
(* (sin k_m) (* (/ (pow t 2.0) l) (/ t l)))))
(/
2.0
(*
(* 2.0 k_m)
(pow (* (cbrt (sin k_m)) (/ t (pow (cbrt l) 2.0))) 3.0))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 9.2e-145) {
tmp = (2.0 / (t * pow(k_m, 2.0))) * ((pow(l, 2.0) * cos(k_m)) / pow(sin(k_m), 2.0));
} else if (t <= 1.18e+154) {
tmp = 2.0 / ((tan(k_m) * (2.0 + pow((k_m / t), 2.0))) * (sin(k_m) * ((pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 / ((2.0 * k_m) * pow((cbrt(sin(k_m)) * (t / pow(cbrt(l), 2.0))), 3.0));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (t <= 9.2e-145) {
tmp = (2.0 / (t * Math.pow(k_m, 2.0))) * ((Math.pow(l, 2.0) * Math.cos(k_m)) / Math.pow(Math.sin(k_m), 2.0));
} else if (t <= 1.18e+154) {
tmp = 2.0 / ((Math.tan(k_m) * (2.0 + Math.pow((k_m / t), 2.0))) * (Math.sin(k_m) * ((Math.pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 / ((2.0 * k_m) * Math.pow((Math.cbrt(Math.sin(k_m)) * (t / Math.pow(Math.cbrt(l), 2.0))), 3.0));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 9.2e-145) tmp = Float64(Float64(2.0 / Float64(t * (k_m ^ 2.0))) * Float64(Float64((l ^ 2.0) * cos(k_m)) / (sin(k_m) ^ 2.0))); elseif (t <= 1.18e+154) tmp = Float64(2.0 / Float64(Float64(tan(k_m) * Float64(2.0 + (Float64(k_m / t) ^ 2.0))) * Float64(sin(k_m) * Float64(Float64((t ^ 2.0) / l) * Float64(t / l))))); else tmp = Float64(2.0 / Float64(Float64(2.0 * k_m) * (Float64(cbrt(sin(k_m)) * Float64(t / (cbrt(l) ^ 2.0))) ^ 3.0))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 9.2e-145], N[(N[(2.0 / N[(t * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.18e+154], N[(2.0 / N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(2.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(2.0 * k$95$m), $MachinePrecision] * N[Power[N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 1/3], $MachinePrecision] * N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 9.2 \cdot 10^{-145}:\\
\;\;\;\;\frac{2}{t \cdot {k\_m}^{2}} \cdot \frac{{\ell}^{2} \cdot \cos k\_m}{{\sin k\_m}^{2}}\\
\mathbf{elif}\;t \leq 1.18 \cdot 10^{+154}:\\
\;\;\;\;\frac{2}{\left(\tan k\_m \cdot \left(2 + {\left(\frac{k\_m}{t}\right)}^{2}\right)\right) \cdot \left(\sin k\_m \cdot \left(\frac{{t}^{2}}{\ell} \cdot \frac{t}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(2 \cdot k\_m\right) \cdot {\left(\sqrt[3]{\sin k\_m} \cdot \frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}}\\
\end{array}
\end{array}
if t < 9.20000000000000028e-145Initial program 51.1%
associate-/r*50.5%
+-commutative50.5%
unpow250.5%
sqr-neg50.5%
distribute-frac-neg250.5%
distribute-frac-neg250.5%
unpow250.5%
+-commutative50.5%
associate-*l*49.2%
associate-*l/49.2%
associate-/r/49.2%
+-commutative49.2%
associate-+r+49.2%
Simplified49.2%
add-cube-cbrt49.2%
pow349.2%
cbrt-prod49.2%
rem-cbrt-cube51.6%
Applied egg-rr51.6%
*-un-lft-identity51.6%
unpow-prod-down49.2%
pow349.2%
add-cube-cbrt49.2%
Applied egg-rr49.2%
*-lft-identity49.2%
*-commutative49.2%
associate-*l*50.4%
Simplified50.4%
Taylor expanded in k around inf 63.3%
associate-*r/63.3%
associate-*r*63.3%
times-frac63.8%
Simplified63.8%
if 9.20000000000000028e-145 < t < 1.18000000000000004e154Initial program 65.1%
Simplified65.1%
unpow365.0%
times-frac86.7%
pow286.7%
Applied egg-rr86.7%
associate-+r+86.7%
metadata-eval86.7%
+-commutative86.7%
distribute-lft-in86.7%
Applied egg-rr86.7%
+-commutative86.7%
distribute-lft-in86.7%
Simplified86.7%
if 1.18000000000000004e154 < t Initial program 50.7%
Simplified50.7%
associate-/r*60.0%
add-cube-cbrt60.0%
*-un-lft-identity60.0%
times-frac60.0%
pow260.0%
cbrt-div60.0%
rem-cbrt-cube60.0%
cbrt-div60.0%
rem-cbrt-cube76.8%
Applied egg-rr76.8%
add-cube-cbrt76.9%
pow376.9%
Applied egg-rr90.3%
*-commutative90.3%
Simplified90.3%
Taylor expanded in k around 0 83.7%
Final simplification71.5%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= t 8.5e-146)
(*
2.0
(* (/ (pow l 2.0) (* t (pow k_m 2.0))) (/ (cos k_m) (pow (sin k_m) 2.0))))
(if (<= t 1.25e+154)
(/
2.0
(*
(* (tan k_m) (+ 2.0 (pow (/ k_m t) 2.0)))
(* (sin k_m) (* (/ (pow t 2.0) l) (/ t l)))))
(/
2.0
(*
(* 2.0 k_m)
(pow (* (cbrt (sin k_m)) (/ t (pow (cbrt l) 2.0))) 3.0))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 8.5e-146) {
tmp = 2.0 * ((pow(l, 2.0) / (t * pow(k_m, 2.0))) * (cos(k_m) / pow(sin(k_m), 2.0)));
} else if (t <= 1.25e+154) {
tmp = 2.0 / ((tan(k_m) * (2.0 + pow((k_m / t), 2.0))) * (sin(k_m) * ((pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 / ((2.0 * k_m) * pow((cbrt(sin(k_m)) * (t / pow(cbrt(l), 2.0))), 3.0));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (t <= 8.5e-146) {
tmp = 2.0 * ((Math.pow(l, 2.0) / (t * Math.pow(k_m, 2.0))) * (Math.cos(k_m) / Math.pow(Math.sin(k_m), 2.0)));
} else if (t <= 1.25e+154) {
tmp = 2.0 / ((Math.tan(k_m) * (2.0 + Math.pow((k_m / t), 2.0))) * (Math.sin(k_m) * ((Math.pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 / ((2.0 * k_m) * Math.pow((Math.cbrt(Math.sin(k_m)) * (t / Math.pow(Math.cbrt(l), 2.0))), 3.0));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 8.5e-146) tmp = Float64(2.0 * Float64(Float64((l ^ 2.0) / Float64(t * (k_m ^ 2.0))) * Float64(cos(k_m) / (sin(k_m) ^ 2.0)))); elseif (t <= 1.25e+154) tmp = Float64(2.0 / Float64(Float64(tan(k_m) * Float64(2.0 + (Float64(k_m / t) ^ 2.0))) * Float64(sin(k_m) * Float64(Float64((t ^ 2.0) / l) * Float64(t / l))))); else tmp = Float64(2.0 / Float64(Float64(2.0 * k_m) * (Float64(cbrt(sin(k_m)) * Float64(t / (cbrt(l) ^ 2.0))) ^ 3.0))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 8.5e-146], N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / N[(t * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k$95$m], $MachinePrecision] / N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.25e+154], N[(2.0 / N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(2.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(2.0 * k$95$m), $MachinePrecision] * N[Power[N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 1/3], $MachinePrecision] * N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 8.5 \cdot 10^{-146}:\\
\;\;\;\;2 \cdot \left(\frac{{\ell}^{2}}{t \cdot {k\_m}^{2}} \cdot \frac{\cos k\_m}{{\sin k\_m}^{2}}\right)\\
\mathbf{elif}\;t \leq 1.25 \cdot 10^{+154}:\\
\;\;\;\;\frac{2}{\left(\tan k\_m \cdot \left(2 + {\left(\frac{k\_m}{t}\right)}^{2}\right)\right) \cdot \left(\sin k\_m \cdot \left(\frac{{t}^{2}}{\ell} \cdot \frac{t}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(2 \cdot k\_m\right) \cdot {\left(\sqrt[3]{\sin k\_m} \cdot \frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}}\\
\end{array}
\end{array}
if t < 8.4999999999999997e-146Initial program 51.1%
associate-/r*50.5%
+-commutative50.5%
unpow250.5%
sqr-neg50.5%
distribute-frac-neg250.5%
distribute-frac-neg250.5%
unpow250.5%
+-commutative50.5%
associate-*l*49.2%
associate-*l/49.2%
associate-/r/49.2%
+-commutative49.2%
associate-+r+49.2%
Simplified49.2%
Taylor expanded in t around 0 63.3%
associate-*r*63.3%
times-frac64.4%
Simplified64.4%
if 8.4999999999999997e-146 < t < 1.25000000000000001e154Initial program 65.1%
Simplified65.1%
unpow365.0%
times-frac86.7%
pow286.7%
Applied egg-rr86.7%
associate-+r+86.7%
metadata-eval86.7%
+-commutative86.7%
distribute-lft-in86.7%
Applied egg-rr86.7%
+-commutative86.7%
distribute-lft-in86.7%
Simplified86.7%
if 1.25000000000000001e154 < t Initial program 50.7%
Simplified50.7%
associate-/r*60.0%
add-cube-cbrt60.0%
*-un-lft-identity60.0%
times-frac60.0%
pow260.0%
cbrt-div60.0%
rem-cbrt-cube60.0%
cbrt-div60.0%
rem-cbrt-cube76.8%
Applied egg-rr76.8%
add-cube-cbrt76.9%
pow376.9%
Applied egg-rr90.3%
*-commutative90.3%
Simplified90.3%
Taylor expanded in k around 0 83.7%
Final simplification71.9%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(let* ((t_1 (/ t (pow (cbrt l) 2.0))))
(if (<= t 2.7e-123)
(/ 2.0 (pow (* t_1 (cbrt (* 2.0 (pow k_m 2.0)))) 3.0))
(if (<= t 1.32e+154)
(/
2.0
(*
(* (tan k_m) (+ 2.0 (pow (/ k_m t) 2.0)))
(* (sin k_m) (* (/ (pow t 2.0) l) (/ t l)))))
(/ 2.0 (* (* 2.0 k_m) (pow (* (cbrt (sin k_m)) t_1) 3.0)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double t_1 = t / pow(cbrt(l), 2.0);
double tmp;
if (t <= 2.7e-123) {
tmp = 2.0 / pow((t_1 * cbrt((2.0 * pow(k_m, 2.0)))), 3.0);
} else if (t <= 1.32e+154) {
tmp = 2.0 / ((tan(k_m) * (2.0 + pow((k_m / t), 2.0))) * (sin(k_m) * ((pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 / ((2.0 * k_m) * pow((cbrt(sin(k_m)) * t_1), 3.0));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double t_1 = t / Math.pow(Math.cbrt(l), 2.0);
double tmp;
if (t <= 2.7e-123) {
tmp = 2.0 / Math.pow((t_1 * Math.cbrt((2.0 * Math.pow(k_m, 2.0)))), 3.0);
} else if (t <= 1.32e+154) {
tmp = 2.0 / ((Math.tan(k_m) * (2.0 + Math.pow((k_m / t), 2.0))) * (Math.sin(k_m) * ((Math.pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 / ((2.0 * k_m) * Math.pow((Math.cbrt(Math.sin(k_m)) * t_1), 3.0));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) t_1 = Float64(t / (cbrt(l) ^ 2.0)) tmp = 0.0 if (t <= 2.7e-123) tmp = Float64(2.0 / (Float64(t_1 * cbrt(Float64(2.0 * (k_m ^ 2.0)))) ^ 3.0)); elseif (t <= 1.32e+154) tmp = Float64(2.0 / Float64(Float64(tan(k_m) * Float64(2.0 + (Float64(k_m / t) ^ 2.0))) * Float64(sin(k_m) * Float64(Float64((t ^ 2.0) / l) * Float64(t / l))))); else tmp = Float64(2.0 / Float64(Float64(2.0 * k_m) * (Float64(cbrt(sin(k_m)) * t_1) ^ 3.0))); end return tmp end
k_m = N[Abs[k], $MachinePrecision]
code[t_, l_, k$95$m_] := Block[{t$95$1 = N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, 2.7e-123], N[(2.0 / N[Power[N[(t$95$1 * N[Power[N[(2.0 * N[Power[k$95$m, 2.0], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.32e+154], N[(2.0 / N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(2.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(2.0 * k$95$m), $MachinePrecision] * N[Power[N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 1/3], $MachinePrecision] * t$95$1), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
t_1 := \frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\\
\mathbf{if}\;t \leq 2.7 \cdot 10^{-123}:\\
\;\;\;\;\frac{2}{{\left(t\_1 \cdot \sqrt[3]{2 \cdot {k\_m}^{2}}\right)}^{3}}\\
\mathbf{elif}\;t \leq 1.32 \cdot 10^{+154}:\\
\;\;\;\;\frac{2}{\left(\tan k\_m \cdot \left(2 + {\left(\frac{k\_m}{t}\right)}^{2}\right)\right) \cdot \left(\sin k\_m \cdot \left(\frac{{t}^{2}}{\ell} \cdot \frac{t}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(2 \cdot k\_m\right) \cdot {\left(\sqrt[3]{\sin k\_m} \cdot t\_1\right)}^{3}}\\
\end{array}
\end{array}
if t < 2.7000000000000001e-123Initial program 51.1%
Simplified57.6%
Taylor expanded in k around 0 60.3%
add-cube-cbrt60.3%
pow360.3%
*-commutative60.3%
cbrt-prod60.3%
associate-/l/54.1%
cbrt-div54.1%
unpow354.1%
add-cbrt-cube61.7%
cbrt-prod70.7%
pow270.7%
Applied egg-rr70.7%
if 2.7000000000000001e-123 < t < 1.31999999999999998e154Initial program 66.1%
Simplified66.2%
unpow366.0%
times-frac87.5%
pow287.5%
Applied egg-rr87.5%
associate-+r+87.6%
metadata-eval87.6%
+-commutative87.6%
distribute-lft-in87.6%
Applied egg-rr87.6%
+-commutative87.6%
distribute-lft-in87.6%
Simplified87.6%
if 1.31999999999999998e154 < t Initial program 50.7%
Simplified50.7%
associate-/r*60.0%
add-cube-cbrt60.0%
*-un-lft-identity60.0%
times-frac60.0%
pow260.0%
cbrt-div60.0%
rem-cbrt-cube60.0%
cbrt-div60.0%
rem-cbrt-cube76.8%
Applied egg-rr76.8%
add-cube-cbrt76.9%
pow376.9%
Applied egg-rr90.3%
*-commutative90.3%
Simplified90.3%
Taylor expanded in k around 0 83.7%
Final simplification75.9%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= t 4.2e-243)
(/
2.0
(* (/ (/ (pow t 3.0) l) l) (* 2.0 (/ (pow (sin k_m) 2.0) (cos k_m)))))
(if (<= t 1.22e-160)
(/ 2.0 (/ (* (/ (pow k_m 3.0) (pow l 2.0)) (* t (sin k_m))) (cos k_m)))
(if (<= t 1.32e+154)
(/
2.0
(*
(* (tan k_m) (+ 2.0 (pow (/ k_m t) 2.0)))
(* (sin k_m) (* (/ (pow t 2.0) l) (/ t l)))))
(/
2.0
(* (* 2.0 k_m) (* (sin k_m) (pow (/ t (pow (cbrt l) 2.0)) 3.0))))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 4.2e-243) {
tmp = 2.0 / (((pow(t, 3.0) / l) / l) * (2.0 * (pow(sin(k_m), 2.0) / cos(k_m))));
} else if (t <= 1.22e-160) {
tmp = 2.0 / (((pow(k_m, 3.0) / pow(l, 2.0)) * (t * sin(k_m))) / cos(k_m));
} else if (t <= 1.32e+154) {
tmp = 2.0 / ((tan(k_m) * (2.0 + pow((k_m / t), 2.0))) * (sin(k_m) * ((pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 / ((2.0 * k_m) * (sin(k_m) * pow((t / pow(cbrt(l), 2.0)), 3.0)));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (t <= 4.2e-243) {
tmp = 2.0 / (((Math.pow(t, 3.0) / l) / l) * (2.0 * (Math.pow(Math.sin(k_m), 2.0) / Math.cos(k_m))));
} else if (t <= 1.22e-160) {
tmp = 2.0 / (((Math.pow(k_m, 3.0) / Math.pow(l, 2.0)) * (t * Math.sin(k_m))) / Math.cos(k_m));
} else if (t <= 1.32e+154) {
tmp = 2.0 / ((Math.tan(k_m) * (2.0 + Math.pow((k_m / t), 2.0))) * (Math.sin(k_m) * ((Math.pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 / ((2.0 * k_m) * (Math.sin(k_m) * Math.pow((t / Math.pow(Math.cbrt(l), 2.0)), 3.0)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 4.2e-243) tmp = Float64(2.0 / Float64(Float64(Float64((t ^ 3.0) / l) / l) * Float64(2.0 * Float64((sin(k_m) ^ 2.0) / cos(k_m))))); elseif (t <= 1.22e-160) tmp = Float64(2.0 / Float64(Float64(Float64((k_m ^ 3.0) / (l ^ 2.0)) * Float64(t * sin(k_m))) / cos(k_m))); elseif (t <= 1.32e+154) tmp = Float64(2.0 / Float64(Float64(tan(k_m) * Float64(2.0 + (Float64(k_m / t) ^ 2.0))) * Float64(sin(k_m) * Float64(Float64((t ^ 2.0) / l) * Float64(t / l))))); else tmp = Float64(2.0 / Float64(Float64(2.0 * k_m) * Float64(sin(k_m) * (Float64(t / (cbrt(l) ^ 2.0)) ^ 3.0)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 4.2e-243], N[(2.0 / N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision] * N[(2.0 * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.22e-160], N[(2.0 / N[(N[(N[(N[Power[k$95$m, 3.0], $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.32e+154], N[(2.0 / N[(N[(N[Tan[k$95$m], $MachinePrecision] * N[(2.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(2.0 * k$95$m), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Power[N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 4.2 \cdot 10^{-243}:\\
\;\;\;\;\frac{2}{\frac{\frac{{t}^{3}}{\ell}}{\ell} \cdot \left(2 \cdot \frac{{\sin k\_m}^{2}}{\cos k\_m}\right)}\\
\mathbf{elif}\;t \leq 1.22 \cdot 10^{-160}:\\
\;\;\;\;\frac{2}{\frac{\frac{{k\_m}^{3}}{{\ell}^{2}} \cdot \left(t \cdot \sin k\_m\right)}{\cos k\_m}}\\
\mathbf{elif}\;t \leq 1.32 \cdot 10^{+154}:\\
\;\;\;\;\frac{2}{\left(\tan k\_m \cdot \left(2 + {\left(\frac{k\_m}{t}\right)}^{2}\right)\right) \cdot \left(\sin k\_m \cdot \left(\frac{{t}^{2}}{\ell} \cdot \frac{t}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(2 \cdot k\_m\right) \cdot \left(\sin k\_m \cdot {\left(\frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}\right)}\\
\end{array}
\end{array}
if t < 4.2000000000000002e-243Initial program 53.3%
Simplified60.1%
Taylor expanded in t around inf 62.1%
if 4.2000000000000002e-243 < t < 1.22000000000000003e-160Initial program 29.4%
Simplified29.4%
Taylor expanded in k around 0 29.4%
Taylor expanded in k around inf 60.8%
times-frac65.6%
associate-*r/65.6%
Simplified65.6%
if 1.22000000000000003e-160 < t < 1.31999999999999998e154Initial program 65.7%
Simplified65.7%
unpow365.6%
times-frac86.9%
pow286.9%
Applied egg-rr86.9%
associate-+r+87.0%
metadata-eval87.0%
+-commutative87.0%
distribute-lft-in87.0%
Applied egg-rr87.0%
+-commutative87.0%
distribute-lft-in87.0%
Simplified87.0%
if 1.31999999999999998e154 < t Initial program 50.7%
Simplified50.7%
associate-/r*60.0%
add-cube-cbrt60.0%
*-un-lft-identity60.0%
times-frac60.0%
pow260.0%
cbrt-div60.0%
rem-cbrt-cube60.0%
cbrt-div60.0%
rem-cbrt-cube76.8%
Applied egg-rr76.8%
add-cube-cbrt76.8%
pow376.7%
frac-times70.2%
unpow270.2%
*-un-lft-identity70.2%
cbrt-div70.2%
add-cbrt-cube80.1%
add-cbrt-cube60.0%
unpow360.0%
cbrt-div60.0%
cbrt-div60.0%
Applied egg-rr80.2%
Taylor expanded in k around 0 73.7%
Final simplification69.7%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= t 9.6e-92)
(/
2.0
(*
(* (sin k_m) (* (/ (pow t 2.0) l) (/ t l)))
(* 2.0 (/ (sin k_m) (cos k_m)))))
(if (<= t 1.9e+92)
(*
(* l (/ 2.0 (* (pow t 3.0) (* (sin k_m) (tan k_m)))))
(/ l (+ 2.0 (pow (/ k_m t) 2.0))))
(/
2.0
(* (* 2.0 k_m) (* (sin k_m) (pow (/ t (pow (cbrt l) 2.0)) 3.0)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (t <= 9.6e-92) {
tmp = 2.0 / ((sin(k_m) * ((pow(t, 2.0) / l) * (t / l))) * (2.0 * (sin(k_m) / cos(k_m))));
} else if (t <= 1.9e+92) {
tmp = (l * (2.0 / (pow(t, 3.0) * (sin(k_m) * tan(k_m))))) * (l / (2.0 + pow((k_m / t), 2.0)));
} else {
tmp = 2.0 / ((2.0 * k_m) * (sin(k_m) * pow((t / pow(cbrt(l), 2.0)), 3.0)));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (t <= 9.6e-92) {
tmp = 2.0 / ((Math.sin(k_m) * ((Math.pow(t, 2.0) / l) * (t / l))) * (2.0 * (Math.sin(k_m) / Math.cos(k_m))));
} else if (t <= 1.9e+92) {
tmp = (l * (2.0 / (Math.pow(t, 3.0) * (Math.sin(k_m) * Math.tan(k_m))))) * (l / (2.0 + Math.pow((k_m / t), 2.0)));
} else {
tmp = 2.0 / ((2.0 * k_m) * (Math.sin(k_m) * Math.pow((t / Math.pow(Math.cbrt(l), 2.0)), 3.0)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (t <= 9.6e-92) tmp = Float64(2.0 / Float64(Float64(sin(k_m) * Float64(Float64((t ^ 2.0) / l) * Float64(t / l))) * Float64(2.0 * Float64(sin(k_m) / cos(k_m))))); elseif (t <= 1.9e+92) tmp = Float64(Float64(l * Float64(2.0 / Float64((t ^ 3.0) * Float64(sin(k_m) * tan(k_m))))) * Float64(l / Float64(2.0 + (Float64(k_m / t) ^ 2.0)))); else tmp = Float64(2.0 / Float64(Float64(2.0 * k_m) * Float64(sin(k_m) * (Float64(t / (cbrt(l) ^ 2.0)) ^ 3.0)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[t, 9.6e-92], N[(2.0 / N[(N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(2.0 * N[(N[Sin[k$95$m], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.9e+92], N[(N[(l * N[(2.0 / N[(N[Power[t, 3.0], $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Tan[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / N[(2.0 + N[Power[N[(k$95$m / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(2.0 * k$95$m), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Power[N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;t \leq 9.6 \cdot 10^{-92}:\\
\;\;\;\;\frac{2}{\left(\sin k\_m \cdot \left(\frac{{t}^{2}}{\ell} \cdot \frac{t}{\ell}\right)\right) \cdot \left(2 \cdot \frac{\sin k\_m}{\cos k\_m}\right)}\\
\mathbf{elif}\;t \leq 1.9 \cdot 10^{+92}:\\
\;\;\;\;\left(\ell \cdot \frac{2}{{t}^{3} \cdot \left(\sin k\_m \cdot \tan k\_m\right)}\right) \cdot \frac{\ell}{2 + {\left(\frac{k\_m}{t}\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(2 \cdot k\_m\right) \cdot \left(\sin k\_m \cdot {\left(\frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}\right)}\\
\end{array}
\end{array}
if t < 9.6000000000000004e-92Initial program 49.9%
Simplified49.9%
unpow349.9%
times-frac62.1%
pow262.1%
Applied egg-rr62.1%
Taylor expanded in t around inf 64.6%
if 9.6000000000000004e-92 < t < 1.9e92Initial program 89.7%
associate-/r*89.6%
+-commutative89.6%
unpow289.6%
sqr-neg89.6%
distribute-frac-neg289.6%
distribute-frac-neg289.6%
unpow289.6%
+-commutative89.6%
associate-*l*89.7%
associate-*l/89.9%
associate-/r/89.6%
+-commutative89.6%
associate-+r+89.7%
Simplified89.7%
associate-*r*89.8%
*-un-lft-identity89.8%
times-frac89.8%
Applied egg-rr89.8%
/-rgt-identity89.8%
*-commutative89.8%
Simplified89.8%
if 1.9e92 < t Initial program 49.1%
Simplified49.1%
associate-/r*57.4%
add-cube-cbrt57.3%
*-un-lft-identity57.3%
times-frac57.3%
pow257.3%
cbrt-div57.4%
rem-cbrt-cube57.4%
cbrt-div57.4%
rem-cbrt-cube84.8%
Applied egg-rr84.8%
add-cube-cbrt84.7%
pow384.6%
frac-times76.3%
unpow276.3%
*-un-lft-identity76.3%
cbrt-div76.4%
add-cbrt-cube86.7%
add-cbrt-cube57.3%
unpow357.3%
cbrt-div57.3%
cbrt-div57.3%
Applied egg-rr86.8%
Taylor expanded in k around 0 76.8%
Final simplification69.7%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 7.8e-17)
(/ 2.0 (* (* 2.0 k_m) (* (sin k_m) (pow (/ t (pow (cbrt l) 2.0)) 3.0))))
(if (<= k_m 2.55e+91)
(/
2.0
(*
(* (sin k_m) (* (/ (pow t 2.0) l) (/ t l)))
(* 2.0 (/ (sin k_m) (cos k_m)))))
(/ 2.0 (/ (* (/ (pow k_m 3.0) (pow l 2.0)) (* t (sin k_m))) (cos k_m))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 7.8e-17) {
tmp = 2.0 / ((2.0 * k_m) * (sin(k_m) * pow((t / pow(cbrt(l), 2.0)), 3.0)));
} else if (k_m <= 2.55e+91) {
tmp = 2.0 / ((sin(k_m) * ((pow(t, 2.0) / l) * (t / l))) * (2.0 * (sin(k_m) / cos(k_m))));
} else {
tmp = 2.0 / (((pow(k_m, 3.0) / pow(l, 2.0)) * (t * sin(k_m))) / cos(k_m));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 7.8e-17) {
tmp = 2.0 / ((2.0 * k_m) * (Math.sin(k_m) * Math.pow((t / Math.pow(Math.cbrt(l), 2.0)), 3.0)));
} else if (k_m <= 2.55e+91) {
tmp = 2.0 / ((Math.sin(k_m) * ((Math.pow(t, 2.0) / l) * (t / l))) * (2.0 * (Math.sin(k_m) / Math.cos(k_m))));
} else {
tmp = 2.0 / (((Math.pow(k_m, 3.0) / Math.pow(l, 2.0)) * (t * Math.sin(k_m))) / Math.cos(k_m));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 7.8e-17) tmp = Float64(2.0 / Float64(Float64(2.0 * k_m) * Float64(sin(k_m) * (Float64(t / (cbrt(l) ^ 2.0)) ^ 3.0)))); elseif (k_m <= 2.55e+91) tmp = Float64(2.0 / Float64(Float64(sin(k_m) * Float64(Float64((t ^ 2.0) / l) * Float64(t / l))) * Float64(2.0 * Float64(sin(k_m) / cos(k_m))))); else tmp = Float64(2.0 / Float64(Float64(Float64((k_m ^ 3.0) / (l ^ 2.0)) * Float64(t * sin(k_m))) / cos(k_m))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 7.8e-17], N[(2.0 / N[(N[(2.0 * k$95$m), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Power[N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2.55e+91], N[(2.0 / N[(N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(2.0 * N[(N[Sin[k$95$m], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(N[Power[k$95$m, 3.0], $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 7.8 \cdot 10^{-17}:\\
\;\;\;\;\frac{2}{\left(2 \cdot k\_m\right) \cdot \left(\sin k\_m \cdot {\left(\frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}\right)}\\
\mathbf{elif}\;k\_m \leq 2.55 \cdot 10^{+91}:\\
\;\;\;\;\frac{2}{\left(\sin k\_m \cdot \left(\frac{{t}^{2}}{\ell} \cdot \frac{t}{\ell}\right)\right) \cdot \left(2 \cdot \frac{\sin k\_m}{\cos k\_m}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\frac{{k\_m}^{3}}{{\ell}^{2}} \cdot \left(t \cdot \sin k\_m\right)}{\cos k\_m}}\\
\end{array}
\end{array}
if k < 7.79999999999999979e-17Initial program 51.6%
Simplified51.6%
associate-/r*58.5%
add-cube-cbrt58.4%
*-un-lft-identity58.4%
times-frac58.4%
pow258.4%
cbrt-div58.5%
rem-cbrt-cube58.5%
cbrt-div58.4%
rem-cbrt-cube71.4%
Applied egg-rr71.4%
add-cube-cbrt71.3%
pow371.3%
frac-times67.7%
unpow267.7%
*-un-lft-identity67.7%
cbrt-div67.6%
add-cbrt-cube72.7%
add-cbrt-cube58.4%
unpow358.4%
cbrt-div58.4%
cbrt-div58.4%
Applied egg-rr72.7%
Taylor expanded in k around 0 69.2%
if 7.79999999999999979e-17 < k < 2.55000000000000007e91Initial program 55.3%
Simplified55.3%
unpow355.2%
times-frac73.8%
pow273.8%
Applied egg-rr73.8%
Taylor expanded in t around inf 69.6%
if 2.55000000000000007e91 < k Initial program 67.1%
Simplified67.1%
Taylor expanded in k around 0 62.5%
Taylor expanded in k around inf 69.0%
times-frac69.0%
associate-*r/69.0%
Simplified69.0%
Final simplification69.2%
k_m = (fabs.f64 k)
(FPCore (t l k_m)
:precision binary64
(if (<= k_m 8.4e-13)
(/ 2.0 (* (* 2.0 k_m) (* (sin k_m) (pow (/ t (pow (cbrt l) 2.0)) 3.0))))
(if (<= k_m 2.25e+83)
(/
2.0
(* (/ (/ (pow t 3.0) l) l) (* 2.0 (/ (pow (sin k_m) 2.0) (cos k_m)))))
(/ (* 2.0 (* (pow l 2.0) (cos k_m))) (* (pow k_m 3.0) (* t (sin k_m)))))))k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 8.4e-13) {
tmp = 2.0 / ((2.0 * k_m) * (sin(k_m) * pow((t / pow(cbrt(l), 2.0)), 3.0)));
} else if (k_m <= 2.25e+83) {
tmp = 2.0 / (((pow(t, 3.0) / l) / l) * (2.0 * (pow(sin(k_m), 2.0) / cos(k_m))));
} else {
tmp = (2.0 * (pow(l, 2.0) * cos(k_m))) / (pow(k_m, 3.0) * (t * sin(k_m)));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 8.4e-13) {
tmp = 2.0 / ((2.0 * k_m) * (Math.sin(k_m) * Math.pow((t / Math.pow(Math.cbrt(l), 2.0)), 3.0)));
} else if (k_m <= 2.25e+83) {
tmp = 2.0 / (((Math.pow(t, 3.0) / l) / l) * (2.0 * (Math.pow(Math.sin(k_m), 2.0) / Math.cos(k_m))));
} else {
tmp = (2.0 * (Math.pow(l, 2.0) * Math.cos(k_m))) / (Math.pow(k_m, 3.0) * (t * Math.sin(k_m)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 8.4e-13) tmp = Float64(2.0 / Float64(Float64(2.0 * k_m) * Float64(sin(k_m) * (Float64(t / (cbrt(l) ^ 2.0)) ^ 3.0)))); elseif (k_m <= 2.25e+83) tmp = Float64(2.0 / Float64(Float64(Float64((t ^ 3.0) / l) / l) * Float64(2.0 * Float64((sin(k_m) ^ 2.0) / cos(k_m))))); else tmp = Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k_m))) / Float64((k_m ^ 3.0) * Float64(t * sin(k_m)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 8.4e-13], N[(2.0 / N[(N[(2.0 * k$95$m), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Power[N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k$95$m, 2.25e+83], N[(2.0 / N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision] * N[(2.0 * N[(N[Power[N[Sin[k$95$m], $MachinePrecision], 2.0], $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Power[k$95$m, 3.0], $MachinePrecision] * N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 8.4 \cdot 10^{-13}:\\
\;\;\;\;\frac{2}{\left(2 \cdot k\_m\right) \cdot \left(\sin k\_m \cdot {\left(\frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}\right)}\\
\mathbf{elif}\;k\_m \leq 2.25 \cdot 10^{+83}:\\
\;\;\;\;\frac{2}{\frac{\frac{{t}^{3}}{\ell}}{\ell} \cdot \left(2 \cdot \frac{{\sin k\_m}^{2}}{\cos k\_m}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\_m\right)}{{k\_m}^{3} \cdot \left(t \cdot \sin k\_m\right)}\\
\end{array}
\end{array}
if k < 8.39999999999999955e-13Initial program 51.6%
Simplified51.6%
associate-/r*58.5%
add-cube-cbrt58.4%
*-un-lft-identity58.4%
times-frac58.4%
pow258.4%
cbrt-div58.5%
rem-cbrt-cube58.5%
cbrt-div58.4%
rem-cbrt-cube71.4%
Applied egg-rr71.4%
add-cube-cbrt71.3%
pow371.3%
frac-times67.7%
unpow267.7%
*-un-lft-identity67.7%
cbrt-div67.6%
add-cbrt-cube72.7%
add-cbrt-cube58.4%
unpow358.4%
cbrt-div58.4%
cbrt-div58.4%
Applied egg-rr72.7%
Taylor expanded in k around 0 69.2%
if 8.39999999999999955e-13 < k < 2.25e83Initial program 57.1%
Simplified73.2%
Taylor expanded in t around inf 68.4%
if 2.25e83 < k Initial program 65.4%
Simplified65.4%
Taylor expanded in k around 0 60.9%
Taylor expanded in k around inf 67.3%
associate-*r/67.3%
Simplified67.3%
Final simplification68.8%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 2.1e+76) (/ 2.0 (* (* 2.0 k_m) (* (sin k_m) (pow (/ t (pow (cbrt l) 2.0)) 3.0)))) (/ (* 2.0 (* (pow l 2.0) (cos k_m))) (* (pow k_m 3.0) (* t (sin k_m))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 2.1e+76) {
tmp = 2.0 / ((2.0 * k_m) * (sin(k_m) * pow((t / pow(cbrt(l), 2.0)), 3.0)));
} else {
tmp = (2.0 * (pow(l, 2.0) * cos(k_m))) / (pow(k_m, 3.0) * (t * sin(k_m)));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 2.1e+76) {
tmp = 2.0 / ((2.0 * k_m) * (Math.sin(k_m) * Math.pow((t / Math.pow(Math.cbrt(l), 2.0)), 3.0)));
} else {
tmp = (2.0 * (Math.pow(l, 2.0) * Math.cos(k_m))) / (Math.pow(k_m, 3.0) * (t * Math.sin(k_m)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 2.1e+76) tmp = Float64(2.0 / Float64(Float64(2.0 * k_m) * Float64(sin(k_m) * (Float64(t / (cbrt(l) ^ 2.0)) ^ 3.0)))); else tmp = Float64(Float64(2.0 * Float64((l ^ 2.0) * cos(k_m))) / Float64((k_m ^ 3.0) * Float64(t * sin(k_m)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 2.1e+76], N[(2.0 / N[(N[(2.0 * k$95$m), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Power[N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] * N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Power[k$95$m, 3.0], $MachinePrecision] * N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 2.1 \cdot 10^{+76}:\\
\;\;\;\;\frac{2}{\left(2 \cdot k\_m\right) \cdot \left(\sin k\_m \cdot {\left(\frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \left({\ell}^{2} \cdot \cos k\_m\right)}{{k\_m}^{3} \cdot \left(t \cdot \sin k\_m\right)}\\
\end{array}
\end{array}
if k < 2.10000000000000007e76Initial program 52.6%
Simplified52.7%
associate-/r*60.8%
add-cube-cbrt60.7%
*-un-lft-identity60.7%
times-frac60.7%
pow260.7%
cbrt-div60.8%
rem-cbrt-cube60.8%
cbrt-div60.7%
rem-cbrt-cube72.4%
Applied egg-rr72.4%
add-cube-cbrt72.3%
pow372.3%
frac-times69.1%
unpow269.1%
*-un-lft-identity69.1%
cbrt-div69.1%
add-cbrt-cube73.5%
add-cbrt-cube60.7%
unpow360.6%
cbrt-div60.7%
cbrt-div60.7%
Applied egg-rr73.5%
Taylor expanded in k around 0 67.6%
if 2.10000000000000007e76 < k Initial program 63.9%
Simplified63.9%
Taylor expanded in k around 0 59.5%
Taylor expanded in k around inf 65.6%
associate-*r/65.6%
Simplified65.6%
Final simplification67.3%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 5.6e+69) (/ 2.0 (* (* 2.0 k_m) (* (sin k_m) (pow (/ t (pow (cbrt l) 2.0)) 3.0)))) (/ 2.0 (/ (* (/ (pow k_m 3.0) (pow l 2.0)) (* t (sin k_m))) (cos k_m)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 5.6e+69) {
tmp = 2.0 / ((2.0 * k_m) * (sin(k_m) * pow((t / pow(cbrt(l), 2.0)), 3.0)));
} else {
tmp = 2.0 / (((pow(k_m, 3.0) / pow(l, 2.0)) * (t * sin(k_m))) / cos(k_m));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 5.6e+69) {
tmp = 2.0 / ((2.0 * k_m) * (Math.sin(k_m) * Math.pow((t / Math.pow(Math.cbrt(l), 2.0)), 3.0)));
} else {
tmp = 2.0 / (((Math.pow(k_m, 3.0) / Math.pow(l, 2.0)) * (t * Math.sin(k_m))) / Math.cos(k_m));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 5.6e+69) tmp = Float64(2.0 / Float64(Float64(2.0 * k_m) * Float64(sin(k_m) * (Float64(t / (cbrt(l) ^ 2.0)) ^ 3.0)))); else tmp = Float64(2.0 / Float64(Float64(Float64((k_m ^ 3.0) / (l ^ 2.0)) * Float64(t * sin(k_m))) / cos(k_m))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 5.6e+69], N[(2.0 / N[(N[(2.0 * k$95$m), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Power[N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(N[Power[k$95$m, 3.0], $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * N[(t * N[Sin[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Cos[k$95$m], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 5.6 \cdot 10^{+69}:\\
\;\;\;\;\frac{2}{\left(2 \cdot k\_m\right) \cdot \left(\sin k\_m \cdot {\left(\frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\frac{{k\_m}^{3}}{{\ell}^{2}} \cdot \left(t \cdot \sin k\_m\right)}{\cos k\_m}}\\
\end{array}
\end{array}
if k < 5.59999999999999964e69Initial program 52.6%
Simplified52.7%
associate-/r*60.8%
add-cube-cbrt60.7%
*-un-lft-identity60.7%
times-frac60.7%
pow260.7%
cbrt-div60.8%
rem-cbrt-cube60.8%
cbrt-div60.7%
rem-cbrt-cube72.4%
Applied egg-rr72.4%
add-cube-cbrt72.3%
pow372.3%
frac-times69.1%
unpow269.1%
*-un-lft-identity69.1%
cbrt-div69.1%
add-cbrt-cube73.5%
add-cbrt-cube60.7%
unpow360.6%
cbrt-div60.7%
cbrt-div60.7%
Applied egg-rr73.5%
Taylor expanded in k around 0 67.6%
if 5.59999999999999964e69 < k Initial program 63.9%
Simplified63.9%
Taylor expanded in k around 0 59.5%
Taylor expanded in k around inf 65.6%
times-frac65.6%
associate-*r/65.6%
Simplified65.6%
Final simplification67.3%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 1.3e+92) (/ 2.0 (* (* 2.0 k_m) (* (sin k_m) (pow (/ t (pow (cbrt l) 2.0)) 3.0)))) (* 2.0 (/ (pow l 2.0) (* t (pow k_m 4.0))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.3e+92) {
tmp = 2.0 / ((2.0 * k_m) * (sin(k_m) * pow((t / pow(cbrt(l), 2.0)), 3.0)));
} else {
tmp = 2.0 * (pow(l, 2.0) / (t * pow(k_m, 4.0)));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.3e+92) {
tmp = 2.0 / ((2.0 * k_m) * (Math.sin(k_m) * Math.pow((t / Math.pow(Math.cbrt(l), 2.0)), 3.0)));
} else {
tmp = 2.0 * (Math.pow(l, 2.0) / (t * Math.pow(k_m, 4.0)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.3e+92) tmp = Float64(2.0 / Float64(Float64(2.0 * k_m) * Float64(sin(k_m) * (Float64(t / (cbrt(l) ^ 2.0)) ^ 3.0)))); else tmp = Float64(2.0 * Float64((l ^ 2.0) / Float64(t * (k_m ^ 4.0)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.3e+92], N[(2.0 / N[(N[(2.0 * k$95$m), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[Power[N[(t / N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] / N[(t * N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.3 \cdot 10^{+92}:\\
\;\;\;\;\frac{2}{\left(2 \cdot k\_m\right) \cdot \left(\sin k\_m \cdot {\left(\frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{3}\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{{\ell}^{2}}{t \cdot {k\_m}^{4}}\\
\end{array}
\end{array}
if k < 1.2999999999999999e92Initial program 52.2%
Simplified52.2%
associate-/r*60.3%
add-cube-cbrt60.1%
*-un-lft-identity60.1%
times-frac60.1%
pow260.1%
cbrt-div60.2%
rem-cbrt-cube60.2%
cbrt-div60.2%
rem-cbrt-cube71.7%
Applied egg-rr71.7%
add-cube-cbrt71.6%
pow371.6%
frac-times68.5%
unpow268.5%
*-un-lft-identity68.5%
cbrt-div68.4%
add-cbrt-cube72.8%
add-cbrt-cube60.1%
unpow360.1%
cbrt-div60.1%
cbrt-div60.2%
Applied egg-rr72.8%
Taylor expanded in k around 0 67.0%
if 1.2999999999999999e92 < k Initial program 67.1%
associate-/r*64.6%
+-commutative64.6%
unpow264.6%
sqr-neg64.6%
distribute-frac-neg264.6%
distribute-frac-neg264.6%
unpow264.6%
+-commutative64.6%
associate-*l*64.6%
associate-*l/64.6%
associate-/r/64.6%
+-commutative64.6%
associate-+r+64.6%
Simplified64.6%
Taylor expanded in k around 0 62.5%
*-commutative62.5%
Simplified62.5%
Taylor expanded in t around 0 69.0%
Final simplification67.3%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 1.3e+92) (/ 2.0 (* (* 2.0 k_m) (* (sin k_m) (* (/ (pow t 2.0) l) (/ t l))))) (* 2.0 (/ (pow l 2.0) (* t (pow k_m 4.0))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.3e+92) {
tmp = 2.0 / ((2.0 * k_m) * (sin(k_m) * ((pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 * (pow(l, 2.0) / (t * pow(k_m, 4.0)));
}
return tmp;
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 1.3d+92) then
tmp = 2.0d0 / ((2.0d0 * k_m) * (sin(k_m) * (((t ** 2.0d0) / l) * (t / l))))
else
tmp = 2.0d0 * ((l ** 2.0d0) / (t * (k_m ** 4.0d0)))
end if
code = tmp
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.3e+92) {
tmp = 2.0 / ((2.0 * k_m) * (Math.sin(k_m) * ((Math.pow(t, 2.0) / l) * (t / l))));
} else {
tmp = 2.0 * (Math.pow(l, 2.0) / (t * Math.pow(k_m, 4.0)));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 1.3e+92: tmp = 2.0 / ((2.0 * k_m) * (math.sin(k_m) * ((math.pow(t, 2.0) / l) * (t / l)))) else: tmp = 2.0 * (math.pow(l, 2.0) / (t * math.pow(k_m, 4.0))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.3e+92) tmp = Float64(2.0 / Float64(Float64(2.0 * k_m) * Float64(sin(k_m) * Float64(Float64((t ^ 2.0) / l) * Float64(t / l))))); else tmp = Float64(2.0 * Float64((l ^ 2.0) / Float64(t * (k_m ^ 4.0)))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 1.3e+92) tmp = 2.0 / ((2.0 * k_m) * (sin(k_m) * (((t ^ 2.0) / l) * (t / l)))); else tmp = 2.0 * ((l ^ 2.0) / (t * (k_m ^ 4.0))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.3e+92], N[(2.0 / N[(N[(2.0 * k$95$m), $MachinePrecision] * N[(N[Sin[k$95$m], $MachinePrecision] * N[(N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] / N[(t * N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.3 \cdot 10^{+92}:\\
\;\;\;\;\frac{2}{\left(2 \cdot k\_m\right) \cdot \left(\sin k\_m \cdot \left(\frac{{t}^{2}}{\ell} \cdot \frac{t}{\ell}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{{\ell}^{2}}{t \cdot {k\_m}^{4}}\\
\end{array}
\end{array}
if k < 1.2999999999999999e92Initial program 52.2%
Simplified52.2%
unpow352.1%
times-frac66.1%
pow266.1%
Applied egg-rr66.1%
Taylor expanded in k around 0 64.0%
if 1.2999999999999999e92 < k Initial program 67.1%
associate-/r*64.6%
+-commutative64.6%
unpow264.6%
sqr-neg64.6%
distribute-frac-neg264.6%
distribute-frac-neg264.6%
unpow264.6%
+-commutative64.6%
associate-*l*64.6%
associate-*l/64.6%
associate-/r/64.6%
+-commutative64.6%
associate-+r+64.6%
Simplified64.6%
Taylor expanded in k around 0 62.5%
*-commutative62.5%
Simplified62.5%
Taylor expanded in t around 0 69.0%
Final simplification64.8%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 1.9e+120) (/ 2.0 (* (/ (pow (/ t (cbrt l)) 3.0) l) (* 2.0 (* k_m k_m)))) (* 2.0 (/ (pow l 2.0) (* t (pow k_m 4.0))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.9e+120) {
tmp = 2.0 / ((pow((t / cbrt(l)), 3.0) / l) * (2.0 * (k_m * k_m)));
} else {
tmp = 2.0 * (pow(l, 2.0) / (t * pow(k_m, 4.0)));
}
return tmp;
}
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.9e+120) {
tmp = 2.0 / ((Math.pow((t / Math.cbrt(l)), 3.0) / l) * (2.0 * (k_m * k_m)));
} else {
tmp = 2.0 * (Math.pow(l, 2.0) / (t * Math.pow(k_m, 4.0)));
}
return tmp;
}
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.9e+120) tmp = Float64(2.0 / Float64(Float64((Float64(t / cbrt(l)) ^ 3.0) / l) * Float64(2.0 * Float64(k_m * k_m)))); else tmp = Float64(2.0 * Float64((l ^ 2.0) / Float64(t * (k_m ^ 4.0)))); end return tmp end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.9e+120], N[(2.0 / N[(N[(N[Power[N[(t / N[Power[l, 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] / l), $MachinePrecision] * N[(2.0 * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] / N[(t * N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.9 \cdot 10^{+120}:\\
\;\;\;\;\frac{2}{\frac{{\left(\frac{t}{\sqrt[3]{\ell}}\right)}^{3}}{\ell} \cdot \left(2 \cdot \left(k\_m \cdot k\_m\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{{\ell}^{2}}{t \cdot {k\_m}^{4}}\\
\end{array}
\end{array}
if k < 1.8999999999999999e120Initial program 53.4%
Simplified58.0%
Taylor expanded in k around 0 56.8%
unpow256.8%
Applied egg-rr56.8%
add-cube-cbrt56.7%
pow356.7%
cbrt-div56.7%
rem-cbrt-cube62.0%
Applied egg-rr62.0%
if 1.8999999999999999e120 < k Initial program 61.2%
associate-/r*58.2%
+-commutative58.2%
unpow258.2%
sqr-neg58.2%
distribute-frac-neg258.2%
distribute-frac-neg258.2%
unpow258.2%
+-commutative58.2%
associate-*l*58.2%
associate-*l/58.2%
associate-/r/58.2%
+-commutative58.2%
associate-+r+58.2%
Simplified58.2%
Taylor expanded in k around 0 58.6%
*-commutative58.6%
Simplified58.6%
Taylor expanded in t around 0 66.2%
Final simplification62.6%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 1.65e+120) (/ 2.0 (* (* 2.0 (* k_m k_m)) (/ (* t (/ (pow t 2.0) l)) l))) (* 2.0 (/ (pow l 2.0) (* t (pow k_m 4.0))))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.65e+120) {
tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t * (pow(t, 2.0) / l)) / l));
} else {
tmp = 2.0 * (pow(l, 2.0) / (t * pow(k_m, 4.0)));
}
return tmp;
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 1.65d+120) then
tmp = 2.0d0 / ((2.0d0 * (k_m * k_m)) * ((t * ((t ** 2.0d0) / l)) / l))
else
tmp = 2.0d0 * ((l ** 2.0d0) / (t * (k_m ** 4.0d0)))
end if
code = tmp
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.65e+120) {
tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t * (Math.pow(t, 2.0) / l)) / l));
} else {
tmp = 2.0 * (Math.pow(l, 2.0) / (t * Math.pow(k_m, 4.0)));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 1.65e+120: tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t * (math.pow(t, 2.0) / l)) / l)) else: tmp = 2.0 * (math.pow(l, 2.0) / (t * math.pow(k_m, 4.0))) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.65e+120) tmp = Float64(2.0 / Float64(Float64(2.0 * Float64(k_m * k_m)) * Float64(Float64(t * Float64((t ^ 2.0) / l)) / l))); else tmp = Float64(2.0 * Float64((l ^ 2.0) / Float64(t * (k_m ^ 4.0)))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 1.65e+120) tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t * ((t ^ 2.0) / l)) / l)); else tmp = 2.0 * ((l ^ 2.0) / (t * (k_m ^ 4.0))); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.65e+120], N[(2.0 / N[(N[(2.0 * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(t * N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[Power[l, 2.0], $MachinePrecision] / N[(t * N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.65 \cdot 10^{+120}:\\
\;\;\;\;\frac{2}{\left(2 \cdot \left(k\_m \cdot k\_m\right)\right) \cdot \frac{t \cdot \frac{{t}^{2}}{\ell}}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{{\ell}^{2}}{t \cdot {k\_m}^{4}}\\
\end{array}
\end{array}
if k < 1.64999999999999995e120Initial program 53.4%
Simplified58.0%
Taylor expanded in k around 0 56.8%
unpow256.8%
Applied egg-rr56.8%
cube-mult56.8%
*-un-lft-identity56.8%
times-frac60.3%
pow260.3%
Applied egg-rr60.3%
if 1.64999999999999995e120 < k Initial program 61.2%
associate-/r*58.2%
+-commutative58.2%
unpow258.2%
sqr-neg58.2%
distribute-frac-neg258.2%
distribute-frac-neg258.2%
unpow258.2%
+-commutative58.2%
associate-*l*58.2%
associate-*l/58.2%
associate-/r/58.2%
+-commutative58.2%
associate-+r+58.2%
Simplified58.2%
Taylor expanded in k around 0 58.6%
*-commutative58.6%
Simplified58.6%
Taylor expanded in t around 0 66.2%
Final simplification61.1%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (if (<= k_m 1.24e+120) (/ 2.0 (* (* 2.0 (* k_m k_m)) (/ (* t (/ (pow t 2.0) l)) l))) (* 2.0 (/ (/ (pow l 2.0) t) (pow k_m 4.0)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.24e+120) {
tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t * (pow(t, 2.0) / l)) / l));
} else {
tmp = 2.0 * ((pow(l, 2.0) / t) / pow(k_m, 4.0));
}
return tmp;
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
real(8) :: tmp
if (k_m <= 1.24d+120) then
tmp = 2.0d0 / ((2.0d0 * (k_m * k_m)) * ((t * ((t ** 2.0d0) / l)) / l))
else
tmp = 2.0d0 * (((l ** 2.0d0) / t) / (k_m ** 4.0d0))
end if
code = tmp
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
double tmp;
if (k_m <= 1.24e+120) {
tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t * (Math.pow(t, 2.0) / l)) / l));
} else {
tmp = 2.0 * ((Math.pow(l, 2.0) / t) / Math.pow(k_m, 4.0));
}
return tmp;
}
k_m = math.fabs(k) def code(t, l, k_m): tmp = 0 if k_m <= 1.24e+120: tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t * (math.pow(t, 2.0) / l)) / l)) else: tmp = 2.0 * ((math.pow(l, 2.0) / t) / math.pow(k_m, 4.0)) return tmp
k_m = abs(k) function code(t, l, k_m) tmp = 0.0 if (k_m <= 1.24e+120) tmp = Float64(2.0 / Float64(Float64(2.0 * Float64(k_m * k_m)) * Float64(Float64(t * Float64((t ^ 2.0) / l)) / l))); else tmp = Float64(2.0 * Float64(Float64((l ^ 2.0) / t) / (k_m ^ 4.0))); end return tmp end
k_m = abs(k); function tmp_2 = code(t, l, k_m) tmp = 0.0; if (k_m <= 1.24e+120) tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t * ((t ^ 2.0) / l)) / l)); else tmp = 2.0 * (((l ^ 2.0) / t) / (k_m ^ 4.0)); end tmp_2 = tmp; end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := If[LessEqual[k$95$m, 1.24e+120], N[(2.0 / N[(N[(2.0 * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(t * N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Power[l, 2.0], $MachinePrecision] / t), $MachinePrecision] / N[Power[k$95$m, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
k_m = \left|k\right|
\\
\begin{array}{l}
\mathbf{if}\;k\_m \leq 1.24 \cdot 10^{+120}:\\
\;\;\;\;\frac{2}{\left(2 \cdot \left(k\_m \cdot k\_m\right)\right) \cdot \frac{t \cdot \frac{{t}^{2}}{\ell}}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{\frac{{\ell}^{2}}{t}}{{k\_m}^{4}}\\
\end{array}
\end{array}
if k < 1.23999999999999998e120Initial program 53.4%
Simplified58.0%
Taylor expanded in k around 0 56.8%
unpow256.8%
Applied egg-rr56.8%
cube-mult56.8%
*-un-lft-identity56.8%
times-frac60.3%
pow260.3%
Applied egg-rr60.3%
if 1.23999999999999998e120 < k Initial program 61.2%
associate-/r*58.2%
+-commutative58.2%
unpow258.2%
sqr-neg58.2%
distribute-frac-neg258.2%
distribute-frac-neg258.2%
unpow258.2%
+-commutative58.2%
associate-*l*58.2%
associate-*l/58.2%
associate-/r/58.2%
+-commutative58.2%
associate-+r+58.2%
Simplified58.2%
Taylor expanded in k around 0 58.6%
*-commutative58.6%
Simplified58.6%
Taylor expanded in t around 0 66.2%
*-commutative66.2%
associate-/r*66.2%
Simplified66.2%
Final simplification61.1%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (/ 2.0 (* (* 2.0 (* k_m k_m)) (/ (* t (/ (pow t 2.0) l)) l))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return 2.0 / ((2.0 * (k_m * k_m)) * ((t * (pow(t, 2.0) / l)) / l));
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = 2.0d0 / ((2.0d0 * (k_m * k_m)) * ((t * ((t ** 2.0d0) / l)) / l))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return 2.0 / ((2.0 * (k_m * k_m)) * ((t * (Math.pow(t, 2.0) / l)) / l));
}
k_m = math.fabs(k) def code(t, l, k_m): return 2.0 / ((2.0 * (k_m * k_m)) * ((t * (math.pow(t, 2.0) / l)) / l))
k_m = abs(k) function code(t, l, k_m) return Float64(2.0 / Float64(Float64(2.0 * Float64(k_m * k_m)) * Float64(Float64(t * Float64((t ^ 2.0) / l)) / l))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = 2.0 / ((2.0 * (k_m * k_m)) * ((t * ((t ^ 2.0) / l)) / l)); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(2.0 / N[(N[(2.0 * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision] * N[(N[(t * N[(N[Power[t, 2.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2}{\left(2 \cdot \left(k\_m \cdot k\_m\right)\right) \cdot \frac{t \cdot \frac{{t}^{2}}{\ell}}{\ell}}
\end{array}
Initial program 54.4%
Simplified58.6%
Taylor expanded in k around 0 56.8%
unpow256.8%
Applied egg-rr56.8%
cube-mult56.8%
*-un-lft-identity56.8%
times-frac59.9%
pow259.9%
Applied egg-rr59.9%
Final simplification59.9%
k_m = (fabs.f64 k) (FPCore (t l k_m) :precision binary64 (/ 2.0 (* (/ (/ (pow t 3.0) l) l) (* 2.0 (* k_m k_m)))))
k_m = fabs(k);
double code(double t, double l, double k_m) {
return 2.0 / (((pow(t, 3.0) / l) / l) * (2.0 * (k_m * k_m)));
}
k_m = abs(k)
real(8) function code(t, l, k_m)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k_m
code = 2.0d0 / ((((t ** 3.0d0) / l) / l) * (2.0d0 * (k_m * k_m)))
end function
k_m = Math.abs(k);
public static double code(double t, double l, double k_m) {
return 2.0 / (((Math.pow(t, 3.0) / l) / l) * (2.0 * (k_m * k_m)));
}
k_m = math.fabs(k) def code(t, l, k_m): return 2.0 / (((math.pow(t, 3.0) / l) / l) * (2.0 * (k_m * k_m)))
k_m = abs(k) function code(t, l, k_m) return Float64(2.0 / Float64(Float64(Float64((t ^ 3.0) / l) / l) * Float64(2.0 * Float64(k_m * k_m)))) end
k_m = abs(k); function tmp = code(t, l, k_m) tmp = 2.0 / ((((t ^ 3.0) / l) / l) * (2.0 * (k_m * k_m))); end
k_m = N[Abs[k], $MachinePrecision] code[t_, l_, k$95$m_] := N[(2.0 / N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision] / l), $MachinePrecision] * N[(2.0 * N[(k$95$m * k$95$m), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
k_m = \left|k\right|
\\
\frac{2}{\frac{\frac{{t}^{3}}{\ell}}{\ell} \cdot \left(2 \cdot \left(k\_m \cdot k\_m\right)\right)}
\end{array}
Initial program 54.4%
Simplified58.6%
Taylor expanded in k around 0 56.8%
unpow256.8%
Applied egg-rr56.8%
herbie shell --seed 2024118
(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))))