Average Error: 48.2 → 16.9
Time: 1.1m
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}\;\ell \le -9.576489350063892240805921733875452410192 \cdot 10^{153}:\\ \;\;\;\;\frac{\frac{2}{\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}}{{\left(\frac{k}{t}\right)}^{2}}\\ \mathbf{else}:\\ \;\;\;\;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)\\ \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}\;\ell \le -9.576489350063892240805921733875452410192 \cdot 10^{153}:\\
\;\;\;\;\frac{\frac{2}{\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}}{{\left(\frac{k}{t}\right)}^{2}}\\

\mathbf{else}:\\
\;\;\;\;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)\\

\end{array}
double f(double t, double l, double k) {
        double r94428 = 2.0;
        double r94429 = t;
        double r94430 = 3.0;
        double r94431 = pow(r94429, r94430);
        double r94432 = l;
        double r94433 = r94432 * r94432;
        double r94434 = r94431 / r94433;
        double r94435 = k;
        double r94436 = sin(r94435);
        double r94437 = r94434 * r94436;
        double r94438 = tan(r94435);
        double r94439 = r94437 * r94438;
        double r94440 = 1.0;
        double r94441 = r94435 / r94429;
        double r94442 = pow(r94441, r94428);
        double r94443 = r94440 + r94442;
        double r94444 = r94443 - r94440;
        double r94445 = r94439 * r94444;
        double r94446 = r94428 / r94445;
        return r94446;
}
↓
double f(double t, double l, double k) {
        double r94447 = l;
        double r94448 = -9.576489350063892e+153;
        bool r94449 = r94447 <= r94448;
        double r94450 = 2.0;
        double r94451 = t;
        double r94452 = cbrt(r94451);
        double r94453 = r94452 * r94452;
        double r94454 = 3.0;
        double r94455 = pow(r94453, r94454);
        double r94456 = r94455 / r94447;
        double r94457 = pow(r94452, r94454);
        double r94458 = r94457 / r94447;
        double r94459 = k;
        double r94460 = sin(r94459);
        double r94461 = r94458 * r94460;
        double r94462 = r94456 * r94461;
        double r94463 = tan(r94459);
        double r94464 = r94462 * r94463;
        double r94465 = r94450 / r94464;
        double r94466 = r94459 / r94451;
        double r94467 = pow(r94466, r94450);
        double r94468 = r94465 / r94467;
        double r94469 = 1.0;
        double r94470 = 2.0;
        double r94471 = r94450 / r94470;
        double r94472 = pow(r94459, r94471);
        double r94473 = r94469 / r94472;
        double r94474 = 1.0;
        double r94475 = pow(r94473, r94474);
        double r94476 = pow(r94451, r94474);
        double r94477 = r94472 * r94476;
        double r94478 = r94469 / r94477;
        double r94479 = pow(r94478, r94474);
        double r94480 = cos(r94459);
        double r94481 = r94480 / r94460;
        double r94482 = pow(r94447, r94470);
        double r94483 = r94482 / r94460;
        double r94484 = r94481 * r94483;
        double r94485 = r94479 * r94484;
        double r94486 = r94475 * r94485;
        double r94487 = r94450 * r94486;
        double r94488 = r94449 ? r94468 : r94487;
        return r94488;
}

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 l < -9.576489350063892e+153

    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. Simplified64.0

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

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

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

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

    7. Applied associate-*l*50.4

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

    if -9.576489350063892e+153 < l

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

      \[\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 18.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-pow18.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*16.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 *-un-lft-identity16.5

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

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

    11. Applied associate-*l*14.4

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

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

    14. Applied times-frac14.1

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

  4. Final simplification16.9

    \[\leadsto \begin{array}{l} \mathbf{if}\;\ell \le -9.576489350063892240805921733875452410192 \cdot 10^{153}:\\ \;\;\;\;\frac{\frac{2}{\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}}{{\left(\frac{k}{t}\right)}^{2}}\\ \mathbf{else}:\\ \;\;\;\;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)\\ \end{array}\]

Reproduce

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