Toniolo and Linder, Equation (10-)

Average Error: 48.6 → 16.0

Time: 29.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 -4.767926568016393126457677186561838999712 \cdot 10^{-144}:\\ \;\;\;\;2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left(\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \left(\cos k \cdot {\ell}^{2}\right)\right) \cdot \frac{1}{{\left(\sin k\right)}^{2}}\right)\right)\\ \mathbf{elif}\;\ell \le 7.032624916942270176139110895833513818529 \cdot 10^{-135}:\\ \;\;\;\;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{else}:\\ \;\;\;\;2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left(\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\cos k}{\sin k}\right) \cdot \frac{{\ell}^{2}}{{\left(\sin k\right)}^{1}}\right)\right)\\ \end{array}\]

C

double f(double t, double l, double k) {
        double r105318 = 2.0;
        double r105319 = t;
        double r105320 = 3.0;
        double r105321 = pow(r105319, r105320);
        double r105322 = l;
        double r105323 = r105322 * r105322;
        double r105324 = r105321 / r105323;
        double r105325 = k;
        double r105326 = sin(r105325);
        double r105327 = r105324 * r105326;
        double r105328 = tan(r105325);
        double r105329 = r105327 * r105328;
        double r105330 = 1.0;
        double r105331 = r105325 / r105319;
        double r105332 = pow(r105331, r105318);
        double r105333 = r105330 + r105332;
        double r105334 = r105333 - r105330;
        double r105335 = r105329 * r105334;
        double r105336 = r105318 / r105335;
        return r105336;
}

↓

double f(double t, double l, double k) {
        double r105337 = l;
        double r105338 = -4.767926568016393e-144;
        bool r105339 = r105337 <= r105338;
        double r105340 = 2.0;
        double r105341 = 1.0;
        double r105342 = k;
        double r105343 = 2.0;
        double r105344 = r105340 / r105343;
        double r105345 = pow(r105342, r105344);
        double r105346 = r105341 / r105345;
        double r105347 = 1.0;
        double r105348 = pow(r105346, r105347);
        double r105349 = t;
        double r105350 = pow(r105349, r105347);
        double r105351 = r105345 * r105350;
        double r105352 = r105341 / r105351;
        double r105353 = pow(r105352, r105347);
        double r105354 = cos(r105342);
        double r105355 = pow(r105337, r105343);
        double r105356 = r105354 * r105355;
        double r105357 = r105353 * r105356;
        double r105358 = sin(r105342);
        double r105359 = pow(r105358, r105343);
        double r105360 = r105341 / r105359;
        double r105361 = r105357 * r105360;
        double r105362 = r105348 * r105361;
        double r105363 = r105340 * r105362;
        double r105364 = 7.03262491694227e-135;
        bool r105365 = r105337 <= r105364;
        double r105366 = r105345 * r105351;
        double r105367 = r105341 / r105366;
        double r105368 = pow(r105367, r105347);
        double r105369 = cbrt(r105358);
        double r105370 = 4.0;
        double r105371 = pow(r105369, r105370);
        double r105372 = r105371 / r105337;
        double r105373 = r105372 / r105337;
        double r105374 = r105354 / r105373;
        double r105375 = pow(r105369, r105343);
        double r105376 = r105374 / r105375;
        double r105377 = r105368 * r105376;
        double r105378 = r105340 * r105377;
        double r105379 = r105354 / r105358;
        double r105380 = r105353 * r105379;
        double r105381 = pow(r105358, r105341);
        double r105382 = r105355 / r105381;
        double r105383 = r105380 * r105382;
        double r105384 = r105348 * r105383;
        double r105385 = r105340 * r105384;
        double r105386 = r105365 ? r105378 : r105385;
        double r105387 = r105339 ? r105363 : r105386;
        return r105387;
}

Error

Bits error versus t

Bits error versus l

Bits error versus k

Derivation

1. Split input into 3 regimes

2. if l < -4.767926568016393e-144

   1. Initial program 49.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 42.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 23.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 23.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* 20.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 *-un-lft-identity 20.0

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

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

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

       \[\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 div-inv 17.3

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

   14. Applied associate-*r* 17.5

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

   if -4.767926568016393e-144 < l < 7.03262491694227e-135

   1. Initial program 46.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 36.2

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

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

      \[\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.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-cube-cbrt 17.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}{\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.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}{{\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 12.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 \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 7.03262491694227e-135 < l

   1. Initial program 50.8

      \[\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 43.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 26.0

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

      \[\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* 22.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 *-un-lft-identity 22.2

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

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

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

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

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

   14. Applied times-frac 19.3

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

   15. Applied associate-*r* 19.1

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

   16. Simplified 19.1

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

4. Final simplification 16.0

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

Reproduce

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