\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{\frac{\pi}{a}}{b}}{b + a}double f(double a, double b) {
double r51441 = atan2(1.0, 0.0);
double r51442 = 2.0;
double r51443 = r51441 / r51442;
double r51444 = 1.0;
double r51445 = b;
double r51446 = r51445 * r51445;
double r51447 = a;
double r51448 = r51447 * r51447;
double r51449 = r51446 - r51448;
double r51450 = r51444 / r51449;
double r51451 = r51443 * r51450;
double r51452 = r51444 / r51447;
double r51453 = r51444 / r51445;
double r51454 = r51452 - r51453;
double r51455 = r51451 * r51454;
return r51455;
}
double f(double a, double b) {
double r51456 = 0.5;
double r51457 = atan2(1.0, 0.0);
double r51458 = a;
double r51459 = r51457 / r51458;
double r51460 = b;
double r51461 = r51459 / r51460;
double r51462 = r51456 * r51461;
double r51463 = r51460 + r51458;
double r51464 = r51462 / r51463;
return r51464;
}



Bits error versus a



Bits error versus b
Results
Initial program 14.3
rmApplied difference-of-squares9.4
Applied *-un-lft-identity9.4
Applied times-frac9.0
Applied associate-*r*9.0
Simplified8.9
rmApplied associate-*l/8.9
Applied associate-*l/0.3
Taylor expanded around 0 0.2
rmApplied associate-/r*0.2
Final simplification0.2
herbie shell --seed 2019326 +o rules:numerics
(FPCore (a b)
:name "NMSE Section 6.1 mentioned, B"
:precision binary64
(* (* (/ PI 2) (/ 1 (- (* b b) (* a a)))) (- (/ 1 a) (/ 1 b))))