src/complex/elem/cl_C_plus.cc

   1 // binary operator +
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_complex.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_C.h"
  13 #include "cl_real.h"
  14
  15 const cl_N operator+ (const cl_N& x, const cl_N& y)
  16 {
  17 // Methode:
  18 // x,y beide reell -> klar.
  19 // x=a, y=b+ci -> (a+b)+ci
  20 // x=a+bi, y=c -> (a+c)+bi
  21 // x=a+bi, y=c+di -> (a+c)+(b+d)i
  22         if (realp(x)) {
  23                 DeclareType(cl_R,x);
  24                 if (realp(y)) {
  25                         DeclareType(cl_R,y);
  26                         return x+y;
  27                 } else {
  28                         DeclareType(cl_C,y);
  29                         // x=a, y=b+ci
  30                         var const cl_R& a = x;
  31                         var const cl_R& b = realpart(y);
  32                         var const cl_R& c = imagpart(y);
  33                         return complex_C(a+b,c);
  34                 }
  35         } else {
  36                 DeclareType(cl_C,x);
  37                 if (realp(y)) {
  38                         DeclareType(cl_R,y);
  39                         // x=a+bi, y=c
  40                         var const cl_R& a = realpart(x);
  41                         var const cl_R& b = imagpart(x);
  42                         var const cl_R& c = y;
  43                         return complex_C(a+c,b);
  44                 } else {
  45                         DeclareType(cl_C,y);
  46                         // x=a+bi, y=c+di
  47                         var const cl_R& a = realpart(x);
  48                         var const cl_R& b = imagpart(x);
  49                         var const cl_R& c = realpart(y);
  50                         var const cl_R& d = imagpart(y);
  51                         return complex(a+c,b+d);
  52                 }
  53         }
  54 }