Average Error: 31.9 → 10.1
Time: 21.3s
Precision: 64
\[\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}\;t \le -3.56310395914880465 \cdot 10^{-34} \lor \neg \left(t \le 6.3604377920357265 \cdot 10^{-52}\right):\\ \;\;\;\;\frac{2}{\frac{\left(\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right)\right) \cdot \sqrt[3]{\sin k}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}\\ \mathbf{else}:\\ \;\;\;\;\frac{2}{\frac{\frac{{k}^{2} \cdot {\left(\sin k\right)}^{2}}{\ell} + 2 \cdot \frac{{t}^{2} \cdot {\left(\sin k\right)}^{2}}{\ell}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}\\ \end{array}\]
\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}\;t \le -3.56310395914880465 \cdot 10^{-34} \lor \neg \left(t \le 6.3604377920357265 \cdot 10^{-52}\right):\\
\;\;\;\;\frac{2}{\frac{\left(\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right)\right) \cdot \sqrt[3]{\sin k}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}\\

\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\frac{{k}^{2} \cdot {\left(\sin k\right)}^{2}}{\ell} + 2 \cdot \frac{{t}^{2} \cdot {\left(\sin k\right)}^{2}}{\ell}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}\\

\end{array}
double f(double t, double l, double k) {
        double r107228 = 2.0;
        double r107229 = t;
        double r107230 = 3.0;
        double r107231 = pow(r107229, r107230);
        double r107232 = l;
        double r107233 = r107232 * r107232;
        double r107234 = r107231 / r107233;
        double r107235 = k;
        double r107236 = sin(r107235);
        double r107237 = r107234 * r107236;
        double r107238 = tan(r107235);
        double r107239 = r107237 * r107238;
        double r107240 = 1.0;
        double r107241 = r107235 / r107229;
        double r107242 = pow(r107241, r107228);
        double r107243 = r107240 + r107242;
        double r107244 = r107243 + r107240;
        double r107245 = r107239 * r107244;
        double r107246 = r107228 / r107245;
        return r107246;
}

double f(double t, double l, double k) {
        double r107247 = t;
        double r107248 = -3.5631039591488047e-34;
        bool r107249 = r107247 <= r107248;
        double r107250 = 6.3604377920357265e-52;
        bool r107251 = r107247 <= r107250;
        double r107252 = !r107251;
        bool r107253 = r107249 || r107252;
        double r107254 = 2.0;
        double r107255 = cbrt(r107247);
        double r107256 = 3.0;
        double r107257 = pow(r107255, r107256);
        double r107258 = l;
        double r107259 = r107257 / r107258;
        double r107260 = k;
        double r107261 = sin(r107260);
        double r107262 = r107259 * r107261;
        double r107263 = r107257 * r107262;
        double r107264 = cbrt(r107261);
        double r107265 = r107264 * r107264;
        double r107266 = r107263 * r107265;
        double r107267 = r107266 * r107264;
        double r107268 = 1.0;
        double r107269 = r107260 / r107247;
        double r107270 = pow(r107269, r107254);
        double r107271 = r107268 + r107270;
        double r107272 = r107271 + r107268;
        double r107273 = r107267 * r107272;
        double r107274 = r107258 / r107257;
        double r107275 = cos(r107260);
        double r107276 = r107274 * r107275;
        double r107277 = r107273 / r107276;
        double r107278 = r107254 / r107277;
        double r107279 = 2.0;
        double r107280 = pow(r107260, r107279);
        double r107281 = pow(r107261, r107279);
        double r107282 = r107280 * r107281;
        double r107283 = r107282 / r107258;
        double r107284 = pow(r107247, r107279);
        double r107285 = r107284 * r107281;
        double r107286 = r107285 / r107258;
        double r107287 = r107254 * r107286;
        double r107288 = r107283 + r107287;
        double r107289 = r107288 / r107276;
        double r107290 = r107254 / r107289;
        double r107291 = r107253 ? r107278 : r107290;
        return r107291;
}

Error

Bits error versus t

Bits error versus l

Bits error versus k

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes
  2. if t < -3.5631039591488047e-34 or 6.3604377920357265e-52 < t

    1. Initial program 21.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. Using strategy rm
    3. Applied add-cube-cbrt22.1

      \[\leadsto \frac{2}{\left(\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\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    4. Applied unpow-prod-down22.1

      \[\leadsto \frac{2}{\left(\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\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    5. Applied times-frac15.8

      \[\leadsto \frac{2}{\left(\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\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    6. Applied associate-*l*13.7

      \[\leadsto \frac{2}{\left(\color{blue}{\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{3}}{\ell} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right)} \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    7. Using strategy rm
    8. Applied unpow-prod-down13.7

      \[\leadsto \frac{2}{\left(\left(\frac{\color{blue}{{\left(\sqrt[3]{t}\right)}^{3} \cdot {\left(\sqrt[3]{t}\right)}^{3}}}{\ell} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    9. Applied associate-/l*8.6

      \[\leadsto \frac{2}{\left(\left(\color{blue}{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    10. Using strategy rm
    11. Applied tan-quot8.6

      \[\leadsto \frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \color{blue}{\frac{\sin k}{\cos k}}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    12. Applied associate-*l/7.4

      \[\leadsto \frac{2}{\left(\color{blue}{\frac{{\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}} \cdot \frac{\sin k}{\cos k}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    13. Applied frac-times5.1

      \[\leadsto \frac{2}{\color{blue}{\frac{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \sin k}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}} \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    14. Applied associate-*l/4.7

      \[\leadsto \frac{2}{\color{blue}{\frac{\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \sin k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}}\]
    15. Using strategy rm
    16. Applied add-cube-cbrt4.7

      \[\leadsto \frac{2}{\frac{\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \color{blue}{\left(\left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right) \cdot \sqrt[3]{\sin k}\right)}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}\]
    17. Applied associate-*r*4.7

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

    if -3.5631039591488047e-34 < t < 6.3604377920357265e-52

    1. Initial program 54.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. Using strategy rm
    3. Applied add-cube-cbrt54.8

      \[\leadsto \frac{2}{\left(\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\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    4. Applied unpow-prod-down54.8

      \[\leadsto \frac{2}{\left(\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\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    5. Applied times-frac46.5

      \[\leadsto \frac{2}{\left(\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\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    6. Applied associate-*l*45.6

      \[\leadsto \frac{2}{\left(\color{blue}{\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{3}}{\ell} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right)} \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    7. Using strategy rm
    8. Applied unpow-prod-down45.6

      \[\leadsto \frac{2}{\left(\left(\frac{\color{blue}{{\left(\sqrt[3]{t}\right)}^{3} \cdot {\left(\sqrt[3]{t}\right)}^{3}}}{\ell} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    9. Applied associate-/l*38.7

      \[\leadsto \frac{2}{\left(\left(\color{blue}{\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    10. Using strategy rm
    11. Applied tan-quot38.7

      \[\leadsto \frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \color{blue}{\frac{\sin k}{\cos k}}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    12. Applied associate-*l/38.7

      \[\leadsto \frac{2}{\left(\color{blue}{\frac{{\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}} \cdot \frac{\sin k}{\cos k}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    13. Applied frac-times40.2

      \[\leadsto \frac{2}{\color{blue}{\frac{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \sin k}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}} \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\]
    14. Applied associate-*l/36.4

      \[\leadsto \frac{2}{\color{blue}{\frac{\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \sin k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}}\]
    15. Taylor expanded around inf 22.5

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

    \[\leadsto \begin{array}{l} \mathbf{if}\;t \le -3.56310395914880465 \cdot 10^{-34} \lor \neg \left(t \le 6.3604377920357265 \cdot 10^{-52}\right):\\ \;\;\;\;\frac{2}{\frac{\left(\left(\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \left(\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}\right)\right) \cdot \sqrt[3]{\sin k}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}\\ \mathbf{else}:\\ \;\;\;\;\frac{2}{\frac{\frac{{k}^{2} \cdot {\left(\sin k\right)}^{2}}{\ell} + 2 \cdot \frac{{t}^{2} \cdot {\left(\sin k\right)}^{2}}{\ell}}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}} \cdot \cos k}}\\ \end{array}\]

Reproduce

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