re \cdot im + im \cdot re
re \cdot \left(im + im\right)
double f(double re, double im) {
double r8421 = re;
double r8422 = im;
double r8423 = r8421 * r8422;
double r8424 = r8422 * r8421;
double r8425 = r8423 + r8424;
return r8425;
}
double f(double re, double im) {
double r8426 = re;
double r8427 = im;
double r8428 = r8427 + r8427;
double r8429 = r8426 * r8428;
return r8429;
}



Bits error versus re



Bits error versus im
Results
Initial program 0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019151 +o rules:numerics
(FPCore (re im)
:name "math.square on complex, imaginary part"
(+ (* re im) (* im re)))