Toniolo and Linder, Equation (10-)

Average Error: 48.4 → 13.2

Time: 1.4m

Precision: 64

Math

\[\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 -3.309694515652526141729414587502190015342 \cdot 10^{-65} \lor \neg \left(\ell \le 3.29157748652316182286011727101396689325 \cdot 10^{146}\right):\\ \;\;\;\;2 \cdot \left(\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{1}{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}}{{\left(\sqrt[3]{1}\right)}^{2}}\right) \cdot \frac{\frac{\cos k}{\frac{1}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\\ \mathbf{else}:\\ \;\;\;\;2 \cdot \left({\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{\sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\right)\\ \end{array}\]

C

double f(double t, double l, double k) {
        double r305664 = 2.0;
        double r305665 = t;
        double r305666 = 3.0;
        double r305667 = pow(r305665, r305666);
        double r305668 = l;
        double r305669 = r305668 * r305668;
        double r305670 = r305667 / r305669;
        double r305671 = k;
        double r305672 = sin(r305671);
        double r305673 = r305670 * r305672;
        double r305674 = tan(r305671);
        double r305675 = r305673 * r305674;
        double r305676 = 1.0;
        double r305677 = r305671 / r305665;
        double r305678 = pow(r305677, r305664);
        double r305679 = r305676 + r305678;
        double r305680 = r305679 - r305676;
        double r305681 = r305675 * r305680;
        double r305682 = r305664 / r305681;
        return r305682;
}

↓

double f(double t, double l, double k) {
        double r305683 = l;
        double r305684 = -3.309694515652526e-65;
        bool r305685 = r305683 <= r305684;
        double r305686 = 3.291577486523162e+146;
        bool r305687 = r305683 <= r305686;
        double r305688 = !r305687;
        bool r305689 = r305685 || r305688;
        double r305690 = 2.0;
        double r305691 = 1.0;
        double r305692 = k;
        double r305693 = 2.0;
        double r305694 = r305690 / r305693;
        double r305695 = pow(r305692, r305694);
        double r305696 = t;
        double r305697 = 1.0;
        double r305698 = pow(r305696, r305697);
        double r305699 = r305695 * r305698;
        double r305700 = r305695 * r305699;
        double r305701 = r305691 / r305700;
        double r305702 = pow(r305701, r305697);
        double r305703 = sin(r305692);
        double r305704 = cbrt(r305703);
        double r305705 = 4.0;
        double r305706 = pow(r305704, r305705);
        double r305707 = r305706 / r305683;
        double r305708 = r305691 / r305707;
        double r305709 = cbrt(r305691);
        double r305710 = pow(r305709, r305693);
        double r305711 = r305708 / r305710;
        double r305712 = r305702 * r305711;
        double r305713 = cos(r305692);
        double r305714 = r305691 / r305683;
        double r305715 = r305713 / r305714;
        double r305716 = pow(r305704, r305693);
        double r305717 = r305715 / r305716;
        double r305718 = r305712 * r305717;
        double r305719 = r305690 * r305718;
        double r305720 = r305709 * r305709;
        double r305721 = r305720 / r305695;
        double r305722 = pow(r305721, r305697);
        double r305723 = r305709 / r305699;
        double r305724 = pow(r305723, r305697);
        double r305725 = r305707 / r305683;
        double r305726 = r305713 / r305725;
        double r305727 = r305726 / r305716;
        double r305728 = r305724 * r305727;
        double r305729 = r305722 * r305728;
        double r305730 = r305690 * r305729;
        double r305731 = r305689 ? r305719 : r305730;
        return r305731;
}

Error

Bits error versus t

Bits error versus l

Bits error versus k

Derivation

1. Split input into 2 regimes

2. if l < -3.309694515652526e-65 or 3.291577486523162e+146 < l

   1. Initial program 54.5

      \[\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)}\]

   2. Simplified 50.3

      \[\leadsto \color{blue}{\frac{2 \cdot \left(\ell \cdot \ell\right)}{\left({\left(\frac{k}{t}\right)}^{2} \cdot \left({t}^{3} \cdot \tan k\right)\right) \cdot \sin k}}\]

   3. Taylor expanded around inf 38.5

      \[\leadsto \color{blue}{2 \cdot \left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)}\]
   4. Using strategy rm

   5. Applied sqr-pow 38.5

      \[\leadsto 2 \cdot \left({\left(\frac{1}{\color{blue}{\left({k}^{\left(\frac{2}{2}\right)} \cdot {k}^{\left(\frac{2}{2}\right)}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)\]

   6. Applied associate-*l* 35.8

      \[\leadsto 2 \cdot \left({\left(\frac{1}{\color{blue}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)\]
   7. Using strategy rm

   8. Applied add-cube-cbrt 36.1

      \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\color{blue}{\left(\left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right) \cdot \sqrt[3]{\sin k}\right)}}^{2}}\right)\]

   9. Applied unpow-prod-down 36.1

      \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{\color{blue}{{\left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right)}^{2} \cdot {\left(\sqrt[3]{\sin k}\right)}^{2}}}\right)\]

   10. Applied associate-/r* 36.1

       \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \color{blue}{\frac{\frac{\cos k \cdot {\ell}^{2}}{{\left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right)}^{2}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}}\right)\]

   11. Simplified 36.1

       \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\color{blue}{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\]
   12. Using strategy rm

   13. Applied *-un-lft-identity 36.1

       \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}{{\left(\sqrt[3]{\color{blue}{1 \cdot \sin k}}\right)}^{2}}\right)\]

   14. Applied cbrt-prod 36.1

       \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}{{\color{blue}{\left(\sqrt[3]{1} \cdot \sqrt[3]{\sin k}\right)}}^{2}}\right)\]

   15. Applied unpow-prod-down 36.1

       \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}{\color{blue}{{\left(\sqrt[3]{1}\right)}^{2} \cdot {\left(\sqrt[3]{\sin k}\right)}^{2}}}\right)\]

   16. Applied div-inv 36.1

       \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{\cos k}{\color{blue}{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell} \cdot \frac{1}{\ell}}}}{{\left(\sqrt[3]{1}\right)}^{2} \cdot {\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\]

   17. Applied *-un-lft-identity 36.1

       \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{\color{blue}{1 \cdot \cos k}}{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell} \cdot \frac{1}{\ell}}}{{\left(\sqrt[3]{1}\right)}^{2} \cdot {\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\]

   18. Applied times-frac 36.1

       \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\color{blue}{\frac{1}{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}} \cdot \frac{\cos k}{\frac{1}{\ell}}}}{{\left(\sqrt[3]{1}\right)}^{2} \cdot {\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\]

   19. Applied times-frac 36.1

       \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \color{blue}{\left(\frac{\frac{1}{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}}{{\left(\sqrt[3]{1}\right)}^{2}} \cdot \frac{\frac{\cos k}{\frac{1}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)}\right)\]

   20. Applied associate-*r* 24.4

       \[\leadsto 2 \cdot \color{blue}{\left(\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{1}{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}}{{\left(\sqrt[3]{1}\right)}^{2}}\right) \cdot \frac{\frac{\cos k}{\frac{1}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)}\]

   if -3.309694515652526e-65 < l < 3.291577486523162e+146

   1. Initial program 45.5

      \[\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)}\]

   2. Simplified 36.7

      \[\leadsto \color{blue}{\frac{2 \cdot \left(\ell \cdot \ell\right)}{\left({\left(\frac{k}{t}\right)}^{2} \cdot \left({t}^{3} \cdot \tan k\right)\right) \cdot \sin k}}\]

   3. Taylor expanded around inf 14.9

      \[\leadsto \color{blue}{2 \cdot \left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)}\]
   4. Using strategy rm

   5. Applied sqr-pow 14.9

      \[\leadsto 2 \cdot \left({\left(\frac{1}{\color{blue}{\left({k}^{\left(\frac{2}{2}\right)} \cdot {k}^{\left(\frac{2}{2}\right)}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)\]

   6. Applied associate-*l* 13.0

      \[\leadsto 2 \cdot \left({\left(\frac{1}{\color{blue}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)\]
   7. Using strategy rm

   8. Applied add-cube-cbrt 13.3

      \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\color{blue}{\left(\left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right) \cdot \sqrt[3]{\sin k}\right)}}^{2}}\right)\]

   9. Applied unpow-prod-down 13.3

      \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{\color{blue}{{\left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right)}^{2} \cdot {\left(\sqrt[3]{\sin k}\right)}^{2}}}\right)\]

   10. Applied associate-/r* 12.8

       \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \color{blue}{\frac{\frac{\cos k \cdot {\ell}^{2}}{{\left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right)}^{2}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}}\right)\]

   11. Simplified 10.1

       \[\leadsto 2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\color{blue}{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\]
   12. Using strategy rm

   13. Applied add-cube-cbrt 10.1

       \[\leadsto 2 \cdot \left({\left(\frac{\color{blue}{\left(\sqrt[3]{1} \cdot \sqrt[3]{1}\right) \cdot \sqrt[3]{1}}}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\]

   14. Applied times-frac 10.0

       \[\leadsto 2 \cdot \left({\color{blue}{\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)}} \cdot \frac{\sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}}^{1} \cdot \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\]

   15. Applied unpow-prod-down 10.0

       \[\leadsto 2 \cdot \left(\color{blue}{\left({\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot {\left(\frac{\sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1}\right)} \cdot \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\]

   16. Applied associate-*l* 7.9

       \[\leadsto 2 \cdot \color{blue}{\left({\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{\sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\right)}\]
3. Recombined 2 regimes into one program.

4. Final simplification 13.2

   \[\leadsto \begin{array}{l} \mathbf{if}\;\ell \le -3.309694515652526141729414587502190015342 \cdot 10^{-65} \lor \neg \left(\ell \le 3.29157748652316182286011727101396689325 \cdot 10^{146}\right):\\ \;\;\;\;2 \cdot \left(\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{1}{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}}{{\left(\sqrt[3]{1}\right)}^{2}}\right) \cdot \frac{\frac{\cos k}{\frac{1}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\\ \mathbf{else}:\\ \;\;\;\;2 \cdot \left({\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{\sqrt[3]{1}}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\right)\\ \end{array}\]

Reproduce

herbie shell --seed 2019346 
(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))))