\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{\pi}{2} \cdot \frac{1}{a \cdot b}}{b + a}double f(double a, double b) {
double r45771 = atan2(1.0, 0.0);
double r45772 = 2.0;
double r45773 = r45771 / r45772;
double r45774 = 1.0;
double r45775 = b;
double r45776 = r45775 * r45775;
double r45777 = a;
double r45778 = r45777 * r45777;
double r45779 = r45776 - r45778;
double r45780 = r45774 / r45779;
double r45781 = r45773 * r45780;
double r45782 = r45774 / r45777;
double r45783 = r45774 / r45775;
double r45784 = r45782 - r45783;
double r45785 = r45781 * r45784;
return r45785;
}
double f(double a, double b) {
double r45786 = atan2(1.0, 0.0);
double r45787 = 2.0;
double r45788 = r45786 / r45787;
double r45789 = 1.0;
double r45790 = a;
double r45791 = b;
double r45792 = r45790 * r45791;
double r45793 = r45789 / r45792;
double r45794 = r45788 * r45793;
double r45795 = r45791 + r45790;
double r45796 = r45794 / r45795;
return r45796;
}



Bits error versus a



Bits error versus b
Results
Initial program 14.7
rmApplied difference-of-squares9.7
Applied *-un-lft-identity9.7
Applied times-frac9.3
Applied associate-*r*9.3
Simplified9.2
rmApplied associate-*l*0.3
rmApplied associate-*l/0.3
Taylor expanded around 0 0.3
Final simplification0.3
herbie shell --seed 2020002
(FPCore (a b)
:name "NMSE Section 6.1 mentioned, B"
:precision binary64
(* (* (/ PI 2) (/ 1 (- (* b b) (* a a)))) (- (/ 1 a) (/ 1 b))))