re \cdot re + im \cdot im
re \cdot re + im \cdot im
(FPCore modulus_sqr (re im) :precision binary64 (+ (* re re) (* im im)))
(FPCore modulus_sqr (re im) :precision binary64 (+ (* re re) (* im im)))
double modulus_sqr(double re, double im) {
return (re * re) + (im * im);
}
double modulus_sqr(double re, double im) {
return (re * re) + (im * im);
}



Bits error versus re



Bits error versus im
Results
Initial program 0.0
Final simplification0.0
herbie shell --seed 2021207
(FPCore modulus_sqr (re im)
:name "math.abs on complex (squared)"
:precision binary64
(+ (* re re) (* im im)))