\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 r72583 = atan2(1.0, 0.0);
double r72584 = 2.0;
double r72585 = r72583 / r72584;
double r72586 = 1.0;
double r72587 = b;
double r72588 = r72587 * r72587;
double r72589 = a;
double r72590 = r72589 * r72589;
double r72591 = r72588 - r72590;
double r72592 = r72586 / r72591;
double r72593 = r72585 * r72592;
double r72594 = r72586 / r72589;
double r72595 = r72586 / r72587;
double r72596 = r72594 - r72595;
double r72597 = r72593 * r72596;
return r72597;
}
double f(double a, double b) {
double r72598 = 0.5;
double r72599 = atan2(1.0, 0.0);
double r72600 = a;
double r72601 = r72599 / r72600;
double r72602 = b;
double r72603 = r72601 / r72602;
double r72604 = r72598 * r72603;
double r72605 = r72602 + r72600;
double r72606 = r72604 / r72605;
return r72606;
}



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))))