src/complex/transcendental/cl_C_cis.cc

   1 // cis().
   2
   3 // General includes.
   4 #include "base/cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cln/complex.h"
   8
   9
  10 // Implementation.
  11
  12 #include "complex/cl_C.h"
  13 #include "cln/real.h"
  14
  15 namespace cln {
  16
  17 const cl_N cis (const cl_N& x)
  18 {
  19 // Methode:
  20 // x reell -> (complex (cos x) (sin x))
  21 // x = a+bi -> (complex (* (exp (- b)) (cos a)) (* (exp (- b)) (sin a)))
  22         if (realp(x)) {
  23                 DeclareType(cl_R,x);
  24                 var cos_sin_t trig = cos_sin(x);
  25                 return complex(trig.cos, trig.sin);
  26         } else {
  27                 DeclareType(cl_C,x);
  28                 // x=a+bi
  29                 var const cl_R& a = realpart(x);
  30                 var const cl_R& b = imagpart(x);
  31                 var cos_sin_t trig_a = cos_sin(a); // cos(a), sin(a) errechnen
  32                 var cl_R exp_minusb = exp(-b); // (exp (- b))
  33                 return complex(exp_minusb*trig_a.cos, // (* (exp (- b)) (cos a))
  34                                exp_minusb*trig_a.sin); // (* (exp (- b)) (sin a))
  35         }
  36 }
  37
  38 }  // namespace cln