re \cdot im + im \cdot re
im \cdot \left(re + re\right)
double f(double re, double im) {
double r98 = re;
double r99 = im;
double r100 = r98 * r99;
double r101 = r99 * r98;
double r102 = r100 + r101;
return r102;
}
double f(double re, double im) {
double r103 = im;
double r104 = re;
double r105 = r104 + r104;
double r106 = r103 * r105;
return r106;
}



Bits error versus re



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