\left(\frac{\pi}{2} \cdot \frac{1}{b \cdot b - a \cdot a}\right) \cdot \left(\frac{1}{a} - \frac{1}{b}\right)\frac{0.5 \cdot \frac{\pi}{b \cdot a}}{b + a}(FPCore (a b) :precision binary64 (* (* (/ PI 2.0) (/ 1.0 (- (* b b) (* a a)))) (- (/ 1.0 a) (/ 1.0 b))))
(FPCore (a b) :precision binary64 (/ (* 0.5 (/ PI (* b a))) (+ b a)))
double code(double a, double b) {
return ((((double) M_PI) / 2.0) * (1.0 / ((b * b) - (a * a)))) * ((1.0 / a) - (1.0 / b));
}
double code(double a, double b) {
return (0.5 * (((double) M_PI) / (b * a))) / (b + a);
}














Bits error versus a














Bits error versus b
Results
| Alternative 1 | |
|---|---|
| Error | 0.3 |
| Cost | 832 |
| Alternative 2 | |
|---|---|
| Error | 0.3 |
| Cost | 1088 |
| Alternative 3 | |
|---|---|
| Error | 0.3 |
| Cost | 1216 |
| Alternative 4 | |
|---|---|
| Error | 0.3 |
| Cost | 1216 |
| Alternative 5 | |
|---|---|
| Error | 0.3 |
| Cost | 1216 |
| Alternative 6 | |
|---|---|
| Error | 0.3 |
| Cost | 1344 |
| Alternative 7 | |
|---|---|
| Error | 0.5 |
| Cost | 14144 |
| Alternative 8 | |
|---|---|
| Error | 0.5 |
| Cost | 14272 |
| Alternative 9 | |
|---|---|
| Error | 0.7 |
| Cost | 14272 |
| Alternative 10 | |
|---|---|
| Error | 0.8 |
| Cost | 22592 |
| Alternative 11 | |
|---|---|
| Error | 0.9 |
| Cost | 40384 |

Initial program 14.5
Simplified14.5
rmApplied difference-of-squares_binary64_7299.7
Applied *-un-lft-identity_binary64_7609.7
Applied times-frac_binary64_7669.2
Applied associate-*l*_binary64_7010.3
Simplified0.3
rmApplied associate-*l/_binary64_7030.3
Simplified0.3
Taylor expanded around 0 0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2021043
(FPCore (a b)
:name "NMSE Section 6.1 mentioned, B"
:precision binary64
(* (* (/ PI 2.0) (/ 1.0 (- (* b b) (* a a)))) (- (/ 1.0 a) (/ 1.0 b))))