\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}{\mathsf{fma}\left(2, {\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}}, {\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}}\right)}\\
\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 r125173 = 2.0;
double r125174 = t;
double r125175 = 3.0;
double r125176 = pow(r125174, r125175);
double r125177 = l;
double r125178 = r125177 * r125177;
double r125179 = r125176 / r125178;
double r125180 = k;
double r125181 = sin(r125180);
double r125182 = r125179 * r125181;
double r125183 = tan(r125180);
double r125184 = r125182 * r125183;
double r125185 = 1.0;
double r125186 = r125180 / r125174;
double r125187 = pow(r125186, r125173);
double r125188 = r125185 + r125187;
double r125189 = r125188 + r125185;
double r125190 = r125184 * r125189;
double r125191 = r125173 / r125190;
return r125191;
}
double f(double t, double l, double k) {
double r125192 = l;
double r125193 = -2.4842396960054388e+123;
bool r125194 = r125192 <= r125193;
double r125195 = 2.0;
double r125196 = t;
double r125197 = cbrt(r125196);
double r125198 = 3.0;
double r125199 = pow(r125197, r125198);
double r125200 = 0.3333333333333333;
double r125201 = r125200 * r125198;
double r125202 = pow(r125196, r125201);
double r125203 = r125202 / r125192;
double r125204 = r125199 * r125203;
double r125205 = r125204 * r125203;
double r125206 = k;
double r125207 = sin(r125206);
double r125208 = r125205 * r125207;
double r125209 = tan(r125206);
double r125210 = r125208 * r125209;
double r125211 = 1.0;
double r125212 = r125206 / r125196;
double r125213 = pow(r125212, r125195);
double r125214 = r125211 + r125213;
double r125215 = r125214 + r125211;
double r125216 = r125210 * r125215;
double r125217 = r125195 / r125216;
double r125218 = -25217102.05555651;
bool r125219 = r125192 <= r125218;
double r125220 = 1.0;
double r125221 = -1.0;
double r125222 = pow(r125221, r125198);
double r125223 = r125220 / r125222;
double r125224 = pow(r125223, r125211);
double r125225 = cbrt(r125221);
double r125226 = 9.0;
double r125227 = pow(r125225, r125226);
double r125228 = 3.0;
double r125229 = pow(r125196, r125228);
double r125230 = 2.0;
double r125231 = pow(r125207, r125230);
double r125232 = r125229 * r125231;
double r125233 = r125227 * r125232;
double r125234 = cos(r125206);
double r125235 = pow(r125192, r125230);
double r125236 = r125234 * r125235;
double r125237 = r125233 / r125236;
double r125238 = r125224 * r125237;
double r125239 = pow(r125206, r125230);
double r125240 = r125239 * r125196;
double r125241 = r125231 * r125240;
double r125242 = r125227 * r125241;
double r125243 = r125242 / r125236;
double r125244 = r125224 * r125243;
double r125245 = fma(r125195, r125238, r125244);
double r125246 = r125195 / r125245;
double r125247 = 2.2970494286072416e-139;
bool r125248 = r125192 <= r125247;
double r125249 = r125199 / r125192;
double r125250 = r125199 * r125249;
double r125251 = r125203 * r125207;
double r125252 = r125250 * r125251;
double r125253 = r125252 * r125209;
double r125254 = r125253 * r125215;
double r125255 = r125195 / r125254;
double r125256 = 5.548140519913771e+121;
bool r125257 = r125192 <= r125256;
double r125258 = pow(r125221, r125195);
double r125259 = r125220 / r125258;
double r125260 = pow(r125259, r125211);
double r125261 = 6.0;
double r125262 = pow(r125225, r125261);
double r125263 = r125262 * r125232;
double r125264 = r125263 / r125236;
double r125265 = r125260 * r125264;
double r125266 = r125195 * r125265;
double r125267 = r125239 * r125231;
double r125268 = r125196 * r125267;
double r125269 = r125268 / r125236;
double r125270 = r125224 * r125269;
double r125271 = r125266 - r125270;
double r125272 = r125195 / r125271;
double r125273 = r125257 ? r125272 : r125217;
double r125274 = r125248 ? r125255 : r125273;
double r125275 = r125219 ? r125246 : r125274;
double r125276 = r125194 ? r125217 : r125275;
return r125276;
}



Bits error versus t



Bits error versus l



Bits error versus k
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
Simplified25.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 +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))))