\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 -2.484239696005438759151946975339159010231 \cdot 10^{123}:\\
\;\;\;\;\frac{2}{\left(\left(\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\ell}\right) \cdot \frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\ell}\right) \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;\ell \le -25217102.0555565096437931060791015625:\\
\;\;\;\;\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{elif}\;\ell \le 2.297049428607241564178931690682164825056 \cdot 10^{-139}:\\
\;\;\;\;\frac{2}{\left(\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \left(\frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;\ell \le 5.548140519913770723669052513231581727514 \cdot 10^{121}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot {\ell}^{2}} - {\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}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\ell}\right) \cdot \frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\ell}\right) \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\end{array}double f(double t, double l, double k) {
double r113918 = 2.0;
double r113919 = t;
double r113920 = 3.0;
double r113921 = pow(r113919, r113920);
double r113922 = l;
double r113923 = r113922 * r113922;
double r113924 = r113921 / r113923;
double r113925 = k;
double r113926 = sin(r113925);
double r113927 = r113924 * r113926;
double r113928 = tan(r113925);
double r113929 = r113927 * r113928;
double r113930 = 1.0;
double r113931 = r113925 / r113919;
double r113932 = pow(r113931, r113918);
double r113933 = r113930 + r113932;
double r113934 = r113933 + r113930;
double r113935 = r113929 * r113934;
double r113936 = r113918 / r113935;
return r113936;
}
double f(double t, double l, double k) {
double r113937 = l;
double r113938 = -2.4842396960054388e+123;
bool r113939 = r113937 <= r113938;
double r113940 = 2.0;
double r113941 = t;
double r113942 = cbrt(r113941);
double r113943 = 3.0;
double r113944 = pow(r113942, r113943);
double r113945 = 0.3333333333333333;
double r113946 = r113945 * r113943;
double r113947 = pow(r113941, r113946);
double r113948 = r113947 / r113937;
double r113949 = r113944 * r113948;
double r113950 = r113949 * r113948;
double r113951 = k;
double r113952 = sin(r113951);
double r113953 = r113950 * r113952;
double r113954 = tan(r113951);
double r113955 = r113953 * r113954;
double r113956 = 1.0;
double r113957 = r113951 / r113941;
double r113958 = pow(r113957, r113940);
double r113959 = r113956 + r113958;
double r113960 = r113959 + r113956;
double r113961 = r113955 * r113960;
double r113962 = r113940 / r113961;
double r113963 = -25217102.05555651;
bool r113964 = r113937 <= r113963;
double r113965 = 1.0;
double r113966 = -1.0;
double r113967 = pow(r113966, r113943);
double r113968 = r113965 / r113967;
double r113969 = pow(r113968, r113956);
double r113970 = cbrt(r113966);
double r113971 = 9.0;
double r113972 = pow(r113970, r113971);
double r113973 = 3.0;
double r113974 = pow(r113941, r113973);
double r113975 = 2.0;
double r113976 = pow(r113952, r113975);
double r113977 = r113974 * r113976;
double r113978 = r113972 * r113977;
double r113979 = cos(r113951);
double r113980 = pow(r113937, r113975);
double r113981 = r113979 * r113980;
double r113982 = r113978 / r113981;
double r113983 = r113969 * r113982;
double r113984 = r113940 * r113983;
double r113985 = pow(r113951, r113975);
double r113986 = r113985 * r113941;
double r113987 = r113976 * r113986;
double r113988 = r113972 * r113987;
double r113989 = r113988 / r113981;
double r113990 = r113969 * r113989;
double r113991 = r113984 + r113990;
double r113992 = r113940 / r113991;
double r113993 = 2.2970494286072416e-139;
bool r113994 = r113937 <= r113993;
double r113995 = r113944 / r113937;
double r113996 = r113944 * r113995;
double r113997 = r113948 * r113952;
double r113998 = r113996 * r113997;
double r113999 = r113998 * r113954;
double r114000 = r113999 * r113960;
double r114001 = r113940 / r114000;
double r114002 = 5.548140519913771e+121;
bool r114003 = r113937 <= r114002;
double r114004 = r113977 / r113981;
double r114005 = r113940 * r114004;
double r114006 = r113985 * r113976;
double r114007 = r113941 * r114006;
double r114008 = r114007 / r113981;
double r114009 = r113969 * r114008;
double r114010 = r114005 - r114009;
double r114011 = r113940 / r114010;
double r114012 = r114003 ? r114011 : r113962;
double r114013 = r113994 ? r114001 : r114012;
double r114014 = r113964 ? r113992 : r114013;
double r114015 = r113939 ? r113962 : r114014;
return r114015;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < -2.4842396960054388e+123 or 5.548140519913771e+121 < l Initial program 59.0
rmApplied add-cube-cbrt59.1
Applied unpow-prod-down59.1
Applied times-frac41.7
rmApplied *-un-lft-identity41.7
Applied unpow-prod-down41.7
Applied times-frac27.4
Simplified27.4
rmApplied pow1/347.0
Applied pow-pow27.2
rmApplied pow1/347.0
Applied pow-pow27.0
if -2.4842396960054388e+123 < l < -25217102.05555651Initial program 31.9
rmApplied add-cube-cbrt32.2
Applied unpow-prod-down32.2
Applied times-frac28.4
Taylor expanded around -inf 25.5
if -25217102.05555651 < l < 2.2970494286072416e-139Initial program 23.0
rmApplied add-cube-cbrt23.1
Applied unpow-prod-down23.1
Applied times-frac17.7
rmApplied *-un-lft-identity17.7
Applied unpow-prod-down17.7
Applied times-frac14.2
Simplified14.2
rmApplied pow1/339.1
Applied pow-pow14.1
rmApplied associate-*l*11.5
if 2.2970494286072416e-139 < l < 5.548140519913771e+121Initial program 26.5
Taylor expanded around -inf 17.9
Final simplification17.5
herbie shell --seed 2020001
(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))))