\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 r44472 = atan2(1.0, 0.0);
double r44473 = 2.0;
double r44474 = r44472 / r44473;
double r44475 = 1.0;
double r44476 = b;
double r44477 = r44476 * r44476;
double r44478 = a;
double r44479 = r44478 * r44478;
double r44480 = r44477 - r44479;
double r44481 = r44475 / r44480;
double r44482 = r44474 * r44481;
double r44483 = r44475 / r44478;
double r44484 = r44475 / r44476;
double r44485 = r44483 - r44484;
double r44486 = r44482 * r44485;
return r44486;
}
double f(double a, double b) {
double r44487 = 0.5;
double r44488 = atan2(1.0, 0.0);
double r44489 = a;
double r44490 = r44488 / r44489;
double r44491 = b;
double r44492 = r44490 / r44491;
double r44493 = r44487 * r44492;
double r44494 = r44491 + r44489;
double r44495 = r44493 / r44494;
return r44495;
}



Bits error versus a



Bits error versus b
Results
Initial program 14.9
rmApplied difference-of-squares9.8
Applied *-un-lft-identity9.8
Applied times-frac9.3
Applied associate-*r*9.3
Simplified9.2
rmApplied associate-*l/9.2
Applied associate-*l/0.3
Taylor expanded around 0 0.2
rmApplied associate-/r*0.3
Final simplification0.3
herbie shell --seed 2020025 +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))))