\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 r116481 = 2.0;
double r116482 = t;
double r116483 = 3.0;
double r116484 = pow(r116482, r116483);
double r116485 = l;
double r116486 = r116485 * r116485;
double r116487 = r116484 / r116486;
double r116488 = k;
double r116489 = sin(r116488);
double r116490 = r116487 * r116489;
double r116491 = tan(r116488);
double r116492 = r116490 * r116491;
double r116493 = 1.0;
double r116494 = r116488 / r116482;
double r116495 = pow(r116494, r116481);
double r116496 = r116493 + r116495;
double r116497 = r116496 + r116493;
double r116498 = r116492 * r116497;
double r116499 = r116481 / r116498;
return r116499;
}
double f(double t, double l, double k) {
double r116500 = l;
double r116501 = -6.624617014004735e+24;
bool r116502 = r116500 <= r116501;
double r116503 = 2.0;
double r116504 = t;
double r116505 = cbrt(r116504);
double r116506 = r116505 * r116505;
double r116507 = 3.0;
double r116508 = 2.0;
double r116509 = r116507 / r116508;
double r116510 = pow(r116506, r116509);
double r116511 = r116510 / r116500;
double r116512 = r116510 * r116511;
double r116513 = pow(r116505, r116507);
double r116514 = r116513 / r116500;
double r116515 = k;
double r116516 = sin(r116515);
double r116517 = r116514 * r116516;
double r116518 = r116512 * r116517;
double r116519 = cbrt(r116518);
double r116520 = r116519 * r116519;
double r116521 = r116520 * r116519;
double r116522 = tan(r116515);
double r116523 = r116521 * r116522;
double r116524 = 1.0;
double r116525 = r116515 / r116504;
double r116526 = pow(r116525, r116503);
double r116527 = r116524 + r116526;
double r116528 = r116527 + r116524;
double r116529 = r116523 * r116528;
double r116530 = r116503 / r116529;
double r116531 = -1.4374109414254249e-122;
bool r116532 = r116500 <= r116531;
double r116533 = 1.0;
double r116534 = -1.0;
double r116535 = pow(r116534, r116507);
double r116536 = r116533 / r116535;
double r116537 = pow(r116536, r116524);
double r116538 = cbrt(r116534);
double r116539 = 9.0;
double r116540 = pow(r116538, r116539);
double r116541 = 3.0;
double r116542 = pow(r116504, r116541);
double r116543 = pow(r116516, r116508);
double r116544 = r116542 * r116543;
double r116545 = r116540 * r116544;
double r116546 = cos(r116515);
double r116547 = pow(r116500, r116508);
double r116548 = r116546 * r116547;
double r116549 = r116545 / r116548;
double r116550 = r116537 * r116549;
double r116551 = r116503 * r116550;
double r116552 = pow(r116515, r116508);
double r116553 = r116552 * r116504;
double r116554 = r116543 * r116553;
double r116555 = r116540 * r116554;
double r116556 = r116555 / r116548;
double r116557 = r116537 * r116556;
double r116558 = r116551 + r116557;
double r116559 = r116503 / r116558;
double r116560 = 3.3841965472717056e-77;
bool r116561 = r116500 <= r116560;
double r116562 = pow(r116504, r116524);
double r116563 = pow(r116562, r116524);
double r116564 = r116516 / r116500;
double r116565 = r116563 * r116564;
double r116566 = r116512 * r116565;
double r116567 = r116566 * r116522;
double r116568 = r116567 * r116528;
double r116569 = r116503 / r116568;
double r116570 = 2.6298142138033434e+112;
bool r116571 = r116500 <= r116570;
double r116572 = r116504 * r116543;
double r116573 = r116552 * r116572;
double r116574 = r116573 / r116548;
double r116575 = r116544 / r116548;
double r116576 = r116503 * r116575;
double r116577 = r116574 + r116576;
double r116578 = r116503 / r116577;
double r116579 = r116500 / r116513;
double r116580 = r116513 / r116579;
double r116581 = r116580 * r116514;
double r116582 = r116581 * r116516;
double r116583 = r116582 * r116522;
double r116584 = r116583 * r116528;
double r116585 = r116503 / r116584;
double r116586 = r116571 ? r116578 : r116585;
double r116587 = r116561 ? r116569 : r116586;
double r116588 = r116532 ? r116559 : r116587;
double r116589 = r116502 ? r116530 : r116588;
return r116589;
}



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