\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}\;k \le -1.932296528705668984693986702426810627746 \cdot 10^{146}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{\sqrt[3]{t}} \cdot \sqrt[3]{\sqrt[3]{t}}\right)}^{3}}{1} \cdot \left(\frac{{\left(\sqrt[3]{\sqrt[3]{t}}\right)}^{3}}{\ell} \cdot \sin k\right)\right)\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;k \le -2.465489873106864392898978142261571099549 \cdot 10^{69}:\\
\;\;\;\;\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}\;k \le \frac{-4688434469438699}{1.333602886575970854815970358144251559429 \cdot 10^{241}}:\\
\;\;\;\;\frac{2}{\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right)\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;k \le \frac{1066040133643155}{2147483648}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell}} \cdot \left(\frac{{t}^{\left(\frac{3}{3}\right)}}{\ell} \cdot \sin k\right)\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;k \le 2.402994973162251042133283519939024529294 \cdot 10^{139}:\\
\;\;\;\;\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}\;k \le 3.506733478446760175749683528136847068451 \cdot 10^{202}:\\
\;\;\;\;\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell}}\right) \cdot \left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right) \cdot \tan k\right)\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\
\mathbf{elif}\;k \le 7.387345331331856574424040723290335999974 \cdot 10^{224}:\\
\;\;\;\;\frac{2}{-\left(2 \cdot \left({\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{{\left(\sqrt[3]{-1}\right)}^{6} \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{t \cdot \left({k}^{2} \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}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right)\right) \cdot \tan k\right) \cdot \left(\sqrt[3]{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1} \cdot \sqrt[3]{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}\right)\right) \cdot \sqrt[3]{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}\\
\end{array}double f(double t, double l, double k) {
double r122563 = 2.0;
double r122564 = t;
double r122565 = 3.0;
double r122566 = pow(r122564, r122565);
double r122567 = l;
double r122568 = r122567 * r122567;
double r122569 = r122566 / r122568;
double r122570 = k;
double r122571 = sin(r122570);
double r122572 = r122569 * r122571;
double r122573 = tan(r122570);
double r122574 = r122572 * r122573;
double r122575 = 1.0;
double r122576 = r122570 / r122564;
double r122577 = pow(r122576, r122563);
double r122578 = r122575 + r122577;
double r122579 = r122578 + r122575;
double r122580 = r122574 * r122579;
double r122581 = r122563 / r122580;
return r122581;
}
double f(double t, double l, double k) {
double r122582 = k;
double r122583 = -1.932296528705669e+146;
bool r122584 = r122582 <= r122583;
double r122585 = 2.0;
double r122586 = t;
double r122587 = cbrt(r122586);
double r122588 = 3.0;
double r122589 = pow(r122587, r122588);
double r122590 = l;
double r122591 = cbrt(r122590);
double r122592 = r122591 * r122591;
double r122593 = r122589 / r122592;
double r122594 = r122589 / r122591;
double r122595 = cbrt(r122587);
double r122596 = r122595 * r122595;
double r122597 = pow(r122596, r122588);
double r122598 = 1.0;
double r122599 = r122597 / r122598;
double r122600 = pow(r122595, r122588);
double r122601 = r122600 / r122590;
double r122602 = sin(r122582);
double r122603 = r122601 * r122602;
double r122604 = r122599 * r122603;
double r122605 = r122594 * r122604;
double r122606 = r122593 * r122605;
double r122607 = tan(r122582);
double r122608 = r122606 * r122607;
double r122609 = 1.0;
double r122610 = r122582 / r122586;
double r122611 = pow(r122610, r122585);
double r122612 = r122609 + r122611;
double r122613 = r122612 + r122609;
double r122614 = r122608 * r122613;
double r122615 = r122585 / r122614;
double r122616 = -2.4654898731068644e+69;
bool r122617 = r122582 <= r122616;
double r122618 = -1.0;
double r122619 = pow(r122618, r122588);
double r122620 = r122598 / r122619;
double r122621 = pow(r122620, r122609);
double r122622 = cbrt(r122618);
double r122623 = 9.0;
double r122624 = pow(r122622, r122623);
double r122625 = 3.0;
double r122626 = pow(r122586, r122625);
double r122627 = 2.0;
double r122628 = pow(r122602, r122627);
double r122629 = r122626 * r122628;
double r122630 = r122624 * r122629;
double r122631 = cos(r122582);
double r122632 = pow(r122590, r122627);
double r122633 = r122631 * r122632;
double r122634 = r122630 / r122633;
double r122635 = r122621 * r122634;
double r122636 = r122585 * r122635;
double r122637 = pow(r122582, r122627);
double r122638 = r122637 * r122586;
double r122639 = r122628 * r122638;
double r122640 = r122624 * r122639;
double r122641 = r122640 / r122633;
double r122642 = r122621 * r122641;
double r122643 = r122636 + r122642;
double r122644 = r122585 / r122643;
double r122645 = -4688434469438699.0;
double r122646 = 1.333602886575971e+241;
double r122647 = r122645 / r122646;
bool r122648 = r122582 <= r122647;
double r122649 = r122589 / r122590;
double r122650 = r122649 * r122602;
double r122651 = r122594 * r122650;
double r122652 = r122651 * r122607;
double r122653 = r122593 * r122652;
double r122654 = r122653 * r122613;
double r122655 = r122585 / r122654;
double r122656 = 1066040133643155.0;
double r122657 = 2147483648.0;
double r122658 = r122656 / r122657;
bool r122659 = r122582 <= r122658;
double r122660 = r122588 / r122625;
double r122661 = pow(r122586, r122660);
double r122662 = r122661 / r122590;
double r122663 = r122662 * r122602;
double r122664 = r122594 * r122663;
double r122665 = r122593 * r122664;
double r122666 = r122665 * r122607;
double r122667 = r122666 * r122613;
double r122668 = r122585 / r122667;
double r122669 = 2.402994973162251e+139;
bool r122670 = r122582 <= r122669;
double r122671 = 3.50673347844676e+202;
bool r122672 = r122582 <= r122671;
double r122673 = r122593 * r122594;
double r122674 = r122650 * r122607;
double r122675 = r122673 * r122674;
double r122676 = r122675 * r122613;
double r122677 = r122585 / r122676;
double r122678 = 7.3873453313318566e+224;
bool r122679 = r122582 <= r122678;
double r122680 = 6.0;
double r122681 = pow(r122622, r122680);
double r122682 = r122681 * r122629;
double r122683 = r122682 / r122633;
double r122684 = r122621 * r122683;
double r122685 = r122585 * r122684;
double r122686 = r122637 * r122628;
double r122687 = r122586 * r122686;
double r122688 = r122687 / r122633;
double r122689 = r122621 * r122688;
double r122690 = r122685 + r122689;
double r122691 = -r122690;
double r122692 = r122585 / r122691;
double r122693 = r122593 * r122651;
double r122694 = r122693 * r122607;
double r122695 = cbrt(r122613);
double r122696 = r122695 * r122695;
double r122697 = r122694 * r122696;
double r122698 = r122697 * r122695;
double r122699 = r122585 / r122698;
double r122700 = r122679 ? r122692 : r122699;
double r122701 = r122672 ? r122677 : r122700;
double r122702 = r122670 ? r122644 : r122701;
double r122703 = r122659 ? r122668 : r122702;
double r122704 = r122648 ? r122655 : r122703;
double r122705 = r122617 ? r122644 : r122704;
double r122706 = r122584 ? r122615 : r122705;
return r122706;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -1.932296528705669e+146Initial program 34.9
rmApplied add-cube-cbrt34.9
Applied unpow-prod-down34.9
Applied times-frac27.8
Applied associate-*l*27.8
rmApplied add-cube-cbrt27.8
Applied unpow-prod-down27.8
Applied times-frac22.4
rmApplied associate-*l*22.4
rmApplied *-un-lft-identity22.4
Applied add-cube-cbrt22.4
Applied unpow-prod-down22.4
Applied times-frac22.4
Applied associate-*l*22.4
if -1.932296528705669e+146 < k < -2.4654898731068644e+69 or 496413.62095398596 < k < 2.402994973162251e+139Initial program 31.9
rmApplied add-cube-cbrt32.0
Applied unpow-prod-down32.0
Applied times-frac24.1
Applied associate-*l*24.1
Taylor expanded around -inf 18.8
if -2.4654898731068644e+69 < k < -3.515615117987835e-226Initial program 30.4
rmApplied add-cube-cbrt30.6
Applied unpow-prod-down30.6
Applied times-frac23.0
Applied associate-*l*20.3
rmApplied add-cube-cbrt20.3
Applied unpow-prod-down20.3
Applied times-frac13.4
rmApplied associate-*l*12.8
rmApplied associate-*l*9.9
if -3.515615117987835e-226 < k < 496413.62095398596Initial program 34.2
rmApplied add-cube-cbrt34.4
Applied unpow-prod-down34.4
Applied times-frac27.6
Applied associate-*l*21.4
rmApplied add-cube-cbrt21.4
Applied unpow-prod-down21.4
Applied times-frac16.3
rmApplied associate-*l*13.0
rmApplied pow1/339.1
Applied pow-pow12.8
Simplified12.8
if 2.402994973162251e+139 < k < 3.50673347844676e+202Initial program 30.6
rmApplied add-cube-cbrt30.7
Applied unpow-prod-down30.7
Applied times-frac24.6
Applied associate-*l*24.6
rmApplied add-cube-cbrt24.6
Applied unpow-prod-down24.6
Applied times-frac19.4
rmApplied associate-*l*19.4
if 3.50673347844676e+202 < k < 7.3873453313318566e+224Initial program 39.1
rmApplied add-cube-cbrt39.1
Applied unpow-prod-down39.1
Applied times-frac33.8
Applied associate-*l*33.8
rmApplied add-cube-cbrt33.8
Applied unpow-prod-down33.8
Applied times-frac27.5
Taylor expanded around -inf 33.3
if 7.3873453313318566e+224 < k Initial program 32.0
rmApplied add-cube-cbrt32.0
Applied unpow-prod-down32.0
Applied times-frac26.2
Applied associate-*l*26.2
rmApplied add-cube-cbrt26.2
Applied unpow-prod-down26.2
Applied times-frac21.3
rmApplied associate-*l*21.3
rmApplied add-cube-cbrt21.3
Applied associate-*r*21.3
Final simplification16.9
herbie shell --seed 2019303
(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))))