re \cdot im + im \cdot re
re \cdot \left(im + im\right)
double f(double re, double im) {
double r8251 = re;
double r8252 = im;
double r8253 = r8251 * r8252;
double r8254 = r8252 * r8251;
double r8255 = r8253 + r8254;
return r8255;
}
double f(double re, double im) {
double r8256 = re;
double r8257 = im;
double r8258 = r8257 + r8257;
double r8259 = r8256 * r8258;
return r8259;
}



Bits error versus re



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