re \cdot im + im \cdot re
re \cdot \left(im + im\right)
double f(double re, double im) {
double r8441 = re;
double r8442 = im;
double r8443 = r8441 * r8442;
double r8444 = r8442 * r8441;
double r8445 = r8443 + r8444;
return r8445;
}
double f(double re, double im) {
double r8446 = re;
double r8447 = im;
double r8448 = r8447 + r8447;
double r8449 = r8446 * r8448;
return r8449;
}



Bits error versus re



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