\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 r7212957 = a;
double r7212958 = r7212957 * r7212957;
double r7212959 = b;
double r7212960 = r7212959 * r7212959;
double r7212961 = r7212958 + r7212960;
double r7212962 = 2.0;
double r7212963 = pow(r7212961, r7212962);
double r7212964 = 4.0;
double r7212965 = r7212964 * r7212960;
double r7212966 = r7212963 + r7212965;
double r7212967 = 1.0;
double r7212968 = r7212966 - r7212967;
return r7212968;
}
double f(double a, double b) {
double r7212969 = 4.0;
double r7212970 = b;
double r7212971 = r7212970 * r7212970;
double r7212972 = r7212969 * r7212971;
double r7212973 = a;
double r7212974 = r7212973 * r7212973;
double r7212975 = r7212971 + r7212974;
double r7212976 = sqrt(r7212975);
double r7212977 = pow(r7212976, r7212969);
double r7212978 = r7212972 + r7212977;
double r7212979 = 1.0;
double r7212980 = r7212978 - r7212979;
return r7212980;
}



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 2019143
(FPCore (a b)
:name "Bouland and Aaronson, Equation (26)"
(- (+ (pow (+ (* a a) (* b b)) 2) (* 4 (* b b))) 1))