re \cdot im + im \cdot re
re \cdot \left(im + im\right)
double f(double re, double im) {
double r8449 = re;
double r8450 = im;
double r8451 = r8449 * r8450;
double r8452 = r8450 * r8449;
double r8453 = r8451 + r8452;
return r8453;
}
double f(double re, double im) {
double r8454 = re;
double r8455 = im;
double r8456 = r8455 + r8455;
double r8457 = r8454 * r8456;
return r8457;
}



Bits error versus re



Bits error versus im
Results
Initial program 0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019104 +o rules:numerics
(FPCore (re im)
:name "math.square on complex, imaginary part"
(+ (* re im) (* im re)))