\left({\left(a \cdot a + b \cdot b\right)}^{2} + 4 \cdot \left(b \cdot b\right)\right) - 1\left({\left(\sqrt{b \cdot b + a \cdot a}\right)}^{4} + -1\right) - -4 \cdot \left(b \cdot b\right)double f(double a, double b) {
double r9132024 = a;
double r9132025 = r9132024 * r9132024;
double r9132026 = b;
double r9132027 = r9132026 * r9132026;
double r9132028 = r9132025 + r9132027;
double r9132029 = 2.0;
double r9132030 = pow(r9132028, r9132029);
double r9132031 = 4.0;
double r9132032 = r9132031 * r9132027;
double r9132033 = r9132030 + r9132032;
double r9132034 = 1.0;
double r9132035 = r9132033 - r9132034;
return r9132035;
}
double f(double a, double b) {
double r9132036 = b;
double r9132037 = r9132036 * r9132036;
double r9132038 = a;
double r9132039 = r9132038 * r9132038;
double r9132040 = r9132037 + r9132039;
double r9132041 = sqrt(r9132040);
double r9132042 = 4.0;
double r9132043 = pow(r9132041, r9132042);
double r9132044 = -1.0;
double r9132045 = r9132043 + r9132044;
double r9132046 = -4.0;
double r9132047 = r9132046 * r9132037;
double r9132048 = r9132045 - r9132047;
return r9132048;
}



Bits error versus a



Bits error versus b
Results
Initial program 0.2
Simplified0.2
rmApplied add-sqr-sqrt0.2
Applied associate-*r*0.1
rmApplied pow10.1
Applied add-sqr-sqrt0.1
Applied pow30.1
Applied pow-prod-up0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019163
(FPCore (a b)
:name "Bouland and Aaronson, Equation (26)"
(- (+ (pow (+ (* a a) (* b b)) 2) (* 4 (* b b))) 1))