\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 \left({\left(\frac{1}{{-1}^{2}}\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}}}\\
\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 r115990 = 2.0;
double r115991 = t;
double r115992 = 3.0;
double r115993 = pow(r115991, r115992);
double r115994 = l;
double r115995 = r115994 * r115994;
double r115996 = r115993 / r115995;
double r115997 = k;
double r115998 = sin(r115997);
double r115999 = r115996 * r115998;
double r116000 = tan(r115997);
double r116001 = r115999 * r116000;
double r116002 = 1.0;
double r116003 = r115997 / r115991;
double r116004 = pow(r116003, r115990);
double r116005 = r116002 + r116004;
double r116006 = r116005 + r116002;
double r116007 = r116001 * r116006;
double r116008 = r115990 / r116007;
return r116008;
}
double f(double t, double l, double k) {
double r116009 = l;
double r116010 = -2.4842396960054388e+123;
bool r116011 = r116009 <= r116010;
double r116012 = 2.0;
double r116013 = t;
double r116014 = cbrt(r116013);
double r116015 = 3.0;
double r116016 = pow(r116014, r116015);
double r116017 = 0.3333333333333333;
double r116018 = r116017 * r116015;
double r116019 = pow(r116013, r116018);
double r116020 = r116019 / r116009;
double r116021 = r116016 * r116020;
double r116022 = r116021 * r116020;
double r116023 = k;
double r116024 = sin(r116023);
double r116025 = r116022 * r116024;
double r116026 = tan(r116023);
double r116027 = r116025 * r116026;
double r116028 = 1.0;
double r116029 = r116023 / r116013;
double r116030 = pow(r116029, r116012);
double r116031 = r116028 + r116030;
double r116032 = r116031 + r116028;
double r116033 = r116027 * r116032;
double r116034 = r116012 / r116033;
double r116035 = -25217102.05555651;
bool r116036 = r116009 <= r116035;
double r116037 = 1.0;
double r116038 = -1.0;
double r116039 = pow(r116038, r116015);
double r116040 = r116037 / r116039;
double r116041 = pow(r116040, r116028);
double r116042 = cbrt(r116038);
double r116043 = 9.0;
double r116044 = pow(r116042, r116043);
double r116045 = 3.0;
double r116046 = pow(r116013, r116045);
double r116047 = 2.0;
double r116048 = pow(r116024, r116047);
double r116049 = r116046 * r116048;
double r116050 = r116044 * r116049;
double r116051 = cos(r116023);
double r116052 = pow(r116009, r116047);
double r116053 = r116051 * r116052;
double r116054 = r116050 / r116053;
double r116055 = r116041 * r116054;
double r116056 = r116012 * r116055;
double r116057 = pow(r116023, r116047);
double r116058 = r116057 * r116013;
double r116059 = r116048 * r116058;
double r116060 = r116044 * r116059;
double r116061 = r116060 / r116053;
double r116062 = r116041 * r116061;
double r116063 = r116056 + r116062;
double r116064 = r116012 / r116063;
double r116065 = 2.2970494286072416e-139;
bool r116066 = r116009 <= r116065;
double r116067 = r116016 / r116009;
double r116068 = r116016 * r116067;
double r116069 = r116020 * r116024;
double r116070 = r116068 * r116069;
double r116071 = r116070 * r116026;
double r116072 = r116071 * r116032;
double r116073 = r116012 / r116072;
double r116074 = 5.548140519913771e+121;
bool r116075 = r116009 <= r116074;
double r116076 = pow(r116038, r116012);
double r116077 = r116037 / r116076;
double r116078 = pow(r116077, r116028);
double r116079 = 6.0;
double r116080 = pow(r116042, r116079);
double r116081 = r116080 * r116049;
double r116082 = r116081 / r116053;
double r116083 = r116078 * r116082;
double r116084 = r116012 * r116083;
double r116085 = r116057 * r116048;
double r116086 = r116013 * r116085;
double r116087 = r116086 / r116053;
double r116088 = r116041 * r116087;
double r116089 = r116084 - r116088;
double r116090 = r116012 / r116089;
double r116091 = r116075 ? r116090 : r116034;
double r116092 = r116066 ? r116073 : r116091;
double r116093 = r116036 ? r116064 : r116092;
double r116094 = r116011 ? r116034 : r116093;
return r116094;
}



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
rmApplied add-cube-cbrt26.7
Applied unpow-prod-down26.7
Applied times-frac24.2
rmApplied *-un-lft-identity24.2
Applied unpow-prod-down24.2
Applied times-frac23.7
Simplified23.7
rmApplied pow1/343.9
Applied pow-pow23.6
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))))