\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 r122391 = 2.0;
double r122392 = t;
double r122393 = 3.0;
double r122394 = pow(r122392, r122393);
double r122395 = l;
double r122396 = r122395 * r122395;
double r122397 = r122394 / r122396;
double r122398 = k;
double r122399 = sin(r122398);
double r122400 = r122397 * r122399;
double r122401 = tan(r122398);
double r122402 = r122400 * r122401;
double r122403 = 1.0;
double r122404 = r122398 / r122392;
double r122405 = pow(r122404, r122391);
double r122406 = r122403 + r122405;
double r122407 = r122406 + r122403;
double r122408 = r122402 * r122407;
double r122409 = r122391 / r122408;
return r122409;
}
double f(double t, double l, double k) {
double r122410 = l;
double r122411 = -2.4842396960054388e+123;
bool r122412 = r122410 <= r122411;
double r122413 = 2.0;
double r122414 = t;
double r122415 = cbrt(r122414);
double r122416 = 3.0;
double r122417 = pow(r122415, r122416);
double r122418 = 0.3333333333333333;
double r122419 = r122418 * r122416;
double r122420 = pow(r122414, r122419);
double r122421 = r122420 / r122410;
double r122422 = r122417 * r122421;
double r122423 = r122422 * r122421;
double r122424 = k;
double r122425 = sin(r122424);
double r122426 = r122423 * r122425;
double r122427 = tan(r122424);
double r122428 = r122426 * r122427;
double r122429 = 1.0;
double r122430 = r122424 / r122414;
double r122431 = pow(r122430, r122413);
double r122432 = r122429 + r122431;
double r122433 = r122432 + r122429;
double r122434 = r122428 * r122433;
double r122435 = r122413 / r122434;
double r122436 = -25217102.05555651;
bool r122437 = r122410 <= r122436;
double r122438 = 1.0;
double r122439 = -1.0;
double r122440 = pow(r122439, r122416);
double r122441 = r122438 / r122440;
double r122442 = pow(r122441, r122429);
double r122443 = cbrt(r122439);
double r122444 = 9.0;
double r122445 = pow(r122443, r122444);
double r122446 = 3.0;
double r122447 = pow(r122414, r122446);
double r122448 = 2.0;
double r122449 = pow(r122425, r122448);
double r122450 = r122447 * r122449;
double r122451 = r122445 * r122450;
double r122452 = cos(r122424);
double r122453 = pow(r122410, r122448);
double r122454 = r122452 * r122453;
double r122455 = r122451 / r122454;
double r122456 = r122442 * r122455;
double r122457 = pow(r122424, r122448);
double r122458 = r122457 * r122414;
double r122459 = r122449 * r122458;
double r122460 = r122445 * r122459;
double r122461 = r122460 / r122454;
double r122462 = r122442 * r122461;
double r122463 = fma(r122413, r122456, r122462);
double r122464 = r122413 / r122463;
double r122465 = 2.2970494286072416e-139;
bool r122466 = r122410 <= r122465;
double r122467 = r122417 / r122410;
double r122468 = r122417 * r122467;
double r122469 = r122421 * r122425;
double r122470 = r122468 * r122469;
double r122471 = r122470 * r122427;
double r122472 = r122471 * r122433;
double r122473 = r122413 / r122472;
double r122474 = 5.548140519913771e+121;
bool r122475 = r122410 <= r122474;
double r122476 = pow(r122439, r122413);
double r122477 = r122438 / r122476;
double r122478 = pow(r122477, r122429);
double r122479 = 6.0;
double r122480 = pow(r122443, r122479);
double r122481 = r122480 * r122450;
double r122482 = r122481 / r122454;
double r122483 = r122478 * r122482;
double r122484 = r122413 * r122483;
double r122485 = r122457 * r122449;
double r122486 = r122414 * r122485;
double r122487 = r122486 / r122454;
double r122488 = r122442 * r122487;
double r122489 = r122484 - r122488;
double r122490 = r122413 / r122489;
double r122491 = r122475 ? r122490 : r122435;
double r122492 = r122466 ? r122473 : r122491;
double r122493 = r122437 ? r122464 : r122492;
double r122494 = r122412 ? r122435 : r122493;
return r122494;
}



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))))