re \cdot im + im \cdot re
re \cdot \left(im + im\right)
double f(double re, double im) {
double r2101 = re;
double r2102 = im;
double r2103 = r2101 * r2102;
double r2104 = r2102 * r2101;
double r2105 = r2103 + r2104;
return r2105;
}
double f(double re, double im) {
double r2106 = re;
double r2107 = im;
double r2108 = r2107 + r2107;
double r2109 = r2106 * r2108;
return r2109;
}



Bits error versus re



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