Average Error: 19.5 → 1.1
Time: 6.3s
Precision: 64
\[c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\]
↓
\[\left(\left|\frac{\frac{\sqrt[3]{A}}{\sqrt[3]{V}}}{\sqrt[3]{\ell}}\right| \cdot c0\right) \cdot \sqrt{\frac{\frac{\sqrt[3]{A}}{\sqrt[3]{V}}}{\sqrt[3]{\ell}}}\]
c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}
↓
\left(\left|\frac{\frac{\sqrt[3]{A}}{\sqrt[3]{V}}}{\sqrt[3]{\ell}}\right| \cdot c0\right) \cdot \sqrt{\frac{\frac{\sqrt[3]{A}}{\sqrt[3]{V}}}{\sqrt[3]{\ell}}}
double f(double c0, double A, double V, double l) {
        double r163622 = c0;
        double r163623 = A;
        double r163624 = V;
        double r163625 = l;
        double r163626 = r163624 * r163625;
        double r163627 = r163623 / r163626;
        double r163628 = sqrt(r163627);
        double r163629 = r163622 * r163628;
        return r163629;
}
↓
double f(double c0, double A, double V, double l) {
        double r163630 = A;
        double r163631 = cbrt(r163630);
        double r163632 = V;
        double r163633 = cbrt(r163632);
        double r163634 = r163631 / r163633;
        double r163635 = l;
        double r163636 = cbrt(r163635);
        double r163637 = r163634 / r163636;
        double r163638 = fabs(r163637);
        double r163639 = c0;
        double r163640 = r163638 * r163639;
        double r163641 = sqrt(r163637);
        double r163642 = r163640 * r163641;
        return r163642;
}

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. Initial program 19.5

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

  3. Applied associate-/r*19.6

    \[\leadsto c0 \cdot \sqrt{\color{blue}{\frac{\frac{A}{V}}{\ell}}}\]
  4. Using strategy rm

  5. Applied add-cube-cbrt19.9

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

  6. Applied add-cube-cbrt20.1

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

  7. Applied add-cube-cbrt20.1

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

  8. Applied times-frac20.1

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

  9. Applied times-frac15.7

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

  10. Applied sqrt-prod7.2

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

  11. Applied associate-*r*7.2

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

  12. Simplified1.1

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

  13. Final simplification1.1

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

Reproduce

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