e^{\log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.re - \tan^{-1}_* \frac{x.im}{x.re} \cdot y.im} \cdot \cos \left(\log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im + \tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right)e^{\log \left(\mathsf{hypot}\left(x.re, x.im\right)\right) \cdot y.re - \tan^{-1}_* \frac{x.im}{x.re} \cdot y.im} \cdot \cos \left(\log \left(\mathsf{hypot}\left(x.re, x.im\right)\right) \cdot y.im + \left(\sqrt[3]{\tan^{-1}_* \frac{x.im}{x.re}} \cdot \left(\left(\sqrt[3]{\sqrt[3]{\tan^{-1}_* \frac{x.im}{x.re}}} \cdot \sqrt[3]{\sqrt[3]{\tan^{-1}_* \frac{x.im}{x.re}}}\right) \cdot \sqrt[3]{\sqrt[3]{\tan^{-1}_* \frac{x.im}{x.re}}}\right)\right) \cdot \left(\sqrt[3]{\tan^{-1}_* \frac{x.im}{x.re}} \cdot y.re\right)\right)double f(double x_re, double x_im, double y_re, double y_im) {
double r13690 = x_re;
double r13691 = r13690 * r13690;
double r13692 = x_im;
double r13693 = r13692 * r13692;
double r13694 = r13691 + r13693;
double r13695 = sqrt(r13694);
double r13696 = log(r13695);
double r13697 = y_re;
double r13698 = r13696 * r13697;
double r13699 = atan2(r13692, r13690);
double r13700 = y_im;
double r13701 = r13699 * r13700;
double r13702 = r13698 - r13701;
double r13703 = exp(r13702);
double r13704 = r13696 * r13700;
double r13705 = r13699 * r13697;
double r13706 = r13704 + r13705;
double r13707 = cos(r13706);
double r13708 = r13703 * r13707;
return r13708;
}
double f(double x_re, double x_im, double y_re, double y_im) {
double r13709 = x_re;
double r13710 = x_im;
double r13711 = hypot(r13709, r13710);
double r13712 = log(r13711);
double r13713 = y_re;
double r13714 = r13712 * r13713;
double r13715 = atan2(r13710, r13709);
double r13716 = y_im;
double r13717 = r13715 * r13716;
double r13718 = r13714 - r13717;
double r13719 = exp(r13718);
double r13720 = r13712 * r13716;
double r13721 = cbrt(r13715);
double r13722 = cbrt(r13721);
double r13723 = r13722 * r13722;
double r13724 = r13723 * r13722;
double r13725 = r13721 * r13724;
double r13726 = r13721 * r13713;
double r13727 = r13725 * r13726;
double r13728 = r13720 + r13727;
double r13729 = cos(r13728);
double r13730 = r13719 * r13729;
return r13730;
}



Bits error versus x.re



Bits error versus x.im



Bits error versus y.re



Bits error versus y.im
Results
Initial program 33.2
rmApplied hypot-def19.3
rmApplied hypot-def2.9
rmApplied add-cube-cbrt2.9
Applied associate-*l*2.9
rmApplied add-cube-cbrt2.9
Final simplification2.9
herbie shell --seed 2020047 +o rules:numerics
(FPCore (x.re x.im y.re y.im)
:name "powComplex, real part"
:precision binary64
(* (exp (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im))) (cos (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re)))))