4 #include "base/cl_sysdep.h"
7 #include "cln/complex.h"
12 #include "complex/cl_C.h"
17 const cl_N exp (const cl_N& x)
21 // x = a+bi -> (exp a) mit (cos b) + i (sin b) multiplizieren:
22 // (complex (* (exp a) (cos b)) (* (exp a) (sin b)))
29 var const cl_R& a = realpart(x);
30 var const cl_R& b = imagpart(x);
31 var cos_sin_t unitvec = cos_sin(b); // (cos b), (sin b)
32 // Da b nicht = Fixnum 0 ist, ist auch sin(b) nicht = Fixnum 0.
33 var cl_R exp_a = exp(a); // (exp a)
34 return complex_C(exp_a * unitvec.cos, exp_a * unitvec.sin);