\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 r66309159 = a;
double r66309160 = r66309159 * r66309159;
double r66309161 = b;
double r66309162 = r66309161 * r66309161;
double r66309163 = r66309160 + r66309162;
double r66309164 = 2.0;
double r66309165 = pow(r66309163, r66309164);
double r66309166 = 4.0;
double r66309167 = r66309166 * r66309162;
double r66309168 = r66309165 + r66309167;
double r66309169 = 1.0;
double r66309170 = r66309168 - r66309169;
return r66309170;
}
double f(double a, double b) {
double r66309171 = 4.0;
double r66309172 = b;
double r66309173 = r66309172 * r66309172;
double r66309174 = r66309171 * r66309173;
double r66309175 = a;
double r66309176 = r66309175 * r66309175;
double r66309177 = r66309173 + r66309176;
double r66309178 = sqrt(r66309177);
double r66309179 = pow(r66309178, r66309171);
double r66309180 = r66309174 + r66309179;
double r66309181 = 1.0;
double r66309182 = r66309180 - r66309181;
return r66309182;
}



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