Toniolo and Linder, Equation (10-)

Average Error: 48.6 → 16.2

Time: 28.4s

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 9.9737193034107518 \cdot 10^{-278}:\\ \;\;\;\;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]{\sin k}\right)}^{2}}\right)\\ \mathbf{elif}\;\ell \cdot \ell \le 3.40947516322032403 \cdot 10^{289}:\\ \;\;\;\;2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{\frac{1}{{k}^{\left(\frac{2}{2}\right)}}}{{t}^{1}}\right)}^{1} \cdot \frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)\right)\\ \mathbf{else}:\\ \;\;\;\;2 \cdot \log \left({\left(e^{{\left(\frac{1}{{t}^{1} \cdot {k}^{2}}\right)}^{1}}\right)}^{\left(\frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)}\right)\\ \end{array}\]

C

double f(double t, double l, double k) {
        double r97012 = 2.0;
        double r97013 = t;
        double r97014 = 3.0;
        double r97015 = pow(r97013, r97014);
        double r97016 = l;
        double r97017 = r97016 * r97016;
        double r97018 = r97015 / r97017;
        double r97019 = k;
        double r97020 = sin(r97019);
        double r97021 = r97018 * r97020;
        double r97022 = tan(r97019);
        double r97023 = r97021 * r97022;
        double r97024 = 1.0;
        double r97025 = r97019 / r97013;
        double r97026 = pow(r97025, r97012);
        double r97027 = r97024 + r97026;
        double r97028 = r97027 - r97024;
        double r97029 = r97023 * r97028;
        double r97030 = r97012 / r97029;
        return r97030;
}

↓

double f(double t, double l, double k) {
        double r97031 = l;
        double r97032 = r97031 * r97031;
        double r97033 = 9.973719303410752e-278;
        bool r97034 = r97032 <= r97033;
        double r97035 = 2.0;
        double r97036 = 1.0;
        double r97037 = k;
        double r97038 = 2.0;
        double r97039 = r97035 / r97038;
        double r97040 = pow(r97037, r97039);
        double r97041 = t;
        double r97042 = 1.0;
        double r97043 = pow(r97041, r97042);
        double r97044 = r97040 * r97043;
        double r97045 = r97040 * r97044;
        double r97046 = r97036 / r97045;
        double r97047 = pow(r97046, r97042);
        double r97048 = cos(r97037);
        double r97049 = sin(r97037);
        double r97050 = cbrt(r97049);
        double r97051 = 4.0;
        double r97052 = pow(r97050, r97051);
        double r97053 = r97052 / r97031;
        double r97054 = r97053 / r97031;
        double r97055 = r97048 / r97054;
        double r97056 = pow(r97050, r97038);
        double r97057 = r97055 / r97056;
        double r97058 = r97047 * r97057;
        double r97059 = r97035 * r97058;
        double r97060 = 3.409475163220324e+289;
        bool r97061 = r97032 <= r97060;
        double r97062 = r97036 / r97040;
        double r97063 = pow(r97062, r97042);
        double r97064 = r97062 / r97043;
        double r97065 = pow(r97064, r97042);
        double r97066 = pow(r97031, r97038);
        double r97067 = r97048 * r97066;
        double r97068 = pow(r97049, r97038);
        double r97069 = r97067 / r97068;
        double r97070 = r97065 * r97069;
        double r97071 = r97063 * r97070;
        double r97072 = r97035 * r97071;
        double r97073 = pow(r97037, r97035);
        double r97074 = r97043 * r97073;
        double r97075 = r97036 / r97074;
        double r97076 = pow(r97075, r97042);
        double r97077 = exp(r97076);
        double r97078 = pow(r97077, r97069);
        double r97079 = log(r97078);
        double r97080 = r97035 * r97079;
        double r97081 = r97061 ? r97072 : r97080;
        double r97082 = r97034 ? r97059 : r97081;
        return r97082;
}

Error

Bits error versus t

Bits error versus l

Bits error versus k

Derivation

1. Split input into 3 regimes

2. if (* l l) < 9.973719303410752e-278

   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 18.6

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

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

      \[\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 18.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 \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 18.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 \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* 18.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}{\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.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)\]

   if 9.973719303410752e-278 < (* l l) < 3.409475163220324e+289

   1. Initial program 45.0

      \[\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 35.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 10.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 10.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* 6.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 *-un-lft-identity 6.5

      \[\leadsto 2 \cdot \left({\left(\frac{\color{blue}{1 \cdot 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 6.3

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

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

       \[\leadsto 2 \cdot \color{blue}{\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{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 associate-/r* 3.1

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

   if 3.409475163220324e+289 < (* l l)

   1. Initial program 63.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 62.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 61.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 61.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* 61.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-log-exp 63.3

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

   9. Simplified 58.6

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

4. Final simplification 16.2

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

Reproduce

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