\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 r5420253 = a;
double r5420254 = r5420253 * r5420253;
double r5420255 = b;
double r5420256 = r5420255 * r5420255;
double r5420257 = r5420254 + r5420256;
double r5420258 = 2.0;
double r5420259 = pow(r5420257, r5420258);
double r5420260 = 4.0;
double r5420261 = r5420260 * r5420256;
double r5420262 = r5420259 + r5420261;
double r5420263 = 1.0;
double r5420264 = r5420262 - r5420263;
return r5420264;
}
double f(double a, double b) {
double r5420265 = 4.0;
double r5420266 = b;
double r5420267 = r5420266 * r5420266;
double r5420268 = r5420265 * r5420267;
double r5420269 = a;
double r5420270 = r5420269 * r5420269;
double r5420271 = r5420267 + r5420270;
double r5420272 = sqrt(r5420271);
double r5420273 = pow(r5420272, r5420265);
double r5420274 = r5420268 + r5420273;
double r5420275 = 1.0;
double r5420276 = r5420274 - r5420275;
return r5420276;
}



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