Toniolo and Linder, Equation (10-)

Average Error: 47.9 → 11.7

Time: 29.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}\;k \le -1.3484728286962528 \cdot 10^{169}:\\ \;\;\;\;2 \cdot \left({\left(\frac{\sqrt{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)\\ \mathbf{elif}\;k \le -1.5790227602943146 \cdot 10^{-138}:\\ \;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\sin k \cdot \frac{\sin k}{\ell}}\\ \mathbf{elif}\;k \le 1.2653785684425997 \cdot 10^{-118}:\\ \;\;\;\;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}{\sin k} \cdot \ell}{\frac{\sin k}{\ell}}\right)\\ \mathbf{elif}\;k \le 2.25324371902540199 \cdot 10^{110}:\\ \;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\sin k \cdot \frac{\sin k}{\ell}}\\ \mathbf{else}:\\ \;\;\;\;2 \cdot \left({\left(\frac{\sqrt{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)\\ \end{array}\]

C

double f(double t, double l, double k) {
        double r91185 = 2.0;
        double r91186 = t;
        double r91187 = 3.0;
        double r91188 = pow(r91186, r91187);
        double r91189 = l;
        double r91190 = r91189 * r91189;
        double r91191 = r91188 / r91190;
        double r91192 = k;
        double r91193 = sin(r91192);
        double r91194 = r91191 * r91193;
        double r91195 = tan(r91192);
        double r91196 = r91194 * r91195;
        double r91197 = 1.0;
        double r91198 = r91192 / r91186;
        double r91199 = pow(r91198, r91185);
        double r91200 = r91197 + r91199;
        double r91201 = r91200 - r91197;
        double r91202 = r91196 * r91201;
        double r91203 = r91185 / r91202;
        return r91203;
}

↓

double f(double t, double l, double k) {
        double r91204 = k;
        double r91205 = -1.3484728286962528e+169;
        bool r91206 = r91204 <= r91205;
        double r91207 = 2.0;
        double r91208 = 1.0;
        double r91209 = sqrt(r91208);
        double r91210 = 2.0;
        double r91211 = r91207 / r91210;
        double r91212 = pow(r91204, r91211);
        double r91213 = r91209 / r91212;
        double r91214 = 1.0;
        double r91215 = pow(r91213, r91214);
        double r91216 = t;
        double r91217 = pow(r91216, r91214);
        double r91218 = r91212 * r91217;
        double r91219 = r91208 / r91218;
        double r91220 = pow(r91219, r91214);
        double r91221 = cos(r91204);
        double r91222 = l;
        double r91223 = pow(r91222, r91210);
        double r91224 = r91221 * r91223;
        double r91225 = sin(r91204);
        double r91226 = pow(r91225, r91210);
        double r91227 = r91224 / r91226;
        double r91228 = r91220 * r91227;
        double r91229 = r91215 * r91228;
        double r91230 = r91207 * r91229;
        double r91231 = -1.5790227602943146e-138;
        bool r91232 = r91204 <= r91231;
        double r91233 = pow(r91204, r91207);
        double r91234 = r91233 * r91217;
        double r91235 = r91208 / r91234;
        double r91236 = pow(r91235, r91214);
        double r91237 = r91236 * r91222;
        double r91238 = r91237 * r91221;
        double r91239 = r91225 / r91222;
        double r91240 = r91225 * r91239;
        double r91241 = r91238 / r91240;
        double r91242 = r91207 * r91241;
        double r91243 = 1.2653785684425997e-118;
        bool r91244 = r91204 <= r91243;
        double r91245 = r91212 * r91218;
        double r91246 = r91208 / r91245;
        double r91247 = pow(r91246, r91214);
        double r91248 = r91221 / r91225;
        double r91249 = r91248 * r91222;
        double r91250 = r91249 / r91239;
        double r91251 = r91247 * r91250;
        double r91252 = r91207 * r91251;
        double r91253 = 2.253243719025402e+110;
        bool r91254 = r91204 <= r91253;
        double r91255 = r91254 ? r91242 : r91230;
        double r91256 = r91244 ? r91252 : r91255;
        double r91257 = r91232 ? r91242 : r91256;
        double r91258 = r91206 ? r91230 : r91257;
        return r91258;
}

Error

Bits error versus t

Bits error versus l

Bits error versus k

Derivation

1. Split input into 3 regimes

2. if k < -1.3484728286962528e+169 or 2.253243719025402e+110 < k

   1. Initial program 40.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.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 23.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 23.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* 18.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 18.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 18.6

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

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

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

       \[\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)\]

   if -1.3484728286962528e+169 < k < -1.5790227602943146e-138 or 1.2653785684425997e-118 < k < 2.253243719025402e+110

   1. Initial program 53.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 42.0

      \[\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.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 17.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* 16.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-sqr-sqrt 40.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{\cos k \cdot {\ell}^{2}}{{\color{blue}{\left(\sqrt{\sin k} \cdot \sqrt{\sin k}\right)}}^{2}}\right)\]

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

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

   14. Applied frac-times 13.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 \color{blue}{\frac{\cos k \cdot \ell}{\sin k \cdot \frac{\sin k}{\ell}}}\right)\]

   15. Applied associate-*r/ 7.1

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

   16. Simplified 7.7

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

   if -1.5790227602943146e-138 < k < 1.2653785684425997e-118

   1. Initial program 64.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 64.0

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

   9. Applied unpow-prod-down 57.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 \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 55.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 55.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 29.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}{\frac{\sin k}{\ell}}}\right)\right)\]
   13. Using strategy rm

   14. Applied associate-*r/ 16.5

       \[\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}{\sin k} \cdot \ell}{\frac{\sin k}{\ell}}}\right)\]
3. Recombined 3 regimes into one program.

4. Final simplification 11.7

   \[\leadsto \begin{array}{l} \mathbf{if}\;k \le -1.3484728286962528 \cdot 10^{169}:\\ \;\;\;\;2 \cdot \left({\left(\frac{\sqrt{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)\\ \mathbf{elif}\;k \le -1.5790227602943146 \cdot 10^{-138}:\\ \;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\sin k \cdot \frac{\sin k}{\ell}}\\ \mathbf{elif}\;k \le 1.2653785684425997 \cdot 10^{-118}:\\ \;\;\;\;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}{\sin k} \cdot \ell}{\frac{\sin k}{\ell}}\right)\\ \mathbf{elif}\;k \le 2.25324371902540199 \cdot 10^{110}:\\ \;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\sin k \cdot \frac{\sin k}{\ell}}\\ \mathbf{else}:\\ \;\;\;\;2 \cdot \left({\left(\frac{\sqrt{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)\\ \end{array}\]

Reproduce

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