re \cdot im + im \cdot re
re \cdot \left(im + im\right)
double f(double re, double im) {
double r8406 = re;
double r8407 = im;
double r8408 = r8406 * r8407;
double r8409 = r8407 * r8406;
double r8410 = r8408 + r8409;
return r8410;
}
double f(double re, double im) {
double r8411 = re;
double r8412 = im;
double r8413 = r8412 + r8412;
double r8414 = r8411 * r8413;
return r8414;
}



Bits error versus re



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