re \cdot im + im \cdot re
re \cdot \left(im + im\right)
double f(double re, double im) {
double r8247 = re;
double r8248 = im;
double r8249 = r8247 * r8248;
double r8250 = r8248 * r8247;
double r8251 = r8249 + r8250;
return r8251;
}
double f(double re, double im) {
double r8252 = re;
double r8253 = im;
double r8254 = r8253 + r8253;
double r8255 = r8252 * r8254;
return r8255;
}



Bits error versus re



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