7 #include "cl_complex.h"
16 // Wert und Branch Cuts nach der Formel CLTL2, S. 307/312/313:
17 // arctan(z) = (log(1+iz)-log(1-iz)) / 2i
18 // Sei z=x+iy, errechne u+iv = artanh(-y+ix) wie oben, Ergebnis v-iu.
19 // Real- und Imaginärteil des Ergebnisses sind Floats, außer wenn z reell oder
22 inline const cl_C_R _atan (const cl_N& z)
29 return atanh(-imagpart(z),realpart(z));
33 const cl_N atan (const cl_N& z)
35 var cl_C_R u_v = _atan(z);
36 var cl_R& u = u_v.realpart;
37 var cl_R& v = u_v.imagpart;
38 return complex(v,-u); // v-iu