\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{1 \cdot \frac{\left(\left(b - a\right) \cdot 1\right) \cdot \pi}{b + a}}{2 \cdot \left(a \cdot b\right)}}{b - a}double f(double a, double b) {
double r2572661 = atan2(1.0, 0.0);
double r2572662 = 2.0;
double r2572663 = r2572661 / r2572662;
double r2572664 = 1.0;
double r2572665 = b;
double r2572666 = r2572665 * r2572665;
double r2572667 = a;
double r2572668 = r2572667 * r2572667;
double r2572669 = r2572666 - r2572668;
double r2572670 = r2572664 / r2572669;
double r2572671 = r2572663 * r2572670;
double r2572672 = r2572664 / r2572667;
double r2572673 = r2572664 / r2572665;
double r2572674 = r2572672 - r2572673;
double r2572675 = r2572671 * r2572674;
return r2572675;
}
double f(double a, double b) {
double r2572676 = 1.0;
double r2572677 = b;
double r2572678 = a;
double r2572679 = r2572677 - r2572678;
double r2572680 = r2572679 * r2572676;
double r2572681 = atan2(1.0, 0.0);
double r2572682 = r2572680 * r2572681;
double r2572683 = r2572677 + r2572678;
double r2572684 = r2572682 / r2572683;
double r2572685 = r2572676 * r2572684;
double r2572686 = 2.0;
double r2572687 = r2572678 * r2572677;
double r2572688 = r2572686 * r2572687;
double r2572689 = r2572685 / r2572688;
double r2572690 = r2572689 / r2572679;
return r2572690;
}



Bits error versus a



Bits error versus b
Results
Initial program 14.1
rmApplied difference-of-squares9.2
Applied *-un-lft-identity9.2
Applied times-frac8.8
Applied associate-*r*8.8
Simplified8.8
rmApplied associate-*r/8.7
Applied associate-*l/0.3
rmApplied frac-sub0.4
Applied associate-*l/0.4
Applied frac-times0.3
Simplified0.3
rmApplied associate-*r/0.5
Final simplification0.5
herbie shell --seed 2019172
(FPCore (a b)
:name "NMSE Section 6.1 mentioned, B"
(* (* (/ PI 2.0) (/ 1.0 (- (* b b) (* a a)))) (- (/ 1.0 a) (/ 1.0 b))))