\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}{\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 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}{\mathsf{fma}\left(2, \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot {\ell}^{2}}, \frac{{k}^{2} \cdot \left(t \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}}\right)}\\
\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 r125465 = 2.0;
double r125466 = t;
double r125467 = 3.0;
double r125468 = pow(r125466, r125467);
double r125469 = l;
double r125470 = r125469 * r125469;
double r125471 = r125468 / r125470;
double r125472 = k;
double r125473 = sin(r125472);
double r125474 = r125471 * r125473;
double r125475 = tan(r125472);
double r125476 = r125474 * r125475;
double r125477 = 1.0;
double r125478 = r125472 / r125466;
double r125479 = pow(r125478, r125465);
double r125480 = r125477 + r125479;
double r125481 = r125480 + r125477;
double r125482 = r125476 * r125481;
double r125483 = r125465 / r125482;
return r125483;
}
double f(double t, double l, double k) {
double r125484 = l;
double r125485 = -6.624617014004735e+24;
bool r125486 = r125484 <= r125485;
double r125487 = 2.0;
double r125488 = t;
double r125489 = cbrt(r125488);
double r125490 = r125489 * r125489;
double r125491 = 3.0;
double r125492 = 2.0;
double r125493 = r125491 / r125492;
double r125494 = pow(r125490, r125493);
double r125495 = r125494 / r125484;
double r125496 = r125494 * r125495;
double r125497 = pow(r125489, r125491);
double r125498 = r125497 / r125484;
double r125499 = k;
double r125500 = sin(r125499);
double r125501 = r125498 * r125500;
double r125502 = r125496 * r125501;
double r125503 = cbrt(r125502);
double r125504 = r125503 * r125503;
double r125505 = r125504 * r125503;
double r125506 = tan(r125499);
double r125507 = r125505 * r125506;
double r125508 = 1.0;
double r125509 = r125499 / r125488;
double r125510 = pow(r125509, r125487);
double r125511 = r125508 + r125510;
double r125512 = r125511 + r125508;
double r125513 = r125507 * r125512;
double r125514 = r125487 / r125513;
double r125515 = -1.4374109414254249e-122;
bool r125516 = r125484 <= r125515;
double r125517 = 1.0;
double r125518 = -1.0;
double r125519 = pow(r125518, r125491);
double r125520 = r125517 / r125519;
double r125521 = pow(r125520, r125508);
double r125522 = cbrt(r125518);
double r125523 = 9.0;
double r125524 = pow(r125522, r125523);
double r125525 = 3.0;
double r125526 = pow(r125488, r125525);
double r125527 = pow(r125500, r125492);
double r125528 = r125526 * r125527;
double r125529 = r125524 * r125528;
double r125530 = cos(r125499);
double r125531 = pow(r125484, r125492);
double r125532 = r125530 * r125531;
double r125533 = r125529 / r125532;
double r125534 = r125521 * r125533;
double r125535 = pow(r125499, r125492);
double r125536 = r125535 * r125488;
double r125537 = r125527 * r125536;
double r125538 = r125524 * r125537;
double r125539 = r125538 / r125532;
double r125540 = r125521 * r125539;
double r125541 = fma(r125487, r125534, r125540);
double r125542 = r125487 / r125541;
double r125543 = 3.3841965472717056e-77;
bool r125544 = r125484 <= r125543;
double r125545 = pow(r125488, r125508);
double r125546 = pow(r125545, r125508);
double r125547 = r125500 / r125484;
double r125548 = r125546 * r125547;
double r125549 = r125496 * r125548;
double r125550 = r125549 * r125506;
double r125551 = r125550 * r125512;
double r125552 = r125487 / r125551;
double r125553 = 2.6298142138033434e+112;
bool r125554 = r125484 <= r125553;
double r125555 = r125528 / r125532;
double r125556 = r125488 * r125527;
double r125557 = r125535 * r125556;
double r125558 = r125557 / r125532;
double r125559 = fma(r125487, r125555, r125558);
double r125560 = r125487 / r125559;
double r125561 = r125484 / r125497;
double r125562 = r125497 / r125561;
double r125563 = r125562 * r125498;
double r125564 = r125563 * r125500;
double r125565 = r125564 * r125506;
double r125566 = r125565 * r125512;
double r125567 = r125487 / r125566;
double r125568 = r125554 ? r125560 : r125567;
double r125569 = r125544 ? r125552 : r125568;
double r125570 = r125516 ? r125542 : r125569;
double r125571 = r125486 ? r125514 : r125570;
return r125571;
}



Bits error versus t



Bits error versus l



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