Toniolo and Linder, Equation (10-)

Average Error: 48.2 → 13.0

Time: 25.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}\;k \le -4.65561929088396835 \cdot 10^{155}:\\ \;\;\;\;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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \frac{\ell}{\frac{\left|\sin k\right|}{\ell}}\right)\right)\right)\\ \mathbf{elif}\;k \le -1.9372850856630931 \cdot 10^{-128}:\\ \;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\left|\sin k\right| \cdot \frac{\left|\sin k\right|}{\ell}}\\ \mathbf{elif}\;k \le 3.08249264848917646 \cdot 10^{-154}:\\ \;\;\;\;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}{\left|\sin k\right|} \cdot \ell\right) \cdot \frac{1}{\frac{\left|\sin k\right|}{\ell}}\right)\right)\\ \mathbf{else}:\\ \;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\left|\sin k\right| \cdot \frac{\left|\sin k\right|}{\ell}}\\ \end{array}\]

C

double f(double t, double l, double k) {
        double r85098 = 2.0;
        double r85099 = t;
        double r85100 = 3.0;
        double r85101 = pow(r85099, r85100);
        double r85102 = l;
        double r85103 = r85102 * r85102;
        double r85104 = r85101 / r85103;
        double r85105 = k;
        double r85106 = sin(r85105);
        double r85107 = r85104 * r85106;
        double r85108 = tan(r85105);
        double r85109 = r85107 * r85108;
        double r85110 = 1.0;
        double r85111 = r85105 / r85099;
        double r85112 = pow(r85111, r85098);
        double r85113 = r85110 + r85112;
        double r85114 = r85113 - r85110;
        double r85115 = r85109 * r85114;
        double r85116 = r85098 / r85115;
        return r85116;
}

↓

double f(double t, double l, double k) {
        double r85117 = k;
        double r85118 = -4.655619290883968e+155;
        bool r85119 = r85117 <= r85118;
        double r85120 = 2.0;
        double r85121 = 1.0;
        double r85122 = sqrt(r85121);
        double r85123 = 2.0;
        double r85124 = r85120 / r85123;
        double r85125 = pow(r85117, r85124);
        double r85126 = r85122 / r85125;
        double r85127 = 1.0;
        double r85128 = pow(r85126, r85127);
        double r85129 = t;
        double r85130 = pow(r85129, r85127);
        double r85131 = r85125 * r85130;
        double r85132 = r85121 / r85131;
        double r85133 = pow(r85132, r85127);
        double r85134 = cos(r85117);
        double r85135 = sin(r85117);
        double r85136 = fabs(r85135);
        double r85137 = r85134 / r85136;
        double r85138 = l;
        double r85139 = r85136 / r85138;
        double r85140 = r85138 / r85139;
        double r85141 = r85137 * r85140;
        double r85142 = r85133 * r85141;
        double r85143 = r85128 * r85142;
        double r85144 = r85120 * r85143;
        double r85145 = -1.937285085663093e-128;
        bool r85146 = r85117 <= r85145;
        double r85147 = pow(r85117, r85120);
        double r85148 = r85147 * r85130;
        double r85149 = r85121 / r85148;
        double r85150 = pow(r85149, r85127);
        double r85151 = r85150 * r85138;
        double r85152 = r85151 * r85134;
        double r85153 = r85136 * r85139;
        double r85154 = r85152 / r85153;
        double r85155 = r85120 * r85154;
        double r85156 = 3.0824926484891765e-154;
        bool r85157 = r85117 <= r85156;
        double r85158 = r85125 * r85131;
        double r85159 = r85121 / r85158;
        double r85160 = pow(r85159, r85127);
        double r85161 = r85137 * r85138;
        double r85162 = r85121 / r85139;
        double r85163 = r85161 * r85162;
        double r85164 = r85160 * r85163;
        double r85165 = r85120 * r85164;
        double r85166 = r85157 ? r85165 : r85155;
        double r85167 = r85146 ? r85155 : r85166;
        double r85168 = r85119 ? r85144 : r85167;
        return r85168;
}

Error

Bits error versus t

Bits error versus l

Bits error versus k

Derivation

1. Split input into 3 regimes

2. if k < -4.655619290883968e+155

   1. Initial program 39.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 34.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 24.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 24.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 add-sqr-sqrt 20.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}{\sqrt{{\left(\sin k\right)}^{2}} \cdot \sqrt{{\left(\sin k\right)}^{2}}}}\right)\]

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

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

   11. Simplified 20.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}{\left|\sin k\right|} \cdot \color{blue}{\frac{\ell}{\frac{\left|\sin k\right|}{\ell}}}\right)\right)\]
   12. Using strategy rm

   13. Applied add-sqr-sqrt 20.1

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

   14. Applied times-frac 19.8

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

   15. Applied unpow-prod-down 19.8

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

   16. Applied associate-*l* 16.4

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

   17. Simplified 16.4

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

   if -4.655619290883968e+155 < k < -1.937285085663093e-128 or 3.0824926484891765e-154 < k

   1. Initial program 49.7

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

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

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

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

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

   11. Simplified 16.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}{\left|\sin k\right|} \cdot \color{blue}{\frac{\ell}{\frac{\left|\sin k\right|}{\ell}}}\right)\right)\]
   12. Using strategy rm

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

   14. Applied associate-*r/ 9.8

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

   15. Simplified 12.0

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

   if -1.937285085663093e-128 < k < 3.0824926484891765e-154

   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 56.5

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

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

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

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

   11. Simplified 28.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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \color{blue}{\frac{\ell}{\frac{\left|\sin k\right|}{\ell}}}\right)\right)\]
   12. Using strategy rm

   13. Applied div-inv 28.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 \left(\frac{\cos k}{\left|\sin k\right|} \cdot \color{blue}{\left(\ell \cdot \frac{1}{\frac{\left|\sin k\right|}{\ell}}\right)}\right)\right)\]

   14. Applied associate-*r* 15.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(\left(\frac{\cos k}{\left|\sin k\right|} \cdot \ell\right) \cdot \frac{1}{\frac{\left|\sin k\right|}{\ell}}\right)}\right)\]
3. Recombined 3 regimes into one program.

4. Final simplification 13.0

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

Reproduce

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