double code(double d, double h, double l, double M, double D) {
return ((double) (((double) (((double) pow(((double) (d / h)), ((double) (1.0 / 2.0)))) * ((double) pow(((double) (d / l)), ((double) (1.0 / 2.0)))))) * ((double) (1.0 - ((double) (((double) (((double) (1.0 / 2.0)) * ((double) pow(((double) (((double) (M * D)) / ((double) (2.0 * d)))), 2.0)))) * ((double) (h / l))))))));
}
double code(double d, double h, double l, double M, double D) {
return ((double) (((double) (((double) pow(((double) (((double) (((double) cbrt(d)) / ((double) cbrt(h)))) * ((double) (((double) cbrt(d)) / ((double) cbrt(h)))))), ((double) (1.0 / 2.0)))) * ((double) pow(((double) (((double) cbrt(d)) / ((double) cbrt(h)))), ((double) (1.0 / 2.0)))))) * ((double) (((double) (((double) pow(((double) (1.0 / ((double) (((double) cbrt(l)) * ((double) cbrt(l)))))), ((double) (1.0 / 2.0)))) * ((double) pow(((double) (d / ((double) cbrt(l)))), ((double) (1.0 / 2.0)))))) * ((double) (1.0 - ((double) (h * ((double) (((double) (1.0 / 2.0)) * ((double) (((double) pow(((double) (M * ((double) (D / ((double) (d * 2.0)))))), ((double) (2.0 / 2.0)))) * ((double) (((double) (((double) pow(M, ((double) (2.0 / 2.0)))) / ((double) (((double) cbrt(l)) * ((double) cbrt(l)))))) * ((double) (((double) pow(((double) (D / ((double) (d * 2.0)))), ((double) (2.0 / 2.0)))) / ((double) cbrt(l))))))))))))))))));
}



Bits error versus d



Bits error versus h



Bits error versus l



Bits error versus M



Bits error versus D
Results
Initial program 26.3
Simplified25.8
rmApplied add-cube-cbrt26.1
Applied add-cube-cbrt26.2
Applied times-frac26.2
Applied unpow-prod-down19.9
Simplified19.9
rmApplied add-cube-cbrt20.0
Applied *-un-lft-identity20.0
Applied times-frac20.0
Applied unpow-prod-down15.0
rmApplied *-un-lft-identity15.0
Applied sqr-pow15.0
Applied times-frac12.7
Simplified12.7
Simplified12.7
rmApplied add-cube-cbrt12.8
Applied unpow-prod-down12.8
Applied times-frac13.5
Final simplification13.5
herbie shell --seed 2020192
(FPCore (d h l M D)
:name "Henrywood and Agarwal, Equation (12)"
:precision binary64
(* (* (pow (/ d h) (/ 1.0 2.0)) (pow (/ d l) (/ 1.0 2.0))) (- 1.0 (* (* (/ 1.0 2.0) (pow (/ (* M D) (* 2.0 d)) 2.0)) (/ h l)))))