\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}\;\ell \le 5.4119406433808711 \cdot 10^{-63}:\\
\;\;\;\;\frac{2}{\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)}{\frac{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}{\tan k}} \cdot \frac{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}{\frac{\sqrt[3]{\ell}}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}}\\
\mathbf{elif}\;\ell \le 1297920887946261000:\\
\;\;\;\;\frac{2}{2 \cdot \left({\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{9} \cdot \left({t}^{3} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}}\right) + {\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{9} \cdot \left({\left(\sin k\right)}^{2} \cdot \left({k}^{2} \cdot t\right)\right)}{\cos k \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\left(\sqrt[3]{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)} \cdot \sqrt[3]{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\right) \cdot \sqrt[3]{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}}{\frac{\ell}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}}\\
\end{array}double f(double t, double l, double k) {
double r141731 = 2.0;
double r141732 = t;
double r141733 = 3.0;
double r141734 = pow(r141732, r141733);
double r141735 = l;
double r141736 = r141735 * r141735;
double r141737 = r141734 / r141736;
double r141738 = k;
double r141739 = sin(r141738);
double r141740 = r141737 * r141739;
double r141741 = tan(r141738);
double r141742 = r141740 * r141741;
double r141743 = 1.0;
double r141744 = r141738 / r141732;
double r141745 = pow(r141744, r141731);
double r141746 = r141743 + r141745;
double r141747 = r141746 + r141743;
double r141748 = r141742 * r141747;
double r141749 = r141731 / r141748;
return r141749;
}
double f(double t, double l, double k) {
double r141750 = l;
double r141751 = 5.411940643380871e-63;
bool r141752 = r141750 <= r141751;
double r141753 = 2.0;
double r141754 = t;
double r141755 = cbrt(r141754);
double r141756 = r141755 * r141755;
double r141757 = 3.0;
double r141758 = 2.0;
double r141759 = r141757 / r141758;
double r141760 = pow(r141756, r141759);
double r141761 = pow(r141755, r141757);
double r141762 = r141761 / r141750;
double r141763 = k;
double r141764 = sin(r141763);
double r141765 = r141762 * r141764;
double r141766 = r141760 * r141765;
double r141767 = cbrt(r141750);
double r141768 = r141767 * r141767;
double r141769 = tan(r141763);
double r141770 = r141768 / r141769;
double r141771 = r141766 / r141770;
double r141772 = 1.0;
double r141773 = r141763 / r141754;
double r141774 = pow(r141773, r141753);
double r141775 = r141772 + r141774;
double r141776 = r141775 + r141772;
double r141777 = r141767 / r141760;
double r141778 = r141776 / r141777;
double r141779 = r141771 * r141778;
double r141780 = r141753 / r141779;
double r141781 = 1.297920887946261e+18;
bool r141782 = r141750 <= r141781;
double r141783 = 1.0;
double r141784 = -1.0;
double r141785 = pow(r141784, r141757);
double r141786 = r141783 / r141785;
double r141787 = pow(r141786, r141772);
double r141788 = cbrt(r141784);
double r141789 = 9.0;
double r141790 = pow(r141788, r141789);
double r141791 = 3.0;
double r141792 = pow(r141754, r141791);
double r141793 = pow(r141764, r141758);
double r141794 = r141792 * r141793;
double r141795 = r141790 * r141794;
double r141796 = cos(r141763);
double r141797 = pow(r141750, r141758);
double r141798 = r141796 * r141797;
double r141799 = r141795 / r141798;
double r141800 = r141787 * r141799;
double r141801 = r141753 * r141800;
double r141802 = pow(r141763, r141758);
double r141803 = r141802 * r141754;
double r141804 = r141793 * r141803;
double r141805 = r141790 * r141804;
double r141806 = r141805 / r141798;
double r141807 = r141787 * r141806;
double r141808 = r141801 + r141807;
double r141809 = r141753 / r141808;
double r141810 = r141766 * r141769;
double r141811 = r141810 * r141776;
double r141812 = cbrt(r141811);
double r141813 = r141812 * r141812;
double r141814 = r141813 * r141812;
double r141815 = r141750 / r141760;
double r141816 = r141814 / r141815;
double r141817 = r141753 / r141816;
double r141818 = r141782 ? r141809 : r141817;
double r141819 = r141752 ? r141780 : r141818;
return r141819;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < 5.411940643380871e-63Initial program 29.6
rmApplied add-cube-cbrt29.7
Applied unpow-prod-down29.7
Applied times-frac23.1
Applied associate-*l*20.8
rmApplied sqr-pow20.8
Applied associate-/l*15.3
rmApplied associate-*l/14.5
Applied associate-*l/14.0
Applied associate-*l/12.9
rmApplied *-un-lft-identity12.9
Applied add-cube-cbrt13.0
Applied times-frac13.0
Applied times-frac13.7
Simplified12.6
if 5.411940643380871e-63 < l < 1.297920887946261e+18Initial program 22.9
rmApplied add-cube-cbrt23.3
Applied unpow-prod-down23.3
Applied times-frac21.9
Applied associate-*l*20.7
Taylor expanded around -inf 12.9
if 1.297920887946261e+18 < l Initial program 48.3
rmApplied add-cube-cbrt48.4
Applied unpow-prod-down48.4
Applied times-frac36.1
Applied associate-*l*35.9
rmApplied sqr-pow35.9
Applied associate-/l*25.9
rmApplied associate-*l/24.1
Applied associate-*l/22.6
Applied associate-*l/20.9
rmApplied add-cube-cbrt20.9
Final simplification14.3
herbie shell --seed 2020083
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10+)"
:precision binary64
(/ 2 (* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (+ (+ 1 (pow (/ k t) 2)) 1))))