re \cdot im + im \cdot re
re \cdot \left(im + im\right)
double f(double re, double im) {
double r198841 = re;
double r198842 = im;
double r198843 = r198841 * r198842;
double r198844 = r198842 * r198841;
double r198845 = r198843 + r198844;
return r198845;
}
double f(double re, double im) {
double r198846 = re;
double r198847 = im;
double r198848 = r198847 + r198847;
double r198849 = r198846 * r198848;
return r198849;
}



Bits error versus re



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