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) {
double VAR;
if ((k <= 7.430325758544677e+150)) {
VAR = ((double) (((double) (a / ((double) sqrt(((double) (1.0 + ((double) (k * ((double) (k + 10.0)))))))))) * ((double) (((double) pow(k, m)) / ((double) sqrt(((double) (1.0 + ((double) (k * ((double) (k + 10.0))))))))))));
} else {
VAR = ((double) (((double) (((double) (((double) pow(((double) exp(m)), ((double) log(k)))) / k)) * ((double) (a / k)))) + ((double) (((double) (((double) (((double) pow(((double) exp(m)), ((double) log(k)))) / k)) * ((double) (a / k)))) * ((double) (((double) (99.0 / ((double) (k * k)))) - ((double) (10.0 / k))))))));
}
return VAR;
}



Bits error versus a



Bits error versus k



Bits error versus m
Results
if k < 7.43032575854467656e150Initial program 0.1
rmApplied add-sqr-sqrt0.2
Applied times-frac0.2
Simplified0.2
Simplified0.2
if 7.43032575854467656e150 < k Initial program 9.5
rmApplied add-sqr-sqrt9.5
Applied times-frac9.5
Simplified9.5
Simplified9.5
Taylor expanded around inf 9.5
Simplified0.1
Final simplification0.2
herbie shell --seed 2020181
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1.0 (* 10.0 k)) (* k k))))