c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\begin{array}{l}
\mathbf{if}\;V \cdot \ell \le -1.62535671440953908 \cdot 10^{68}:\\
\;\;\;\;\left(c0 \cdot \frac{\sqrt[3]{\sqrt{\sqrt[3]{A} \cdot \left(\frac{\sqrt[3]{A}}{\sqrt[3]{V}} \cdot \frac{\sqrt[3]{A}}{\ell}\right)}} \cdot \sqrt[3]{\sqrt{\sqrt[3]{A} \cdot \left(\frac{\sqrt[3]{A}}{\sqrt[3]{V}} \cdot \frac{\sqrt[3]{A}}{\ell}\right)}}}{\sqrt{\left|\sqrt[3]{V}\right|}}\right) \cdot \frac{\sqrt[3]{\sqrt{\sqrt[3]{A} \cdot \left(\frac{\sqrt[3]{A}}{\sqrt[3]{V}} \cdot \frac{\sqrt[3]{A}}{\ell}\right)}}}{\sqrt{\left|\sqrt[3]{V}\right|}}\\
\mathbf{elif}\;V \cdot \ell \le -5.05709688187036769 \cdot 10^{-237}:\\
\;\;\;\;c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}\\
\mathbf{elif}\;V \cdot \ell \le 0.0:\\
\;\;\;\;\left(c0 \cdot \frac{\sqrt[3]{\sqrt{\sqrt[3]{A} \cdot \left(\frac{\sqrt[3]{A}}{\sqrt[3]{V}} \cdot \frac{\sqrt[3]{A}}{\ell}\right)}} \cdot \sqrt[3]{\sqrt{\sqrt[3]{A} \cdot \left(\frac{\sqrt[3]{A}}{\sqrt[3]{V}} \cdot \frac{\sqrt[3]{A}}{\ell}\right)}}}{\sqrt{\left|\sqrt[3]{V}\right|}}\right) \cdot \frac{\sqrt[3]{\sqrt{\sqrt[3]{A} \cdot \left(\frac{\sqrt[3]{A}}{\sqrt[3]{V}} \cdot \frac{\sqrt[3]{A}}{\ell}\right)}}}{\sqrt{\left|\sqrt[3]{V}\right|}}\\
\mathbf{else}:\\
\;\;\;\;c0 \cdot \frac{\frac{\sqrt{A}}{\sqrt{\sqrt[3]{V} \cdot \ell}}}{\left|\sqrt[3]{V}\right|}\\
\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)) <= -1.625356714409539e+68)) {
VAR = ((double) (((double) (c0 * ((double) (((double) (((double) cbrt(((double) sqrt(((double) (((double) cbrt(A)) * ((double) (((double) (((double) cbrt(A)) / ((double) cbrt(V)))) * ((double) (((double) cbrt(A)) / l)))))))))) * ((double) cbrt(((double) sqrt(((double) (((double) cbrt(A)) * ((double) (((double) (((double) cbrt(A)) / ((double) cbrt(V)))) * ((double) (((double) cbrt(A)) / l)))))))))))) / ((double) sqrt(((double) fabs(((double) cbrt(V)))))))))) * ((double) (((double) cbrt(((double) sqrt(((double) (((double) cbrt(A)) * ((double) (((double) (((double) cbrt(A)) / ((double) cbrt(V)))) * ((double) (((double) cbrt(A)) / l)))))))))) / ((double) sqrt(((double) fabs(((double) cbrt(V))))))))));
} else {
double VAR_1;
if ((((double) (V * l)) <= -5.057096881870368e-237)) {
VAR_1 = ((double) (c0 * ((double) sqrt(((double) (A / ((double) (V * l))))))));
} else {
double VAR_2;
if ((((double) (V * l)) <= 0.0)) {
VAR_2 = ((double) (((double) (c0 * ((double) (((double) (((double) cbrt(((double) sqrt(((double) (((double) cbrt(A)) * ((double) (((double) (((double) cbrt(A)) / ((double) cbrt(V)))) * ((double) (((double) cbrt(A)) / l)))))))))) * ((double) cbrt(((double) sqrt(((double) (((double) cbrt(A)) * ((double) (((double) (((double) cbrt(A)) / ((double) cbrt(V)))) * ((double) (((double) cbrt(A)) / l)))))))))))) / ((double) sqrt(((double) fabs(((double) cbrt(V)))))))))) * ((double) (((double) cbrt(((double) sqrt(((double) (((double) cbrt(A)) * ((double) (((double) (((double) cbrt(A)) / ((double) cbrt(V)))) * ((double) (((double) cbrt(A)) / l)))))))))) / ((double) sqrt(((double) fabs(((double) cbrt(V))))))))));
} else {
VAR_2 = ((double) (c0 * ((double) (((double) (((double) sqrt(A)) / ((double) sqrt(((double) (((double) cbrt(V)) * l)))))) / ((double) fabs(((double) cbrt(V))))))));
}
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) < -1.625356714409539e+68 or -5.057096881870368e-237 < (* V l) < 0.0Initial program 32.2
rmApplied add-cube-cbrt32.4
Applied times-frac22.7
rmApplied add-cube-cbrt22.8
Applied times-frac22.8
Applied associate-*l*22.3
rmApplied associate-*l/22.8
Applied sqrt-div14.7
Simplified14.7
rmApplied add-sqr-sqrt14.7
Applied add-cube-cbrt15.0
Applied times-frac15.0
Applied associate-*r*15.0
if -1.625356714409539e+68 < (* V l) < -5.057096881870368e-237Initial program 7.1
if 0.0 < (* V l) Initial program 14.9
rmApplied add-cube-cbrt15.3
Applied times-frac17.1
rmApplied add-cube-cbrt17.2
Applied times-frac17.2
Applied associate-*l*13.8
rmApplied associate-*l/15.3
Applied sqrt-div10.8
Simplified10.8
rmApplied frac-times11.6
Applied associate-*r/11.6
Applied sqrt-div3.1
Simplified2.9
Final simplification7.9
herbie shell --seed 2020128
(FPCore (c0 A V l)
:name "Henrywood and Agarwal, Equation (3)"
:precision binary64
(* c0 (sqrt (/ A (* V l)))))