\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}\;t \le -2.71319654436438530491844512119381433671 \cdot 10^{-119}:\\
\;\;\;\;\frac{2}{\frac{\frac{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\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)}{\frac{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}{\sqrt[3]{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}} \cdot \sqrt[3]{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}}}{\frac{\sqrt[3]{\ell}}{\sqrt[3]{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}}}\\
\mathbf{elif}\;t \le -3.396358842890294350312657745725957851569 \cdot 10^{-234}:\\
\;\;\;\;\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}\;t \le -1.889101273676244999235720045944996713013 \cdot 10^{-261}:\\
\;\;\;\;\frac{2}{\frac{\frac{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\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)}{\frac{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}{\sqrt[3]{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}} \cdot \sqrt[3]{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}}}{\frac{\sqrt[3]{\ell}}{\sqrt[3]{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}}}\\
\mathbf{elif}\;t \le -8.381505698861730984029302269680564346164 \cdot 10^{-301}:\\
\;\;\;\;\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}\;t \le 2.007838435386157428168197010709776023626 \cdot 10^{-184}:\\
\;\;\;\;\frac{2}{\frac{\frac{{k}^{4}}{\ell} + 2 \cdot \frac{{k}^{2} \cdot {t}^{2}}{\ell}}{\frac{\ell}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}{\frac{\ell}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}}\\
\end{array}double f(double t, double l, double k) {
double r124463 = 2.0;
double r124464 = t;
double r124465 = 3.0;
double r124466 = pow(r124464, r124465);
double r124467 = l;
double r124468 = r124467 * r124467;
double r124469 = r124466 / r124468;
double r124470 = k;
double r124471 = sin(r124470);
double r124472 = r124469 * r124471;
double r124473 = tan(r124470);
double r124474 = r124472 * r124473;
double r124475 = 1.0;
double r124476 = r124470 / r124464;
double r124477 = pow(r124476, r124463);
double r124478 = r124475 + r124477;
double r124479 = r124478 + r124475;
double r124480 = r124474 * r124479;
double r124481 = r124463 / r124480;
return r124481;
}
double f(double t, double l, double k) {
double r124482 = t;
double r124483 = -2.7131965443643853e-119;
bool r124484 = r124482 <= r124483;
double r124485 = 2.0;
double r124486 = cbrt(r124482);
double r124487 = r124486 * r124486;
double r124488 = 3.0;
double r124489 = 2.0;
double r124490 = r124488 / r124489;
double r124491 = pow(r124487, r124490);
double r124492 = pow(r124486, r124488);
double r124493 = l;
double r124494 = r124492 / r124493;
double r124495 = k;
double r124496 = sin(r124495);
double r124497 = r124494 * r124496;
double r124498 = r124491 * r124497;
double r124499 = tan(r124495);
double r124500 = r124498 * r124499;
double r124501 = 1.0;
double r124502 = r124495 / r124482;
double r124503 = pow(r124502, r124485);
double r124504 = r124501 + r124503;
double r124505 = r124504 + r124501;
double r124506 = r124500 * r124505;
double r124507 = cbrt(r124493);
double r124508 = r124507 * r124507;
double r124509 = cbrt(r124491);
double r124510 = r124509 * r124509;
double r124511 = r124508 / r124510;
double r124512 = r124506 / r124511;
double r124513 = r124507 / r124509;
double r124514 = r124512 / r124513;
double r124515 = r124485 / r124514;
double r124516 = -3.3963588428902944e-234;
bool r124517 = r124482 <= r124516;
double r124518 = 1.0;
double r124519 = -1.0;
double r124520 = pow(r124519, r124488);
double r124521 = r124518 / r124520;
double r124522 = pow(r124521, r124501);
double r124523 = cbrt(r124519);
double r124524 = 9.0;
double r124525 = pow(r124523, r124524);
double r124526 = 3.0;
double r124527 = pow(r124482, r124526);
double r124528 = pow(r124496, r124489);
double r124529 = r124527 * r124528;
double r124530 = r124525 * r124529;
double r124531 = cos(r124495);
double r124532 = pow(r124493, r124489);
double r124533 = r124531 * r124532;
double r124534 = r124530 / r124533;
double r124535 = r124522 * r124534;
double r124536 = r124485 * r124535;
double r124537 = pow(r124495, r124489);
double r124538 = r124537 * r124482;
double r124539 = r124528 * r124538;
double r124540 = r124525 * r124539;
double r124541 = r124540 / r124533;
double r124542 = r124522 * r124541;
double r124543 = r124536 + r124542;
double r124544 = r124485 / r124543;
double r124545 = -1.889101273676245e-261;
bool r124546 = r124482 <= r124545;
double r124547 = -8.381505698861731e-301;
bool r124548 = r124482 <= r124547;
double r124549 = 2.0078384353861574e-184;
bool r124550 = r124482 <= r124549;
double r124551 = 4.0;
double r124552 = pow(r124495, r124551);
double r124553 = r124552 / r124493;
double r124554 = pow(r124482, r124489);
double r124555 = r124537 * r124554;
double r124556 = r124555 / r124493;
double r124557 = r124485 * r124556;
double r124558 = r124553 + r124557;
double r124559 = r124493 / r124491;
double r124560 = r124558 / r124559;
double r124561 = r124485 / r124560;
double r124562 = r124499 * r124505;
double r124563 = r124498 * r124562;
double r124564 = r124563 / r124559;
double r124565 = r124485 / r124564;
double r124566 = r124550 ? r124561 : r124565;
double r124567 = r124548 ? r124544 : r124566;
double r124568 = r124546 ? r124515 : r124567;
double r124569 = r124517 ? r124544 : r124568;
double r124570 = r124484 ? r124515 : r124569;
return r124570;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -2.7131965443643853e-119 or -3.3963588428902944e-234 < t < -1.889101273676245e-261Initial program 25.0
rmApplied add-cube-cbrt25.2
Applied unpow-prod-down25.2
Applied times-frac19.1
Applied associate-*l*17.0
rmApplied sqr-pow17.0
Applied associate-/l*11.9
rmApplied associate-*l/10.9
Applied associate-*l/9.4
Applied associate-*l/8.4
rmApplied add-cube-cbrt8.4
Applied add-cube-cbrt8.4
Applied times-frac8.4
Applied associate-/r*8.4
if -2.7131965443643853e-119 < t < -3.3963588428902944e-234 or -1.889101273676245e-261 < t < -8.381505698861731e-301Initial program 64.0
rmApplied add-cube-cbrt64.0
Applied unpow-prod-down64.0
Applied times-frac54.1
Applied associate-*l*54.1
Taylor expanded around -inf 41.0
if -8.381505698861731e-301 < t < 2.0078384353861574e-184Initial program 64.0
rmApplied add-cube-cbrt64.0
Applied unpow-prod-down64.0
Applied times-frac64.0
Applied associate-*l*64.0
rmApplied sqr-pow64.0
Applied associate-/l*54.3
rmApplied associate-*l/54.3
Applied associate-*l/54.6
Applied associate-*l/52.6
Taylor expanded around 0 52.1
if 2.0078384353861574e-184 < t Initial program 26.9
rmApplied add-cube-cbrt27.1
Applied unpow-prod-down27.1
Applied times-frac19.3
Applied associate-*l*17.5
rmApplied sqr-pow17.5
Applied associate-/l*12.8
rmApplied associate-*l/11.9
Applied associate-*l/10.6
Applied associate-*l/9.4
rmApplied associate-*l*9.2
Final simplification14.8
herbie shell --seed 2019351
(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))))