\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}\;k \le -9.39222433932268557 \cdot 10^{-95}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\sqrt{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1} \cdot {\left(\sqrt[3]{t}\right)}^{3}\right) \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot {\left(\sin k\right)}^{2}\right)\right) \cdot \sqrt{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}\\
\mathbf{elif}\;k \le 2.12710102744016808 \cdot 10^{-292}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{\sqrt[3]{t} \cdot \sqrt[3]{t}}\right)}^{3}}{\sqrt[3]{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}} \cdot \sqrt[3]{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}} \cdot \left(\frac{{\left(\sqrt[3]{\sqrt[3]{t}}\right)}^{3}}{\sqrt[3]{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\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) \cdot \sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k}\right)\right) \cdot \sin k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}\\
\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 ((k <= -9.392224339322686e-95)) {
VAR = ((double) (2.0 / ((double) (((double) (((double) (((double) (((double) sqrt(((double) (((double) (1.0 + ((double) pow(((double) (k / t)), 2.0)))) + 1.0)))) * ((double) pow(((double) cbrt(t)), 3.0)))) * ((double) (((double) (((double) pow(((double) cbrt(t)), 3.0)) / l)) * ((double) pow(((double) sin(k)), 2.0)))))) * ((double) sqrt(((double) (((double) (1.0 + ((double) pow(((double) (k / t)), 2.0)))) + 1.0)))))) / ((double) (((double) (l / ((double) pow(((double) cbrt(t)), 3.0)))) * ((double) cos(k))))))));
} else {
double VAR_1;
if ((k <= 2.127101027440168e-292)) {
VAR_1 = ((double) (2.0 / ((double) (((double) (((double) (((double) (((double) pow(((double) cbrt(((double) (((double) cbrt(t)) * ((double) cbrt(t)))))), 3.0)) / ((double) (((double) cbrt(((double) (l / ((double) pow(((double) cbrt(t)), 3.0)))))) * ((double) cbrt(((double) (l / ((double) pow(((double) cbrt(t)), 3.0)))))))))) * ((double) (((double) (((double) pow(((double) cbrt(((double) cbrt(t)))), 3.0)) / ((double) cbrt(((double) (l / ((double) pow(((double) cbrt(t)), 3.0)))))))) * ((double) (((double) (((double) pow(((double) cbrt(t)), 3.0)) / l)) * ((double) sin(k)))))))) * ((double) tan(k)))) * ((double) (((double) (1.0 + ((double) pow(((double) (k / t)), 2.0)))) + 1.0))))));
} else {
VAR_1 = ((double) (2.0 / ((double) (((double) (((double) (((double) (((double) pow(((double) cbrt(t)), 3.0)) * ((double) (((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) cbrt(((double) (((double) (((double) pow(((double) cbrt(t)), 3.0)) / l)) * ((double) sin(k)))))))))) * ((double) sin(k)))) * ((double) (((double) (1.0 + ((double) pow(((double) (k / t)), 2.0)))) + 1.0)))) / ((double) (((double) (l / ((double) pow(((double) cbrt(t)), 3.0)))) * ((double) cos(k))))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -9.39222433932268557e-95Initial program 32.7
rmApplied add-cube-cbrt32.8
Applied unpow-prod-down32.8
Applied times-frac24.6
Applied associate-*l*24.3
rmApplied unpow-prod-down24.3
Applied associate-/l*19.2
rmApplied tan-quot19.2
Applied associate-*l/19.1
Applied frac-times18.7
Applied associate-*l/16.6
rmApplied add-sqr-sqrt16.6
Applied associate-*r*16.6
Simplified14.5
if -9.39222433932268557e-95 < k < 2.12710102744016808e-292Initial program 36.3
rmApplied add-cube-cbrt36.4
Applied unpow-prod-down36.4
Applied times-frac31.0
Applied associate-*l*23.2
rmApplied unpow-prod-down23.2
Applied associate-/l*16.9
rmApplied add-cube-cbrt16.9
Applied add-cube-cbrt17.0
Applied cbrt-prod17.0
Applied unpow-prod-down17.0
Applied times-frac17.0
Applied associate-*l*13.0
if 2.12710102744016808e-292 < k Initial program 32.3
rmApplied add-cube-cbrt32.4
Applied unpow-prod-down32.4
Applied times-frac25.4
Applied associate-*l*24.0
rmApplied unpow-prod-down24.0
Applied associate-/l*18.2
rmApplied tan-quot18.2
Applied associate-*l/17.6
Applied frac-times16.4
Applied associate-*l/14.8
rmApplied add-cube-cbrt14.8
Final simplification14.5
herbie shell --seed 2020150
(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))))