Toniolo and Linder, Equation (10-)

Average Error: 47.9 → 16.2

Time: 51.0s

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 \frac{-1031436374413907}{3.685510180489786476798393145496356338786 \cdot 10^{165}} \lor \neg \left(\ell \le \frac{1273299905791625}{1.207667975942893232717297465916246045093 \cdot 10^{170}}\right):\\ \;\;\;\;2 \cdot \left({\left(\frac{\sqrt{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{\sqrt{1}}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\frac{\cos k}{\frac{\sin k}{{\ell}^{2}}}}{\sin k}\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{2 \cdot {\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1}}{\frac{\frac{\frac{{\left({\left(\sqrt[3]{\sin k}\right)}^{2}\right)}^{3}}{\ell}}{\ell}}{\cos k}}\\ \end{array}\]

C

double f(double t, double l, double k) {
        double r82724 = 2.0;
        double r82725 = t;
        double r82726 = 3.0;
        double r82727 = pow(r82725, r82726);
        double r82728 = l;
        double r82729 = r82728 * r82728;
        double r82730 = r82727 / r82729;
        double r82731 = k;
        double r82732 = sin(r82731);
        double r82733 = r82730 * r82732;
        double r82734 = tan(r82731);
        double r82735 = r82733 * r82734;
        double r82736 = 1.0;
        double r82737 = r82731 / r82725;
        double r82738 = pow(r82737, r82724);
        double r82739 = r82736 + r82738;
        double r82740 = r82739 - r82736;
        double r82741 = r82735 * r82740;
        double r82742 = r82724 / r82741;
        return r82742;
}

↓

double f(double t, double l, double k) {
        double r82743 = l;
        double r82744 = -1031436374413907.0;
        double r82745 = 3.6855101804897865e+165;
        double r82746 = r82744 / r82745;
        bool r82747 = r82743 <= r82746;
        double r82748 = 1273299905791625.0;
        double r82749 = 1.2076679759428932e+170;
        double r82750 = r82748 / r82749;
        bool r82751 = r82743 <= r82750;
        double r82752 = !r82751;
        bool r82753 = r82747 || r82752;
        double r82754 = 2.0;
        double r82755 = 1.0;
        double r82756 = sqrt(r82755);
        double r82757 = k;
        double r82758 = 2.0;
        double r82759 = r82754 / r82758;
        double r82760 = pow(r82757, r82759);
        double r82761 = r82756 / r82760;
        double r82762 = 1.0;
        double r82763 = pow(r82761, r82762);
        double r82764 = t;
        double r82765 = pow(r82764, r82762);
        double r82766 = r82760 * r82765;
        double r82767 = r82756 / r82766;
        double r82768 = pow(r82767, r82762);
        double r82769 = cos(r82757);
        double r82770 = sin(r82757);
        double r82771 = pow(r82743, r82758);
        double r82772 = r82770 / r82771;
        double r82773 = r82769 / r82772;
        double r82774 = r82773 / r82770;
        double r82775 = r82768 * r82774;
        double r82776 = r82763 * r82775;
        double r82777 = r82754 * r82776;
        double r82778 = r82760 * r82766;
        double r82779 = r82755 / r82778;
        double r82780 = pow(r82779, r82762);
        double r82781 = r82754 * r82780;
        double r82782 = cbrt(r82770);
        double r82783 = pow(r82782, r82758);
        double r82784 = 3.0;
        double r82785 = pow(r82783, r82784);
        double r82786 = r82785 / r82743;
        double r82787 = r82786 / r82743;
        double r82788 = r82787 / r82769;
        double r82789 = r82781 / r82788;
        double r82790 = r82753 ? r82777 : r82789;
        return r82790;
}

Error

Bits error versus t

Bits error versus l

Bits error versus k

Derivation

1. Split input into 3 regimes

2. if l < -2.7986257638741188e-151

   1. Initial program 48.1

      \[\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 40.9

      \[\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 22.8

      \[\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 22.8

      \[\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* 19.2

      \[\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-sqr-sqrt 19.2

      \[\leadsto 2 \cdot \left({\left(\frac{\color{blue}{\sqrt{1} \cdot \sqrt{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}}{{\left(\sin k\right)}^{2}}\right)\]

   9. Applied times-frac 18.9

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

   10. Applied unpow-prod-down 18.9

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

   11. Applied associate-*l* 16.9

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

   13. Applied unpow2 16.9

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

   14. Applied times-frac 16.8

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

   15. Applied associate-*r* 16.7

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

   16. Simplified 16.7

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

   if -2.7986257638741188e-151 < l < 1.0543460049915535e-155

   1. Initial program 46.1

      \[\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 37.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 19.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 19.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* 19.5

      \[\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 19.5

      \[\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 19.5

      \[\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* 19.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 \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 13.4

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

   if 1.0543460049915535e-155 < l

   1. Initial program 49.4

      \[\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 41.8

      \[\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 22.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 22.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* 19.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-sqr-sqrt 19.8

      \[\leadsto 2 \cdot \left({\left(\frac{\color{blue}{\sqrt{1} \cdot \sqrt{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}}{{\left(\sin k\right)}^{2}}\right)\]

   9. Applied times-frac 19.6

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

   10. Applied unpow-prod-down 19.6

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

   11. Applied associate-*l* 17.8

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

   13. Applied unpow2 17.8

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

   14. Applied associate-/r* 17.1

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

   15. Simplified 17.1

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

4. Final simplification 16.2

   \[\leadsto \begin{array}{l} \mathbf{if}\;\ell \le \frac{-1031436374413907}{3.685510180489786476798393145496356338786 \cdot 10^{165}} \lor \neg \left(\ell \le \frac{1273299905791625}{1.207667975942893232717297465916246045093 \cdot 10^{170}}\right):\\ \;\;\;\;2 \cdot \left({\left(\frac{\sqrt{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{\sqrt{1}}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\frac{\cos k}{\frac{\sin k}{{\ell}^{2}}}}{\sin k}\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{2 \cdot {\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1}}{\frac{\frac{\frac{{\left({\left(\sqrt[3]{\sin k}\right)}^{2}\right)}^{3}}{\ell}}{\ell}}{\cos k}}\\ \end{array}\]

Reproduce

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))))