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) (((double) pow(((double) (M * ((double) (D / ((double) (d * 2.0)))))), ((double) (2.0 / 2.0)))) / ((double) (((double) cbrt(l)) * ((double) cbrt(l)))))) * ((double) (((double) pow(((double) (M * ((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.8
Simplified26.3
rmApplied add-cube-cbrt26.6
Applied *-un-lft-identity26.6
Applied times-frac26.6
Applied unpow-prod-down21.4
rmApplied add-cube-cbrt21.6
Applied add-cube-cbrt21.6
Applied times-frac21.6
Applied unpow-prod-down15.5
Simplified15.5
rmApplied add-cube-cbrt15.6
Applied sqr-pow15.6
Applied times-frac13.2
Final simplification13.2
herbie shell --seed 2020191
(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)))))