double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return ((double) (((double) (((double) (x_46_im * y_46_re)) - ((double) (x_46_re * y_46_im)))) / ((double) (((double) (y_46_re * y_46_re)) + ((double) (y_46_im * y_46_im))))));
}
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return ((double) (((double) (((double) (x_46_im / ((double) sqrt(((double) (((double) (y_46_re * y_46_re)) + ((double) (y_46_im * y_46_im)))))))) * ((double) (y_46_re / ((double) sqrt(((double) (((double) (y_46_re * y_46_re)) + ((double) (y_46_im * y_46_im)))))))))) - ((double) (x_46_re * ((double) (y_46_im / ((double) (((double) (y_46_re * y_46_re)) + ((double) (y_46_im * y_46_im))))))))));
}



Bits error versus x.re



Bits error versus x.im



Bits error versus y.re



Bits error versus y.im
Results
Initial program 26.1
rmApplied div-sub26.1
Simplified24.9
Simplified23.3
rmApplied add-sqr-sqrt23.3
Applied *-un-lft-identity23.3
Applied times-frac23.3
Applied associate-*r*23.4
Simplified23.4
Final simplification23.4
herbie shell --seed 2020180
(FPCore (x.re x.im y.re y.im)
:name "_divideComplex, imaginary part"
:precision binary64
(/ (- (* x.im y.re) (* x.re y.im)) (+ (* y.re y.re) (* y.im y.im))))