
(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 23 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}
(FPCore (t l k)
:precision binary64
(if (or (<= t -4.35e+71) (not (<= t 31000000000.0)))
(/
2.0
(*
(pow (/ t (* (cbrt (/ l (sin k))) (cbrt l))) 3.0)
(* (tan k) (+ 2.0 (pow (/ k t) 2.0)))))
(/ (/ 2.0 (* t (/ k l))) (* (pow (sin k) 2.0) (/ (/ k l) (cos k))))))
double code(double t, double l, double k) {
double tmp;
if ((t <= -4.35e+71) || !(t <= 31000000000.0)) {
tmp = 2.0 / (pow((t / (cbrt((l / sin(k))) * cbrt(l))), 3.0) * (tan(k) * (2.0 + pow((k / t), 2.0))));
} else {
tmp = (2.0 / (t * (k / l))) / (pow(sin(k), 2.0) * ((k / l) / cos(k)));
}
return tmp;
}
public static double code(double t, double l, double k) {
double tmp;
if ((t <= -4.35e+71) || !(t <= 31000000000.0)) {
tmp = 2.0 / (Math.pow((t / (Math.cbrt((l / Math.sin(k))) * Math.cbrt(l))), 3.0) * (Math.tan(k) * (2.0 + Math.pow((k / t), 2.0))));
} else {
tmp = (2.0 / (t * (k / l))) / (Math.pow(Math.sin(k), 2.0) * ((k / l) / Math.cos(k)));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if ((t <= -4.35e+71) || !(t <= 31000000000.0)) tmp = Float64(2.0 / Float64((Float64(t / Float64(cbrt(Float64(l / sin(k))) * cbrt(l))) ^ 3.0) * Float64(tan(k) * Float64(2.0 + (Float64(k / t) ^ 2.0))))); else tmp = Float64(Float64(2.0 / Float64(t * Float64(k / l))) / Float64((sin(k) ^ 2.0) * Float64(Float64(k / l) / cos(k)))); end return tmp end
code[t_, l_, k_] := If[Or[LessEqual[t, -4.35e+71], N[Not[LessEqual[t, 31000000000.0]], $MachinePrecision]], N[(2.0 / N[(N[Power[N[(t / N[(N[Power[N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] * N[Power[l, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] * N[(N[(k / l), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -4.35 \cdot 10^{+71} \lor \neg \left(t \leq 31000000000\right):\\
\;\;\;\;\frac{2}{{\left(\frac{t}{\sqrt[3]{\frac{\ell}{\sin k}} \cdot \sqrt[3]{\ell}}\right)}^{3} \cdot \left(\tan k \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{t \cdot \frac{k}{\ell}}}{{\sin k}^{2} \cdot \frac{\frac{k}{\ell}}{\cos k}}\\
\end{array}
\end{array}
if t < -4.3499999999999999e71 or 3.1e10 < t Initial program 56.7%
associate-/r/55.5%
add-cube-cbrt55.5%
pow355.5%
cbrt-div55.4%
rem-cbrt-cube68.7%
Applied egg-rr68.7%
distribute-lft-in68.7%
associate-/l*68.7%
associate-/l*68.7%
Applied egg-rr68.7%
distribute-lft-out68.7%
+-commutative68.7%
associate-*r*68.7%
associate-/r/68.7%
associate-+r+68.7%
metadata-eval68.7%
Simplified68.7%
cbrt-prod85.2%
Applied egg-rr85.2%
if -4.3499999999999999e71 < t < 3.1e10Initial program 50.2%
Taylor expanded in t around 0 75.5%
associate-/l*74.9%
unpow274.9%
unpow274.9%
associate-*r*74.9%
times-frac84.8%
Simplified84.8%
Taylor expanded in k around inf 75.5%
associate-/l*74.9%
unpow274.9%
unpow274.9%
associate-*r*74.9%
times-frac84.8%
times-frac93.1%
associate-/r/93.2%
*-commutative93.2%
associate-/r/95.6%
Simplified95.6%
Taylor expanded in k around 0 95.0%
div-inv95.0%
*-commutative95.0%
associate-/l*93.2%
Applied egg-rr93.2%
associate-*r/93.2%
metadata-eval93.2%
associate-/r*93.2%
associate-/r/95.7%
*-commutative95.7%
associate-/r*95.7%
Simplified95.7%
Final simplification90.9%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (pow (/ k t) 2.0)))
(if (<=
(/
2.0
(*
(* (tan k) (* (sin k) (/ (pow t 3.0) (* l l))))
(+ 1.0 (+ t_1 1.0))))
1e+296)
(/
2.0
(* (* (tan k) (+ 2.0 t_1)) (pow (/ t (cbrt (* l (/ l (sin k))))) 3.0)))
(/ (/ 2.0 (* t (/ k l))) (* (pow (sin k) 2.0) (/ (/ k l) (cos k)))))))
double code(double t, double l, double k) {
double t_1 = pow((k / t), 2.0);
double tmp;
if ((2.0 / ((tan(k) * (sin(k) * (pow(t, 3.0) / (l * l)))) * (1.0 + (t_1 + 1.0)))) <= 1e+296) {
tmp = 2.0 / ((tan(k) * (2.0 + t_1)) * pow((t / cbrt((l * (l / sin(k))))), 3.0));
} else {
tmp = (2.0 / (t * (k / l))) / (pow(sin(k), 2.0) * ((k / l) / cos(k)));
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = Math.pow((k / t), 2.0);
double tmp;
if ((2.0 / ((Math.tan(k) * (Math.sin(k) * (Math.pow(t, 3.0) / (l * l)))) * (1.0 + (t_1 + 1.0)))) <= 1e+296) {
tmp = 2.0 / ((Math.tan(k) * (2.0 + t_1)) * Math.pow((t / Math.cbrt((l * (l / Math.sin(k))))), 3.0));
} else {
tmp = (2.0 / (t * (k / l))) / (Math.pow(Math.sin(k), 2.0) * ((k / l) / Math.cos(k)));
}
return tmp;
}
function code(t, l, k) t_1 = Float64(k / t) ^ 2.0 tmp = 0.0 if (Float64(2.0 / Float64(Float64(tan(k) * Float64(sin(k) * Float64((t ^ 3.0) / Float64(l * l)))) * Float64(1.0 + Float64(t_1 + 1.0)))) <= 1e+296) tmp = Float64(2.0 / Float64(Float64(tan(k) * Float64(2.0 + t_1)) * (Float64(t / cbrt(Float64(l * Float64(l / sin(k))))) ^ 3.0))); else tmp = Float64(Float64(2.0 / Float64(t * Float64(k / l))) / Float64((sin(k) ^ 2.0) * Float64(Float64(k / l) / cos(k)))); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(t$95$1 + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1e+296], N[(2.0 / N[(N[(N[Tan[k], $MachinePrecision] * N[(2.0 + t$95$1), $MachinePrecision]), $MachinePrecision] * N[Power[N[(t / N[Power[N[(l * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] * N[(N[(k / l), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {\left(\frac{k}{t}\right)}^{2}\\
\mathbf{if}\;\frac{2}{\left(\tan k \cdot \left(\sin k \cdot \frac{{t}^{3}}{\ell \cdot \ell}\right)\right) \cdot \left(1 + \left(t_1 + 1\right)\right)} \leq 10^{+296}:\\
\;\;\;\;\frac{2}{\left(\tan k \cdot \left(2 + t_1\right)\right) \cdot {\left(\frac{t}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{t \cdot \frac{k}{\ell}}}{{\sin k}^{2} \cdot \frac{\frac{k}{\ell}}{\cos k}}\\
\end{array}
\end{array}
if (/.f64 2 (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t 3) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 1 (pow.f64 (/.f64 k t) 2)) 1))) < 9.99999999999999981e295Initial program 81.1%
associate-/r/80.1%
add-cube-cbrt80.0%
pow380.0%
cbrt-div80.0%
rem-cbrt-cube87.0%
Applied egg-rr87.0%
distribute-lft-in87.0%
associate-/l*87.0%
associate-/l*88.7%
Applied egg-rr88.7%
distribute-lft-out88.7%
+-commutative88.7%
associate-*r*88.7%
associate-/r/88.7%
associate-+r+88.7%
metadata-eval88.7%
Simplified88.7%
if 9.99999999999999981e295 < (/.f64 2 (*.f64 (*.f64 (*.f64 (/.f64 (pow.f64 t 3) (*.f64 l l)) (sin.f64 k)) (tan.f64 k)) (+.f64 (+.f64 1 (pow.f64 (/.f64 k t) 2)) 1))) Initial program 20.4%
Taylor expanded in t around 0 56.8%
associate-/l*55.3%
unpow255.3%
unpow255.3%
associate-*r*55.3%
times-frac67.8%
Simplified67.8%
Taylor expanded in k around inf 56.8%
associate-/l*55.3%
unpow255.3%
unpow255.3%
associate-*r*55.3%
times-frac67.8%
times-frac80.2%
associate-/r/80.2%
*-commutative80.2%
associate-/r/83.2%
Simplified83.2%
Taylor expanded in k around 0 81.7%
div-inv81.7%
*-commutative81.7%
associate-/l*80.2%
Applied egg-rr80.2%
associate-*r/80.2%
metadata-eval80.2%
associate-/r*80.4%
associate-/r/83.4%
*-commutative83.4%
associate-/r*83.4%
Simplified83.4%
Final simplification86.3%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (* l (/ l (* k (* k (pow t 3.0)))))))
(if (<= t -7.5e+180)
(cbrt (* t_1 (* t_1 t_1)))
(if (<= t 64000000000.0)
(/ 2.0 (* (/ k l) (* t (* (pow (sin k) 2.0) (/ (/ k l) (cos k))))))
(/
2.0
(*
(+ 2.0 (pow (/ k t) 2.0))
(* (tan k) (* (sin k) (pow (/ (pow t 1.5) l) 2.0)))))))))
double code(double t, double l, double k) {
double t_1 = l * (l / (k * (k * pow(t, 3.0))));
double tmp;
if (t <= -7.5e+180) {
tmp = cbrt((t_1 * (t_1 * t_1)));
} else if (t <= 64000000000.0) {
tmp = 2.0 / ((k / l) * (t * (pow(sin(k), 2.0) * ((k / l) / cos(k)))));
} else {
tmp = 2.0 / ((2.0 + pow((k / t), 2.0)) * (tan(k) * (sin(k) * pow((pow(t, 1.5) / l), 2.0))));
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = l * (l / (k * (k * Math.pow(t, 3.0))));
double tmp;
if (t <= -7.5e+180) {
tmp = Math.cbrt((t_1 * (t_1 * t_1)));
} else if (t <= 64000000000.0) {
tmp = 2.0 / ((k / l) * (t * (Math.pow(Math.sin(k), 2.0) * ((k / l) / Math.cos(k)))));
} else {
tmp = 2.0 / ((2.0 + Math.pow((k / t), 2.0)) * (Math.tan(k) * (Math.sin(k) * Math.pow((Math.pow(t, 1.5) / l), 2.0))));
}
return tmp;
}
function code(t, l, k) t_1 = Float64(l * Float64(l / Float64(k * Float64(k * (t ^ 3.0))))) tmp = 0.0 if (t <= -7.5e+180) tmp = cbrt(Float64(t_1 * Float64(t_1 * t_1))); elseif (t <= 64000000000.0) tmp = Float64(2.0 / Float64(Float64(k / l) * Float64(t * Float64((sin(k) ^ 2.0) * Float64(Float64(k / l) / cos(k)))))); else tmp = Float64(2.0 / Float64(Float64(2.0 + (Float64(k / t) ^ 2.0)) * Float64(tan(k) * Float64(sin(k) * (Float64((t ^ 1.5) / l) ^ 2.0))))); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[(l * N[(l / N[(k * N[(k * N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -7.5e+180], N[Power[N[(t$95$1 * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision], If[LessEqual[t, 64000000000.0], N[(2.0 / N[(N[(k / l), $MachinePrecision] * N[(t * N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] * N[(N[(k / l), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(N[Sin[k], $MachinePrecision] * N[Power[N[(N[Power[t, 1.5], $MachinePrecision] / l), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \ell \cdot \frac{\ell}{k \cdot \left(k \cdot {t}^{3}\right)}\\
\mathbf{if}\;t \leq -7.5 \cdot 10^{+180}:\\
\;\;\;\;\sqrt[3]{t_1 \cdot \left(t_1 \cdot t_1\right)}\\
\mathbf{elif}\;t \leq 64000000000:\\
\;\;\;\;\frac{2}{\frac{k}{\ell} \cdot \left(t \cdot \left({\sin k}^{2} \cdot \frac{\frac{k}{\ell}}{\cos k}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(2 + {\left(\frac{k}{t}\right)}^{2}\right) \cdot \left(\tan k \cdot \left(\sin k \cdot {\left(\frac{{t}^{1.5}}{\ell}\right)}^{2}\right)\right)}\\
\end{array}
\end{array}
if t < -7.5000000000000003e180Initial program 46.2%
associate-/r*46.2%
associate-*l*40.4%
sqr-neg40.4%
associate-*l*46.2%
*-commutative46.2%
sqr-neg46.2%
associate-*l/46.2%
associate-*r/46.2%
associate-/r/46.2%
Simplified46.2%
Taylor expanded in k around 0 40.4%
unpow240.4%
*-commutative40.4%
unpow240.4%
Simplified40.4%
add-cbrt-cube40.4%
associate-/l*40.4%
*-commutative40.4%
associate-/l*40.4%
*-commutative40.4%
associate-/l*60.8%
*-commutative60.8%
Applied egg-rr60.8%
associate-*l*60.8%
associate-/r/60.8%
associate-*l*60.8%
associate-/r/60.8%
associate-*l*60.8%
associate-/r/60.8%
associate-*l*66.8%
Simplified66.8%
if -7.5000000000000003e180 < t < 6.4e10Initial program 49.3%
Taylor expanded in t around 0 71.7%
associate-/l*71.2%
unpow271.2%
unpow271.2%
associate-*r*71.2%
times-frac79.9%
Simplified79.9%
Taylor expanded in k around inf 71.7%
associate-/l*71.2%
unpow271.2%
unpow271.2%
associate-*r*71.2%
times-frac79.9%
times-frac87.9%
associate-/r/88.0%
*-commutative88.0%
associate-/r/90.0%
Simplified90.0%
div-inv90.0%
*-commutative90.0%
Applied egg-rr90.0%
associate-*r/90.0%
metadata-eval90.0%
associate-*l*90.1%
*-commutative90.1%
associate-/r*90.0%
Simplified90.0%
if 6.4e10 < t Initial program 64.8%
distribute-lft-in64.8%
associate-*l*57.1%
add-sqr-sqrt57.1%
pow257.1%
sqrt-div57.1%
sqrt-pow157.2%
metadata-eval57.2%
sqrt-prod33.0%
add-sqr-sqrt57.2%
Applied egg-rr70.7%
distribute-lft-out70.7%
+-commutative70.7%
associate-+r+70.7%
metadata-eval70.7%
associate-*l*70.7%
unpow270.7%
times-frac57.0%
pow-sqr56.9%
metadata-eval56.9%
unpow256.9%
associate-*r*56.9%
Simplified84.0%
Final simplification86.7%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (pow (sin k) 2.0)))
(if (<= k 3.9e-39)
(/
2.0
(*
(pow (/ t (cbrt (* l (/ l (sin k))))) 3.0)
(/ 2.0 (/ (cos k) (sin k)))))
(if (<= k 5e+223)
(/ 2.0 (* (* t_1 (/ k (* l (cos k)))) (/ (* t k) l)))
(/ (/ 2.0 (* t (/ k l))) (* t_1 (/ (/ k l) (cos k))))))))
double code(double t, double l, double k) {
double t_1 = pow(sin(k), 2.0);
double tmp;
if (k <= 3.9e-39) {
tmp = 2.0 / (pow((t / cbrt((l * (l / sin(k))))), 3.0) * (2.0 / (cos(k) / sin(k))));
} else if (k <= 5e+223) {
tmp = 2.0 / ((t_1 * (k / (l * cos(k)))) * ((t * k) / l));
} else {
tmp = (2.0 / (t * (k / l))) / (t_1 * ((k / l) / cos(k)));
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = Math.pow(Math.sin(k), 2.0);
double tmp;
if (k <= 3.9e-39) {
tmp = 2.0 / (Math.pow((t / Math.cbrt((l * (l / Math.sin(k))))), 3.0) * (2.0 / (Math.cos(k) / Math.sin(k))));
} else if (k <= 5e+223) {
tmp = 2.0 / ((t_1 * (k / (l * Math.cos(k)))) * ((t * k) / l));
} else {
tmp = (2.0 / (t * (k / l))) / (t_1 * ((k / l) / Math.cos(k)));
}
return tmp;
}
function code(t, l, k) t_1 = sin(k) ^ 2.0 tmp = 0.0 if (k <= 3.9e-39) tmp = Float64(2.0 / Float64((Float64(t / cbrt(Float64(l * Float64(l / sin(k))))) ^ 3.0) * Float64(2.0 / Float64(cos(k) / sin(k))))); elseif (k <= 5e+223) tmp = Float64(2.0 / Float64(Float64(t_1 * Float64(k / Float64(l * cos(k)))) * Float64(Float64(t * k) / l))); else tmp = Float64(Float64(2.0 / Float64(t * Float64(k / l))) / Float64(t_1 * Float64(Float64(k / l) / cos(k)))); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[k, 3.9e-39], N[(2.0 / N[(N[Power[N[(t / N[Power[N[(l * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(2.0 / N[(N[Cos[k], $MachinePrecision] / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 5e+223], N[(2.0 / N[(N[(t$95$1 * N[(k / N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(t * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$1 * N[(N[(k / l), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {\sin k}^{2}\\
\mathbf{if}\;k \leq 3.9 \cdot 10^{-39}:\\
\;\;\;\;\frac{2}{{\left(\frac{t}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3} \cdot \frac{2}{\frac{\cos k}{\sin k}}}\\
\mathbf{elif}\;k \leq 5 \cdot 10^{+223}:\\
\;\;\;\;\frac{2}{\left(t_1 \cdot \frac{k}{\ell \cdot \cos k}\right) \cdot \frac{t \cdot k}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{t \cdot \frac{k}{\ell}}}{t_1 \cdot \frac{\frac{k}{\ell}}{\cos k}}\\
\end{array}
\end{array}
if k < 3.9000000000000003e-39Initial program 52.5%
associate-/r/51.8%
add-cube-cbrt51.7%
pow351.7%
cbrt-div51.7%
rem-cbrt-cube61.3%
Applied egg-rr61.3%
distribute-lft-in61.3%
associate-/l*61.3%
associate-/l*66.1%
Applied egg-rr66.1%
distribute-lft-out66.1%
+-commutative66.1%
associate-*r*66.1%
associate-/r/66.1%
associate-+r+66.1%
metadata-eval66.1%
Simplified66.1%
Taylor expanded in t around inf 63.4%
associate-*r/63.4%
associate-/l*63.4%
Simplified63.4%
if 3.9000000000000003e-39 < k < 4.99999999999999985e223Initial program 57.5%
Taylor expanded in t around 0 82.3%
associate-/l*80.4%
unpow280.4%
unpow280.4%
associate-*r*80.4%
times-frac84.0%
Simplified84.0%
Taylor expanded in k around inf 82.3%
associate-/l*80.4%
unpow280.4%
unpow280.4%
associate-*r*80.4%
times-frac84.0%
times-frac87.9%
associate-/r/87.9%
*-commutative87.9%
associate-/r/89.8%
Simplified89.8%
Taylor expanded in k around 0 93.8%
if 4.99999999999999985e223 < k Initial program 46.2%
Taylor expanded in t around 0 62.3%
associate-/l*62.0%
unpow262.0%
unpow262.0%
associate-*r*62.0%
times-frac62.4%
Simplified62.4%
Taylor expanded in k around inf 62.3%
associate-/l*62.0%
unpow262.0%
unpow262.0%
associate-*r*62.0%
times-frac62.4%
times-frac99.1%
associate-/r/99.4%
*-commutative99.4%
associate-/r/99.7%
Simplified99.7%
Taylor expanded in k around 0 77.9%
div-inv77.9%
*-commutative77.9%
associate-/l*99.4%
Applied egg-rr99.4%
associate-*r/99.4%
metadata-eval99.4%
associate-/r*99.4%
associate-/r/99.4%
*-commutative99.4%
associate-/r*99.5%
Simplified99.5%
Final simplification71.0%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (* (pow (sin k) 2.0) (/ k (* l (cos k))))))
(if (<= k 1.95e-40)
(/ 2.0 (* (pow (/ t (cbrt (* l (/ l (sin k))))) 3.0) (* 2.0 k)))
(if (<= k 9.5e+218)
(/ 2.0 (* t_1 (/ (* t k) l)))
(/ 2.0 (* (* t (/ k l)) t_1))))))
double code(double t, double l, double k) {
double t_1 = pow(sin(k), 2.0) * (k / (l * cos(k)));
double tmp;
if (k <= 1.95e-40) {
tmp = 2.0 / (pow((t / cbrt((l * (l / sin(k))))), 3.0) * (2.0 * k));
} else if (k <= 9.5e+218) {
tmp = 2.0 / (t_1 * ((t * k) / l));
} else {
tmp = 2.0 / ((t * (k / l)) * t_1);
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = Math.pow(Math.sin(k), 2.0) * (k / (l * Math.cos(k)));
double tmp;
if (k <= 1.95e-40) {
tmp = 2.0 / (Math.pow((t / Math.cbrt((l * (l / Math.sin(k))))), 3.0) * (2.0 * k));
} else if (k <= 9.5e+218) {
tmp = 2.0 / (t_1 * ((t * k) / l));
} else {
tmp = 2.0 / ((t * (k / l)) * t_1);
}
return tmp;
}
function code(t, l, k) t_1 = Float64((sin(k) ^ 2.0) * Float64(k / Float64(l * cos(k)))) tmp = 0.0 if (k <= 1.95e-40) tmp = Float64(2.0 / Float64((Float64(t / cbrt(Float64(l * Float64(l / sin(k))))) ^ 3.0) * Float64(2.0 * k))); elseif (k <= 9.5e+218) tmp = Float64(2.0 / Float64(t_1 * Float64(Float64(t * k) / l))); else tmp = Float64(2.0 / Float64(Float64(t * Float64(k / l)) * t_1)); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] * N[(k / N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, 1.95e-40], N[(2.0 / N[(N[Power[N[(t / N[Power[N[(l * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(2.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 9.5e+218], N[(2.0 / N[(t$95$1 * N[(N[(t * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {\sin k}^{2} \cdot \frac{k}{\ell \cdot \cos k}\\
\mathbf{if}\;k \leq 1.95 \cdot 10^{-40}:\\
\;\;\;\;\frac{2}{{\left(\frac{t}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3} \cdot \left(2 \cdot k\right)}\\
\mathbf{elif}\;k \leq 9.5 \cdot 10^{+218}:\\
\;\;\;\;\frac{2}{t_1 \cdot \frac{t \cdot k}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(t \cdot \frac{k}{\ell}\right) \cdot t_1}\\
\end{array}
\end{array}
if k < 1.9499999999999999e-40Initial program 52.5%
associate-/r/51.8%
add-cube-cbrt51.7%
pow351.7%
cbrt-div51.7%
rem-cbrt-cube61.3%
Applied egg-rr61.3%
distribute-lft-in61.3%
associate-/l*61.3%
associate-/l*66.1%
Applied egg-rr66.1%
distribute-lft-out66.1%
+-commutative66.1%
associate-*r*66.1%
associate-/r/66.1%
associate-+r+66.1%
metadata-eval66.1%
Simplified66.1%
Taylor expanded in k around 0 63.2%
*-commutative63.2%
Simplified63.2%
if 1.9499999999999999e-40 < k < 9.4999999999999999e218Initial program 57.5%
Taylor expanded in t around 0 82.3%
associate-/l*80.4%
unpow280.4%
unpow280.4%
associate-*r*80.4%
times-frac84.0%
Simplified84.0%
Taylor expanded in k around inf 82.3%
associate-/l*80.4%
unpow280.4%
unpow280.4%
associate-*r*80.4%
times-frac84.0%
times-frac87.9%
associate-/r/87.9%
*-commutative87.9%
associate-/r/89.8%
Simplified89.8%
Taylor expanded in k around 0 93.8%
if 9.4999999999999999e218 < k Initial program 46.2%
Taylor expanded in t around 0 62.3%
associate-/l*62.0%
unpow262.0%
unpow262.0%
associate-*r*62.0%
times-frac62.4%
Simplified62.4%
Taylor expanded in k around inf 62.3%
associate-/l*62.0%
unpow262.0%
unpow262.0%
associate-*r*62.0%
times-frac62.4%
times-frac99.1%
associate-/r/99.4%
*-commutative99.4%
associate-/r/99.7%
Simplified99.7%
Final simplification70.9%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (pow (sin k) 2.0)))
(if (<= k 1.45e-40)
(/ 2.0 (* (pow (/ t (cbrt (* l (/ l (sin k))))) 3.0) (* 2.0 k)))
(if (<= k 4e+218)
(/ 2.0 (* (* t_1 (/ k (* l (cos k)))) (/ (* t k) l)))
(/ 2.0 (* (/ k l) (* t (* t_1 (/ (/ k l) (cos k))))))))))
double code(double t, double l, double k) {
double t_1 = pow(sin(k), 2.0);
double tmp;
if (k <= 1.45e-40) {
tmp = 2.0 / (pow((t / cbrt((l * (l / sin(k))))), 3.0) * (2.0 * k));
} else if (k <= 4e+218) {
tmp = 2.0 / ((t_1 * (k / (l * cos(k)))) * ((t * k) / l));
} else {
tmp = 2.0 / ((k / l) * (t * (t_1 * ((k / l) / cos(k)))));
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = Math.pow(Math.sin(k), 2.0);
double tmp;
if (k <= 1.45e-40) {
tmp = 2.0 / (Math.pow((t / Math.cbrt((l * (l / Math.sin(k))))), 3.0) * (2.0 * k));
} else if (k <= 4e+218) {
tmp = 2.0 / ((t_1 * (k / (l * Math.cos(k)))) * ((t * k) / l));
} else {
tmp = 2.0 / ((k / l) * (t * (t_1 * ((k / l) / Math.cos(k)))));
}
return tmp;
}
function code(t, l, k) t_1 = sin(k) ^ 2.0 tmp = 0.0 if (k <= 1.45e-40) tmp = Float64(2.0 / Float64((Float64(t / cbrt(Float64(l * Float64(l / sin(k))))) ^ 3.0) * Float64(2.0 * k))); elseif (k <= 4e+218) tmp = Float64(2.0 / Float64(Float64(t_1 * Float64(k / Float64(l * cos(k)))) * Float64(Float64(t * k) / l))); else tmp = Float64(2.0 / Float64(Float64(k / l) * Float64(t * Float64(t_1 * Float64(Float64(k / l) / cos(k)))))); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[k, 1.45e-40], N[(2.0 / N[(N[Power[N[(t / N[Power[N[(l * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(2.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 4e+218], N[(2.0 / N[(N[(t$95$1 * N[(k / N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(t * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(k / l), $MachinePrecision] * N[(t * N[(t$95$1 * N[(N[(k / l), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {\sin k}^{2}\\
\mathbf{if}\;k \leq 1.45 \cdot 10^{-40}:\\
\;\;\;\;\frac{2}{{\left(\frac{t}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3} \cdot \left(2 \cdot k\right)}\\
\mathbf{elif}\;k \leq 4 \cdot 10^{+218}:\\
\;\;\;\;\frac{2}{\left(t_1 \cdot \frac{k}{\ell \cdot \cos k}\right) \cdot \frac{t \cdot k}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{k}{\ell} \cdot \left(t \cdot \left(t_1 \cdot \frac{\frac{k}{\ell}}{\cos k}\right)\right)}\\
\end{array}
\end{array}
if k < 1.4499999999999999e-40Initial program 52.5%
associate-/r/51.8%
add-cube-cbrt51.7%
pow351.7%
cbrt-div51.7%
rem-cbrt-cube61.3%
Applied egg-rr61.3%
distribute-lft-in61.3%
associate-/l*61.3%
associate-/l*66.1%
Applied egg-rr66.1%
distribute-lft-out66.1%
+-commutative66.1%
associate-*r*66.1%
associate-/r/66.1%
associate-+r+66.1%
metadata-eval66.1%
Simplified66.1%
Taylor expanded in k around 0 63.2%
*-commutative63.2%
Simplified63.2%
if 1.4499999999999999e-40 < k < 4.00000000000000033e218Initial program 57.5%
Taylor expanded in t around 0 82.3%
associate-/l*80.4%
unpow280.4%
unpow280.4%
associate-*r*80.4%
times-frac84.0%
Simplified84.0%
Taylor expanded in k around inf 82.3%
associate-/l*80.4%
unpow280.4%
unpow280.4%
associate-*r*80.4%
times-frac84.0%
times-frac87.9%
associate-/r/87.9%
*-commutative87.9%
associate-/r/89.8%
Simplified89.8%
Taylor expanded in k around 0 93.8%
if 4.00000000000000033e218 < k Initial program 46.2%
Taylor expanded in t around 0 62.3%
associate-/l*62.0%
unpow262.0%
unpow262.0%
associate-*r*62.0%
times-frac62.4%
Simplified62.4%
Taylor expanded in k around inf 62.3%
associate-/l*62.0%
unpow262.0%
unpow262.0%
associate-*r*62.0%
times-frac62.4%
times-frac99.1%
associate-/r/99.4%
*-commutative99.4%
associate-/r/99.7%
Simplified99.7%
div-inv99.7%
*-commutative99.7%
Applied egg-rr99.7%
associate-*r/99.7%
metadata-eval99.7%
associate-*l*99.6%
*-commutative99.6%
associate-/r*99.6%
Simplified99.6%
Final simplification70.9%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (pow (sin k) 2.0)) (t_2 (* l (cos k))))
(if (<= k 1.4e-40)
(/ 2.0 (* (pow (/ t (cbrt (* l (/ l (sin k))))) 3.0) (* 2.0 k)))
(if (<= k 2.45e+218)
(/ 2.0 (* (* t_1 (/ k t_2)) (/ (* t k) l)))
(/ 2.0 (* (* t (/ k l)) (/ (* k t_1) t_2)))))))
double code(double t, double l, double k) {
double t_1 = pow(sin(k), 2.0);
double t_2 = l * cos(k);
double tmp;
if (k <= 1.4e-40) {
tmp = 2.0 / (pow((t / cbrt((l * (l / sin(k))))), 3.0) * (2.0 * k));
} else if (k <= 2.45e+218) {
tmp = 2.0 / ((t_1 * (k / t_2)) * ((t * k) / l));
} else {
tmp = 2.0 / ((t * (k / l)) * ((k * t_1) / t_2));
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = Math.pow(Math.sin(k), 2.0);
double t_2 = l * Math.cos(k);
double tmp;
if (k <= 1.4e-40) {
tmp = 2.0 / (Math.pow((t / Math.cbrt((l * (l / Math.sin(k))))), 3.0) * (2.0 * k));
} else if (k <= 2.45e+218) {
tmp = 2.0 / ((t_1 * (k / t_2)) * ((t * k) / l));
} else {
tmp = 2.0 / ((t * (k / l)) * ((k * t_1) / t_2));
}
return tmp;
}
function code(t, l, k) t_1 = sin(k) ^ 2.0 t_2 = Float64(l * cos(k)) tmp = 0.0 if (k <= 1.4e-40) tmp = Float64(2.0 / Float64((Float64(t / cbrt(Float64(l * Float64(l / sin(k))))) ^ 3.0) * Float64(2.0 * k))); elseif (k <= 2.45e+218) tmp = Float64(2.0 / Float64(Float64(t_1 * Float64(k / t_2)) * Float64(Float64(t * k) / l))); else tmp = Float64(2.0 / Float64(Float64(t * Float64(k / l)) * Float64(Float64(k * t_1) / t_2))); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, 1.4e-40], N[(2.0 / N[(N[Power[N[(t / N[Power[N[(l * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(2.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 2.45e+218], N[(2.0 / N[(N[(t$95$1 * N[(k / t$95$2), $MachinePrecision]), $MachinePrecision] * N[(N[(t * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision] * N[(N[(k * t$95$1), $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {\sin k}^{2}\\
t_2 := \ell \cdot \cos k\\
\mathbf{if}\;k \leq 1.4 \cdot 10^{-40}:\\
\;\;\;\;\frac{2}{{\left(\frac{t}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3} \cdot \left(2 \cdot k\right)}\\
\mathbf{elif}\;k \leq 2.45 \cdot 10^{+218}:\\
\;\;\;\;\frac{2}{\left(t_1 \cdot \frac{k}{t_2}\right) \cdot \frac{t \cdot k}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(t \cdot \frac{k}{\ell}\right) \cdot \frac{k \cdot t_1}{t_2}}\\
\end{array}
\end{array}
if k < 1.4e-40Initial program 52.5%
associate-/r/51.8%
add-cube-cbrt51.7%
pow351.7%
cbrt-div51.7%
rem-cbrt-cube61.3%
Applied egg-rr61.3%
distribute-lft-in61.3%
associate-/l*61.3%
associate-/l*66.1%
Applied egg-rr66.1%
distribute-lft-out66.1%
+-commutative66.1%
associate-*r*66.1%
associate-/r/66.1%
associate-+r+66.1%
metadata-eval66.1%
Simplified66.1%
Taylor expanded in k around 0 63.2%
*-commutative63.2%
Simplified63.2%
if 1.4e-40 < k < 2.4499999999999999e218Initial program 57.5%
Taylor expanded in t around 0 82.3%
associate-/l*80.4%
unpow280.4%
unpow280.4%
associate-*r*80.4%
times-frac84.0%
Simplified84.0%
Taylor expanded in k around inf 82.3%
associate-/l*80.4%
unpow280.4%
unpow280.4%
associate-*r*80.4%
times-frac84.0%
times-frac87.9%
associate-/r/87.9%
*-commutative87.9%
associate-/r/89.8%
Simplified89.8%
Taylor expanded in k around 0 93.8%
if 2.4499999999999999e218 < k Initial program 46.2%
Taylor expanded in t around 0 62.3%
associate-/l*62.0%
unpow262.0%
unpow262.0%
associate-*r*62.0%
times-frac62.4%
Simplified62.4%
Taylor expanded in k around inf 62.3%
associate-/l*62.0%
unpow262.0%
unpow262.0%
associate-*r*62.0%
times-frac62.4%
times-frac99.1%
associate-/r/99.4%
*-commutative99.4%
associate-/r/99.7%
Simplified99.7%
Taylor expanded in k around inf 99.7%
Final simplification70.9%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (pow (sin k) 2.0)))
(if (<= k 1.8e-41)
(/ 2.0 (* (pow (/ t (cbrt (* l (/ l (sin k))))) 3.0) (* 2.0 k)))
(if (<= k 2.3e+218)
(/ 2.0 (* (* t_1 (/ k (* l (cos k)))) (/ (* t k) l)))
(/ (/ 2.0 (* t (/ k l))) (* t_1 (/ (/ k l) (cos k))))))))
double code(double t, double l, double k) {
double t_1 = pow(sin(k), 2.0);
double tmp;
if (k <= 1.8e-41) {
tmp = 2.0 / (pow((t / cbrt((l * (l / sin(k))))), 3.0) * (2.0 * k));
} else if (k <= 2.3e+218) {
tmp = 2.0 / ((t_1 * (k / (l * cos(k)))) * ((t * k) / l));
} else {
tmp = (2.0 / (t * (k / l))) / (t_1 * ((k / l) / cos(k)));
}
return tmp;
}
public static double code(double t, double l, double k) {
double t_1 = Math.pow(Math.sin(k), 2.0);
double tmp;
if (k <= 1.8e-41) {
tmp = 2.0 / (Math.pow((t / Math.cbrt((l * (l / Math.sin(k))))), 3.0) * (2.0 * k));
} else if (k <= 2.3e+218) {
tmp = 2.0 / ((t_1 * (k / (l * Math.cos(k)))) * ((t * k) / l));
} else {
tmp = (2.0 / (t * (k / l))) / (t_1 * ((k / l) / Math.cos(k)));
}
return tmp;
}
function code(t, l, k) t_1 = sin(k) ^ 2.0 tmp = 0.0 if (k <= 1.8e-41) tmp = Float64(2.0 / Float64((Float64(t / cbrt(Float64(l * Float64(l / sin(k))))) ^ 3.0) * Float64(2.0 * k))); elseif (k <= 2.3e+218) tmp = Float64(2.0 / Float64(Float64(t_1 * Float64(k / Float64(l * cos(k)))) * Float64(Float64(t * k) / l))); else tmp = Float64(Float64(2.0 / Float64(t * Float64(k / l))) / Float64(t_1 * Float64(Float64(k / l) / cos(k)))); end return tmp end
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[k, 1.8e-41], N[(2.0 / N[(N[Power[N[(t / N[Power[N[(l * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(2.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 2.3e+218], N[(2.0 / N[(N[(t$95$1 * N[(k / N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(t * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$1 * N[(N[(k / l), $MachinePrecision] / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := {\sin k}^{2}\\
\mathbf{if}\;k \leq 1.8 \cdot 10^{-41}:\\
\;\;\;\;\frac{2}{{\left(\frac{t}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3} \cdot \left(2 \cdot k\right)}\\
\mathbf{elif}\;k \leq 2.3 \cdot 10^{+218}:\\
\;\;\;\;\frac{2}{\left(t_1 \cdot \frac{k}{\ell \cdot \cos k}\right) \cdot \frac{t \cdot k}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{t \cdot \frac{k}{\ell}}}{t_1 \cdot \frac{\frac{k}{\ell}}{\cos k}}\\
\end{array}
\end{array}
if k < 1.8e-41Initial program 52.5%
associate-/r/51.8%
add-cube-cbrt51.7%
pow351.7%
cbrt-div51.7%
rem-cbrt-cube61.3%
Applied egg-rr61.3%
distribute-lft-in61.3%
associate-/l*61.3%
associate-/l*66.1%
Applied egg-rr66.1%
distribute-lft-out66.1%
+-commutative66.1%
associate-*r*66.1%
associate-/r/66.1%
associate-+r+66.1%
metadata-eval66.1%
Simplified66.1%
Taylor expanded in k around 0 63.2%
*-commutative63.2%
Simplified63.2%
if 1.8e-41 < k < 2.3000000000000001e218Initial program 57.5%
Taylor expanded in t around 0 82.3%
associate-/l*80.4%
unpow280.4%
unpow280.4%
associate-*r*80.4%
times-frac84.0%
Simplified84.0%
Taylor expanded in k around inf 82.3%
associate-/l*80.4%
unpow280.4%
unpow280.4%
associate-*r*80.4%
times-frac84.0%
times-frac87.9%
associate-/r/87.9%
*-commutative87.9%
associate-/r/89.8%
Simplified89.8%
Taylor expanded in k around 0 93.8%
if 2.3000000000000001e218 < k Initial program 46.2%
Taylor expanded in t around 0 62.3%
associate-/l*62.0%
unpow262.0%
unpow262.0%
associate-*r*62.0%
times-frac62.4%
Simplified62.4%
Taylor expanded in k around inf 62.3%
associate-/l*62.0%
unpow262.0%
unpow262.0%
associate-*r*62.0%
times-frac62.4%
times-frac99.1%
associate-/r/99.4%
*-commutative99.4%
associate-/r/99.7%
Simplified99.7%
Taylor expanded in k around 0 77.9%
div-inv77.9%
*-commutative77.9%
associate-/l*99.4%
Applied egg-rr99.4%
associate-*r/99.4%
metadata-eval99.4%
associate-/r*99.4%
associate-/r/99.4%
*-commutative99.4%
associate-/r*99.5%
Simplified99.5%
Final simplification70.9%
(FPCore (t l k) :precision binary64 (if (<= k 1.35e-13) (/ 2.0 (* (pow (/ t (cbrt (* l (/ l (sin k))))) 3.0) (* 2.0 k))) (* 2.0 (* (/ (cos k) (* k k)) (/ (* l l) (* t (pow (sin k) 2.0)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 1.35e-13) {
tmp = 2.0 / (pow((t / cbrt((l * (l / sin(k))))), 3.0) * (2.0 * k));
} else {
tmp = 2.0 * ((cos(k) / (k * k)) * ((l * l) / (t * pow(sin(k), 2.0))));
}
return tmp;
}
public static double code(double t, double l, double k) {
double tmp;
if (k <= 1.35e-13) {
tmp = 2.0 / (Math.pow((t / Math.cbrt((l * (l / Math.sin(k))))), 3.0) * (2.0 * k));
} else {
tmp = 2.0 * ((Math.cos(k) / (k * k)) * ((l * l) / (t * Math.pow(Math.sin(k), 2.0))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (k <= 1.35e-13) tmp = Float64(2.0 / Float64((Float64(t / cbrt(Float64(l * Float64(l / sin(k))))) ^ 3.0) * Float64(2.0 * k))); else tmp = Float64(2.0 * Float64(Float64(cos(k) / Float64(k * k)) * Float64(Float64(l * l) / Float64(t * (sin(k) ^ 2.0))))); end return tmp end
code[t_, l_, k_] := If[LessEqual[k, 1.35e-13], N[(2.0 / N[(N[Power[N[(t / N[Power[N[(l * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(2.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(l * l), $MachinePrecision] / N[(t * N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.35 \cdot 10^{-13}:\\
\;\;\;\;\frac{2}{{\left(\frac{t}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3} \cdot \left(2 \cdot k\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\frac{\cos k}{k \cdot k} \cdot \frac{\ell \cdot \ell}{t \cdot {\sin k}^{2}}\right)\\
\end{array}
\end{array}
if k < 1.35000000000000005e-13Initial program 53.6%
associate-/r/53.0%
add-cube-cbrt52.9%
pow352.9%
cbrt-div52.9%
rem-cbrt-cube62.2%
Applied egg-rr62.2%
distribute-lft-in62.2%
associate-/l*62.2%
associate-/l*66.8%
Applied egg-rr66.8%
distribute-lft-out66.8%
+-commutative66.8%
associate-*r*66.8%
associate-/r/66.8%
associate-+r+66.8%
metadata-eval66.8%
Simplified66.8%
Taylor expanded in k around 0 63.6%
*-commutative63.6%
Simplified63.6%
if 1.35000000000000005e-13 < k Initial program 51.3%
Taylor expanded in t around 0 77.0%
associate-/l*75.3%
unpow275.3%
unpow275.3%
associate-*r*75.3%
times-frac76.9%
Simplified76.9%
Taylor expanded in k around inf 77.0%
times-frac75.3%
unpow275.3%
unpow275.3%
*-commutative75.3%
Simplified75.3%
Final simplification66.1%
(FPCore (t l k) :precision binary64 (if (<= k 1.15e-13) (/ 2.0 (* (pow (/ t (cbrt (* l (/ l (sin k))))) 3.0) (* 2.0 k))) (* 2.0 (/ (* (cos k) (* l l)) (* (* k k) (* t (pow (sin k) 2.0)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 1.15e-13) {
tmp = 2.0 / (pow((t / cbrt((l * (l / sin(k))))), 3.0) * (2.0 * k));
} else {
tmp = 2.0 * ((cos(k) * (l * l)) / ((k * k) * (t * pow(sin(k), 2.0))));
}
return tmp;
}
public static double code(double t, double l, double k) {
double tmp;
if (k <= 1.15e-13) {
tmp = 2.0 / (Math.pow((t / Math.cbrt((l * (l / Math.sin(k))))), 3.0) * (2.0 * k));
} else {
tmp = 2.0 * ((Math.cos(k) * (l * l)) / ((k * k) * (t * Math.pow(Math.sin(k), 2.0))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (k <= 1.15e-13) tmp = Float64(2.0 / Float64((Float64(t / cbrt(Float64(l * Float64(l / sin(k))))) ^ 3.0) * Float64(2.0 * k))); else tmp = Float64(2.0 * Float64(Float64(cos(k) * Float64(l * l)) / Float64(Float64(k * k) * Float64(t * (sin(k) ^ 2.0))))); end return tmp end
code[t_, l_, k_] := If[LessEqual[k, 1.15e-13], N[(2.0 / N[(N[Power[N[(t / N[Power[N[(l * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(2.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[(l * l), $MachinePrecision]), $MachinePrecision] / N[(N[(k * k), $MachinePrecision] * N[(t * N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.15 \cdot 10^{-13}:\\
\;\;\;\;\frac{2}{{\left(\frac{t}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3} \cdot \left(2 \cdot k\right)}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{\cos k \cdot \left(\ell \cdot \ell\right)}{\left(k \cdot k\right) \cdot \left(t \cdot {\sin k}^{2}\right)}\\
\end{array}
\end{array}
if k < 1.1499999999999999e-13Initial program 53.6%
associate-/r/53.0%
add-cube-cbrt52.9%
pow352.9%
cbrt-div52.9%
rem-cbrt-cube62.2%
Applied egg-rr62.2%
distribute-lft-in62.2%
associate-/l*62.2%
associate-/l*66.8%
Applied egg-rr66.8%
distribute-lft-out66.8%
+-commutative66.8%
associate-*r*66.8%
associate-/r/66.8%
associate-+r+66.8%
metadata-eval66.8%
Simplified66.8%
Taylor expanded in k around 0 63.6%
*-commutative63.6%
Simplified63.6%
if 1.1499999999999999e-13 < k Initial program 51.3%
associate-/r*51.3%
associate-*l*51.3%
sqr-neg51.3%
associate-*l*51.3%
*-commutative51.3%
sqr-neg51.3%
associate-*l/51.3%
associate-*r/51.3%
associate-/r/51.4%
Simplified51.4%
Taylor expanded in k around inf 77.0%
*-commutative77.0%
unpow277.0%
unpow277.0%
*-commutative77.0%
Simplified77.0%
Final simplification66.4%
(FPCore (t l k) :precision binary64 (if (<= k 1.2e-40) (/ 2.0 (* (pow (/ t (cbrt (* l (/ l (sin k))))) 3.0) (* 2.0 k))) (* (/ 2.0 (* k k)) (* (/ (* l (cos k)) (pow (sin k) 2.0)) (/ l t)))))
double code(double t, double l, double k) {
double tmp;
if (k <= 1.2e-40) {
tmp = 2.0 / (pow((t / cbrt((l * (l / sin(k))))), 3.0) * (2.0 * k));
} else {
tmp = (2.0 / (k * k)) * (((l * cos(k)) / pow(sin(k), 2.0)) * (l / t));
}
return tmp;
}
public static double code(double t, double l, double k) {
double tmp;
if (k <= 1.2e-40) {
tmp = 2.0 / (Math.pow((t / Math.cbrt((l * (l / Math.sin(k))))), 3.0) * (2.0 * k));
} else {
tmp = (2.0 / (k * k)) * (((l * Math.cos(k)) / Math.pow(Math.sin(k), 2.0)) * (l / t));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (k <= 1.2e-40) tmp = Float64(2.0 / Float64((Float64(t / cbrt(Float64(l * Float64(l / sin(k))))) ^ 3.0) * Float64(2.0 * k))); else tmp = Float64(Float64(2.0 / Float64(k * k)) * Float64(Float64(Float64(l * cos(k)) / (sin(k) ^ 2.0)) * Float64(l / t))); end return tmp end
code[t_, l_, k_] := If[LessEqual[k, 1.2e-40], N[(2.0 / N[(N[Power[N[(t / N[Power[N[(l * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(2.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.2 \cdot 10^{-40}:\\
\;\;\;\;\frac{2}{{\left(\frac{t}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3} \cdot \left(2 \cdot k\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{k \cdot k} \cdot \left(\frac{\ell \cdot \cos k}{{\sin k}^{2}} \cdot \frac{\ell}{t}\right)\\
\end{array}
\end{array}
if k < 1.19999999999999996e-40Initial program 52.5%
associate-/r/51.8%
add-cube-cbrt51.7%
pow351.7%
cbrt-div51.7%
rem-cbrt-cube61.3%
Applied egg-rr61.3%
distribute-lft-in61.3%
associate-/l*61.3%
associate-/l*66.1%
Applied egg-rr66.1%
distribute-lft-out66.1%
+-commutative66.1%
associate-*r*66.1%
associate-/r/66.1%
associate-+r+66.1%
metadata-eval66.1%
Simplified66.1%
Taylor expanded in k around 0 63.2%
*-commutative63.2%
Simplified63.2%
if 1.19999999999999996e-40 < k Initial program 55.1%
associate-/r*55.1%
associate-*l*55.2%
sqr-neg55.2%
associate-*l*55.1%
*-commutative55.1%
sqr-neg55.1%
associate-*l/55.2%
associate-*r/55.2%
associate-/r/53.6%
Simplified53.6%
Taylor expanded in k around inf 78.1%
associate-*r/78.1%
times-frac76.5%
unpow276.5%
unpow276.5%
associate-*r*76.5%
times-frac79.4%
Simplified79.4%
Final simplification67.1%
(FPCore (t l k) :precision binary64 (if (<= k 4.2e-39) (/ 2.0 (* (pow (/ t (cbrt (* l (/ l (sin k))))) 3.0) (* 2.0 k))) (/ 2.0 (* (* t (/ k l)) (* (pow (sin k) 2.0) (/ k (* l (cos k))))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 4.2e-39) {
tmp = 2.0 / (pow((t / cbrt((l * (l / sin(k))))), 3.0) * (2.0 * k));
} else {
tmp = 2.0 / ((t * (k / l)) * (pow(sin(k), 2.0) * (k / (l * cos(k)))));
}
return tmp;
}
public static double code(double t, double l, double k) {
double tmp;
if (k <= 4.2e-39) {
tmp = 2.0 / (Math.pow((t / Math.cbrt((l * (l / Math.sin(k))))), 3.0) * (2.0 * k));
} else {
tmp = 2.0 / ((t * (k / l)) * (Math.pow(Math.sin(k), 2.0) * (k / (l * Math.cos(k)))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (k <= 4.2e-39) tmp = Float64(2.0 / Float64((Float64(t / cbrt(Float64(l * Float64(l / sin(k))))) ^ 3.0) * Float64(2.0 * k))); else tmp = Float64(2.0 / Float64(Float64(t * Float64(k / l)) * Float64((sin(k) ^ 2.0) * Float64(k / Float64(l * cos(k)))))); end return tmp end
code[t_, l_, k_] := If[LessEqual[k, 4.2e-39], N[(2.0 / N[(N[Power[N[(t / N[Power[N[(l * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(2.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] * N[(k / N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 4.2 \cdot 10^{-39}:\\
\;\;\;\;\frac{2}{{\left(\frac{t}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3} \cdot \left(2 \cdot k\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(t \cdot \frac{k}{\ell}\right) \cdot \left({\sin k}^{2} \cdot \frac{k}{\ell \cdot \cos k}\right)}\\
\end{array}
\end{array}
if k < 4.19999999999999987e-39Initial program 52.5%
associate-/r/51.8%
add-cube-cbrt51.7%
pow351.7%
cbrt-div51.7%
rem-cbrt-cube61.3%
Applied egg-rr61.3%
distribute-lft-in61.3%
associate-/l*61.3%
associate-/l*66.1%
Applied egg-rr66.1%
distribute-lft-out66.1%
+-commutative66.1%
associate-*r*66.1%
associate-/r/66.1%
associate-+r+66.1%
metadata-eval66.1%
Simplified66.1%
Taylor expanded in k around 0 63.2%
*-commutative63.2%
Simplified63.2%
if 4.19999999999999987e-39 < k Initial program 55.1%
Taylor expanded in t around 0 78.1%
associate-/l*76.6%
unpow276.6%
unpow276.6%
associate-*r*76.6%
times-frac79.5%
Simplified79.5%
Taylor expanded in k around inf 78.1%
associate-/l*76.6%
unpow276.6%
unpow276.6%
associate-*r*76.6%
times-frac79.5%
times-frac90.3%
associate-/r/90.3%
*-commutative90.3%
associate-/r/91.8%
Simplified91.8%
Final simplification70.1%
(FPCore (t l k) :precision binary64 (if (<= k 3.3e-41) (/ 2.0 (* (pow (/ t (cbrt (* l (/ l (sin k))))) 3.0) (* 2.0 k))) (/ 2.0 (* (* (/ t l) (/ (* k k) l)) (* k (/ k (cos k)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 3.3e-41) {
tmp = 2.0 / (pow((t / cbrt((l * (l / sin(k))))), 3.0) * (2.0 * k));
} else {
tmp = 2.0 / (((t / l) * ((k * k) / l)) * (k * (k / cos(k))));
}
return tmp;
}
public static double code(double t, double l, double k) {
double tmp;
if (k <= 3.3e-41) {
tmp = 2.0 / (Math.pow((t / Math.cbrt((l * (l / Math.sin(k))))), 3.0) * (2.0 * k));
} else {
tmp = 2.0 / (((t / l) * ((k * k) / l)) * (k * (k / Math.cos(k))));
}
return tmp;
}
function code(t, l, k) tmp = 0.0 if (k <= 3.3e-41) tmp = Float64(2.0 / Float64((Float64(t / cbrt(Float64(l * Float64(l / sin(k))))) ^ 3.0) * Float64(2.0 * k))); else tmp = Float64(2.0 / Float64(Float64(Float64(t / l) * Float64(Float64(k * k) / l)) * Float64(k * Float64(k / cos(k))))); end return tmp end
code[t_, l_, k_] := If[LessEqual[k, 3.3e-41], N[(2.0 / N[(N[Power[N[(t / N[Power[N[(l * N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(2.0 * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(t / l), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(k * N[(k / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 3.3 \cdot 10^{-41}:\\
\;\;\;\;\frac{2}{{\left(\frac{t}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3} \cdot \left(2 \cdot k\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\frac{t}{\ell} \cdot \frac{k \cdot k}{\ell}\right) \cdot \left(k \cdot \frac{k}{\cos k}\right)}\\
\end{array}
\end{array}
if k < 3.30000000000000024e-41Initial program 52.5%
associate-/r/51.8%
add-cube-cbrt51.7%
pow351.7%
cbrt-div51.7%
rem-cbrt-cube61.3%
Applied egg-rr61.3%
distribute-lft-in61.3%
associate-/l*61.3%
associate-/l*66.1%
Applied egg-rr66.1%
distribute-lft-out66.1%
+-commutative66.1%
associate-*r*66.1%
associate-/r/66.1%
associate-+r+66.1%
metadata-eval66.1%
Simplified66.1%
Taylor expanded in k around 0 63.2%
*-commutative63.2%
Simplified63.2%
if 3.30000000000000024e-41 < k Initial program 55.1%
Taylor expanded in t around 0 78.1%
times-frac76.6%
unpow276.6%
*-commutative76.6%
unpow276.6%
Simplified76.6%
Taylor expanded in k around 0 63.2%
*-commutative63.2%
unpow263.2%
times-frac63.5%
unpow263.5%
Simplified63.5%
pow163.5%
*-commutative63.5%
frac-times63.2%
associate-/l*63.2%
Applied egg-rr63.2%
unpow163.2%
times-frac63.5%
associate-/r/63.5%
Simplified63.5%
Final simplification63.3%
(FPCore (t l k)
:precision binary64
(if (<= t -4.8e+206)
(* l (/ l (* (pow t 3.0) (* k k))))
(if (<= t 4e-67)
(/ 2.0 (* (* (/ t l) (/ (* k k) l)) (* k (/ k (cos k)))))
(if (<= t 9.6e+111)
(/
(* 2.0 (/ (* (/ l k) (/ l k)) (pow t 3.0)))
(+ 2.0 (pow (/ k t) 2.0)))
(if (<= t 1.4e+181)
(/ 2.0 (/ (* k k) (log (exp (* (/ l t) (/ l (* k k)))))))
(/ 2.0 (* 2.0 (* (* k (/ k l)) (/ (pow t 3.0) l)))))))))
double code(double t, double l, double k) {
double tmp;
if (t <= -4.8e+206) {
tmp = l * (l / (pow(t, 3.0) * (k * k)));
} else if (t <= 4e-67) {
tmp = 2.0 / (((t / l) * ((k * k) / l)) * (k * (k / cos(k))));
} else if (t <= 9.6e+111) {
tmp = (2.0 * (((l / k) * (l / k)) / pow(t, 3.0))) / (2.0 + pow((k / t), 2.0));
} else if (t <= 1.4e+181) {
tmp = 2.0 / ((k * k) / log(exp(((l / t) * (l / (k * k))))));
} else {
tmp = 2.0 / (2.0 * ((k * (k / l)) * (pow(t, 3.0) / l)));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t <= (-4.8d+206)) then
tmp = l * (l / ((t ** 3.0d0) * (k * k)))
else if (t <= 4d-67) then
tmp = 2.0d0 / (((t / l) * ((k * k) / l)) * (k * (k / cos(k))))
else if (t <= 9.6d+111) then
tmp = (2.0d0 * (((l / k) * (l / k)) / (t ** 3.0d0))) / (2.0d0 + ((k / t) ** 2.0d0))
else if (t <= 1.4d+181) then
tmp = 2.0d0 / ((k * k) / log(exp(((l / t) * (l / (k * k))))))
else
tmp = 2.0d0 / (2.0d0 * ((k * (k / l)) * ((t ** 3.0d0) / l)))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (t <= -4.8e+206) {
tmp = l * (l / (Math.pow(t, 3.0) * (k * k)));
} else if (t <= 4e-67) {
tmp = 2.0 / (((t / l) * ((k * k) / l)) * (k * (k / Math.cos(k))));
} else if (t <= 9.6e+111) {
tmp = (2.0 * (((l / k) * (l / k)) / Math.pow(t, 3.0))) / (2.0 + Math.pow((k / t), 2.0));
} else if (t <= 1.4e+181) {
tmp = 2.0 / ((k * k) / Math.log(Math.exp(((l / t) * (l / (k * k))))));
} else {
tmp = 2.0 / (2.0 * ((k * (k / l)) * (Math.pow(t, 3.0) / l)));
}
return tmp;
}
def code(t, l, k): tmp = 0 if t <= -4.8e+206: tmp = l * (l / (math.pow(t, 3.0) * (k * k))) elif t <= 4e-67: tmp = 2.0 / (((t / l) * ((k * k) / l)) * (k * (k / math.cos(k)))) elif t <= 9.6e+111: tmp = (2.0 * (((l / k) * (l / k)) / math.pow(t, 3.0))) / (2.0 + math.pow((k / t), 2.0)) elif t <= 1.4e+181: tmp = 2.0 / ((k * k) / math.log(math.exp(((l / t) * (l / (k * k)))))) else: tmp = 2.0 / (2.0 * ((k * (k / l)) * (math.pow(t, 3.0) / l))) return tmp
function code(t, l, k) tmp = 0.0 if (t <= -4.8e+206) tmp = Float64(l * Float64(l / Float64((t ^ 3.0) * Float64(k * k)))); elseif (t <= 4e-67) tmp = Float64(2.0 / Float64(Float64(Float64(t / l) * Float64(Float64(k * k) / l)) * Float64(k * Float64(k / cos(k))))); elseif (t <= 9.6e+111) tmp = Float64(Float64(2.0 * Float64(Float64(Float64(l / k) * Float64(l / k)) / (t ^ 3.0))) / Float64(2.0 + (Float64(k / t) ^ 2.0))); elseif (t <= 1.4e+181) tmp = Float64(2.0 / Float64(Float64(k * k) / log(exp(Float64(Float64(l / t) * Float64(l / Float64(k * k))))))); else tmp = Float64(2.0 / Float64(2.0 * Float64(Float64(k * Float64(k / l)) * Float64((t ^ 3.0) / l)))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (t <= -4.8e+206) tmp = l * (l / ((t ^ 3.0) * (k * k))); elseif (t <= 4e-67) tmp = 2.0 / (((t / l) * ((k * k) / l)) * (k * (k / cos(k)))); elseif (t <= 9.6e+111) tmp = (2.0 * (((l / k) * (l / k)) / (t ^ 3.0))) / (2.0 + ((k / t) ^ 2.0)); elseif (t <= 1.4e+181) tmp = 2.0 / ((k * k) / log(exp(((l / t) * (l / (k * k)))))); else tmp = 2.0 / (2.0 * ((k * (k / l)) * ((t ^ 3.0) / l))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[t, -4.8e+206], N[(l * N[(l / N[(N[Power[t, 3.0], $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 4e-67], N[(2.0 / N[(N[(N[(t / l), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(k * N[(k / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 9.6e+111], N[(N[(2.0 * N[(N[(N[(l / k), $MachinePrecision] * N[(l / k), $MachinePrecision]), $MachinePrecision] / N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.4e+181], N[(2.0 / N[(N[(k * k), $MachinePrecision] / N[Log[N[Exp[N[(N[(l / t), $MachinePrecision] * N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(2.0 * N[(N[(k * N[(k / l), $MachinePrecision]), $MachinePrecision] * N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -4.8 \cdot 10^{+206}:\\
\;\;\;\;\ell \cdot \frac{\ell}{{t}^{3} \cdot \left(k \cdot k\right)}\\
\mathbf{elif}\;t \leq 4 \cdot 10^{-67}:\\
\;\;\;\;\frac{2}{\left(\frac{t}{\ell} \cdot \frac{k \cdot k}{\ell}\right) \cdot \left(k \cdot \frac{k}{\cos k}\right)}\\
\mathbf{elif}\;t \leq 9.6 \cdot 10^{+111}:\\
\;\;\;\;\frac{2 \cdot \frac{\frac{\ell}{k} \cdot \frac{\ell}{k}}{{t}^{3}}}{2 + {\left(\frac{k}{t}\right)}^{2}}\\
\mathbf{elif}\;t \leq 1.4 \cdot 10^{+181}:\\
\;\;\;\;\frac{2}{\frac{k \cdot k}{\log \left(e^{\frac{\ell}{t} \cdot \frac{\ell}{k \cdot k}}\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\left(k \cdot \frac{k}{\ell}\right) \cdot \frac{{t}^{3}}{\ell}\right)}\\
\end{array}
\end{array}
if t < -4.7999999999999999e206Initial program 42.6%
associate-/r*42.6%
associate-*l*35.8%
sqr-neg35.8%
associate-*l*42.6%
*-commutative42.6%
sqr-neg42.6%
associate-*l/42.6%
associate-*r/42.6%
associate-/r/42.6%
Simplified42.6%
Taylor expanded in k around 0 35.8%
unpow235.8%
*-commutative35.8%
unpow235.8%
Simplified35.8%
*-un-lft-identity35.8%
associate-/l*59.6%
*-commutative59.6%
Applied egg-rr59.6%
associate-/r/59.6%
Applied egg-rr59.6%
if -4.7999999999999999e206 < t < 3.99999999999999977e-67Initial program 46.4%
Taylor expanded in t around 0 70.3%
times-frac70.4%
unpow270.4%
*-commutative70.4%
unpow270.4%
Simplified70.4%
Taylor expanded in k around 0 59.1%
*-commutative59.1%
unpow259.1%
times-frac68.1%
unpow268.1%
Simplified68.1%
pow168.1%
*-commutative68.1%
frac-times59.1%
associate-/l*59.1%
Applied egg-rr59.1%
unpow159.1%
times-frac68.1%
associate-/r/68.1%
Simplified68.1%
if 3.99999999999999977e-67 < t < 9.60000000000000023e111Initial program 80.2%
associate-/r*80.3%
associate-*l*73.6%
sqr-neg73.6%
associate-*l*80.3%
*-commutative80.3%
sqr-neg80.3%
associate-*l/83.5%
associate-*r/83.7%
associate-/r/83.7%
Simplified83.7%
Taylor expanded in k around 0 68.3%
associate-/r*65.1%
unpow265.1%
unpow265.1%
times-frac78.5%
Simplified78.5%
if 9.60000000000000023e111 < t < 1.39999999999999992e181Initial program 40.0%
Taylor expanded in t around 0 57.1%
associate-/l*57.1%
unpow257.1%
unpow257.1%
associate-*r*57.1%
times-frac57.4%
Simplified57.4%
Taylor expanded in k around 0 57.4%
unpow257.4%
Simplified57.4%
add-log-exp62.8%
Applied egg-rr62.8%
if 1.39999999999999992e181 < t Initial program 73.1%
associate-/r/72.8%
add-cube-cbrt72.8%
pow372.8%
cbrt-div72.8%
rem-cbrt-cube73.5%
Applied egg-rr73.5%
Taylor expanded in k around 0 63.6%
unpow263.6%
unpow263.6%
times-frac66.9%
associate-*l/76.0%
Simplified76.0%
Final simplification69.4%
(FPCore (t l k)
:precision binary64
(if (<= t -1.76e+208)
(* l (/ l (* (pow t 3.0) (* k k))))
(if (<= t 8.6e-35)
(/ 2.0 (* (* (/ t l) (/ (* k k) l)) (* k (/ k (cos k)))))
(if (<= t 1.6e+112)
(/ (* (/ l k) (/ l k)) (pow t 3.0))
(if (<= t 1.4e+181)
(/ 2.0 (/ (* k k) (log (exp (* (/ l t) (/ l (* k k)))))))
(/ 2.0 (* 2.0 (* (* k (/ k l)) (/ (pow t 3.0) l)))))))))
double code(double t, double l, double k) {
double tmp;
if (t <= -1.76e+208) {
tmp = l * (l / (pow(t, 3.0) * (k * k)));
} else if (t <= 8.6e-35) {
tmp = 2.0 / (((t / l) * ((k * k) / l)) * (k * (k / cos(k))));
} else if (t <= 1.6e+112) {
tmp = ((l / k) * (l / k)) / pow(t, 3.0);
} else if (t <= 1.4e+181) {
tmp = 2.0 / ((k * k) / log(exp(((l / t) * (l / (k * k))))));
} else {
tmp = 2.0 / (2.0 * ((k * (k / l)) * (pow(t, 3.0) / l)));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t <= (-1.76d+208)) then
tmp = l * (l / ((t ** 3.0d0) * (k * k)))
else if (t <= 8.6d-35) then
tmp = 2.0d0 / (((t / l) * ((k * k) / l)) * (k * (k / cos(k))))
else if (t <= 1.6d+112) then
tmp = ((l / k) * (l / k)) / (t ** 3.0d0)
else if (t <= 1.4d+181) then
tmp = 2.0d0 / ((k * k) / log(exp(((l / t) * (l / (k * k))))))
else
tmp = 2.0d0 / (2.0d0 * ((k * (k / l)) * ((t ** 3.0d0) / l)))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (t <= -1.76e+208) {
tmp = l * (l / (Math.pow(t, 3.0) * (k * k)));
} else if (t <= 8.6e-35) {
tmp = 2.0 / (((t / l) * ((k * k) / l)) * (k * (k / Math.cos(k))));
} else if (t <= 1.6e+112) {
tmp = ((l / k) * (l / k)) / Math.pow(t, 3.0);
} else if (t <= 1.4e+181) {
tmp = 2.0 / ((k * k) / Math.log(Math.exp(((l / t) * (l / (k * k))))));
} else {
tmp = 2.0 / (2.0 * ((k * (k / l)) * (Math.pow(t, 3.0) / l)));
}
return tmp;
}
def code(t, l, k): tmp = 0 if t <= -1.76e+208: tmp = l * (l / (math.pow(t, 3.0) * (k * k))) elif t <= 8.6e-35: tmp = 2.0 / (((t / l) * ((k * k) / l)) * (k * (k / math.cos(k)))) elif t <= 1.6e+112: tmp = ((l / k) * (l / k)) / math.pow(t, 3.0) elif t <= 1.4e+181: tmp = 2.0 / ((k * k) / math.log(math.exp(((l / t) * (l / (k * k)))))) else: tmp = 2.0 / (2.0 * ((k * (k / l)) * (math.pow(t, 3.0) / l))) return tmp
function code(t, l, k) tmp = 0.0 if (t <= -1.76e+208) tmp = Float64(l * Float64(l / Float64((t ^ 3.0) * Float64(k * k)))); elseif (t <= 8.6e-35) tmp = Float64(2.0 / Float64(Float64(Float64(t / l) * Float64(Float64(k * k) / l)) * Float64(k * Float64(k / cos(k))))); elseif (t <= 1.6e+112) tmp = Float64(Float64(Float64(l / k) * Float64(l / k)) / (t ^ 3.0)); elseif (t <= 1.4e+181) tmp = Float64(2.0 / Float64(Float64(k * k) / log(exp(Float64(Float64(l / t) * Float64(l / Float64(k * k))))))); else tmp = Float64(2.0 / Float64(2.0 * Float64(Float64(k * Float64(k / l)) * Float64((t ^ 3.0) / l)))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (t <= -1.76e+208) tmp = l * (l / ((t ^ 3.0) * (k * k))); elseif (t <= 8.6e-35) tmp = 2.0 / (((t / l) * ((k * k) / l)) * (k * (k / cos(k)))); elseif (t <= 1.6e+112) tmp = ((l / k) * (l / k)) / (t ^ 3.0); elseif (t <= 1.4e+181) tmp = 2.0 / ((k * k) / log(exp(((l / t) * (l / (k * k)))))); else tmp = 2.0 / (2.0 * ((k * (k / l)) * ((t ^ 3.0) / l))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[t, -1.76e+208], N[(l * N[(l / N[(N[Power[t, 3.0], $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 8.6e-35], N[(2.0 / N[(N[(N[(t / l), $MachinePrecision] * N[(N[(k * k), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision] * N[(k * N[(k / N[Cos[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.6e+112], N[(N[(N[(l / k), $MachinePrecision] * N[(l / k), $MachinePrecision]), $MachinePrecision] / N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.4e+181], N[(2.0 / N[(N[(k * k), $MachinePrecision] / N[Log[N[Exp[N[(N[(l / t), $MachinePrecision] * N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(2.0 * N[(N[(k * N[(k / l), $MachinePrecision]), $MachinePrecision] * N[(N[Power[t, 3.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.76 \cdot 10^{+208}:\\
\;\;\;\;\ell \cdot \frac{\ell}{{t}^{3} \cdot \left(k \cdot k\right)}\\
\mathbf{elif}\;t \leq 8.6 \cdot 10^{-35}:\\
\;\;\;\;\frac{2}{\left(\frac{t}{\ell} \cdot \frac{k \cdot k}{\ell}\right) \cdot \left(k \cdot \frac{k}{\cos k}\right)}\\
\mathbf{elif}\;t \leq 1.6 \cdot 10^{+112}:\\
\;\;\;\;\frac{\frac{\ell}{k} \cdot \frac{\ell}{k}}{{t}^{3}}\\
\mathbf{elif}\;t \leq 1.4 \cdot 10^{+181}:\\
\;\;\;\;\frac{2}{\frac{k \cdot k}{\log \left(e^{\frac{\ell}{t} \cdot \frac{\ell}{k \cdot k}}\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{2 \cdot \left(\left(k \cdot \frac{k}{\ell}\right) \cdot \frac{{t}^{3}}{\ell}\right)}\\
\end{array}
\end{array}
if t < -1.76000000000000007e208Initial program 42.6%
associate-/r*42.6%
associate-*l*35.8%
sqr-neg35.8%
associate-*l*42.6%
*-commutative42.6%
sqr-neg42.6%
associate-*l/42.6%
associate-*r/42.6%
associate-/r/42.6%
Simplified42.6%
Taylor expanded in k around 0 35.8%
unpow235.8%
*-commutative35.8%
unpow235.8%
Simplified35.8%
*-un-lft-identity35.8%
associate-/l*59.6%
*-commutative59.6%
Applied egg-rr59.6%
associate-/r/59.6%
Applied egg-rr59.6%
if -1.76000000000000007e208 < t < 8.6000000000000004e-35Initial program 48.4%
Taylor expanded in t around 0 71.4%
times-frac71.5%
unpow271.5%
*-commutative71.5%
unpow271.5%
Simplified71.5%
Taylor expanded in k around 0 60.2%
*-commutative60.2%
unpow260.2%
times-frac68.8%
unpow268.8%
Simplified68.8%
pow168.8%
*-commutative68.8%
frac-times60.2%
associate-/l*60.2%
Applied egg-rr60.2%
unpow160.2%
times-frac68.8%
associate-/r/68.8%
Simplified68.8%
if 8.6000000000000004e-35 < t < 1.59999999999999993e112Initial program 75.4%
associate-/r*75.4%
associate-*l*67.0%
sqr-neg67.0%
associate-*l*75.4%
*-commutative75.4%
sqr-neg75.4%
associate-*l/79.5%
associate-*r/79.7%
associate-/r/79.6%
Simplified79.6%
Taylor expanded in k around 0 63.5%
associate-/r*59.5%
unpow259.5%
unpow259.5%
times-frac76.1%
Simplified76.1%
if 1.59999999999999993e112 < t < 1.39999999999999992e181Initial program 40.0%
Taylor expanded in t around 0 57.1%
associate-/l*57.1%
unpow257.1%
unpow257.1%
associate-*r*57.1%
times-frac57.4%
Simplified57.4%
Taylor expanded in k around 0 57.4%
unpow257.4%
Simplified57.4%
add-log-exp62.8%
Applied egg-rr62.8%
if 1.39999999999999992e181 < t Initial program 73.1%
associate-/r/72.8%
add-cube-cbrt72.8%
pow372.8%
cbrt-div72.8%
rem-cbrt-cube73.5%
Applied egg-rr73.5%
Taylor expanded in k around 0 63.6%
unpow263.6%
unpow263.6%
times-frac66.9%
associate-*l/76.0%
Simplified76.0%
Final simplification69.4%
(FPCore (t l k)
:precision binary64
(let* ((t_1 (/ l (* k k))))
(if (<= l 4.3e-50)
(/ 2.0 (/ (* k k) (* (/ l t) (+ t_1 (* l -0.16666666666666666)))))
(if (<= l 5.8e+207)
(* t_1 (/ l (pow t 3.0)))
(* (/ 2.0 (* k k)) (* (/ l t) t_1))))))
double code(double t, double l, double k) {
double t_1 = l / (k * k);
double tmp;
if (l <= 4.3e-50) {
tmp = 2.0 / ((k * k) / ((l / t) * (t_1 + (l * -0.16666666666666666))));
} else if (l <= 5.8e+207) {
tmp = t_1 * (l / pow(t, 3.0));
} else {
tmp = (2.0 / (k * k)) * ((l / t) * t_1);
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = l / (k * k)
if (l <= 4.3d-50) then
tmp = 2.0d0 / ((k * k) / ((l / t) * (t_1 + (l * (-0.16666666666666666d0)))))
else if (l <= 5.8d+207) then
tmp = t_1 * (l / (t ** 3.0d0))
else
tmp = (2.0d0 / (k * k)) * ((l / t) * t_1)
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double t_1 = l / (k * k);
double tmp;
if (l <= 4.3e-50) {
tmp = 2.0 / ((k * k) / ((l / t) * (t_1 + (l * -0.16666666666666666))));
} else if (l <= 5.8e+207) {
tmp = t_1 * (l / Math.pow(t, 3.0));
} else {
tmp = (2.0 / (k * k)) * ((l / t) * t_1);
}
return tmp;
}
def code(t, l, k): t_1 = l / (k * k) tmp = 0 if l <= 4.3e-50: tmp = 2.0 / ((k * k) / ((l / t) * (t_1 + (l * -0.16666666666666666)))) elif l <= 5.8e+207: tmp = t_1 * (l / math.pow(t, 3.0)) else: tmp = (2.0 / (k * k)) * ((l / t) * t_1) return tmp
function code(t, l, k) t_1 = Float64(l / Float64(k * k)) tmp = 0.0 if (l <= 4.3e-50) tmp = Float64(2.0 / Float64(Float64(k * k) / Float64(Float64(l / t) * Float64(t_1 + Float64(l * -0.16666666666666666))))); elseif (l <= 5.8e+207) tmp = Float64(t_1 * Float64(l / (t ^ 3.0))); else tmp = Float64(Float64(2.0 / Float64(k * k)) * Float64(Float64(l / t) * t_1)); end return tmp end
function tmp_2 = code(t, l, k) t_1 = l / (k * k); tmp = 0.0; if (l <= 4.3e-50) tmp = 2.0 / ((k * k) / ((l / t) * (t_1 + (l * -0.16666666666666666)))); elseif (l <= 5.8e+207) tmp = t_1 * (l / (t ^ 3.0)); else tmp = (2.0 / (k * k)) * ((l / t) * t_1); end tmp_2 = tmp; end
code[t_, l_, k_] := Block[{t$95$1 = N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[l, 4.3e-50], N[(2.0 / N[(N[(k * k), $MachinePrecision] / N[(N[(l / t), $MachinePrecision] * N[(t$95$1 + N[(l * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[l, 5.8e+207], N[(t$95$1 * N[(l / N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(l / t), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{\ell}{k \cdot k}\\
\mathbf{if}\;\ell \leq 4.3 \cdot 10^{-50}:\\
\;\;\;\;\frac{2}{\frac{k \cdot k}{\frac{\ell}{t} \cdot \left(t_1 + \ell \cdot -0.16666666666666666\right)}}\\
\mathbf{elif}\;\ell \leq 5.8 \cdot 10^{+207}:\\
\;\;\;\;t_1 \cdot \frac{\ell}{{t}^{3}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{k \cdot k} \cdot \left(\frac{\ell}{t} \cdot t_1\right)\\
\end{array}
\end{array}
if l < 4.29999999999999997e-50Initial program 55.4%
Taylor expanded in t around 0 60.9%
associate-/l*62.1%
unpow262.1%
unpow262.1%
associate-*r*62.1%
times-frac70.1%
Simplified70.1%
Taylor expanded in k around 0 64.3%
associate--l+64.3%
unpow264.3%
distribute-rgt-out--64.3%
metadata-eval64.3%
Simplified64.3%
if 4.29999999999999997e-50 < l < 5.79999999999999994e207Initial program 52.6%
associate-/r*50.7%
associate-*l*50.2%
sqr-neg50.2%
associate-*l*50.7%
*-commutative50.7%
sqr-neg50.7%
associate-*l/50.7%
associate-*r/50.7%
associate-/r/50.8%
Simplified50.8%
Taylor expanded in k around 0 47.3%
unpow247.3%
*-commutative47.3%
unpow247.3%
Simplified47.3%
times-frac47.5%
Applied egg-rr47.5%
if 5.79999999999999994e207 < l Initial program 34.2%
Taylor expanded in t around 0 67.8%
associate-/l*67.8%
unpow267.8%
unpow267.8%
associate-*r*67.8%
times-frac72.8%
Simplified72.8%
Taylor expanded in k around 0 72.6%
unpow272.6%
Simplified72.6%
associate-/r/72.6%
Applied egg-rr72.6%
Final simplification61.7%
(FPCore (t l k)
:precision binary64
(if (<= k 1.15e-177)
(/ (* (/ l k) (/ l k)) (pow t 3.0))
(if (<= k 5.1e+95)
(* l (/ l (* (pow t 3.0) (* k k))))
(/ (- l) (* k (* t (/ k l)))))))
double code(double t, double l, double k) {
double tmp;
if (k <= 1.15e-177) {
tmp = ((l / k) * (l / k)) / pow(t, 3.0);
} else if (k <= 5.1e+95) {
tmp = l * (l / (pow(t, 3.0) * (k * k)));
} else {
tmp = -l / (k * (t * (k / l)));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (k <= 1.15d-177) then
tmp = ((l / k) * (l / k)) / (t ** 3.0d0)
else if (k <= 5.1d+95) then
tmp = l * (l / ((t ** 3.0d0) * (k * k)))
else
tmp = -l / (k * (t * (k / l)))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (k <= 1.15e-177) {
tmp = ((l / k) * (l / k)) / Math.pow(t, 3.0);
} else if (k <= 5.1e+95) {
tmp = l * (l / (Math.pow(t, 3.0) * (k * k)));
} else {
tmp = -l / (k * (t * (k / l)));
}
return tmp;
}
def code(t, l, k): tmp = 0 if k <= 1.15e-177: tmp = ((l / k) * (l / k)) / math.pow(t, 3.0) elif k <= 5.1e+95: tmp = l * (l / (math.pow(t, 3.0) * (k * k))) else: tmp = -l / (k * (t * (k / l))) return tmp
function code(t, l, k) tmp = 0.0 if (k <= 1.15e-177) tmp = Float64(Float64(Float64(l / k) * Float64(l / k)) / (t ^ 3.0)); elseif (k <= 5.1e+95) tmp = Float64(l * Float64(l / Float64((t ^ 3.0) * Float64(k * k)))); else tmp = Float64(Float64(-l) / Float64(k * Float64(t * Float64(k / l)))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (k <= 1.15e-177) tmp = ((l / k) * (l / k)) / (t ^ 3.0); elseif (k <= 5.1e+95) tmp = l * (l / ((t ^ 3.0) * (k * k))); else tmp = -l / (k * (t * (k / l))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[k, 1.15e-177], N[(N[(N[(l / k), $MachinePrecision] * N[(l / k), $MachinePrecision]), $MachinePrecision] / N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 5.1e+95], N[(l * N[(l / N[(N[Power[t, 3.0], $MachinePrecision] * N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[((-l) / N[(k * N[(t * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;k \leq 1.15 \cdot 10^{-177}:\\
\;\;\;\;\frac{\frac{\ell}{k} \cdot \frac{\ell}{k}}{{t}^{3}}\\
\mathbf{elif}\;k \leq 5.1 \cdot 10^{+95}:\\
\;\;\;\;\ell \cdot \frac{\ell}{{t}^{3} \cdot \left(k \cdot k\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{-\ell}{k \cdot \left(t \cdot \frac{k}{\ell}\right)}\\
\end{array}
\end{array}
if k < 1.15000000000000011e-177Initial program 52.2%
associate-/r*51.6%
associate-*l*47.8%
sqr-neg47.8%
associate-*l*51.6%
*-commutative51.6%
sqr-neg51.6%
associate-*l/52.1%
associate-*r/52.2%
associate-/r/52.2%
Simplified52.2%
Taylor expanded in k around 0 49.8%
associate-/r*49.1%
unpow249.1%
unpow249.1%
times-frac61.6%
Simplified61.6%
if 1.15000000000000011e-177 < k < 5.10000000000000003e95Initial program 59.2%
associate-/r*59.3%
associate-*l*59.1%
sqr-neg59.1%
associate-*l*59.3%
*-commutative59.3%
sqr-neg59.3%
associate-*l/59.3%
associate-*r/59.3%
associate-/r/57.4%
Simplified57.4%
Taylor expanded in k around 0 61.2%
unpow261.2%
*-commutative61.2%
unpow261.2%
Simplified61.2%
*-un-lft-identity61.2%
associate-/l*65.6%
*-commutative65.6%
Applied egg-rr65.6%
associate-/r/65.6%
Applied egg-rr65.6%
if 5.10000000000000003e95 < k Initial program 49.0%
Taylor expanded in t around 0 78.8%
times-frac78.8%
unpow278.8%
*-commutative78.8%
unpow278.8%
Simplified78.8%
Taylor expanded in k around 0 58.6%
*-commutative58.6%
unpow258.6%
times-frac58.9%
unpow258.9%
Simplified58.9%
Taylor expanded in k around 0 59.3%
fma-def59.3%
unpow259.3%
unpow259.3%
associate-/l*59.3%
associate-*r/59.3%
unpow259.3%
*-commutative59.3%
Simplified59.3%
Taylor expanded in k around inf 59.6%
associate-*r/59.6%
times-frac59.4%
unpow259.4%
associate-/l*59.5%
times-frac60.3%
unpow260.3%
associate-*r*60.9%
mul-1-neg60.9%
associate-*r/60.9%
associate-*l/60.9%
*-commutative60.9%
Simplified60.9%
Final simplification62.3%
(FPCore (t l k) :precision binary64 (if (<= t 5.1e-36) (/ 2.0 (/ (* k k) (* (/ l t) (/ l (* k k))))) (/ (* l l) (* k (* k (pow t 3.0))))))
double code(double t, double l, double k) {
double tmp;
if (t <= 5.1e-36) {
tmp = 2.0 / ((k * k) / ((l / t) * (l / (k * k))));
} else {
tmp = (l * l) / (k * (k * pow(t, 3.0)));
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: tmp
if (t <= 5.1d-36) then
tmp = 2.0d0 / ((k * k) / ((l / t) * (l / (k * k))))
else
tmp = (l * l) / (k * (k * (t ** 3.0d0)))
end if
code = tmp
end function
public static double code(double t, double l, double k) {
double tmp;
if (t <= 5.1e-36) {
tmp = 2.0 / ((k * k) / ((l / t) * (l / (k * k))));
} else {
tmp = (l * l) / (k * (k * Math.pow(t, 3.0)));
}
return tmp;
}
def code(t, l, k): tmp = 0 if t <= 5.1e-36: tmp = 2.0 / ((k * k) / ((l / t) * (l / (k * k)))) else: tmp = (l * l) / (k * (k * math.pow(t, 3.0))) return tmp
function code(t, l, k) tmp = 0.0 if (t <= 5.1e-36) tmp = Float64(2.0 / Float64(Float64(k * k) / Float64(Float64(l / t) * Float64(l / Float64(k * k))))); else tmp = Float64(Float64(l * l) / Float64(k * Float64(k * (t ^ 3.0)))); end return tmp end
function tmp_2 = code(t, l, k) tmp = 0.0; if (t <= 5.1e-36) tmp = 2.0 / ((k * k) / ((l / t) * (l / (k * k)))); else tmp = (l * l) / (k * (k * (t ^ 3.0))); end tmp_2 = tmp; end
code[t_, l_, k_] := If[LessEqual[t, 5.1e-36], N[(2.0 / N[(N[(k * k), $MachinePrecision] / N[(N[(l / t), $MachinePrecision] * N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(l * l), $MachinePrecision] / N[(k * N[(k * N[Power[t, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 5.1 \cdot 10^{-36}:\\
\;\;\;\;\frac{2}{\frac{k \cdot k}{\frac{\ell}{t} \cdot \frac{\ell}{k \cdot k}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\ell \cdot \ell}{k \cdot \left(k \cdot {t}^{3}\right)}\\
\end{array}
\end{array}
if t < 5.09999999999999973e-36Initial program 47.8%
Taylor expanded in t around 0 66.5%
associate-/l*67.1%
unpow267.1%
unpow267.1%
associate-*r*67.1%
times-frac75.9%
Simplified75.9%
Taylor expanded in k around 0 64.4%
unpow264.4%
Simplified64.4%
if 5.09999999999999973e-36 < t Initial program 65.9%
associate-/r*65.9%
associate-*l*58.9%
sqr-neg58.9%
associate-*l*65.9%
*-commutative65.9%
sqr-neg65.9%
associate-*l/67.2%
associate-*r/67.3%
associate-/r/67.2%
Simplified67.2%
Taylor expanded in k around 0 57.7%
unpow257.7%
*-commutative57.7%
unpow257.7%
Simplified57.7%
Taylor expanded in l around 0 57.7%
unpow257.7%
unpow257.7%
associate-*l*66.1%
Simplified66.1%
Final simplification64.9%
(FPCore (t l k) :precision binary64 (* (/ 2.0 (* k k)) (* (/ l t) (/ l (* k k)))))
double code(double t, double l, double k) {
return (2.0 / (k * k)) * ((l / t) * (l / (k * k)));
}
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 / (k * k)) * ((l / t) * (l / (k * k)))
end function
public static double code(double t, double l, double k) {
return (2.0 / (k * k)) * ((l / t) * (l / (k * k)));
}
def code(t, l, k): return (2.0 / (k * k)) * ((l / t) * (l / (k * k)))
function code(t, l, k) return Float64(Float64(2.0 / Float64(k * k)) * Float64(Float64(l / t) * Float64(l / Float64(k * k)))) end
function tmp = code(t, l, k) tmp = (2.0 / (k * k)) * ((l / t) * (l / (k * k))); end
code[t_, l_, k_] := N[(N[(2.0 / N[(k * k), $MachinePrecision]), $MachinePrecision] * N[(N[(l / t), $MachinePrecision] * N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{k \cdot k} \cdot \left(\frac{\ell}{t} \cdot \frac{\ell}{k \cdot k}\right)
\end{array}
Initial program 53.1%
Taylor expanded in t around 0 63.0%
associate-/l*63.1%
unpow263.1%
unpow263.1%
associate-*r*63.1%
times-frac69.3%
Simplified69.3%
Taylor expanded in k around 0 60.9%
unpow260.9%
Simplified60.9%
associate-/r/60.5%
Applied egg-rr60.5%
Final simplification60.5%
(FPCore (t l k) :precision binary64 (* (/ l (* (* k k) (/ t l))) (/ (/ 2.0 k) k)))
double code(double t, double l, double k) {
return (l / ((k * k) * (t / l))) * ((2.0 / k) / k);
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = (l / ((k * k) * (t / l))) * ((2.0d0 / k) / k)
end function
public static double code(double t, double l, double k) {
return (l / ((k * k) * (t / l))) * ((2.0 / k) / k);
}
def code(t, l, k): return (l / ((k * k) * (t / l))) * ((2.0 / k) / k)
function code(t, l, k) return Float64(Float64(l / Float64(Float64(k * k) * Float64(t / l))) * Float64(Float64(2.0 / k) / k)) end
function tmp = code(t, l, k) tmp = (l / ((k * k) * (t / l))) * ((2.0 / k) / k); end
code[t_, l_, k_] := N[(N[(l / N[(N[(k * k), $MachinePrecision] * N[(t / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(2.0 / k), $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell}{\left(k \cdot k\right) \cdot \frac{t}{\ell}} \cdot \frac{\frac{2}{k}}{k}
\end{array}
Initial program 53.1%
Taylor expanded in t around 0 63.0%
associate-/l*63.1%
unpow263.1%
unpow263.1%
associate-*r*63.1%
times-frac69.3%
Simplified69.3%
Taylor expanded in k around 0 60.9%
unpow260.9%
Simplified60.9%
*-un-lft-identity60.9%
associate-/r/60.5%
Applied egg-rr60.5%
*-lft-identity60.5%
*-commutative60.5%
associate-/r*60.5%
associate-*r/59.8%
associate-/l*60.5%
associate-/r*60.5%
Simplified60.5%
Final simplification60.5%
(FPCore (t l k) :precision binary64 (/ 2.0 (* (/ k (/ l (* k k))) (/ k (/ l t)))))
double code(double t, double l, double k) {
return 2.0 / ((k / (l / (k * k))) * (k / (l / t)));
}
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 / ((k / (l / (k * k))) * (k / (l / t)))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((k / (l / (k * k))) * (k / (l / t)));
}
def code(t, l, k): return 2.0 / ((k / (l / (k * k))) * (k / (l / t)))
function code(t, l, k) return Float64(2.0 / Float64(Float64(k / Float64(l / Float64(k * k))) * Float64(k / Float64(l / t)))) end
function tmp = code(t, l, k) tmp = 2.0 / ((k / (l / (k * k))) * (k / (l / t))); end
code[t_, l_, k_] := N[(2.0 / N[(N[(k / N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(k / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\frac{k}{\frac{\ell}{k \cdot k}} \cdot \frac{k}{\frac{\ell}{t}}}
\end{array}
Initial program 53.1%
Taylor expanded in t around 0 63.0%
associate-/l*63.1%
unpow263.1%
unpow263.1%
associate-*r*63.1%
times-frac69.3%
Simplified69.3%
Taylor expanded in k around 0 60.9%
unpow260.9%
Simplified60.9%
times-frac60.6%
Applied egg-rr60.6%
Final simplification60.6%
(FPCore (t l k) :precision binary64 (/ 2.0 (/ (* k k) (* (/ l t) (/ l (* k k))))))
double code(double t, double l, double k) {
return 2.0 / ((k * k) / ((l / t) * (l / (k * k))));
}
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 / ((k * k) / ((l / t) * (l / (k * k))))
end function
public static double code(double t, double l, double k) {
return 2.0 / ((k * k) / ((l / t) * (l / (k * k))));
}
def code(t, l, k): return 2.0 / ((k * k) / ((l / t) * (l / (k * k))))
function code(t, l, k) return Float64(2.0 / Float64(Float64(k * k) / Float64(Float64(l / t) * Float64(l / Float64(k * k))))) end
function tmp = code(t, l, k) tmp = 2.0 / ((k * k) / ((l / t) * (l / (k * k)))); end
code[t_, l_, k_] := N[(2.0 / N[(N[(k * k), $MachinePrecision] / N[(N[(l / t), $MachinePrecision] * N[(l / N[(k * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{2}{\frac{k \cdot k}{\frac{\ell}{t} \cdot \frac{\ell}{k \cdot k}}}
\end{array}
Initial program 53.1%
Taylor expanded in t around 0 63.0%
associate-/l*63.1%
unpow263.1%
unpow263.1%
associate-*r*63.1%
times-frac69.3%
Simplified69.3%
Taylor expanded in k around 0 60.9%
unpow260.9%
Simplified60.9%
Final simplification60.9%
(FPCore (t l k) :precision binary64 (/ (* l (- l)) (* k (* t k))))
double code(double t, double l, double k) {
return (l * -l) / (k * (t * k));
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = (l * -l) / (k * (t * k))
end function
public static double code(double t, double l, double k) {
return (l * -l) / (k * (t * k));
}
def code(t, l, k): return (l * -l) / (k * (t * k))
function code(t, l, k) return Float64(Float64(l * Float64(-l)) / Float64(k * Float64(t * k))) end
function tmp = code(t, l, k) tmp = (l * -l) / (k * (t * k)); end
code[t_, l_, k_] := N[(N[(l * (-l)), $MachinePrecision] / N[(k * N[(t * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\ell \cdot \left(-\ell\right)}{k \cdot \left(t \cdot k\right)}
\end{array}
Initial program 53.1%
Taylor expanded in t around 0 63.0%
times-frac63.1%
unpow263.1%
*-commutative63.1%
unpow263.1%
Simplified63.1%
Taylor expanded in k around 0 56.2%
*-commutative56.2%
unpow256.2%
times-frac61.6%
unpow261.6%
Simplified61.6%
Taylor expanded in k around 0 30.7%
fma-def30.7%
unpow230.7%
unpow230.7%
associate-/l*30.7%
associate-*r/30.7%
unpow230.7%
*-commutative30.7%
Simplified30.7%
Taylor expanded in k around inf 31.4%
mul-1-neg31.4%
unpow231.4%
unpow231.4%
associate-*l*32.6%
Simplified32.6%
Final simplification32.6%
herbie shell --seed 2023279
(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))))