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



Bits error versus re



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