\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)}\begin{array}{l}
\mathbf{if}\;t \le -1.3078668923076751 \cdot 10^{-190}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}} \cdot \left(\sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k} \cdot \left(\sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k} \cdot \sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k}\right)\right)\right) \cdot \tan k\right) \cdot \left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right)}\\
\mathbf{elif}\;t \le 1.3337804818741451 \cdot 10^{-151}:\\
\;\;\;\;0\\
\mathbf{elif}\;t \le 2.67098617363212593 \cdot 10^{169}:\\
\;\;\;\;\frac{2}{\left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right) \cdot \left(\tan k \cdot \left(\frac{{t}^{\left(\frac{3}{2}\right)}}{\ell} \cdot \left(\sin k \cdot \frac{{t}^{\left(\frac{3}{2}\right)}}{\ell}\right)\right)\right)}\\
\mathbf{elif}\;t \le 1.45585033432551279 \cdot 10^{227}:\\
\;\;\;\;\frac{2}{\frac{\left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right) \cdot \left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\tan k \cdot \left(\sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k} \cdot \left(\sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k} \cdot \sqrt[3]{{\left(\sqrt[3]{t}\right)}^{3} \cdot \sin k}\right)\right)\right)\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \sqrt[3]{\ell}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right) \cdot \left(\tan k \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}} \cdot \left(\frac{{\left(\sqrt[3]{\sqrt{t}}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\sin k \cdot \frac{{\left(\sqrt[3]{\sqrt{t}}\right)}^{3}}{\sqrt[3]{\ell}}\right)\right)\right)\right)}\\
\end{array}double code(double t, double l, double k) {
return ((double) (2.0 / ((double) (((double) (((double) (((double) (((double) pow(t, 3.0)) / ((double) (l * l)))) * ((double) sin(k)))) * ((double) tan(k)))) * ((double) (((double) (1.0 + ((double) pow(((double) (k / t)), 2.0)))) + 1.0))))));
}
double code(double t, double l, double k) {
double VAR;
if ((t <= -1.3078668923076751e-190)) {
VAR = ((double) (2.0 / ((double) (((double) (((double) (((double) (((double) pow(((double) cbrt(t)), 3.0)) / ((double) (l / ((double) pow(((double) cbrt(t)), 3.0)))))) * ((double) (((double) cbrt(((double) (((double) (((double) pow(((double) cbrt(t)), 3.0)) / l)) * ((double) sin(k)))))) * ((double) (((double) cbrt(((double) (((double) (((double) pow(((double) cbrt(t)), 3.0)) / l)) * ((double) sin(k)))))) * ((double) cbrt(((double) (((double) (((double) pow(((double) cbrt(t)), 3.0)) / l)) * ((double) sin(k)))))))))))) * ((double) tan(k)))) * ((double) (1.0 + ((double) (1.0 + ((double) pow(((double) (k / t)), 2.0))))))))));
} else {
double VAR_1;
if ((t <= 1.3337804818741451e-151)) {
VAR_1 = 0.0;
} else {
double VAR_2;
if ((t <= 2.670986173632126e+169)) {
VAR_2 = ((double) (2.0 / ((double) (((double) (1.0 + ((double) (1.0 + ((double) pow(((double) (k / t)), 2.0)))))) * ((double) (((double) tan(k)) * ((double) (((double) (((double) pow(t, ((double) (3.0 / 2.0)))) / l)) * ((double) (((double) sin(k)) * ((double) (((double) pow(t, ((double) (3.0 / 2.0)))) / l))))))))))));
} else {
double VAR_3;
if ((t <= 1.4558503343255128e+227)) {
VAR_3 = ((double) (2.0 / ((double) (((double) (((double) (1.0 + ((double) (1.0 + ((double) pow(((double) (k / t)), 2.0)))))) * ((double) (((double) pow(((double) cbrt(t)), 3.0)) * ((double) (((double) tan(k)) * ((double) (((double) cbrt(((double) (((double) (((double) pow(((double) cbrt(t)), 3.0)) / l)) * ((double) sin(k)))))) * ((double) (((double) cbrt(((double) (((double) (((double) pow(((double) cbrt(t)), 3.0)) / l)) * ((double) sin(k)))))) * ((double) cbrt(((double) (((double) pow(((double) cbrt(t)), 3.0)) * ((double) sin(k)))))))))))))))) / ((double) (((double) (l / ((double) pow(((double) cbrt(t)), 3.0)))) * ((double) cbrt(l))))))));
} else {
VAR_3 = ((double) (2.0 / ((double) (((double) (1.0 + ((double) (1.0 + ((double) pow(((double) (k / t)), 2.0)))))) * ((double) (((double) tan(k)) * ((double) (((double) (((double) pow(((double) cbrt(t)), 3.0)) / ((double) (l / ((double) pow(((double) cbrt(t)), 3.0)))))) * ((double) (((double) (((double) pow(((double) cbrt(((double) sqrt(t)))), 3.0)) / ((double) (((double) cbrt(l)) * ((double) cbrt(l)))))) * ((double) (((double) sin(k)) * ((double) (((double) pow(((double) cbrt(((double) sqrt(t)))), 3.0)) / ((double) cbrt(l))))))))))))))));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -1.3078668923076751e-190Initial program 27.7
rmApplied add-cube-cbrt27.8
Applied unpow-prod-down27.8
Applied times-frac19.7
Applied associate-*l*17.9
rmApplied unpow-prod-down17.9
Applied associate-/l*13.0
rmApplied add-cube-cbrt13.1
if -1.3078668923076751e-190 < t < 1.3337804818741451e-151Initial program 64.0
rmApplied add-cube-cbrt64.0
Applied unpow-prod-down64.0
Applied times-frac63.0
Applied associate-*l*63.0
rmApplied unpow-prod-down63.0
Applied associate-/l*51.8
Taylor expanded around inf 41.6
if 1.3337804818741451e-151 < t < 2.67098617363212593e169Initial program 29.6
rmApplied sqr-pow29.6
Applied times-frac17.5
Applied associate-*l*13.7
if 2.67098617363212593e169 < t < 1.45585033432551279e227Initial program 22.6
rmApplied add-cube-cbrt22.6
Applied unpow-prod-down22.6
Applied times-frac21.6
Applied associate-*l*21.6
rmApplied unpow-prod-down21.6
Applied associate-/l*7.0
rmApplied add-cube-cbrt7.1
rmApplied associate-*l/7.1
Applied cbrt-div7.1
Applied associate-*r/7.0
Applied frac-times5.3
Applied associate-*l/2.0
Applied associate-*l/2.7
Simplified6.2
if 1.45585033432551279e227 < t Initial program 18.6
rmApplied add-cube-cbrt18.6
Applied unpow-prod-down18.6
Applied times-frac13.6
Applied associate-*l*13.6
rmApplied unpow-prod-down13.6
Applied associate-/l*7.1
rmApplied add-cube-cbrt7.1
Applied add-sqr-sqrt7.1
Applied cbrt-prod7.1
Applied unpow-prod-down7.1
Applied times-frac7.1
Applied associate-*l*7.1
Simplified7.1
Final simplification16.5
herbie shell --seed 2020179
(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))))