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 r166047 = c0;
        double r166048 = A;
        double r166049 = V;
        double r166050 = l;
        double r166051 = r166049 * r166050;
        double r166052 = r166048 / r166051;
        double r166053 = sqrt(r166052);
        double r166054 = r166047 * r166053;
        return r166054;
}
↓
double f(double c0, double A, double V, double l) {
        double r166055 = A;
        double r166056 = cbrt(r166055);
        double r166057 = V;
        double r166058 = cbrt(r166057);
        double r166059 = r166056 / r166058;
        double r166060 = l;
        double r166061 = cbrt(r166060);
        double r166062 = r166059 / r166061;
        double r166063 = fabs(r166062);
        double r166064 = c0;
        double r166065 = r166063 * r166064;
        double r166066 = sqrt(r166062);
        double r166067 = r166065 * r166066;
        return r166067;
}

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)))))