\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 r7192966 = a;
double r7192967 = r7192966 * r7192966;
double r7192968 = b;
double r7192969 = r7192968 * r7192968;
double r7192970 = r7192967 + r7192969;
double r7192971 = 2.0;
double r7192972 = pow(r7192970, r7192971);
double r7192973 = 4.0;
double r7192974 = r7192973 * r7192969;
double r7192975 = r7192972 + r7192974;
double r7192976 = 1.0;
double r7192977 = r7192975 - r7192976;
return r7192977;
}
double f(double a, double b) {
double r7192978 = 4.0;
double r7192979 = b;
double r7192980 = r7192979 * r7192979;
double r7192981 = r7192978 * r7192980;
double r7192982 = a;
double r7192983 = r7192982 * r7192982;
double r7192984 = r7192980 + r7192983;
double r7192985 = sqrt(r7192984);
double r7192986 = pow(r7192985, r7192978);
double r7192987 = r7192981 + r7192986;
double r7192988 = 1.0;
double r7192989 = r7192987 - r7192988;
return r7192989;
}



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