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



Bits error versus a



Bits error versus k



Bits error versus m
Results
Initial program 2.0
Simplified1.9
Final simplification1.9
herbie shell --seed 2020338
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))