Toniolo and Linder, Equation (10-)

Average Error: 48.3 → 16.0

Time: 27.1s

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 4.18838592989303695132438800625329553817 \cdot 10^{-248}:\\ \;\;\;\;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 5.603600301497313251444735398660715198036 \cdot 10^{283}:\\ \;\;\;\;2 \cdot \left(\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) \cdot {\left(\frac{\sqrt{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1}\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 r92106 = 2.0;
        double r92107 = t;
        double r92108 = 3.0;
        double r92109 = pow(r92107, r92108);
        double r92110 = l;
        double r92111 = r92110 * r92110;
        double r92112 = r92109 / r92111;
        double r92113 = k;
        double r92114 = sin(r92113);
        double r92115 = r92112 * r92114;
        double r92116 = tan(r92113);
        double r92117 = r92115 * r92116;
        double r92118 = 1.0;
        double r92119 = r92113 / r92107;
        double r92120 = pow(r92119, r92106);
        double r92121 = r92118 + r92120;
        double r92122 = r92121 - r92118;
        double r92123 = r92117 * r92122;
        double r92124 = r92106 / r92123;
        return r92124;
}

↓

double f(double t, double l, double k) {
        double r92125 = l;
        double r92126 = r92125 * r92125;
        double r92127 = 4.188385929893037e-248;
        bool r92128 = r92126 <= r92127;
        double r92129 = 2.0;
        double r92130 = 1.0;
        double r92131 = k;
        double r92132 = 2.0;
        double r92133 = r92129 / r92132;
        double r92134 = pow(r92131, r92133);
        double r92135 = t;
        double r92136 = 1.0;
        double r92137 = pow(r92135, r92136);
        double r92138 = r92134 * r92137;
        double r92139 = r92134 * r92138;
        double r92140 = r92130 / r92139;
        double r92141 = pow(r92140, r92136);
        double r92142 = cos(r92131);
        double r92143 = sin(r92131);
        double r92144 = cbrt(r92143);
        double r92145 = 4.0;
        double r92146 = pow(r92144, r92145);
        double r92147 = r92146 / r92125;
        double r92148 = r92147 / r92125;
        double r92149 = r92142 / r92148;
        double r92150 = pow(r92144, r92132);
        double r92151 = r92149 / r92150;
        double r92152 = r92141 * r92151;
        double r92153 = r92129 * r92152;
        double r92154 = 5.603600301497313e+283;
        bool r92155 = r92126 <= r92154;
        double r92156 = r92130 / r92138;
        double r92157 = pow(r92156, r92136);
        double r92158 = pow(r92125, r92132);
        double r92159 = r92142 * r92158;
        double r92160 = pow(r92143, r92132);
        double r92161 = r92159 / r92160;
        double r92162 = r92157 * r92161;
        double r92163 = sqrt(r92130);
        double r92164 = r92163 / r92134;
        double r92165 = pow(r92164, r92136);
        double r92166 = r92162 * r92165;
        double r92167 = r92129 * r92166;
        double r92168 = pow(r92131, r92129);
        double r92169 = r92137 * r92168;
        double r92170 = r92130 / r92169;
        double r92171 = pow(r92170, r92136);
        double r92172 = exp(r92171);
        double r92173 = pow(r92172, r92161);
        double r92174 = log(r92173);
        double r92175 = r92129 * r92174;
        double r92176 = r92155 ? r92167 : r92175;
        double r92177 = r92128 ? r92153 : r92176;
        return r92177;
}

Error

Bits error versus t

Bits error versus l

Bits error versus k

Derivation

1. Split input into 3 regimes

2. if (* l l) < 4.188385929893037e-248

   1. Initial program 46.3

      \[\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.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 17.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 17.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* 17.9

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

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

      \[\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* 17.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 \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 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 \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 4.188385929893037e-248 < (* l l) < 5.603600301497313e+283

   1. Initial program 44.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 34.5

      \[\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.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 11.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* 6.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 6.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 6.4

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

       \[\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* 2.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. Simplified 2.9

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

   14. Applied *-commutative 2.9

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

   if 5.603600301497313e+283 < (* l l)

   1. Initial program 63.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 62.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 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.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-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.8

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

   \[\leadsto \begin{array}{l} \mathbf{if}\;\ell \cdot \ell \le 4.18838592989303695132438800625329553817 \cdot 10^{-248}:\\ \;\;\;\;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 5.603600301497313251444735398660715198036 \cdot 10^{283}:\\ \;\;\;\;2 \cdot \left(\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) \cdot {\left(\frac{\sqrt{1}}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1}\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 2019354 
(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))))