re \cdot re - im \cdot im
\left(im + re\right) \cdot \left(re - im\right)
(FPCore (re im) :precision binary64 (- (* re re) (* im im)))
(FPCore (re im) :precision binary64 (* (+ im re) (- re im)))
double code(double re, double im) {
return (re * re) - (im * im);
}
double code(double re, double im) {
return (im + re) * (re - im);
}



Bits error versus re



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