\left(\frac{\pi}{2} \cdot \frac{1}{b \cdot b - a \cdot a}\right) \cdot \left(\frac{1}{a} - \frac{1}{b}\right)\frac{\frac{\frac{\pi}{2}}{a + b}}{a \cdot \left(b - a\right)} - \frac{\frac{\frac{\frac{\pi}{2}}{a + b}}{b}}{b - a}double f(double a, double b) {
double r1683370 = atan2(1.0, 0.0);
double r1683371 = 2.0;
double r1683372 = r1683370 / r1683371;
double r1683373 = 1.0;
double r1683374 = b;
double r1683375 = r1683374 * r1683374;
double r1683376 = a;
double r1683377 = r1683376 * r1683376;
double r1683378 = r1683375 - r1683377;
double r1683379 = r1683373 / r1683378;
double r1683380 = r1683372 * r1683379;
double r1683381 = r1683373 / r1683376;
double r1683382 = r1683373 / r1683374;
double r1683383 = r1683381 - r1683382;
double r1683384 = r1683380 * r1683383;
return r1683384;
}
double f(double a, double b) {
double r1683385 = atan2(1.0, 0.0);
double r1683386 = 2.0;
double r1683387 = r1683385 / r1683386;
double r1683388 = a;
double r1683389 = b;
double r1683390 = r1683388 + r1683389;
double r1683391 = r1683387 / r1683390;
double r1683392 = r1683389 - r1683388;
double r1683393 = r1683388 * r1683392;
double r1683394 = r1683391 / r1683393;
double r1683395 = r1683391 / r1683389;
double r1683396 = r1683395 / r1683392;
double r1683397 = r1683394 - r1683396;
return r1683397;
}



Bits error versus a



Bits error versus b
Results
Initial program 14.3
Simplified9.0
rmApplied div-inv9.1
Applied associate-/l*4.6
rmApplied associate-/l/0.3
rmApplied associate-/r/0.2
Simplified0.2
rmApplied associate-/r*0.2
Final simplification0.2
herbie shell --seed 2019153 +o rules:numerics
(FPCore (a b)
:name "NMSE Section 6.1 mentioned, B"
(* (* (/ PI 2) (/ 1 (- (* b b) (* a a)))) (- (/ 1 a) (/ 1 b))))