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) (y_46_re * ((double) (x_46_im / ((double) sqrt(((double) (((double) (y_46_re * y_46_re)) + ((double) (y_46_im * y_46_im)))))))))) - ((double) (y_46_im * ((double) (x_46_re / ((double) sqrt(((double) (((double) (y_46_re * y_46_re)) + ((double) (y_46_im * y_46_im)))))))))))) / ((double) sqrt(((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 25.3
rmApplied add-sqr-sqrt25.3
Applied associate-/r*25.3
rmApplied div-sub25.3
Simplified24.2
Simplified22.6
Final simplification22.6
herbie shell --seed 2020192
(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))))