c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\begin{array}{l}
\mathbf{if}\;V \cdot \ell \le -5.8266415200593072 \cdot 10^{-268}:\\
\;\;\;\;\left(c0 \cdot \sqrt{\sqrt{\frac{A}{V \cdot \ell}}}\right) \cdot \sqrt{\left|\sqrt[3]{\frac{A}{V \cdot \ell}}\right| \cdot \sqrt{\sqrt[3]{\sqrt[3]{\frac{A}{V \cdot \ell}} \cdot \sqrt[3]{\frac{A}{V \cdot \ell}}} \cdot \sqrt[3]{\sqrt[3]{\frac{A}{V \cdot \ell}}}}}\\
\mathbf{elif}\;V \cdot \ell \le 5.1533916046317492 \cdot 10^{-288}:\\
\;\;\;\;\left(c0 \cdot \sqrt{\frac{1}{V}}\right) \cdot \sqrt{\frac{A}{\ell}}\\
\mathbf{elif}\;V \cdot \ell \le 1.0536369811328654 \cdot 10^{297}:\\
\;\;\;\;c0 \cdot \frac{\sqrt{A}}{\sqrt{V \cdot \ell}}\\
\mathbf{else}:\\
\;\;\;\;c0 \cdot \sqrt{\frac{\frac{A}{V}}{\ell}}\\
\end{array}double code(double c0, double A, double V, double l) {
return ((double) (c0 * ((double) sqrt(((double) (A / ((double) (V * l))))))));
}
double code(double c0, double A, double V, double l) {
double VAR;
if ((((double) (V * l)) <= -5.826641520059307e-268)) {
VAR = ((double) (((double) (c0 * ((double) sqrt(((double) sqrt(((double) (A / ((double) (V * l)))))))))) * ((double) sqrt(((double) (((double) fabs(((double) cbrt(((double) (A / ((double) (V * l)))))))) * ((double) sqrt(((double) (((double) cbrt(((double) (((double) cbrt(((double) (A / ((double) (V * l)))))) * ((double) cbrt(((double) (A / ((double) (V * l)))))))))) * ((double) cbrt(((double) cbrt(((double) (A / ((double) (V * l))))))))))))))))));
} else {
double VAR_1;
if ((((double) (V * l)) <= 5.153391604631749e-288)) {
VAR_1 = ((double) (((double) (c0 * ((double) sqrt(((double) (1.0 / V)))))) * ((double) sqrt(((double) (A / l))))));
} else {
double VAR_2;
if ((((double) (V * l)) <= 1.0536369811328654e+297)) {
VAR_2 = ((double) (c0 * ((double) (((double) sqrt(A)) / ((double) sqrt(((double) (V * l))))))));
} else {
VAR_2 = ((double) (c0 * ((double) sqrt(((double) (((double) (A / V)) / l))))));
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus c0



Bits error versus A



Bits error versus V



Bits error versus l
Results
if (* V l) < -5.826641520059307e-268Initial program 14.5
rmApplied add-sqr-sqrt14.5
Applied sqrt-prod14.7
Applied associate-*r*14.7
rmApplied add-cube-cbrt14.7
Applied sqrt-prod14.7
Simplified14.7
rmApplied add-cube-cbrt14.7
Applied cbrt-prod14.7
if -5.826641520059307e-268 < (* V l) < 5.153391604631749e-288Initial program 53.7
rmApplied *-un-lft-identity53.7
Applied times-frac35.4
Applied sqrt-prod39.9
Applied associate-*r*40.7
if 5.153391604631749e-288 < (* V l) < 1.0536369811328654e+297Initial program 9.3
rmApplied sqrt-div0.4
if 1.0536369811328654e+297 < (* V l) Initial program 38.0
rmApplied associate-/r*21.1
Final simplification12.9
herbie shell --seed 2020148
(FPCore (c0 A V l)
:name "Henrywood and Agarwal, Equation (3)"
:precision binary64
(* c0 (sqrt (/ A (* V l)))))