Average Error: 1.1 → 1.1
Time: 12.7s
Precision: 64
\[\frac{\left(\left(x.im \cdot y.re\right) - \left(x.re \cdot y.im\right)\right)}{\left(\frac{\left(y.re \cdot y.re\right)}{\left(y.im \cdot y.im\right)}\right)}\]
↓
\[\frac{\frac{x.im \cdot y.re + x.re \cdot y.im}{\frac{x.im \cdot y.re + x.re \cdot y.im}{x.im \cdot y.re - x.re \cdot y.im}}}{y.re \cdot y.re + y.im \cdot y.im}\]
\frac{\left(\left(x.im \cdot y.re\right) - \left(x.re \cdot y.im\right)\right)}{\left(\frac{\left(y.re \cdot y.re\right)}{\left(y.im \cdot y.im\right)}\right)}
↓
\frac{\frac{x.im \cdot y.re + x.re \cdot y.im}{\frac{x.im \cdot y.re + x.re \cdot y.im}{x.im \cdot y.re - x.re \cdot y.im}}}{y.re \cdot y.re + y.im \cdot y.im}
double f(double x_re, double x_im, double y_re, double y_im) {
        double r1474893 = x_im;
        double r1474894 = y_re;
        double r1474895 = r1474893 * r1474894;
        double r1474896 = x_re;
        double r1474897 = y_im;
        double r1474898 = r1474896 * r1474897;
        double r1474899 = r1474895 - r1474898;
        double r1474900 = r1474894 * r1474894;
        double r1474901 = r1474897 * r1474897;
        double r1474902 = r1474900 + r1474901;
        double r1474903 = r1474899 / r1474902;
        return r1474903;
}
↓
double f(double x_re, double x_im, double y_re, double y_im) {
        double r1474904 = x_im;
        double r1474905 = y_re;
        double r1474906 = r1474904 * r1474905;
        double r1474907 = x_re;
        double r1474908 = y_im;
        double r1474909 = r1474907 * r1474908;
        double r1474910 = r1474906 + r1474909;
        double r1474911 = r1474906 - r1474909;
        double r1474912 = r1474910 / r1474911;
        double r1474913 = r1474910 / r1474912;
        double r1474914 = r1474905 * r1474905;
        double r1474915 = r1474908 * r1474908;
        double r1474916 = r1474914 + r1474915;
        double r1474917 = r1474913 / r1474916;
        return r1474917;
}

Error

Bits error versus x.re

Bits error versus x.im

Bits error versus y.re

Bits error versus y.im

Derivation

  1. Initial program 1.1

    \[\frac{\left(\left(x.im \cdot y.re\right) - \left(x.re \cdot y.im\right)\right)}{\left(\frac{\left(y.re \cdot y.re\right)}{\left(y.im \cdot y.im\right)}\right)}\]
  2. Using strategy rm

  3. Applied p16-flip--2.1

    \[\leadsto \frac{\color{blue}{\left(\frac{\left(\left(\left(x.im \cdot y.re\right) \cdot \left(x.im \cdot y.re\right)\right) - \left(\left(x.re \cdot y.im\right) \cdot \left(x.re \cdot y.im\right)\right)\right)}{\left(\frac{\left(x.im \cdot y.re\right)}{\left(x.re \cdot y.im\right)}\right)}\right)}}{\left(\frac{\left(y.re \cdot y.re\right)}{\left(y.im \cdot y.im\right)}\right)}\]
  4. Using strategy rm

  5. Applied difference-of-squares2.0

    \[\leadsto \frac{\left(\frac{\color{blue}{\left(\left(\frac{\left(x.im \cdot y.re\right)}{\left(x.re \cdot y.im\right)}\right) \cdot \left(\left(x.im \cdot y.re\right) - \left(x.re \cdot y.im\right)\right)\right)}}{\left(\frac{\left(x.im \cdot y.re\right)}{\left(x.re \cdot y.im\right)}\right)}\right)}{\left(\frac{\left(y.re \cdot y.re\right)}{\left(y.im \cdot y.im\right)}\right)}\]

  6. Applied associate-/l*1.1

    \[\leadsto \frac{\color{blue}{\left(\frac{\left(\frac{\left(x.im \cdot y.re\right)}{\left(x.re \cdot y.im\right)}\right)}{\left(\frac{\left(\frac{\left(x.im \cdot y.re\right)}{\left(x.re \cdot y.im\right)}\right)}{\left(\left(x.im \cdot y.re\right) - \left(x.re \cdot y.im\right)\right)}\right)}\right)}}{\left(\frac{\left(y.re \cdot y.re\right)}{\left(y.im \cdot y.im\right)}\right)}\]

  7. Final simplification1.1

    \[\leadsto \frac{\frac{x.im \cdot y.re + x.re \cdot y.im}{\frac{x.im \cdot y.re + x.re \cdot y.im}{x.im \cdot y.re - x.re \cdot y.im}}}{y.re \cdot y.re + y.im \cdot y.im}\]

Reproduce

herbie shell --seed 2019128 
(FPCore (x.re x.im y.re y.im)
  :name "_divideComplex, imaginary part"
  (/.p16 (-.p16 (*.p16 x.im y.re) (*.p16 x.re y.im)) (+.p16 (*.p16 y.re y.re) (*.p16 y.im y.im))))