Average Error: 32.1 → 13.6
Time: 25.6s
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.610783428824226 \cdot 10^{-248} \lor \neg \left(t \le 2.60350273525137057 \cdot 10^{-209}\right):\\ \;\;\;\;\frac{2}{\frac{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \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} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}} \cdot \cos k}}\\ \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 -5.610783428824226 \cdot 10^{-248} \lor \neg \left(t \le 2.60350273525137057 \cdot 10^{-209}\right):\\
\;\;\;\;\frac{2}{\frac{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \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} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}} \cdot \cos k}}\\

\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 r133228 = 2.0;
        double r133229 = t;
        double r133230 = 3.0;
        double r133231 = pow(r133229, r133230);
        double r133232 = l;
        double r133233 = r133232 * r133232;
        double r133234 = r133231 / r133233;
        double r133235 = k;
        double r133236 = sin(r133235);
        double r133237 = r133234 * r133236;
        double r133238 = tan(r133235);
        double r133239 = r133237 * r133238;
        double r133240 = 1.0;
        double r133241 = r133235 / r133229;
        double r133242 = pow(r133241, r133228);
        double r133243 = r133240 + r133242;
        double r133244 = r133243 + r133240;
        double r133245 = r133239 * r133244;
        double r133246 = r133228 / r133245;
        return r133246;
}
↓
double f(double t, double l, double k) {
        double r133247 = t;
        double r133248 = -5.610783428824226e-248;
        bool r133249 = r133247 <= r133248;
        double r133250 = 2.6035027352513706e-209;
        bool r133251 = r133247 <= r133250;
        double r133252 = !r133251;
        bool r133253 = r133249 || r133252;
        double r133254 = 2.0;
        double r133255 = cbrt(r133247);
        double r133256 = r133255 * r133255;
        double r133257 = 3.0;
        double r133258 = 2.0;
        double r133259 = r133257 / r133258;
        double r133260 = pow(r133256, r133259);
        double r133261 = pow(r133255, r133257);
        double r133262 = l;
        double r133263 = r133261 / r133262;
        double r133264 = k;
        double r133265 = sin(r133264);
        double r133266 = r133263 * r133265;
        double r133267 = r133260 * r133266;
        double r133268 = r133267 * r133265;
        double r133269 = 1.0;
        double r133270 = r133264 / r133247;
        double r133271 = pow(r133270, r133254);
        double r133272 = r133269 + r133271;
        double r133273 = r133272 + r133269;
        double r133274 = r133268 * r133273;
        double r133275 = r133262 / r133260;
        double r133276 = cos(r133264);
        double r133277 = r133275 * r133276;
        double r133278 = r133274 / r133277;
        double r133279 = r133254 / r133278;
        double r133280 = pow(r133264, r133258);
        double r133281 = pow(r133265, r133258);
        double r133282 = r133247 * r133281;
        double r133283 = r133280 * r133282;
        double r133284 = pow(r133262, r133258);
        double r133285 = r133276 * r133284;
        double r133286 = r133283 / r133285;
        double r133287 = 3.0;
        double r133288 = pow(r133247, r133287);
        double r133289 = r133288 * r133281;
        double r133290 = r133289 / r133285;
        double r133291 = r133254 * r133290;
        double r133292 = r133286 + r133291;
        double r133293 = r133254 / r133292;
        double r133294 = r133253 ? r133279 : r133293;
        return r133294;
}

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.610783428824226e-248 or 2.6035027352513706e-209 < t

    1. Initial program 29.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-cbrt29.6

      \[\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-down29.6

      \[\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-frac22.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*20.5

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

      \[\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)}}}{\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*14.8

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

    11. Applied tan-quot14.8

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

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

    13. Applied frac-times12.5

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

    14. Applied associate-*l/11.1

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

    if -5.610783428824226e-248 < t < 2.6035027352513706e-209

    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 42.9

      \[\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 simplification13.6

    \[\leadsto \begin{array}{l} \mathbf{if}\;t \le -5.610783428824226 \cdot 10^{-248} \lor \neg \left(t \le 2.60350273525137057 \cdot 10^{-209}\right):\\ \;\;\;\;\frac{2}{\frac{\left(\left({\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)} \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} \cdot \sqrt[3]{t}\right)}^{\left(\frac{3}{2}\right)}} \cdot \cos k}}\\ \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 2020045 
(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))))