Average Error: 32.7 → 17.5
Time: 24.4s
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 -5.481507593861878460145816755067202944483 \cdot 10^{-263} \lor \neg \left(t \le 6.565726687888966078976606764192118335339 \cdot 10^{-111}\right):\\ \;\;\;\;\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\sqrt[3]{\ell}}\right) \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}\\ \mathbf{else}:\\ \;\;\;\;\frac{2}{2 \cdot \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot {\ell}^{2}} - {\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{t \cdot \left({k}^{2} \cdot {\left(\sin k\right)}^{2}\right)}{\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 -5.481507593861878460145816755067202944483 \cdot 10^{-263} \lor \neg \left(t \le 6.565726687888966078976606764192118335339 \cdot 10^{-111}\right):\\
\;\;\;\;\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{{t}^{\left(\frac{1}{3} \cdot 3\right)}}{\sqrt[3]{\ell}}\right) \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}\\

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

\end{array}
double f(double t, double l, double k) {
        double r117325 = 2.0;
        double r117326 = t;
        double r117327 = 3.0;
        double r117328 = pow(r117326, r117327);
        double r117329 = l;
        double r117330 = r117329 * r117329;
        double r117331 = r117328 / r117330;
        double r117332 = k;
        double r117333 = sin(r117332);
        double r117334 = r117331 * r117333;
        double r117335 = tan(r117332);
        double r117336 = r117334 * r117335;
        double r117337 = 1.0;
        double r117338 = r117332 / r117326;
        double r117339 = pow(r117338, r117325);
        double r117340 = r117337 + r117339;
        double r117341 = r117340 + r117337;
        double r117342 = r117336 * r117341;
        double r117343 = r117325 / r117342;
        return r117343;
}

double f(double t, double l, double k) {
        double r117344 = t;
        double r117345 = -5.4815075938618785e-263;
        bool r117346 = r117344 <= r117345;
        double r117347 = 6.565726687888966e-111;
        bool r117348 = r117344 <= r117347;
        double r117349 = !r117348;
        bool r117350 = r117346 || r117349;
        double r117351 = 2.0;
        double r117352 = cbrt(r117344);
        double r117353 = 3.0;
        double r117354 = pow(r117352, r117353);
        double r117355 = l;
        double r117356 = cbrt(r117355);
        double r117357 = r117356 * r117356;
        double r117358 = r117354 / r117357;
        double r117359 = 0.3333333333333333;
        double r117360 = r117359 * r117353;
        double r117361 = pow(r117344, r117360);
        double r117362 = r117361 / r117356;
        double r117363 = r117358 * r117362;
        double r117364 = r117354 / r117355;
        double r117365 = k;
        double r117366 = sin(r117365);
        double r117367 = r117364 * r117366;
        double r117368 = r117363 * r117367;
        double r117369 = tan(r117365);
        double r117370 = 1.0;
        double r117371 = r117365 / r117344;
        double r117372 = pow(r117371, r117351);
        double r117373 = r117370 + r117372;
        double r117374 = r117373 + r117370;
        double r117375 = r117369 * r117374;
        double r117376 = r117368 * r117375;
        double r117377 = r117351 / r117376;
        double r117378 = 3.0;
        double r117379 = pow(r117344, r117378);
        double r117380 = 2.0;
        double r117381 = pow(r117366, r117380);
        double r117382 = r117379 * r117381;
        double r117383 = cos(r117365);
        double r117384 = pow(r117355, r117380);
        double r117385 = r117383 * r117384;
        double r117386 = r117382 / r117385;
        double r117387 = r117351 * r117386;
        double r117388 = 1.0;
        double r117389 = -1.0;
        double r117390 = pow(r117389, r117353);
        double r117391 = r117388 / r117390;
        double r117392 = pow(r117391, r117370);
        double r117393 = pow(r117365, r117380);
        double r117394 = r117393 * r117381;
        double r117395 = r117344 * r117394;
        double r117396 = r117395 / r117385;
        double r117397 = r117392 * r117396;
        double r117398 = r117387 - r117397;
        double r117399 = r117351 / r117398;
        double r117400 = r117350 ? r117377 : r117399;
        return r117400;
}

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 < -5.4815075938618785e-263 or 6.565726687888966e-111 < t

    1. Initial program 27.6

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

      \[\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.3

      \[\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.3

      \[\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 unpow-prod-down19.3

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

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

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

      \[\leadsto \frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}} \cdot \frac{{\color{blue}{\left({t}^{\frac{1}{3}}\right)}}^{3}}{\sqrt[3]{\ell}}\right) \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}\]
    15. Applied pow-pow13.9

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

    if -5.4815075938618785e-263 < t < 6.565726687888966e-111

    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. Taylor expanded around -inf 39.7

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

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

Reproduce

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