re \cdot re - im \cdot im
re \cdot re - im \cdot im
double code(double re, double im) {
return ((double) (((double) (re * re)) - ((double) (im * im))));
}
double code(double re, double im) {
return ((double) (((double) (re * re)) - ((double) (im * im))));
}



Bits error versus re



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