\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 r69387202 = a;
double r69387203 = r69387202 * r69387202;
double r69387204 = b;
double r69387205 = r69387204 * r69387204;
double r69387206 = r69387203 + r69387205;
double r69387207 = 2.0;
double r69387208 = pow(r69387206, r69387207);
double r69387209 = 4.0;
double r69387210 = r69387209 * r69387205;
double r69387211 = r69387208 + r69387210;
double r69387212 = 1.0;
double r69387213 = r69387211 - r69387212;
return r69387213;
}
double f(double a, double b) {
double r69387214 = 4.0;
double r69387215 = b;
double r69387216 = r69387215 * r69387215;
double r69387217 = r69387214 * r69387216;
double r69387218 = a;
double r69387219 = r69387218 * r69387218;
double r69387220 = r69387216 + r69387219;
double r69387221 = sqrt(r69387220);
double r69387222 = pow(r69387221, r69387214);
double r69387223 = r69387217 + r69387222;
double r69387224 = 1.0;
double r69387225 = r69387223 - r69387224;
return r69387225;
}



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