\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 -1.479648986687081907909669591174294921845 \cdot 10^{163}:\\
\;\;\;\;\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 \tan k}{\frac{\ell}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}} \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;\ell \le \frac{-1743626443581345}{2.879304828507645684998744644919028389677 \cdot 10^{163}}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot {\ell}^{2}} - {\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{elif}\;\ell \le \frac{8201058597573213}{1.405910560794748869628293283651869330897 \cdot 10^{160}}:\\
\;\;\;\;\frac{2}{\frac{\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \sin k}{\frac{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}{1}} \cdot \left(\frac{\tan k}{\frac{\sqrt[3]{\ell}}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}} \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}\\
\mathbf{elif}\;\ell \le 2.196899654104361112665217492985535582092 \cdot 10^{73}:\\
\;\;\;\;\frac{2}{2 \cdot \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot {\ell}^{2}} - {\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}{\frac{\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \sin k}{\sqrt[3]{\frac{\ell}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}} \cdot \sqrt[3]{\frac{\ell}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}} \cdot \left(\frac{\tan k}{\sqrt[3]{\frac{\ell}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}} \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}\\
\end{array}double f(double t, double l, double k) {
double r117534 = 2.0;
double r117535 = t;
double r117536 = 3.0;
double r117537 = pow(r117535, r117536);
double r117538 = l;
double r117539 = r117538 * r117538;
double r117540 = r117537 / r117539;
double r117541 = k;
double r117542 = sin(r117541);
double r117543 = r117540 * r117542;
double r117544 = tan(r117541);
double r117545 = r117543 * r117544;
double r117546 = 1.0;
double r117547 = r117541 / r117535;
double r117548 = pow(r117547, r117534);
double r117549 = r117546 + r117548;
double r117550 = r117549 + r117546;
double r117551 = r117545 * r117550;
double r117552 = r117534 / r117551;
return r117552;
}
double f(double t, double l, double k) {
double r117553 = l;
double r117554 = -1.479648986687082e+163;
bool r117555 = r117553 <= r117554;
double r117556 = 2.0;
double r117557 = t;
double r117558 = cbrt(r117557);
double r117559 = r117558 * r117558;
double r117560 = 3.0;
double r117561 = 2.0;
double r117562 = r117560 / r117561;
double r117563 = pow(r117559, r117562);
double r117564 = pow(r117558, r117560);
double r117565 = r117564 / r117553;
double r117566 = k;
double r117567 = sin(r117566);
double r117568 = r117565 * r117567;
double r117569 = r117563 * r117568;
double r117570 = tan(r117566);
double r117571 = r117569 * r117570;
double r117572 = r117553 / r117563;
double r117573 = r117571 / r117572;
double r117574 = 1.0;
double r117575 = r117566 / r117557;
double r117576 = pow(r117575, r117556);
double r117577 = r117574 + r117576;
double r117578 = r117577 + r117574;
double r117579 = r117573 * r117578;
double r117580 = r117556 / r117579;
double r117581 = -1743626443581345.0;
double r117582 = 2.879304828507646e+163;
double r117583 = r117581 / r117582;
bool r117584 = r117553 <= r117583;
double r117585 = 3.0;
double r117586 = pow(r117557, r117585);
double r117587 = pow(r117567, r117561);
double r117588 = r117586 * r117587;
double r117589 = cos(r117566);
double r117590 = pow(r117553, r117561);
double r117591 = r117589 * r117590;
double r117592 = r117588 / r117591;
double r117593 = r117556 * r117592;
double r117594 = 1.0;
double r117595 = -1.0;
double r117596 = pow(r117595, r117560);
double r117597 = r117594 / r117596;
double r117598 = pow(r117597, r117574);
double r117599 = pow(r117566, r117561);
double r117600 = r117599 * r117587;
double r117601 = r117557 * r117600;
double r117602 = r117601 / r117591;
double r117603 = r117598 * r117602;
double r117604 = r117593 - r117603;
double r117605 = r117556 / r117604;
double r117606 = 8201058597573213.0;
double r117607 = 1.405910560794749e+160;
double r117608 = r117606 / r117607;
bool r117609 = r117553 <= r117608;
double r117610 = r117563 * r117565;
double r117611 = r117610 * r117567;
double r117612 = cbrt(r117553);
double r117613 = r117612 * r117612;
double r117614 = r117613 / r117594;
double r117615 = r117611 / r117614;
double r117616 = r117612 / r117563;
double r117617 = r117570 / r117616;
double r117618 = r117617 * r117578;
double r117619 = r117615 * r117618;
double r117620 = r117556 / r117619;
double r117621 = 2.196899654104361e+73;
bool r117622 = r117553 <= r117621;
double r117623 = cbrt(r117572);
double r117624 = r117623 * r117623;
double r117625 = r117611 / r117624;
double r117626 = r117570 / r117623;
double r117627 = r117626 * r117578;
double r117628 = r117625 * r117627;
double r117629 = r117556 / r117628;
double r117630 = r117622 ? r117605 : r117629;
double r117631 = r117609 ? r117620 : r117630;
double r117632 = r117584 ? r117605 : r117631;
double r117633 = r117555 ? r117580 : r117632;
return r117633;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < -1.479648986687082e+163Initial program 64.0
rmApplied add-cube-cbrt64.0
Applied unpow-prod-down64.0
Applied times-frac43.9
rmApplied sqr-pow43.9
Applied associate-/l*27.4
rmApplied associate-*l/27.4
Applied associate-*l/25.9
Applied associate-*l/25.8
rmApplied associate-*l*22.8
if -1.479648986687082e+163 < l < -6.055720208287474e-149 or 5.833271920894617e-145 < l < 2.196899654104361e+73Initial program 26.7
Taylor expanded around -inf 18.4
if -6.055720208287474e-149 < l < 5.833271920894617e-145Initial program 25.2
rmApplied add-cube-cbrt25.3
Applied unpow-prod-down25.3
Applied times-frac18.7
rmApplied sqr-pow18.7
Applied associate-/l*12.9
rmApplied associate-*l/12.9
Applied associate-*l/9.8
Applied associate-*l/10.5
rmApplied *-un-lft-identity10.5
Applied add-cube-cbrt10.5
Applied times-frac10.5
Applied times-frac8.7
Applied associate-*l*7.5
if 2.196899654104361e+73 < l Initial program 52.0
rmApplied add-cube-cbrt52.1
Applied unpow-prod-down52.1
Applied times-frac38.1
rmApplied sqr-pow38.1
Applied associate-/l*28.0
rmApplied associate-*l/28.0
Applied associate-*l/26.2
Applied associate-*l/26.2
rmApplied add-cube-cbrt26.2
Applied times-frac26.0
Applied associate-*l*25.1
Final simplification16.2
herbie shell --seed 2019304
(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))))