c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\begin{array}{l}
\mathbf{if}\;\frac{A}{V \cdot \ell} \le 4.189537493318518753431675781489977417082 \cdot 10^{-310}:\\
\;\;\;\;\left(c0 \cdot \left|\frac{\sqrt[3]{A}}{\sqrt[3]{\ell}}\right|\right) \cdot \left(\sqrt{\frac{\frac{\sqrt[3]{\sqrt[3]{A} \cdot \sqrt[3]{A}}}{\sqrt[3]{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}}{\sqrt[3]{V} \cdot \sqrt[3]{V}}} \cdot \sqrt{\frac{\frac{\sqrt[3]{\sqrt[3]{A}}}{\sqrt[3]{\sqrt[3]{\ell}}}}{\sqrt[3]{V}}}\right)\\
\mathbf{elif}\;\frac{A}{V \cdot \ell} \le 8.02558393739106920381790045636343785553 \cdot 10^{288}:\\
\;\;\;\;c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(c0 \cdot \left|\frac{\sqrt[3]{A}}{\sqrt[3]{\ell}}\right|\right) \cdot \sqrt{\frac{\sqrt[3]{\frac{\sqrt[3]{A}}{\sqrt[3]{\ell}}} \cdot \sqrt[3]{\frac{\sqrt[3]{A}}{\sqrt[3]{\ell}}}}{\sqrt[3]{V} \cdot \sqrt[3]{V}}}\right) \cdot \sqrt{\frac{\sqrt[3]{\frac{\sqrt[3]{A}}{\sqrt[3]{\ell}}}}{\sqrt[3]{V}}}\\
\end{array}double f(double c0, double A, double V, double l) {
double r101936 = c0;
double r101937 = A;
double r101938 = V;
double r101939 = l;
double r101940 = r101938 * r101939;
double r101941 = r101937 / r101940;
double r101942 = sqrt(r101941);
double r101943 = r101936 * r101942;
return r101943;
}
double f(double c0, double A, double V, double l) {
double r101944 = A;
double r101945 = V;
double r101946 = l;
double r101947 = r101945 * r101946;
double r101948 = r101944 / r101947;
double r101949 = 4.1895374933185e-310;
bool r101950 = r101948 <= r101949;
double r101951 = c0;
double r101952 = cbrt(r101944);
double r101953 = cbrt(r101946);
double r101954 = r101952 / r101953;
double r101955 = fabs(r101954);
double r101956 = r101951 * r101955;
double r101957 = r101952 * r101952;
double r101958 = cbrt(r101957);
double r101959 = r101953 * r101953;
double r101960 = cbrt(r101959);
double r101961 = r101958 / r101960;
double r101962 = cbrt(r101945);
double r101963 = r101962 * r101962;
double r101964 = r101961 / r101963;
double r101965 = sqrt(r101964);
double r101966 = cbrt(r101952);
double r101967 = cbrt(r101953);
double r101968 = r101966 / r101967;
double r101969 = r101968 / r101962;
double r101970 = sqrt(r101969);
double r101971 = r101965 * r101970;
double r101972 = r101956 * r101971;
double r101973 = 8.025583937391069e+288;
bool r101974 = r101948 <= r101973;
double r101975 = sqrt(r101948);
double r101976 = r101951 * r101975;
double r101977 = cbrt(r101954);
double r101978 = r101977 * r101977;
double r101979 = r101978 / r101963;
double r101980 = sqrt(r101979);
double r101981 = r101956 * r101980;
double r101982 = r101977 / r101962;
double r101983 = sqrt(r101982);
double r101984 = r101981 * r101983;
double r101985 = r101974 ? r101976 : r101984;
double r101986 = r101950 ? r101972 : r101985;
return r101986;
}



Bits error versus c0



Bits error versus A



Bits error versus V



Bits error versus l
Results
if (/ A (* V l)) < 4.1895374933185e-310Initial program 41.4
rmApplied *-un-lft-identity41.4
Applied times-frac31.5
rmApplied associate-*l/31.5
Simplified31.5
rmApplied *-un-lft-identity31.5
Applied add-cube-cbrt31.6
Applied add-cube-cbrt31.6
Applied times-frac31.6
Applied times-frac31.3
Applied sqrt-prod11.2
Applied associate-*r*11.5
Simplified8.5
rmApplied add-cube-cbrt8.6
Applied add-cube-cbrt8.6
Applied cbrt-prod8.6
Applied add-cube-cbrt8.6
Applied cbrt-prod8.6
Applied times-frac8.6
Applied times-frac8.6
Applied sqrt-prod1.3
if 4.1895374933185e-310 < (/ A (* V l)) < 8.025583937391069e+288Initial program 0.4
if 8.025583937391069e+288 < (/ A (* V l)) Initial program 59.7
rmApplied *-un-lft-identity59.7
Applied times-frac47.8
rmApplied associate-*l/47.8
Simplified47.8
rmApplied *-un-lft-identity47.8
Applied add-cube-cbrt47.9
Applied add-cube-cbrt48.0
Applied times-frac48.0
Applied times-frac47.1
Applied sqrt-prod17.9
Applied associate-*r*18.1
Simplified14.0
rmApplied add-cube-cbrt14.1
Applied add-cube-cbrt14.2
Applied times-frac14.2
Applied sqrt-prod1.7
Applied associate-*r*1.7
Final simplification0.8
herbie shell --seed 2019323
(FPCore (c0 A V l)
:name "Henrywood and Agarwal, Equation (3)"
:precision binary64
(* c0 (sqrt (/ A (* V l)))))