\left({\left(a \cdot a + b \cdot b\right)}^{2} + 4 \cdot \left(b \cdot b\right)\right) - 1\left(4 \cdot \left(b \cdot b\right) + {\left(\sqrt{b \cdot b + a \cdot a}\right)}^{4}\right) - 1double f(double a, double b) {
double r4559881 = a;
double r4559882 = r4559881 * r4559881;
double r4559883 = b;
double r4559884 = r4559883 * r4559883;
double r4559885 = r4559882 + r4559884;
double r4559886 = 2.0;
double r4559887 = pow(r4559885, r4559886);
double r4559888 = 4.0;
double r4559889 = r4559888 * r4559884;
double r4559890 = r4559887 + r4559889;
double r4559891 = 1.0;
double r4559892 = r4559890 - r4559891;
return r4559892;
}
double f(double a, double b) {
double r4559893 = 4.0;
double r4559894 = b;
double r4559895 = r4559894 * r4559894;
double r4559896 = r4559893 * r4559895;
double r4559897 = a;
double r4559898 = r4559897 * r4559897;
double r4559899 = r4559895 + r4559898;
double r4559900 = sqrt(r4559899);
double r4559901 = pow(r4559900, r4559893);
double r4559902 = r4559896 + r4559901;
double r4559903 = 1.0;
double r4559904 = r4559902 - r4559903;
return r4559904;
}



Bits error versus a



Bits error versus b
Results
Initial program 0.2
Simplified0.2
rmApplied add-sqr-sqrt0.2
Applied associate-*l*0.1
rmApplied add-sqr-sqrt0.1
Applied cube-unmult0.1
Applied pow10.1
Applied pow-prod-up0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019152
(FPCore (a b)
:name "Bouland and Aaronson, Equation (26)"
(- (+ (pow (+ (* a a) (* b b)) 2) (* 4 (* b b))) 1))