Toniolo and Linder, Equation (10-)

Average Error: 48.1 → 12.7

Time: 24.8s

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 \cdot \ell \le 2.6145740768465546 \cdot 10^{-212}:\\ \;\;\;\;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}{\left|\sin k\right|} \cdot \ell}{\frac{\left|\sin k\right|}{\ell}}\right)\\ \mathbf{elif}\;\ell \cdot \ell \le 8.3711679961397521 \cdot 10^{265}:\\ \;\;\;\;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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \frac{\ell}{\frac{\left|\sin k\right|}{\ell}}\right)\right)\right)\\ \mathbf{else}:\\ \;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\left|\sin k\right| \cdot \frac{\left|\sin k\right|}{\ell}}\\ \end{array}\]

C

double f(double t, double l, double k) {
        double r79766 = 2.0;
        double r79767 = t;
        double r79768 = 3.0;
        double r79769 = pow(r79767, r79768);
        double r79770 = l;
        double r79771 = r79770 * r79770;
        double r79772 = r79769 / r79771;
        double r79773 = k;
        double r79774 = sin(r79773);
        double r79775 = r79772 * r79774;
        double r79776 = tan(r79773);
        double r79777 = r79775 * r79776;
        double r79778 = 1.0;
        double r79779 = r79773 / r79767;
        double r79780 = pow(r79779, r79766);
        double r79781 = r79778 + r79780;
        double r79782 = r79781 - r79778;
        double r79783 = r79777 * r79782;
        double r79784 = r79766 / r79783;
        return r79784;
}

↓

double f(double t, double l, double k) {
        double r79785 = l;
        double r79786 = r79785 * r79785;
        double r79787 = 2.6145740768465546e-212;
        bool r79788 = r79786 <= r79787;
        double r79789 = 2.0;
        double r79790 = 1.0;
        double r79791 = k;
        double r79792 = 2.0;
        double r79793 = r79789 / r79792;
        double r79794 = pow(r79791, r79793);
        double r79795 = t;
        double r79796 = 1.0;
        double r79797 = pow(r79795, r79796);
        double r79798 = r79794 * r79797;
        double r79799 = r79794 * r79798;
        double r79800 = r79790 / r79799;
        double r79801 = pow(r79800, r79796);
        double r79802 = cos(r79791);
        double r79803 = sin(r79791);
        double r79804 = fabs(r79803);
        double r79805 = r79802 / r79804;
        double r79806 = r79805 * r79785;
        double r79807 = r79804 / r79785;
        double r79808 = r79806 / r79807;
        double r79809 = r79801 * r79808;
        double r79810 = r79789 * r79809;
        double r79811 = 8.371167996139752e+265;
        bool r79812 = r79786 <= r79811;
        double r79813 = cbrt(r79790);
        double r79814 = r79813 * r79813;
        double r79815 = r79814 / r79794;
        double r79816 = pow(r79815, r79796);
        double r79817 = r79813 / r79798;
        double r79818 = pow(r79817, r79796);
        double r79819 = r79785 / r79807;
        double r79820 = r79805 * r79819;
        double r79821 = r79818 * r79820;
        double r79822 = r79816 * r79821;
        double r79823 = r79789 * r79822;
        double r79824 = pow(r79791, r79789);
        double r79825 = r79824 * r79797;
        double r79826 = r79790 / r79825;
        double r79827 = pow(r79826, r79796);
        double r79828 = r79827 * r79785;
        double r79829 = r79828 * r79802;
        double r79830 = r79804 * r79807;
        double r79831 = r79829 / r79830;
        double r79832 = r79789 * r79831;
        double r79833 = r79812 ? r79823 : r79832;
        double r79834 = r79788 ? r79810 : r79833;
        return r79834;
}

Error

Bits error versus t

Bits error versus l

Bits error versus k

Derivation

1. Split input into 3 regimes

2. if (* l l) < 2.6145740768465546e-212

   1. Initial program 46.2

      \[\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.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 16.3

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

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

      \[\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 16.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}{\sqrt{{\left(\sin k\right)}^{2}} \cdot \sqrt{{\left(\sin k\right)}^{2}}}}\right)\]

   9. Applied times-frac 16.2

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

   10. Simplified 16.2

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

   11. Simplified 11.6

       \[\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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \color{blue}{\frac{\ell}{\frac{\left|\sin k\right|}{\ell}}}\right)\right)\]
   12. Using strategy rm

   13. Applied associate-*r/ 8.6

       \[\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}{\left|\sin k\right|} \cdot \ell}{\frac{\left|\sin k\right|}{\ell}}}\right)\]

   if 2.6145740768465546e-212 < (* l l) < 8.371167996139752e+265

   1. Initial program 43.9

      \[\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 34.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 11.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 11.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* 7.1

      \[\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 7.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}{\sqrt{{\left(\sin k\right)}^{2}} \cdot \sqrt{{\left(\sin k\right)}^{2}}}}\right)\]

   9. Applied times-frac 7.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{\cos k}{\sqrt{{\left(\sin k\right)}^{2}}} \cdot \frac{{\ell}^{2}}{\sqrt{{\left(\sin k\right)}^{2}}}\right)}\right)\]

   10. Simplified 7.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 \left(\color{blue}{\frac{\cos k}{\left|\sin k\right|}} \cdot \frac{{\ell}^{2}}{\sqrt{{\left(\sin k\right)}^{2}}}\right)\right)\]

   11. Simplified 7.0

       \[\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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \color{blue}{\frac{\ell}{\frac{\left|\sin k\right|}{\ell}}}\right)\right)\]
   12. Using strategy rm

   13. Applied add-cube-cbrt 7.0

       \[\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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \frac{\ell}{\frac{\left|\sin k\right|}{\ell}}\right)\right)\]

   14. Applied times-frac 6.4

       \[\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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \frac{\ell}{\frac{\left|\sin k\right|}{\ell}}\right)\right)\]

   15. Applied unpow-prod-down 6.4

       \[\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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \frac{\ell}{\frac{\left|\sin k\right|}{\ell}}\right)\right)\]

   16. Applied associate-*l* 2.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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \frac{\ell}{\frac{\left|\sin k\right|}{\ell}}\right)\right)\right)}\]

   if 8.371167996139752e+265 < (* l l)

   1. Initial program 61.7

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

      \[\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 58.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 58.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* 57.7

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

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

   9. Applied times-frac 57.7

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

   10. Simplified 57.7

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

   11. Simplified 57.7

       \[\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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \color{blue}{\frac{\ell}{\frac{\left|\sin k\right|}{\ell}}}\right)\right)\]
   12. Using strategy rm

   13. Applied frac-times 57.7

       \[\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{\cos k \cdot \ell}{\left|\sin k\right| \cdot \frac{\left|\sin k\right|}{\ell}}}\right)\]

   14. Applied associate-*r/ 36.7

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

   15. Simplified 44.3

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

4. Final simplification 12.7

   \[\leadsto \begin{array}{l} \mathbf{if}\;\ell \cdot \ell \le 2.6145740768465546 \cdot 10^{-212}:\\ \;\;\;\;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}{\left|\sin k\right|} \cdot \ell}{\frac{\left|\sin k\right|}{\ell}}\right)\\ \mathbf{elif}\;\ell \cdot \ell \le 8.3711679961397521 \cdot 10^{265}:\\ \;\;\;\;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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \frac{\ell}{\frac{\left|\sin k\right|}{\ell}}\right)\right)\right)\\ \mathbf{else}:\\ \;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\left|\sin k\right| \cdot \frac{\left|\sin k\right|}{\ell}}\\ \end{array}\]

Reproduce

herbie shell --seed 2020056 +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))))