re \cdot im + im \cdot re
im \cdot \left(re + re\right)
double f(double re, double im) {
double r56470 = re;
double r56471 = im;
double r56472 = r56470 * r56471;
double r56473 = r56471 * r56470;
double r56474 = r56472 + r56473;
return r56474;
}
double f(double re, double im) {
double r56475 = im;
double r56476 = re;
double r56477 = r56476 + r56476;
double r56478 = r56475 * r56477;
return r56478;
}



Bits error versus re



Bits error versus im
Results
Initial program 0.0
Simplified0.0
rmApplied distribute-rgt-out0
Final simplification0
herbie shell --seed 2019146
(FPCore (re im)
:name "math.square on complex, imaginary part"
(+ (* re im) (* im re)))