Extend the exponent range from 32 bits to 64 bits on selected platforms.

[cln.git] / src / complex / transcendental / cl_C_cis.cc

diff --git a/src/complex/transcendental/cl_C_cis.cc b/src/complex/transcendental/cl_C_cis.cc
index 0753a2729aaeccf036a3ddf01c518679c33d78d0..fc29e4313917db810451788fb8c7ee6694dfeb0d 100644 (file)
--- a/src/complex/transcendental/cl_C_cis.cc
+++ b/src/complex/transcendental/cl_C_cis.cc
@@ -4,13 +4,15 @@
 #include "cl_sysdep.h"
 
 // Specification.
-#include "cl_complex.h"
+#include "cln/complex.h"
 
 
 // Implementation.
 
 #include "cl_C.h"
-#include "cl_real.h"
+#include "cln/real.h"
+
+namespace cln {
 
 const cl_N cis (const cl_N& x)
 {
@@ -19,16 +21,18 @@ const cl_N cis (const cl_N& x)
 // x = a+bi -> (complex (* (exp (- b)) (cos a)) (* (exp (- b)) (sin a)))
        if (realp(x)) {
                DeclareType(cl_R,x);
-               var cl_cos_sin_t trig = cl_cos_sin(x);
+               var cos_sin_t trig = cos_sin(x);
                return complex(trig.cos, trig.sin);
        } else {
                DeclareType(cl_C,x);
                // x=a+bi
                var const cl_R& a = realpart(x);
                var const cl_R& b = imagpart(x);
-               var cl_cos_sin_t trig_a = cl_cos_sin(a); // cos(a), sin(a) errechnen
+               var cos_sin_t trig_a = cos_sin(a); // cos(a), sin(a) errechnen
                var cl_R exp_minusb = exp(-b); // (exp (- b))
                return complex(exp_minusb*trig_a.cos, // (* (exp (- b)) (cos a))
                               exp_minusb*trig_a.sin); // (* (exp (- b)) (sin a))
        }
 }
+
+}  // namespace cln