\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 -6624617014004734683512832:\\
\;\;\;\;\frac{2}{\left(\left(\left(\sqrt[3]{\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\ell}\right) \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)} \cdot \sqrt[3]{\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\ell}\right) \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)}\right) \cdot \sqrt[3]{\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\ell}\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)}\\
\mathbf{elif}\;\ell \le -1.43741094142542485967444676002052768528 \cdot 10^{-122}:\\
\;\;\;\;\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 3.384196547271705606214536606391309350479 \cdot 10^{-77}:\\
\;\;\;\;\frac{2}{\left(\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\ell}\right) \cdot \left({\left({t}^{1}\right)}^{1} \cdot \frac{\sin k}{\ell}\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;\ell \le 2.62981421380334335697111505575029823235 \cdot 10^{112}:\\
\;\;\;\;\frac{2}{\frac{{k}^{2} \cdot \left(t \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}} + 2 \cdot \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot {\ell}^{2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\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 r116891 = 2.0;
double r116892 = t;
double r116893 = 3.0;
double r116894 = pow(r116892, r116893);
double r116895 = l;
double r116896 = r116895 * r116895;
double r116897 = r116894 / r116896;
double r116898 = k;
double r116899 = sin(r116898);
double r116900 = r116897 * r116899;
double r116901 = tan(r116898);
double r116902 = r116900 * r116901;
double r116903 = 1.0;
double r116904 = r116898 / r116892;
double r116905 = pow(r116904, r116891);
double r116906 = r116903 + r116905;
double r116907 = r116906 + r116903;
double r116908 = r116902 * r116907;
double r116909 = r116891 / r116908;
return r116909;
}
double f(double t, double l, double k) {
double r116910 = l;
double r116911 = -6.624617014004735e+24;
bool r116912 = r116910 <= r116911;
double r116913 = 2.0;
double r116914 = t;
double r116915 = cbrt(r116914);
double r116916 = r116915 * r116915;
double r116917 = 3.0;
double r116918 = 2.0;
double r116919 = r116917 / r116918;
double r116920 = pow(r116916, r116919);
double r116921 = r116920 / r116910;
double r116922 = r116920 * r116921;
double r116923 = pow(r116915, r116917);
double r116924 = r116923 / r116910;
double r116925 = k;
double r116926 = sin(r116925);
double r116927 = r116924 * r116926;
double r116928 = r116922 * r116927;
double r116929 = cbrt(r116928);
double r116930 = r116929 * r116929;
double r116931 = r116930 * r116929;
double r116932 = tan(r116925);
double r116933 = r116931 * r116932;
double r116934 = 1.0;
double r116935 = r116925 / r116914;
double r116936 = pow(r116935, r116913);
double r116937 = r116934 + r116936;
double r116938 = r116937 + r116934;
double r116939 = r116933 * r116938;
double r116940 = r116913 / r116939;
double r116941 = -1.4374109414254249e-122;
bool r116942 = r116910 <= r116941;
double r116943 = 1.0;
double r116944 = -1.0;
double r116945 = pow(r116944, r116917);
double r116946 = r116943 / r116945;
double r116947 = pow(r116946, r116934);
double r116948 = cbrt(r116944);
double r116949 = 9.0;
double r116950 = pow(r116948, r116949);
double r116951 = 3.0;
double r116952 = pow(r116914, r116951);
double r116953 = pow(r116926, r116918);
double r116954 = r116952 * r116953;
double r116955 = r116950 * r116954;
double r116956 = cos(r116925);
double r116957 = pow(r116910, r116918);
double r116958 = r116956 * r116957;
double r116959 = r116955 / r116958;
double r116960 = r116947 * r116959;
double r116961 = r116913 * r116960;
double r116962 = pow(r116925, r116918);
double r116963 = r116962 * r116914;
double r116964 = r116953 * r116963;
double r116965 = r116950 * r116964;
double r116966 = r116965 / r116958;
double r116967 = r116947 * r116966;
double r116968 = r116961 + r116967;
double r116969 = r116913 / r116968;
double r116970 = 3.3841965472717056e-77;
bool r116971 = r116910 <= r116970;
double r116972 = pow(r116914, r116934);
double r116973 = pow(r116972, r116934);
double r116974 = r116926 / r116910;
double r116975 = r116973 * r116974;
double r116976 = r116922 * r116975;
double r116977 = r116976 * r116932;
double r116978 = r116977 * r116938;
double r116979 = r116913 / r116978;
double r116980 = 2.6298142138033434e+112;
bool r116981 = r116910 <= r116980;
double r116982 = r116914 * r116953;
double r116983 = r116962 * r116982;
double r116984 = r116983 / r116958;
double r116985 = r116954 / r116958;
double r116986 = r116913 * r116985;
double r116987 = r116984 + r116986;
double r116988 = r116913 / r116987;
double r116989 = r116910 / r116923;
double r116990 = r116923 / r116989;
double r116991 = r116990 * r116924;
double r116992 = r116991 * r116926;
double r116993 = r116992 * r116932;
double r116994 = r116993 * r116938;
double r116995 = r116913 / r116994;
double r116996 = r116981 ? r116988 : r116995;
double r116997 = r116971 ? r116979 : r116996;
double r116998 = r116942 ? r116969 : r116997;
double r116999 = r116912 ? r116940 : r116998;
return r116999;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < -6.624617014004735e+24Initial program 46.9
rmApplied add-cube-cbrt47.1
Applied unpow-prod-down47.1
Applied times-frac35.2
rmApplied *-un-lft-identity35.2
Applied sqr-pow35.2
Applied times-frac26.4
Simplified26.4
rmApplied associate-*l*24.3
rmApplied add-cube-cbrt24.3
if -6.624617014004735e+24 < l < -1.4374109414254249e-122Initial program 22.6
rmApplied add-cube-cbrt22.9
Applied unpow-prod-down22.9
Applied times-frac21.5
Taylor expanded around -inf 15.1
if -1.4374109414254249e-122 < l < 3.3841965472717056e-77Initial program 23.4
rmApplied add-cube-cbrt23.4
Applied unpow-prod-down23.4
Applied times-frac17.4
rmApplied *-un-lft-identity17.4
Applied sqr-pow17.4
Applied times-frac12.1
Simplified12.1
rmApplied associate-*l*10.0
Taylor expanded around inf 9.9
if 3.3841965472717056e-77 < l < 2.6298142138033434e+112Initial program 27.2
Taylor expanded around inf 18.1
if 2.6298142138033434e+112 < l Initial program 58.2
rmApplied add-cube-cbrt58.2
Applied unpow-prod-down58.2
Applied times-frac40.8
rmApplied unpow-prod-down40.8
Applied associate-/l*26.5
Final simplification16.5
herbie shell --seed 2019353
(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))))