\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 r7950131 = a;
double r7950132 = r7950131 * r7950131;
double r7950133 = b;
double r7950134 = r7950133 * r7950133;
double r7950135 = r7950132 + r7950134;
double r7950136 = 2.0;
double r7950137 = pow(r7950135, r7950136);
double r7950138 = 4.0;
double r7950139 = r7950138 * r7950134;
double r7950140 = r7950137 + r7950139;
double r7950141 = 1.0;
double r7950142 = r7950140 - r7950141;
return r7950142;
}
double f(double a, double b) {
double r7950143 = 4.0;
double r7950144 = b;
double r7950145 = r7950144 * r7950144;
double r7950146 = r7950143 * r7950145;
double r7950147 = a;
double r7950148 = r7950147 * r7950147;
double r7950149 = r7950145 + r7950148;
double r7950150 = sqrt(r7950149);
double r7950151 = pow(r7950150, r7950143);
double r7950152 = r7950146 + r7950151;
double r7950153 = 1.0;
double r7950154 = r7950152 - r7950153;
return r7950154;
}



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