Average Error: 32.7 → 15.9
Time: 23.8s
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.22413082050212105 \cdot 10^{-15} \lor \neg \left(t \le 4.896278986279809 \cdot 10^{-272}\right):\\ \;\;\;\;\frac{2}{\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right)\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\ \mathbf{else}:\\ \;\;\;\;\frac{2}{\frac{{k}^{2} \cdot \left(t \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot {\ell}^{2}} + 2 \cdot \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot {\ell}^{2}}}\\ \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.22413082050212105 \cdot 10^{-15} \lor \neg \left(t \le 4.896278986279809 \cdot 10^{-272}\right):\\
\;\;\;\;\frac{2}{\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell}} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k\right)\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}\\

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

\end{array}
double f(double t, double l, double k) {
        double r141484 = 2.0;
        double r141485 = t;
        double r141486 = 3.0;
        double r141487 = pow(r141485, r141486);
        double r141488 = l;
        double r141489 = r141488 * r141488;
        double r141490 = r141487 / r141489;
        double r141491 = k;
        double r141492 = sin(r141491);
        double r141493 = r141490 * r141492;
        double r141494 = tan(r141491);
        double r141495 = r141493 * r141494;
        double r141496 = 1.0;
        double r141497 = r141491 / r141485;
        double r141498 = pow(r141497, r141484);
        double r141499 = r141496 + r141498;
        double r141500 = r141499 + r141496;
        double r141501 = r141495 * r141500;
        double r141502 = r141484 / r141501;
        return r141502;
}
↓
double f(double t, double l, double k) {
        double r141503 = t;
        double r141504 = -3.224130820502121e-15;
        bool r141505 = r141503 <= r141504;
        double r141506 = 4.896278986279809e-272;
        bool r141507 = r141503 <= r141506;
        double r141508 = !r141507;
        bool r141509 = r141505 || r141508;
        double r141510 = 2.0;
        double r141511 = cbrt(r141503);
        double r141512 = r141511 * r141511;
        double r141513 = 3.0;
        double r141514 = 2.0;
        double r141515 = r141513 / r141514;
        double r141516 = pow(r141512, r141515);
        double r141517 = l;
        double r141518 = cbrt(r141517);
        double r141519 = r141518 * r141518;
        double r141520 = r141516 / r141519;
        double r141521 = r141516 / r141518;
        double r141522 = pow(r141511, r141513);
        double r141523 = r141522 / r141517;
        double r141524 = k;
        double r141525 = sin(r141524);
        double r141526 = r141523 * r141525;
        double r141527 = r141521 * r141526;
        double r141528 = tan(r141524);
        double r141529 = r141527 * r141528;
        double r141530 = r141520 * r141529;
        double r141531 = 1.0;
        double r141532 = r141524 / r141503;
        double r141533 = pow(r141532, r141510);
        double r141534 = r141531 + r141533;
        double r141535 = r141534 + r141531;
        double r141536 = r141530 * r141535;
        double r141537 = r141510 / r141536;
        double r141538 = pow(r141524, r141514);
        double r141539 = pow(r141525, r141514);
        double r141540 = r141503 * r141539;
        double r141541 = r141538 * r141540;
        double r141542 = cos(r141524);
        double r141543 = pow(r141517, r141514);
        double r141544 = r141542 * r141543;
        double r141545 = r141541 / r141544;
        double r141546 = 3.0;
        double r141547 = pow(r141503, r141546);
        double r141548 = r141547 * r141539;
        double r141549 = r141548 / r141544;
        double r141550 = r141510 * r141549;
        double r141551 = r141545 + r141550;
        double r141552 = r141510 / r141551;
        double r141553 = r141509 ? r141537 : r141552;
        return r141553;
}

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.224130820502121e-15 or 4.896278986279809e-272 < t

    1. Initial program 27.4

      \[\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-cbrt27.5

      \[\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-down27.5

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

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

      \[\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 add-cube-cbrt19.1

      \[\leadsto \frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{3}}{\color{blue}{\left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right) \cdot \sqrt[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 sqr-pow19.1

      \[\leadsto \frac{2}{\left(\left(\frac{\color{blue}{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \cdot {\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}}{\left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right) \cdot \sqrt[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)}\]

    10. Applied times-frac13.3

      \[\leadsto \frac{2}{\left(\left(\color{blue}{\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}}{\sqrt[3]{\ell}}\right)} \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)}\]
    11. Using strategy rm

    12. Applied associate-*l*12.3

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

    14. Applied associate-*l*10.6

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

    if -3.224130820502121e-15 < t < 4.896278986279809e-272

    1. Initial program 53.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. Taylor expanded around inf 36.5

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

  4. Final simplification15.9

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

Reproduce

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