double code(double a, double k, double m) {
return ((double) (((double) (a * ((double) pow(k, m)))) / ((double) (((double) (1.0 + ((double) (10.0 * k)))) + ((double) (k * k))))));
}
double code(double a, double k, double m) {
return ((double) (((double) (a * ((double) pow(k, m)))) / ((double) (1.0 + ((double) (k * ((double) (k + 10.0))))))));
}



Bits error versus a



Bits error versus k



Bits error versus m
Results
Initial program 2.0
Simplified2.0
rmApplied associate-*r/2.0
Final simplification2.0
herbie shell --seed 2020192
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))