src/complex/algebraic/cl_C_abs.cc

   1 // abs().
   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 #include "base/cl_inline.h"
  16 #include "complex/algebraic/cl_C_abs_aux.cc"
  17
  18 namespace cln {
  19
  20 const cl_R abs (const cl_N& x)
  21 {
  22 // Methode:
  23 // Falls x reell: klar
  24 // Falls x=a+bi: sqrt(a^2+b^2)
  25         if (realp(x)) {
  26                 DeclareType(cl_R,x);
  27                 return abs(x);
  28         } else {
  29                 DeclareType(cl_C,x);
  30                 return abs_inline(x);
  31         }
  32 }
  33
  34 }  // namespace cln