Toniolo and Linder, Equation (10-)

Average Error: 48.4 → 12.9

Time: 53.3s

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 36341823884774118312709128287999229952:\\ \;\;\;\;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(\left(\frac{\cos k}{\sin k} \cdot \ell\right) \cdot \frac{\ell}{\sin k}\right)\right)\\ \mathbf{elif}\;\ell \cdot \ell \le 2.702815358468802609888516667751806875166 \cdot 10^{307}:\\ \;\;\;\;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 \left(\frac{\cos k}{\sin k} \cdot \frac{{\ell}^{2}}{\sin k}\right)\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{3}}{\ell} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \sin k\right) \cdot \tan k}}{{\left(\frac{k}{t}\right)}^{2}}\\ \end{array}\]

C

double f(double t, double l, double k) {
        double r91130 = 2.0;
        double r91131 = t;
        double r91132 = 3.0;
        double r91133 = pow(r91131, r91132);
        double r91134 = l;
        double r91135 = r91134 * r91134;
        double r91136 = r91133 / r91135;
        double r91137 = k;
        double r91138 = sin(r91137);
        double r91139 = r91136 * r91138;
        double r91140 = tan(r91137);
        double r91141 = r91139 * r91140;
        double r91142 = 1.0;
        double r91143 = r91137 / r91131;
        double r91144 = pow(r91143, r91130);
        double r91145 = r91142 + r91144;
        double r91146 = r91145 - r91142;
        double r91147 = r91141 * r91146;
        double r91148 = r91130 / r91147;
        return r91148;
}

↓

double f(double t, double l, double k) {
        double r91149 = l;
        double r91150 = r91149 * r91149;
        double r91151 = 3.634182388477412e+37;
        bool r91152 = r91150 <= r91151;
        double r91153 = 2.0;
        double r91154 = 1.0;
        double r91155 = k;
        double r91156 = 2.0;
        double r91157 = r91153 / r91156;
        double r91158 = pow(r91155, r91157);
        double r91159 = t;
        double r91160 = 1.0;
        double r91161 = pow(r91159, r91160);
        double r91162 = r91158 * r91161;
        double r91163 = r91158 * r91162;
        double r91164 = r91154 / r91163;
        double r91165 = pow(r91164, r91160);
        double r91166 = cos(r91155);
        double r91167 = sin(r91155);
        double r91168 = r91166 / r91167;
        double r91169 = r91168 * r91149;
        double r91170 = r91149 / r91167;
        double r91171 = r91169 * r91170;
        double r91172 = r91165 * r91171;
        double r91173 = r91153 * r91172;
        double r91174 = 2.7028153584688026e+307;
        bool r91175 = r91150 <= r91174;
        double r91176 = r91154 / r91158;
        double r91177 = pow(r91176, r91160);
        double r91178 = r91154 / r91162;
        double r91179 = pow(r91178, r91160);
        double r91180 = pow(r91149, r91156);
        double r91181 = r91180 / r91167;
        double r91182 = r91168 * r91181;
        double r91183 = r91179 * r91182;
        double r91184 = r91177 * r91183;
        double r91185 = r91153 * r91184;
        double r91186 = cbrt(r91159);
        double r91187 = r91186 * r91186;
        double r91188 = 3.0;
        double r91189 = pow(r91187, r91188);
        double r91190 = r91189 / r91149;
        double r91191 = pow(r91186, r91188);
        double r91192 = r91191 / r91149;
        double r91193 = r91190 * r91192;
        double r91194 = r91193 * r91167;
        double r91195 = tan(r91155);
        double r91196 = r91194 * r91195;
        double r91197 = r91153 / r91196;
        double r91198 = r91155 / r91159;
        double r91199 = pow(r91198, r91153);
        double r91200 = r91197 / r91199;
        double r91201 = r91175 ? r91185 : r91200;
        double r91202 = r91152 ? r91173 : r91201;
        return r91202;
}

Error

Bits error versus t

Bits error versus l

Bits error versus k

Derivation

1. Split input into 3 regimes

2. if (* l l) < 3.634182388477412e+37

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

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

   3. Taylor expanded around inf 12.4

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

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

   9. Applied unpow-prod-down 37.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}{{\left(\sqrt{\sin k}\right)}^{2} \cdot {\left(\sqrt{\sin k}\right)}^{2}}}\right)\]

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

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

   12. Simplified 10.8

       \[\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}{\sin k} \cdot \color{blue}{\frac{{\ell}^{2}}{\sin k}}\right)\right)\]
   13. Using strategy rm

   14. Applied *-un-lft-identity 10.8

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

   15. Applied unpow2 10.8

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

   16. Applied times-frac 8.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(\frac{\cos k}{\sin k} \cdot \color{blue}{\left(\frac{\ell}{1} \cdot \frac{\ell}{\sin k}\right)}\right)\right)\]

   17. Applied associate-*r* 6.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(\left(\frac{\cos k}{\sin k} \cdot \frac{\ell}{1}\right) \cdot \frac{\ell}{\sin k}\right)}\right)\]

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

   if 3.634182388477412e+37 < (* l l) < 2.7028153584688026e+307

   1. Initial program 48.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 41.1

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

   3. Taylor expanded around inf 21.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 21.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* 15.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 add-sqr-sqrt 39.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 \frac{\cos k \cdot {\ell}^{2}}{{\color{blue}{\left(\sqrt{\sin k} \cdot \sqrt{\sin k}\right)}}^{2}}\right)\]

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

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

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

   12. Simplified 15.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}{\sin k} \cdot \color{blue}{\frac{{\ell}^{2}}{\sin k}}\right)\right)\]
   13. Using strategy rm

   14. Applied *-un-lft-identity 15.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 \left(\frac{\cos k}{\sin k} \cdot \frac{{\ell}^{2}}{\sin k}\right)\right)\]

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

   16. Applied unpow-prod-down 14.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 \left(\frac{\cos k}{\sin k} \cdot \frac{{\ell}^{2}}{\sin k}\right)\right)\]

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

   if 2.7028153584688026e+307 < (* l l)

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

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

   4. Applied add-cube-cbrt 63.9

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

   5. Applied unpow-prod-down 63.9

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

   6. Applied times-frac 48.9

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

4. Final simplification 12.9

   \[\leadsto \begin{array}{l} \mathbf{if}\;\ell \cdot \ell \le 36341823884774118312709128287999229952:\\ \;\;\;\;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(\left(\frac{\cos k}{\sin k} \cdot \ell\right) \cdot \frac{\ell}{\sin k}\right)\right)\\ \mathbf{elif}\;\ell \cdot \ell \le 2.702815358468802609888516667751806875166 \cdot 10^{307}:\\ \;\;\;\;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 \left(\frac{\cos k}{\sin k} \cdot \frac{{\ell}^{2}}{\sin k}\right)\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{3}}{\ell} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \sin k\right) \cdot \tan k}}{{\left(\frac{k}{t}\right)}^{2}}\\ \end{array}\]

Reproduce

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