\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 -5.41720610659310475 \cdot 10^{-101}:\\
\;\;\;\;\frac{2}{\frac{\left(\sqrt[3]{\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)} \cdot \sqrt[3]{\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)}\right) \cdot \sqrt[3]{\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{\ell}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}}\\
\mathbf{elif}\;t \le 1.6736513015160059 \cdot 10^{-53}:\\
\;\;\;\;\frac{2}{\frac{\mathsf{fma}\left({\left(\frac{1}{{-1}^{2}}\right)}^{1}, \frac{{k}^{2} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot \ell}, 2 \cdot \left({\left(\frac{1}{{-1}^{2}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{6} \cdot \left({t}^{2} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot \ell}\right)\right)}{\frac{\ell}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \left(\left(\sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k} \cdot \sqrt[3]{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k}\right) \cdot \sqrt[3]{\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{\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 r115424 = 2.0;
double r115425 = t;
double r115426 = 3.0;
double r115427 = pow(r115425, r115426);
double r115428 = l;
double r115429 = r115428 * r115428;
double r115430 = r115427 / r115429;
double r115431 = k;
double r115432 = sin(r115431);
double r115433 = r115430 * r115432;
double r115434 = tan(r115431);
double r115435 = r115433 * r115434;
double r115436 = 1.0;
double r115437 = r115431 / r115425;
double r115438 = pow(r115437, r115424);
double r115439 = r115436 + r115438;
double r115440 = r115439 + r115436;
double r115441 = r115435 * r115440;
double r115442 = r115424 / r115441;
return r115442;
}
double f(double t, double l, double k) {
double r115443 = t;
double r115444 = -5.417206106593105e-101;
bool r115445 = r115443 <= r115444;
double r115446 = 2.0;
double r115447 = cbrt(r115443);
double r115448 = r115447 * r115447;
double r115449 = 3.0;
double r115450 = 2.0;
double r115451 = r115449 / r115450;
double r115452 = pow(r115448, r115451);
double r115453 = pow(r115447, r115449);
double r115454 = l;
double r115455 = r115453 / r115454;
double r115456 = k;
double r115457 = sin(r115456);
double r115458 = r115455 * r115457;
double r115459 = r115452 * r115458;
double r115460 = tan(r115456);
double r115461 = r115459 * r115460;
double r115462 = 1.0;
double r115463 = r115456 / r115443;
double r115464 = pow(r115463, r115446);
double r115465 = r115462 + r115464;
double r115466 = r115465 + r115462;
double r115467 = r115461 * r115466;
double r115468 = cbrt(r115467);
double r115469 = r115468 * r115468;
double r115470 = r115469 * r115468;
double r115471 = r115454 / r115452;
double r115472 = r115470 / r115471;
double r115473 = r115446 / r115472;
double r115474 = 1.673651301516006e-53;
bool r115475 = r115443 <= r115474;
double r115476 = 1.0;
double r115477 = -1.0;
double r115478 = pow(r115477, r115446);
double r115479 = r115476 / r115478;
double r115480 = pow(r115479, r115462);
double r115481 = pow(r115456, r115450);
double r115482 = pow(r115457, r115450);
double r115483 = r115481 * r115482;
double r115484 = cos(r115456);
double r115485 = r115484 * r115454;
double r115486 = r115483 / r115485;
double r115487 = cbrt(r115477);
double r115488 = 6.0;
double r115489 = pow(r115487, r115488);
double r115490 = pow(r115443, r115450);
double r115491 = r115490 * r115482;
double r115492 = r115489 * r115491;
double r115493 = r115492 / r115485;
double r115494 = r115480 * r115493;
double r115495 = r115446 * r115494;
double r115496 = fma(r115480, r115486, r115495);
double r115497 = r115496 / r115471;
double r115498 = r115446 / r115497;
double r115499 = cbrt(r115458);
double r115500 = r115499 * r115499;
double r115501 = r115500 * r115499;
double r115502 = r115452 * r115501;
double r115503 = r115502 * r115460;
double r115504 = r115503 * r115466;
double r115505 = r115504 / r115471;
double r115506 = r115446 / r115505;
double r115507 = r115475 ? r115498 : r115506;
double r115508 = r115445 ? r115473 : r115507;
return r115508;
}



Bits error versus t



Bits error versus l



Bits error versus k
if t < -5.417206106593105e-101Initial program 22.9
rmApplied add-cube-cbrt23.1
Applied unpow-prod-down23.1
Applied times-frac16.7
Applied associate-*l*14.3
rmApplied sqr-pow14.3
Applied associate-/l*10.0
rmApplied associate-*l/9.1
Applied associate-*l/7.0
Applied associate-*l/6.3
rmApplied add-cube-cbrt6.4
if -5.417206106593105e-101 < t < 1.673651301516006e-53Initial program 59.3
rmApplied add-cube-cbrt59.3
Applied unpow-prod-down59.3
Applied times-frac50.7
Applied associate-*l*50.3
rmApplied sqr-pow50.3
Applied associate-/l*42.1
rmApplied associate-*l/42.1
Applied associate-*l/43.2
Applied associate-*l/39.3
Taylor expanded around -inf 30.2
Simplified30.2
if 1.673651301516006e-53 < t Initial program 22.0
rmApplied add-cube-cbrt22.2
Applied unpow-prod-down22.2
Applied times-frac16.0
Applied associate-*l*13.9
rmApplied sqr-pow13.9
Applied associate-/l*8.8
rmApplied associate-*l/7.4
Applied associate-*l/5.6
Applied associate-*l/5.1
rmApplied add-cube-cbrt5.1
Final simplification12.0
herbie shell --seed 2020047 +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))))