\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 -1.350083645795687900514049661869243401303 \cdot 10^{154}:\\
\;\;\;\;\frac{\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{3}}{\ell} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \sin k\right) \cdot \tan k}}{{\left(\frac{k}{t}\right)}^{2}}\\
\mathbf{elif}\;\ell \le 1.324056689281534162233822277534588328969 \cdot 10^{154}:\\
\;\;\;\;2 \cdot \left(\frac{{\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\ell}^{2} \cdot \cos k\right)}{{\left(\sin k\right)}^{2}} \cdot {\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\left(\left(\frac{{t}^{\left(\frac{3}{2}\right)}}{\ell} \cdot \frac{{t}^{\left(\frac{3}{2}\right)}}{\ell}\right) \cdot \sin k\right) \cdot \tan k}}{{\left(\frac{k}{t}\right)}^{2}}\\
\end{array}double f(double t, double l, double k) {
double r99829 = 2.0;
double r99830 = t;
double r99831 = 3.0;
double r99832 = pow(r99830, r99831);
double r99833 = l;
double r99834 = r99833 * r99833;
double r99835 = r99832 / r99834;
double r99836 = k;
double r99837 = sin(r99836);
double r99838 = r99835 * r99837;
double r99839 = tan(r99836);
double r99840 = r99838 * r99839;
double r99841 = 1.0;
double r99842 = r99836 / r99830;
double r99843 = pow(r99842, r99829);
double r99844 = r99841 + r99843;
double r99845 = r99844 - r99841;
double r99846 = r99840 * r99845;
double r99847 = r99829 / r99846;
return r99847;
}
double f(double t, double l, double k) {
double r99848 = l;
double r99849 = -1.350083645795688e+154;
bool r99850 = r99848 <= r99849;
double r99851 = 2.0;
double r99852 = t;
double r99853 = cbrt(r99852);
double r99854 = r99853 * r99853;
double r99855 = 3.0;
double r99856 = pow(r99854, r99855);
double r99857 = r99856 / r99848;
double r99858 = pow(r99853, r99855);
double r99859 = r99858 / r99848;
double r99860 = r99857 * r99859;
double r99861 = k;
double r99862 = sin(r99861);
double r99863 = r99860 * r99862;
double r99864 = tan(r99861);
double r99865 = r99863 * r99864;
double r99866 = r99851 / r99865;
double r99867 = r99861 / r99852;
double r99868 = pow(r99867, r99851);
double r99869 = r99866 / r99868;
double r99870 = 1.3240566892815342e+154;
bool r99871 = r99848 <= r99870;
double r99872 = 1.0;
double r99873 = 2.0;
double r99874 = r99851 / r99873;
double r99875 = pow(r99861, r99874);
double r99876 = r99872 / r99875;
double r99877 = 1.0;
double r99878 = pow(r99876, r99877);
double r99879 = pow(r99848, r99873);
double r99880 = cos(r99861);
double r99881 = r99879 * r99880;
double r99882 = r99878 * r99881;
double r99883 = pow(r99862, r99873);
double r99884 = r99882 / r99883;
double r99885 = pow(r99852, r99877);
double r99886 = r99875 * r99885;
double r99887 = r99872 / r99886;
double r99888 = pow(r99887, r99877);
double r99889 = r99884 * r99888;
double r99890 = r99851 * r99889;
double r99891 = r99855 / r99873;
double r99892 = pow(r99852, r99891);
double r99893 = r99892 / r99848;
double r99894 = r99893 * r99893;
double r99895 = r99894 * r99862;
double r99896 = r99895 * r99864;
double r99897 = r99851 / r99896;
double r99898 = r99897 / r99868;
double r99899 = r99871 ? r99890 : r99898;
double r99900 = r99850 ? r99869 : r99899;
return r99900;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < -1.350083645795688e+154Initial program 64.0
Simplified64.0
rmApplied add-cube-cbrt64.0
Applied unpow-prod-down64.0
Applied times-frac48.4
if -1.350083645795688e+154 < l < 1.3240566892815342e+154Initial program 45.6
Simplified36.9
Taylor expanded around inf 15.2
rmApplied sqr-pow15.2
Applied associate-*l*12.7
rmApplied add-cube-cbrt12.7
Applied times-frac12.4
Applied unpow-prod-down12.4
Applied associate-*l*10.5
Simplified10.5
rmApplied associate-*r*10.6
Simplified10.6
if 1.3240566892815342e+154 < l Initial program 64.0
Simplified64.0
rmApplied sqr-pow64.0
Applied times-frac53.7
Final simplification16.9
herbie shell --seed 2019326
(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))))