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