\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.361561607826163938590433986165123404617 \cdot 10^{74}:\\
\;\;\;\;\frac{\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell}}\right) \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \sin k}}{\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;\ell \le -3.671469527357900756356766403755497030844 \cdot 10^{-125}:\\
\;\;\;\;\frac{2}{-\left(2 \cdot \left({\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{6} \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{t \cdot \left({k}^{2} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}}\right)}\\
\mathbf{elif}\;\ell \le 1.147150844170901392374916916178894499446 \cdot 10^{-119}:\\
\;\;\;\;\frac{2}{\left(\left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell}}\right) \cdot \frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\ell}\right) \cdot \sin k\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}\\
\mathbf{elif}\;\ell \le 2.547304444790666703643914625432659827925 \cdot 10^{119}:\\
\;\;\;\;\frac{2}{-\left(2 \cdot \left({\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{6} \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{t \cdot \left({k}^{2} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell}}\right) \cdot \frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\ell}\right) \cdot \sin k\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}\\
\end{array}double f(double t, double l, double k) {
double r136826 = 2.0;
double r136827 = t;
double r136828 = 3.0;
double r136829 = pow(r136827, r136828);
double r136830 = l;
double r136831 = r136830 * r136830;
double r136832 = r136829 / r136831;
double r136833 = k;
double r136834 = sin(r136833);
double r136835 = r136832 * r136834;
double r136836 = tan(r136833);
double r136837 = r136835 * r136836;
double r136838 = 1.0;
double r136839 = r136833 / r136827;
double r136840 = pow(r136839, r136826);
double r136841 = r136838 + r136840;
double r136842 = r136841 + r136838;
double r136843 = r136837 * r136842;
double r136844 = r136826 / r136843;
return r136844;
}
double f(double t, double l, double k) {
double r136845 = l;
double r136846 = -1.361561607826164e+74;
bool r136847 = r136845 <= r136846;
double r136848 = 2.0;
double r136849 = t;
double r136850 = cbrt(r136849);
double r136851 = r136850 * r136850;
double r136852 = 3.0;
double r136853 = 2.0;
double r136854 = r136852 / r136853;
double r136855 = pow(r136851, r136854);
double r136856 = cbrt(r136845);
double r136857 = r136856 * r136856;
double r136858 = r136855 / r136857;
double r136859 = r136855 / r136856;
double r136860 = r136858 * r136859;
double r136861 = pow(r136850, r136852);
double r136862 = r136861 / r136845;
double r136863 = r136860 * r136862;
double r136864 = k;
double r136865 = sin(r136864);
double r136866 = r136863 * r136865;
double r136867 = r136848 / r136866;
double r136868 = tan(r136864);
double r136869 = 1.0;
double r136870 = r136864 / r136849;
double r136871 = pow(r136870, r136848);
double r136872 = r136869 + r136871;
double r136873 = r136872 + r136869;
double r136874 = r136868 * r136873;
double r136875 = r136867 / r136874;
double r136876 = -3.671469527357901e-125;
bool r136877 = r136845 <= r136876;
double r136878 = 1.0;
double r136879 = -1.0;
double r136880 = pow(r136879, r136852);
double r136881 = r136878 / r136880;
double r136882 = pow(r136881, r136869);
double r136883 = cbrt(r136879);
double r136884 = 6.0;
double r136885 = pow(r136883, r136884);
double r136886 = 3.0;
double r136887 = pow(r136849, r136886);
double r136888 = pow(r136865, r136853);
double r136889 = r136887 * r136888;
double r136890 = r136885 * r136889;
double r136891 = cos(r136864);
double r136892 = pow(r136845, r136853);
double r136893 = r136891 * r136892;
double r136894 = r136890 / r136893;
double r136895 = r136882 * r136894;
double r136896 = r136848 * r136895;
double r136897 = pow(r136864, r136853);
double r136898 = r136897 * r136888;
double r136899 = r136849 * r136898;
double r136900 = r136899 / r136893;
double r136901 = r136882 * r136900;
double r136902 = r136896 + r136901;
double r136903 = -r136902;
double r136904 = r136848 / r136903;
double r136905 = 1.1471508441709014e-119;
bool r136906 = r136845 <= r136905;
double r136907 = 0.3333333333333333;
double r136908 = r136907 * r136852;
double r136909 = pow(r136849, r136908);
double r136910 = r136909 / r136845;
double r136911 = r136860 * r136910;
double r136912 = r136911 * r136865;
double r136913 = r136912 * r136874;
double r136914 = r136848 / r136913;
double r136915 = 2.5473044447906667e+119;
bool r136916 = r136845 <= r136915;
double r136917 = r136916 ? r136904 : r136914;
double r136918 = r136906 ? r136914 : r136917;
double r136919 = r136877 ? r136904 : r136918;
double r136920 = r136847 ? r136875 : r136919;
return r136920;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < -1.361561607826164e+74Initial program 53.0
rmApplied add-cube-cbrt53.1
Applied unpow-prod-down53.1
Applied times-frac38.5
rmApplied add-cube-cbrt38.5
Applied sqr-pow38.5
Applied times-frac27.9
rmApplied associate-*l*27.9
rmApplied associate-/r*28.1
if -1.361561607826164e+74 < l < -3.671469527357901e-125 or 1.1471508441709014e-119 < l < 2.5473044447906667e+119Initial program 25.5
rmApplied add-cube-cbrt25.8
Applied unpow-prod-down25.8
Applied times-frac23.2
rmApplied add-cube-cbrt23.2
Applied sqr-pow23.2
Applied times-frac22.9
Taylor expanded around -inf 16.8
if -3.671469527357901e-125 < l < 1.1471508441709014e-119 or 2.5473044447906667e+119 < l Initial program 32.2
rmApplied add-cube-cbrt32.2
Applied unpow-prod-down32.2
Applied times-frac23.5
rmApplied add-cube-cbrt23.5
Applied sqr-pow23.5
Applied times-frac16.3
rmApplied associate-*l*15.9
rmApplied pow1/339.8
Applied pow-pow15.9
Final simplification17.9
herbie shell --seed 2019347 +o rules:numerics
(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))))