Average Error: 19.3 → 9.8
Time: 29.0s
Precision: 64
\[c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\]
↓
\[\begin{array}{l} \mathbf{if}\;\ell \cdot V = -\infty:\\ \;\;\;\;\left(\sqrt{\frac{A}{\ell}} \cdot \sqrt{\frac{1}{V}}\right) \cdot c0\\ \mathbf{elif}\;\ell \cdot V \le -5.601891409264343123426103753339848933114 \cdot 10^{-287}:\\ \;\;\;\;\sqrt{\sqrt[3]{\frac{1}{\ell \cdot V}} \cdot \sqrt[3]{A}} \cdot \left(c0 \cdot \frac{\left|\sqrt[3]{A}\right|}{\left|\sqrt[3]{\ell \cdot V}\right|}\right)\\ \mathbf{elif}\;\ell \cdot V \le 2.442774188137666401114116722762846155888 \cdot 10^{-317}:\\ \;\;\;\;\left(\sqrt{\frac{A}{\ell}} \cdot \sqrt{\frac{1}{V}}\right) \cdot c0\\ \mathbf{else}:\\ \;\;\;\;c0 \cdot \frac{\sqrt{A}}{\sqrt{\ell \cdot V}}\\ \end{array}\]
c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}
↓
\begin{array}{l}
\mathbf{if}\;\ell \cdot V = -\infty:\\
\;\;\;\;\left(\sqrt{\frac{A}{\ell}} \cdot \sqrt{\frac{1}{V}}\right) \cdot c0\\

\mathbf{elif}\;\ell \cdot V \le -5.601891409264343123426103753339848933114 \cdot 10^{-287}:\\
\;\;\;\;\sqrt{\sqrt[3]{\frac{1}{\ell \cdot V}} \cdot \sqrt[3]{A}} \cdot \left(c0 \cdot \frac{\left|\sqrt[3]{A}\right|}{\left|\sqrt[3]{\ell \cdot V}\right|}\right)\\

\mathbf{elif}\;\ell \cdot V \le 2.442774188137666401114116722762846155888 \cdot 10^{-317}:\\
\;\;\;\;\left(\sqrt{\frac{A}{\ell}} \cdot \sqrt{\frac{1}{V}}\right) \cdot c0\\

\mathbf{else}:\\
\;\;\;\;c0 \cdot \frac{\sqrt{A}}{\sqrt{\ell \cdot V}}\\

\end{array}
double f(double c0, double A, double V, double l) {
        double r4617300 = c0;
        double r4617301 = A;
        double r4617302 = V;
        double r4617303 = l;
        double r4617304 = r4617302 * r4617303;
        double r4617305 = r4617301 / r4617304;
        double r4617306 = sqrt(r4617305);
        double r4617307 = r4617300 * r4617306;
        return r4617307;
}
↓
double f(double c0, double A, double V, double l) {
        double r4617308 = l;
        double r4617309 = V;
        double r4617310 = r4617308 * r4617309;
        double r4617311 = -inf.0;
        bool r4617312 = r4617310 <= r4617311;
        double r4617313 = A;
        double r4617314 = r4617313 / r4617308;
        double r4617315 = sqrt(r4617314);
        double r4617316 = 1.0;
        double r4617317 = r4617316 / r4617309;
        double r4617318 = sqrt(r4617317);
        double r4617319 = r4617315 * r4617318;
        double r4617320 = c0;
        double r4617321 = r4617319 * r4617320;
        double r4617322 = -5.601891409264343e-287;
        bool r4617323 = r4617310 <= r4617322;
        double r4617324 = r4617316 / r4617310;
        double r4617325 = cbrt(r4617324);
        double r4617326 = cbrt(r4617313);
        double r4617327 = r4617325 * r4617326;
        double r4617328 = sqrt(r4617327);
        double r4617329 = fabs(r4617326);
        double r4617330 = cbrt(r4617310);
        double r4617331 = fabs(r4617330);
        double r4617332 = r4617329 / r4617331;
        double r4617333 = r4617320 * r4617332;
        double r4617334 = r4617328 * r4617333;
        double r4617335 = 2.4427741881377e-317;
        bool r4617336 = r4617310 <= r4617335;
        double r4617337 = sqrt(r4617313);
        double r4617338 = sqrt(r4617310);
        double r4617339 = r4617337 / r4617338;
        double r4617340 = r4617320 * r4617339;
        double r4617341 = r4617336 ? r4617321 : r4617340;
        double r4617342 = r4617323 ? r4617334 : r4617341;
        double r4617343 = r4617312 ? r4617321 : r4617342;
        return r4617343;
}

Error

Bits error versus c0

Bits error versus A

Bits error versus V

Bits error versus l

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 3 regimes

  2. if (* V l) < -inf.0 or -5.601891409264343e-287 < (* V l) < 2.4427741881377e-317

    1. Initial program 53.0

      \[c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\]
    2. Using strategy rm

    3. Applied *-un-lft-identity53.0

      \[\leadsto c0 \cdot \sqrt{\frac{\color{blue}{1 \cdot A}}{V \cdot \ell}}\]

    4. Applied times-frac30.9

      \[\leadsto c0 \cdot \sqrt{\color{blue}{\frac{1}{V} \cdot \frac{A}{\ell}}}\]

    5. Applied sqrt-prod38.2

      \[\leadsto c0 \cdot \color{blue}{\left(\sqrt{\frac{1}{V}} \cdot \sqrt{\frac{A}{\ell}}\right)}\]

    if -inf.0 < (* V l) < -5.601891409264343e-287

    1. Initial program 9.7

      \[c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\]
    2. Using strategy rm

    3. Applied add-cube-cbrt10.2

      \[\leadsto c0 \cdot \sqrt{\color{blue}{\left(\sqrt[3]{\frac{A}{V \cdot \ell}} \cdot \sqrt[3]{\frac{A}{V \cdot \ell}}\right) \cdot \sqrt[3]{\frac{A}{V \cdot \ell}}}}\]

    4. Applied sqrt-prod10.2

      \[\leadsto c0 \cdot \color{blue}{\left(\sqrt{\sqrt[3]{\frac{A}{V \cdot \ell}} \cdot \sqrt[3]{\frac{A}{V \cdot \ell}}} \cdot \sqrt{\sqrt[3]{\frac{A}{V \cdot \ell}}}\right)}\]

    5. Applied associate-*r*10.2

      \[\leadsto \color{blue}{\left(c0 \cdot \sqrt{\sqrt[3]{\frac{A}{V \cdot \ell}} \cdot \sqrt[3]{\frac{A}{V \cdot \ell}}}\right) \cdot \sqrt{\sqrt[3]{\frac{A}{V \cdot \ell}}}}\]
    6. Using strategy rm

    7. Applied div-inv10.2

      \[\leadsto \left(c0 \cdot \sqrt{\sqrt[3]{\frac{A}{V \cdot \ell}} \cdot \sqrt[3]{\frac{A}{V \cdot \ell}}}\right) \cdot \sqrt{\sqrt[3]{\color{blue}{A \cdot \frac{1}{V \cdot \ell}}}}\]

    8. Applied cbrt-prod10.1

      \[\leadsto \left(c0 \cdot \sqrt{\sqrt[3]{\frac{A}{V \cdot \ell}} \cdot \sqrt[3]{\frac{A}{V \cdot \ell}}}\right) \cdot \sqrt{\color{blue}{\sqrt[3]{A} \cdot \sqrt[3]{\frac{1}{V \cdot \ell}}}}\]
    9. Using strategy rm

    10. Applied cbrt-div10.1

      \[\leadsto \left(c0 \cdot \sqrt{\sqrt[3]{\frac{A}{V \cdot \ell}} \cdot \color{blue}{\frac{\sqrt[3]{A}}{\sqrt[3]{V \cdot \ell}}}}\right) \cdot \sqrt{\sqrt[3]{A} \cdot \sqrt[3]{\frac{1}{V \cdot \ell}}}\]

    11. Applied cbrt-div2.6

      \[\leadsto \left(c0 \cdot \sqrt{\color{blue}{\frac{\sqrt[3]{A}}{\sqrt[3]{V \cdot \ell}}} \cdot \frac{\sqrt[3]{A}}{\sqrt[3]{V \cdot \ell}}}\right) \cdot \sqrt{\sqrt[3]{A} \cdot \sqrt[3]{\frac{1}{V \cdot \ell}}}\]

    12. Applied frac-times2.6

      \[\leadsto \left(c0 \cdot \sqrt{\color{blue}{\frac{\sqrt[3]{A} \cdot \sqrt[3]{A}}{\sqrt[3]{V \cdot \ell} \cdot \sqrt[3]{V \cdot \ell}}}}\right) \cdot \sqrt{\sqrt[3]{A} \cdot \sqrt[3]{\frac{1}{V \cdot \ell}}}\]

    13. Applied sqrt-div1.0

      \[\leadsto \left(c0 \cdot \color{blue}{\frac{\sqrt{\sqrt[3]{A} \cdot \sqrt[3]{A}}}{\sqrt{\sqrt[3]{V \cdot \ell} \cdot \sqrt[3]{V \cdot \ell}}}}\right) \cdot \sqrt{\sqrt[3]{A} \cdot \sqrt[3]{\frac{1}{V \cdot \ell}}}\]

    14. Simplified1.0

      \[\leadsto \left(c0 \cdot \frac{\color{blue}{\left|\sqrt[3]{A}\right|}}{\sqrt{\sqrt[3]{V \cdot \ell} \cdot \sqrt[3]{V \cdot \ell}}}\right) \cdot \sqrt{\sqrt[3]{A} \cdot \sqrt[3]{\frac{1}{V \cdot \ell}}}\]

    15. Simplified1.0

      \[\leadsto \left(c0 \cdot \frac{\left|\sqrt[3]{A}\right|}{\color{blue}{\left|\sqrt[3]{V \cdot \ell}\right|}}\right) \cdot \sqrt{\sqrt[3]{A} \cdot \sqrt[3]{\frac{1}{V \cdot \ell}}}\]

    if 2.4427741881377e-317 < (* V l)

    1. Initial program 14.9

      \[c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\]
    2. Using strategy rm

    3. Applied sqrt-div6.6

      \[\leadsto c0 \cdot \color{blue}{\frac{\sqrt{A}}{\sqrt{V \cdot \ell}}}\]
  3. Recombined 3 regimes into one program.

  4. Final simplification9.8

    \[\leadsto \begin{array}{l} \mathbf{if}\;\ell \cdot V = -\infty:\\ \;\;\;\;\left(\sqrt{\frac{A}{\ell}} \cdot \sqrt{\frac{1}{V}}\right) \cdot c0\\ \mathbf{elif}\;\ell \cdot V \le -5.601891409264343123426103753339848933114 \cdot 10^{-287}:\\ \;\;\;\;\sqrt{\sqrt[3]{\frac{1}{\ell \cdot V}} \cdot \sqrt[3]{A}} \cdot \left(c0 \cdot \frac{\left|\sqrt[3]{A}\right|}{\left|\sqrt[3]{\ell \cdot V}\right|}\right)\\ \mathbf{elif}\;\ell \cdot V \le 2.442774188137666401114116722762846155888 \cdot 10^{-317}:\\ \;\;\;\;\left(\sqrt{\frac{A}{\ell}} \cdot \sqrt{\frac{1}{V}}\right) \cdot c0\\ \mathbf{else}:\\ \;\;\;\;c0 \cdot \frac{\sqrt{A}}{\sqrt{\ell \cdot V}}\\ \end{array}\]

Reproduce

herbie shell --seed 2019200 +o rules:numerics
(FPCore (c0 A V l)
  :name "Henrywood and Agarwal, Equation (3)"
  (* c0 (sqrt (/ A (* V l)))))