c0 \cdot \sqrt{\frac{A}{V \cdot \ell}}
\begin{array}{l}
t_0 := \sqrt{\sqrt[3]{V} \cdot \sqrt[3]{V}}\\
\mathbf{if}\;\ell \leq -9.463703661952438 \cdot 10^{-308}:\\
\;\;\;\;c0 \cdot \frac{\frac{\sqrt{A}}{\sqrt{\ell \cdot \sqrt[3]{V}}}}{t_0}\\
\mathbf{else}:\\
\;\;\;\;c0 \cdot \frac{\frac{\sqrt{A \cdot \frac{1}{\sqrt[3]{V}}}}{\sqrt{\ell}}}{t_0}\\
\end{array}
(FPCore (c0 A V l) :precision binary64 (* c0 (sqrt (/ A (* V l)))))
(FPCore (c0 A V l)
:precision binary64
(let* ((t_0 (sqrt (* (cbrt V) (cbrt V)))))
(if (<= l -9.463703661952438e-308)
(* c0 (/ (/ (sqrt A) (sqrt (* l (cbrt V)))) t_0))
(* c0 (/ (/ (sqrt (* A (/ 1.0 (cbrt V)))) (sqrt l)) t_0)))))double code(double c0, double A, double V, double l) {
return c0 * sqrt(A / (V * l));
}
double code(double c0, double A, double V, double l) {
double t_0 = sqrt(cbrt(V) * cbrt(V));
double tmp;
if (l <= -9.463703661952438e-308) {
tmp = c0 * ((sqrt(A) / sqrt(l * cbrt(V))) / t_0);
} else {
tmp = c0 * ((sqrt(A * (1.0 / cbrt(V))) / sqrt(l)) / t_0);
}
return tmp;
}



Bits error versus c0



Bits error versus A



Bits error versus V



Bits error versus l
Results
if l < -9.46370366195243773e-308Initial program 18.0
Applied *-un-lft-identity_binary6418.0
Applied times-frac_binary6419.4
Applied add-cube-cbrt_binary6419.7
Applied *-un-lft-identity_binary6419.7
Applied times-frac_binary6419.7
Applied associate-*l*_binary6419.7
Applied associate-*l/_binary6419.7
Applied sqrt-div_binary6414.2
Applied frac-times_binary6413.2
Applied associate-*r/_binary6413.2
Applied sqrt-div_binary643.9
Simplified3.9
Simplified3.9
if -9.46370366195243773e-308 < l Initial program 19.2
Applied *-un-lft-identity_binary6419.2
Applied times-frac_binary6419.4
Applied add-cube-cbrt_binary6419.7
Applied *-un-lft-identity_binary6419.7
Applied times-frac_binary6419.7
Applied associate-*l*_binary6419.7
Applied associate-*l/_binary6419.7
Applied sqrt-div_binary6414.0
Applied associate-*r/_binary6413.3
Applied associate-*r/_binary6413.3
Applied sqrt-div_binary644.6
Final simplification4.4
herbie shell --seed 2021224
(FPCore (c0 A V l)
:name "Henrywood and Agarwal, Equation (3)"
:precision binary64
(* c0 (sqrt (/ A (* V l)))))