src/complex/transcendental/cl_C_tan.cc

   1 // tan().
   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 tan (const cl_N& x)
  18 {
  19 // Methode:
  20 // x reell -> (/ (sin x) (cos x))
  21 // x = a+bi -> (/ (complex (* (sin a) (cosh b)) (* (cos a) (sinh b)))
  22 //                (complex (* (cos a) (cosh b)) (- (* (sin a) (sinh b)))) )
  23         if (realp(x)) {
  24                 DeclareType(cl_R,x);
  25                 var cos_sin_t trig = cos_sin(x);
  26                 return trig.sin / trig.cos;
  27         } else {
  28                 DeclareType(cl_C,x);
  29                 // x=a+bi
  30                 var const cl_R& a = realpart(x);
  31                 var const cl_R& b = imagpart(x);
  32                 var cosh_sinh_t hyp_b = cosh_sinh(b); // cosh(b), sinh(b) errechnen
  33                 var cos_sin_t trig_a = cos_sin(a); // cos(a), sin(a) errechnen
  34                 return
  35                         complex_C(trig_a.sin * hyp_b.cosh, // sin(a)*cosh(b)
  36                                   trig_a.cos * hyp_b.sinh // cos(a)*sinh(b), nicht Fixnum 0
  37                                  )
  38                       / complex(trig_a.cos * hyp_b.cosh, // cos(a)*cosh(b)
  39                                 - (trig_a.sin * hyp_b.sinh) // -sin(a)*sinh(b)
  40                                );
  41         }
  42 }
  43
  44 }  // namespace cln