re \cdot re - im \cdot im
\left(im + re\right) \cdot \left(re - im\right)
double f(double re, double im) {
double r8151 = re;
double r8152 = r8151 * r8151;
double r8153 = im;
double r8154 = r8153 * r8153;
double r8155 = r8152 - r8154;
return r8155;
}
double f(double re, double im) {
double r8156 = im;
double r8157 = re;
double r8158 = r8156 + r8157;
double r8159 = r8157 - r8156;
double r8160 = r8158 * r8159;
return r8160;
}



Bits error versus re



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