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 \sin \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)\begin{array}{l}
\mathbf{if}\;x.re \leq -3.354388985317229 \cdot 10^{-26}:\\
\;\;\;\;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 \sin \left(y.im \cdot \log \left(-x.re\right) + y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right)\\
\mathbf{elif}\;x.re \leq -1.127312069324031 \cdot 10^{-250}:\\
\;\;\;\;\frac{{\left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right)}^{y.re}}{{\left(e^{\tan^{-1}_* \frac{x.im}{x.re}}\right)}^{y.im}} \cdot \sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re} + \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im\right)\\
\mathbf{elif}\;x.re \leq 3.0314558693474 \cdot 10^{-310}:\\
\;\;\;\;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 \sin \left(y.im \cdot \log \left(-x.re\right) + y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right)\\
\mathbf{elif}\;x.re \leq 2.2739435938255087 \cdot 10^{-107} \lor \neg \left(x.re \leq 4.699363203431519 \cdot 10^{+24}\right):\\
\;\;\;\;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 \sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re} + y.im \cdot \log x.re\right)\\
\mathbf{else}:\\
\;\;\;\;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 \left(\sqrt[3]{\sqrt[3]{\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re} + \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im\right)} \cdot \left(\sqrt[3]{\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re} + \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im\right)} \cdot \sqrt[3]{\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re} + \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im\right)}\right)} \cdot \left(\sqrt[3]{\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re} + \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im\right)} \cdot \sqrt[3]{\sqrt[3]{\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re} + \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im\right)} \cdot \left(\sqrt[3]{\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re} + \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im\right)} \cdot \sqrt[3]{\sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re} + \log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im\right)}\right)}\right)\right)\\
\end{array}(FPCore (x.re x.im y.re y.im)
:precision binary64
(*
(exp
(-
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re)
(* (atan2 x.im x.re) y.im)))
(sin
(+
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im)
(* (atan2 x.im x.re) y.re)))))(FPCore (x.re x.im y.re y.im)
:precision binary64
(if (<= x.re -3.354388985317229e-26)
(*
(exp
(-
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re)
(* (atan2 x.im x.re) y.im)))
(sin (+ (* y.im (log (- x.re))) (* y.re (atan2 x.im x.re)))))
(if (<= x.re -1.127312069324031e-250)
(*
(/
(pow (sqrt (+ (* x.re x.re) (* x.im x.im))) y.re)
(pow (exp (atan2 x.im x.re)) y.im))
(sin
(+
(* y.re (atan2 x.im x.re))
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im))))
(if (<= x.re 3.0314558693474e-310)
(*
(exp
(-
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re)
(* (atan2 x.im x.re) y.im)))
(sin (+ (* y.im (log (- x.re))) (* y.re (atan2 x.im x.re)))))
(if (or (<= x.re 2.2739435938255087e-107)
(not (<= x.re 4.699363203431519e+24)))
(*
(exp
(-
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re)
(* (atan2 x.im x.re) y.im)))
(sin (+ (* y.re (atan2 x.im x.re)) (* y.im (log x.re)))))
(*
(exp
(-
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re)
(* (atan2 x.im x.re) y.im)))
(*
(cbrt
(*
(cbrt
(sin
(+
(* y.re (atan2 x.im x.re))
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im))))
(*
(cbrt
(sin
(+
(* y.re (atan2 x.im x.re))
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im))))
(cbrt
(sin
(+
(* y.re (atan2 x.im x.re))
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im)))))))
(*
(cbrt
(sin
(+
(* y.re (atan2 x.im x.re))
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im))))
(cbrt
(*
(cbrt
(sin
(+
(* y.re (atan2 x.im x.re))
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im))))
(*
(cbrt
(sin
(+
(* y.re (atan2 x.im x.re))
(* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im))))
(cbrt
(sin
(+
(* y.re (atan2 x.im x.re))
(*
(log (sqrt (+ (* x.re x.re) (* x.im x.im))))
y.im)))))))))))))))double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
return ((double) (((double) exp(((double) (((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_re)) - ((double) (((double) atan2(x_46_im, x_46_re)) * y_46_im)))))) * ((double) sin(((double) (((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_im)) + ((double) (((double) atan2(x_46_im, x_46_re)) * y_46_re))))))));
}
double code(double x_46_re, double x_46_im, double y_46_re, double y_46_im) {
double tmp;
if ((x_46_re <= -3.354388985317229e-26)) {
tmp = ((double) (((double) exp(((double) (((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_re)) - ((double) (((double) atan2(x_46_im, x_46_re)) * y_46_im)))))) * ((double) sin(((double) (((double) (y_46_im * ((double) log(((double) -(x_46_re)))))) + ((double) (y_46_re * ((double) atan2(x_46_im, x_46_re))))))))));
} else {
double tmp_1;
if ((x_46_re <= -1.127312069324031e-250)) {
tmp_1 = ((double) ((((double) pow(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))), y_46_re)) / ((double) pow(((double) exp(((double) atan2(x_46_im, x_46_re)))), y_46_im))) * ((double) sin(((double) (((double) (y_46_re * ((double) atan2(x_46_im, x_46_re)))) + ((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_im))))))));
} else {
double tmp_2;
if ((x_46_re <= 3.0314558693474e-310)) {
tmp_2 = ((double) (((double) exp(((double) (((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_re)) - ((double) (((double) atan2(x_46_im, x_46_re)) * y_46_im)))))) * ((double) sin(((double) (((double) (y_46_im * ((double) log(((double) -(x_46_re)))))) + ((double) (y_46_re * ((double) atan2(x_46_im, x_46_re))))))))));
} else {
double tmp_3;
if (((x_46_re <= 2.2739435938255087e-107) || !(x_46_re <= 4.699363203431519e+24))) {
tmp_3 = ((double) (((double) exp(((double) (((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_re)) - ((double) (((double) atan2(x_46_im, x_46_re)) * y_46_im)))))) * ((double) sin(((double) (((double) (y_46_re * ((double) atan2(x_46_im, x_46_re)))) + ((double) (y_46_im * ((double) log(x_46_re))))))))));
} else {
tmp_3 = ((double) (((double) exp(((double) (((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_re)) - ((double) (((double) atan2(x_46_im, x_46_re)) * y_46_im)))))) * ((double) (((double) cbrt(((double) (((double) cbrt(((double) sin(((double) (((double) (y_46_re * ((double) atan2(x_46_im, x_46_re)))) + ((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_im)))))))) * ((double) (((double) cbrt(((double) sin(((double) (((double) (y_46_re * ((double) atan2(x_46_im, x_46_re)))) + ((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_im)))))))) * ((double) cbrt(((double) sin(((double) (((double) (y_46_re * ((double) atan2(x_46_im, x_46_re)))) + ((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_im)))))))))))))) * ((double) (((double) cbrt(((double) sin(((double) (((double) (y_46_re * ((double) atan2(x_46_im, x_46_re)))) + ((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_im)))))))) * ((double) cbrt(((double) (((double) cbrt(((double) sin(((double) (((double) (y_46_re * ((double) atan2(x_46_im, x_46_re)))) + ((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_im)))))))) * ((double) (((double) cbrt(((double) sin(((double) (((double) (y_46_re * ((double) atan2(x_46_im, x_46_re)))) + ((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_im)))))))) * ((double) cbrt(((double) sin(((double) (((double) (y_46_re * ((double) atan2(x_46_im, x_46_re)))) + ((double) (((double) log(((double) sqrt(((double) (((double) (x_46_re * x_46_re)) + ((double) (x_46_im * x_46_im)))))))) * y_46_im))))))))))))))))))));
}
tmp_2 = tmp_3;
}
tmp_1 = tmp_2;
}
tmp = tmp_1;
}
return tmp;
}



Bits error versus x.re



Bits error versus x.im



Bits error versus y.re



Bits error versus y.im
Results
if x.re < -3.3543889853172293e-26 or -1.1273120693240311e-250 < x.re < 3.03145586934742e-310Initial program Error: 37.2 bits
Taylor expanded around -inf Error: 21.7 bits
SimplifiedError: 21.7 bits
if -3.3543889853172293e-26 < x.re < -1.1273120693240311e-250Initial program Error: 21.4 bits
SimplifiedError: 24.9 bits
if 3.03145586934742e-310 < x.re < 2.2739435938255087e-107 or 4.699363203431519e24 < x.re Initial program Error: 39.5 bits
Taylor expanded around inf Error: 26.4 bits
SimplifiedError: 26.4 bits
if 2.2739435938255087e-107 < x.re < 4.699363203431519e24Initial program Error: 19.9 bits
rmApplied add-cube-cbrtError: 20.2 bits
rmApplied add-cube-cbrtError: 20.2 bits
rmApplied add-cube-cbrtError: 20.3 bits
Final simplificationError: 24.0 bits
herbie shell --seed 2020203
(FPCore (x.re x.im y.re y.im)
:name "powComplex, imaginary part"
:precision binary64
(* (exp (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im))) (sin (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re)))))