cumo 0.1.0 → 0.1.1

data/ext/cumo/include/cumo/template.h

@@ -1,146 +1,151 @@
 #ifndef CUMO_TEMPLATE_H
 #define CUMO_TEMPLATE_H
 
-#define INIT_COUNTER( lp, c )                   \
+#define CUMO_INIT_COUNTER( lp, c )                   \
     {   c = (lp)->n[0]; }
 
-#define NDL_CNT(lp) ((lp)->n[0])
-#define NDL_PTR(lp,i) ((lp)->args[i].ptr + (lp)->args[i].iter[0].pos)
-#define NDL_STEP(lp,i) ((lp)->args[i].iter[0].step)
-#define NDL_IDX(lp,i) ((lp)->args[i].iter[0].idx)
-#define NDL_ESZ(lp,i) ((lp)->args[i].elmsz)
-#define NDL_SHAPE(lp,i) ((lp)->args[i].shape)
+#define CUMO_NDL_CNT(lp) ((lp)->n[0])
+#define CUMO_NDL_ARG(lp,i) ((lp)->args[i])
+#define CUMO_NDL_PTR(lp,i) ((lp)->args[i].ptr + (lp)->args[i].iter[0].pos)
+#define CUMO_NDL_STEP(lp,i) ((lp)->args[i].iter[0].step)
+#define CUMO_NDL_IDX(lp,i) ((lp)->args[i].iter[0].idx)
+#define CUMO_NDL_ESZ(lp,i) ((lp)->args[i].elmsz)
+#define CUMO_NDL_SHAPE(lp,i) ((lp)->args[i].shape)
 
-#define INIT_PTR( lp, i, pt, st )                               \
+#define CUMO_NDL_ARG_STEP(arg,idim) ((arg).iter[idim].step)
+#define CUMO_NDL_ARG_IDX(arg,idim) ((arg).iter[idim].idx)
+#define CUMO_NDL_ARG_SHAPE(arg,idim) ((arg).shape[idim])
+
+#define CUMO_INIT_PTR( lp, i, pt, st )                               \
     {                                                           \
         pt = ((lp)->args[i]).ptr + ((lp)->args[i].iter[0]).pos;         \
         st = ((lp)->args[i].iter[0]).step;                              \
     }
 
-#define INIT_PTR_IDX( lp, i, pt, st, id )                       \
+#define CUMO_INIT_PTR_IDX( lp, i, pt, st, id )                       \
     {                                                           \
         pt = ((lp)->args[i]).ptr + ((lp)->args[i].iter[0]).pos;         \
         st = ((lp)->args[i].iter[0]).step;                              \
         id = ((lp)->args[i].iter[0]).idx;                               \
     }
 
-#define INIT_ELMSIZE( lp, i, es )                               \
+#define CUMO_INIT_ELMSIZE( lp, i, es )                               \
     {                                                           \
         es = ((lp)->args[i]).elmsz;                             \
     }
 
-#define INIT_PTR_BIT( lp, i, ad, ps, st )               \
+#define CUMO_INIT_PTR_BIT( lp, i, ad, ps, st )               \
     {                                                   \
         ps = ((lp)->args[i].iter[0]).pos;                       \
-        ad = (BIT_DIGIT*)(((lp)->args[i]).ptr) + ps/NB; \
-        ps %= NB;                                       \
+        ad = (CUMO_BIT_DIGIT*)(((lp)->args[i]).ptr) + ps/CUMO_NB; \
+        ps %= CUMO_NB;                                       \
         st = ((lp)->args[i].iter[0]).step;                      \
     }
 
-#define INIT_PTR_BIT_IDX( lp, i, ad, ps, st, id )       \
+#define CUMO_INIT_PTR_BIT_IDX( lp, i, ad, ps, st, id )       \
     {                                                   \
         ps = ((lp)->args[i].iter[0]).pos;                       \
-        ad = (BIT_DIGIT*)(((lp)->args[i]).ptr) + ps/NB; \
-        ps %= NB;                                       \
+        ad = (CUMO_BIT_DIGIT*)(((lp)->args[i]).ptr) + ps/CUMO_NB; \
+        ps %= CUMO_NB;                                       \
         st = ((lp)->args[i].iter[0]).step;                      \
         id = ((lp)->args[i].iter[0]).idx;                       \
     }
 
-#define GET_DATA( ptr, type, val )                 \
+#define CUMO_GET_DATA( ptr, type, val )                 \
     {                                              \
         val = *(type*)(ptr);                       \
     }
 
-#define SET_DATA( ptr, type, val )                 \
+#define CUMO_SET_DATA( ptr, type, val )                 \
     {                                              \
         *(type*)(ptr) = val;                       \
     }
 
-#define GET_DATA_STRIDE( ptr, step, type, val )    \
+#define CUMO_GET_DATA_STRIDE( ptr, step, type, val )    \
     {                                              \
         val = *(type*)(ptr);                       \
         ptr += step;                               \
     }
 
-#define GET_DATA_INDEX( ptr, idx, type, val )     \
+#define CUMO_GET_DATA_INDEX( ptr, idx, type, val )     \
     {                                           \
         val = *(type*)(ptr + *idx);             \
         idx++;                                  \
     }
 
-#define SET_DATA_STRIDE( ptr, step, type, val ) \
+#define CUMO_SET_DATA_STRIDE( ptr, step, type, val ) \
     {                                           \
         *(type*)(ptr) = val;                    \
         ptr += step;                            \
     }
 
-#define SET_DATA_INDEX( ptr, idx, type, val )   \
+#define CUMO_SET_DATA_INDEX( ptr, idx, type, val )   \
     {                                           \
         *(type*)(ptr + *idx) = val;             \
         idx++;                                  \
     }
 
-#define LOAD_BIT( adr, pos, val )                       \
+#define CUMO_LOAD_BIT( adr, pos, val )                       \
     {                                                   \
-        size_t dig = (pos) / NB;                        \
-        int    bit = (pos) % NB;                        \
-        val = (((BIT_DIGIT*)(adr))[dig]>>(bit)) & 1u;   \
+        size_t dig = (pos) / CUMO_NB;                        \
+        int    bit = (pos) % CUMO_NB;                        \
+        val = (((CUMO_BIT_DIGIT*)(adr))[dig]>>(bit)) & 1u;   \
     }
 
-#define LOAD_BIT_STEP( adr, pos, step, idx, val )       \
+#define CUMO_LOAD_BIT_STEP( adr, pos, step, idx, val )       \
     {                                                   \
         size_t dig; int bit;                            \
         if (idx) {                                      \
-            dig = ((pos) + *(idx)) / NB;                \
-            bit = ((pos) + *(idx)) % NB;                \
+            dig = ((pos) + *(idx)) / CUMO_NB;                \
+            bit = ((pos) + *(idx)) % CUMO_NB;                \
             idx++;                                      \
         } else {                                        \
-            dig = (pos) / NB;                           \
-            bit = (pos) % NB;                           \
+            dig = (pos) / CUMO_NB;                           \
+            bit = (pos) % CUMO_NB;                           \
             pos += step;                                \
         }                                               \
-        val = (((BIT_DIGIT*)(adr))[dig]>>bit) & 1u;     \
+        val = (((CUMO_BIT_DIGIT*)(adr))[dig]>>bit) & 1u;     \
     }
 
-#define STORE_BIT(adr,pos,val)                  \
+#define CUMO_STORE_BIT(adr,pos,val)                  \
     {                                           \
-        size_t dig = (pos) / NB;                \
-        int    bit = (pos) % NB;                \
-        ((BIT_DIGIT*)(adr))[dig] =              \
-            (((BIT_DIGIT*)(adr))[dig] & ~(1u<<(bit))) | ((val)<<(bit)); \
+        size_t dig = (pos) / CUMO_NB;                \
+        int    bit = (pos) % CUMO_NB;                \
+        ((CUMO_BIT_DIGIT*)(adr))[dig] =              \
+            (((CUMO_BIT_DIGIT*)(adr))[dig] & ~(1u<<(bit))) | ((val)<<(bit)); \
     }
 // val -> val&1 ??
 
-#define STORE_BIT_STEP( adr, pos, step, idx, val )\
+#define CUMO_STORE_BIT_STEP( adr, pos, step, idx, val )\
     {                                           \
         size_t dig; int bit;                    \
         if (idx) {                              \
-            dig = ((pos) + *(idx)) / NB;        \
-            bit = ((pos) + *(idx)) % NB;        \
+            dig = ((pos) + *(idx)) / CUMO_NB;        \
+            bit = ((pos) + *(idx)) % CUMO_NB;        \
             idx++;                              \
         } else {                                \
-            dig = (pos) / NB;                   \
-            bit = (pos) % NB;                   \
+            dig = (pos) / CUMO_NB;                   \
+            bit = (pos) % CUMO_NB;                   \
             pos += step;                        \
         }                                       \
-        ((BIT_DIGIT*)(adr))[dig] =              \
-            (((BIT_DIGIT*)(adr))[dig] & ~(1u<<(bit))) | ((val)<<(bit)); \
+        ((CUMO_BIT_DIGIT*)(adr))[dig] =              \
+            (((CUMO_BIT_DIGIT*)(adr))[dig] & ~(1u<<(bit))) | ((val)<<(bit)); \
     }
 // val -> val&1 ??
 
 static inline int
-is_aligned(const void *ptr, const size_t alignment)
+cumo_is_aligned(const void *ptr, const size_t alignment)
 {
     return ((size_t)(ptr) & ((alignment)-1)) == 0;
 }
 
 static inline int
-is_aligned_step(const ssize_t step, const size_t alignment)
+cumo_is_aligned_step(const ssize_t step, const size_t alignment)
 {
     return ((step) & ((alignment)-1)) == 0;
 }
 
-#define SHOW_WARNING_ONCE( c_str ) \
+#define CUMO_SHOW_WARNING_ONCE( c_str ) \
     { \
         static bool show_warning = true; \
         if (show_warning) { \
@@ -149,10 +154,10 @@ is_aligned_step(const ssize_t step, const size_t alignment)
         } \
     }
 
-#define SHOW_SYNCHRONIZE_FIXME_WARNING_ONCE( func_name, type_name ) \
-    SHOW_WARNING_ONCE("Warning: FIXME: Method \"" func_name "\" for dtype \"" type_name "\" synchronizes with CPU.\n")
+#define CUMO_SHOW_SYNCHRONIZE_FIXME_WARNING_ONCE( func_name, type_name ) \
+    CUMO_SHOW_WARNING_ONCE("Warning: FIXME: Method \"" func_name "\" for dtype \"" type_name "\" synchronizes with CPU.\n")
 
-#define SHOW_SYNCHRONIZE_WARNING_ONCE( func_name, type_name ) \
-    SHOW_WARNING_ONCE("Warning: Method \"" func_name "\" for dtype \"" type_name "\" synchronizes with CPU.\n")
+#define CUMO_SHOW_SYNCHRONIZE_WARNING_ONCE( func_name, type_name ) \
+    CUMO_SHOW_WARNING_ONCE("Warning: Method \"" func_name "\" for dtype \"" type_name "\" synchronizes with CPU.\n")
 
 #endif /* ifndef CUMO_TEMPLATE_H */

data/ext/cumo/include/cumo/template_kernel.h

@@ -1,76 +1,76 @@
 #ifndef CUMO_TEMPLATE_KERNEL_H
 #define CUMO_TEMPLATE_KERNEL_H
 
-#define LOAD_BIT( adr, pos, val )                       \
+#define CUMO_LOAD_BIT( adr, pos, val )                       \
     {                                                   \
-        size_t dig = (size_t)(pos) / NB;                \
-        int    bit = (size_t)(pos) % NB;                \
-        val = (((BIT_DIGIT*)(adr))[dig]>>(bit)) & 1u;   \
+        size_t dig = (size_t)(pos) / CUMO_NB;                \
+        int    bit = (size_t)(pos) % CUMO_NB;                \
+        val = (((CUMO_BIT_DIGIT*)(adr))[dig]>>(bit)) & 1u;   \
     }
 
-#define LOAD_BIT_STEP( adr, pos, step, idx, val )       \
+#define CUMO_LOAD_BIT_STEP( adr, pos, step, idx, val )       \
     {                                                   \
         size_t dig; int bit;                            \
         if (idx) {                                      \
-            dig = (size_t)((pos) + *(idx)) / NB;        \
-            bit = (size_t)((pos) + *(idx)) % NB;        \
+            dig = (size_t)((pos) + *(idx)) / CUMO_NB;        \
+            bit = (size_t)((pos) + *(idx)) % CUMO_NB;        \
             idx++;                                      \
         } else {                                        \
-            dig = (size_t)(pos) / NB;                   \
-            bit = (size_t)(pos) % NB;                   \
+            dig = (size_t)(pos) / CUMO_NB;                   \
+            bit = (size_t)(pos) % CUMO_NB;                   \
             pos += step;                                \
         }                                               \
-        val = (((BIT_DIGIT*)(adr))[dig]>>bit) & 1u;     \
+        val = (((CUMO_BIT_DIGIT*)(adr))[dig]>>bit) & 1u;     \
     }
 
-#define STORE_BIT(adr,pos,val)                                     \
+#define CUMO_STORE_BIT(adr,pos,val)                                     \
     {                                                              \
-        size_t dig = (size_t)(pos) / NB;                           \
-        int    bit = (size_t)(pos) % NB;                           \
+        size_t dig = (size_t)(pos) / CUMO_NB;                           \
+        int    bit = (size_t)(pos) % CUMO_NB;                           \
         if (val) {                                                 \
-            atomicOr((BIT_DIGIT*)(adr) + (dig), (val)<<(bit));     \
+            atomicOr((CUMO_BIT_DIGIT*)(adr) + (dig), (val)<<(bit));     \
         } else {                                                   \
-            atomicAnd((BIT_DIGIT*)(adr) + (dig), ~(1u<<(bit)));    \
+            atomicAnd((CUMO_BIT_DIGIT*)(adr) + (dig), ~(1u<<(bit)));    \
         }                                                          \
     }
 // val -> val&1 ??
 
-#define STORE_BIT_STEP( adr, pos, step, idx, val )                 \
+#define CUMO_STORE_BIT_STEP( adr, pos, step, idx, val )                 \
     {                                                              \
         size_t dig; int bit;                                       \
         if (idx) {                                                 \
-            dig = (size_t)((pos) + *(idx)) / NB;                   \
-            bit = (size_t)((pos) + *(idx)) % NB;                   \
+            dig = (size_t)((pos) + *(idx)) / CUMO_NB;                   \
+            bit = (size_t)((pos) + *(idx)) % CUMO_NB;                   \
             idx++;                                                 \
         } else {                                                   \
-            dig = (size_t)(pos) / NB;                              \
-            bit = (size_t)(pos) % NB;                              \
+            dig = (size_t)(pos) / CUMO_NB;                              \
+            bit = (size_t)(pos) % CUMO_NB;                              \
             pos += step;                                           \
         }                                                          \
         if (val) {                                                 \
-            atomicOr((BIT_DIGIT*)(adr) + (dig), (val)<<(bit));     \
+            atomicOr((CUMO_BIT_DIGIT*)(adr) + (dig), (val)<<(bit));     \
         } else {                                                   \
-            atomicAnd((BIT_DIGIT*)(adr) + (dig), ~((1u)<<(bit)));  \
+            atomicAnd((CUMO_BIT_DIGIT*)(adr) + (dig), ~((1u)<<(bit)));  \
         }                                                          \
     }
 // val -> val&1 ??
 
-#define MAX_BLOCK_DIM 128
-#define MAX_GRID_DIM 2147483647 // ref. http://docs.nvidia.com/cuda/cuda-c-programming-guide/index.html#compute-capabilities
+#define CUMO_MAX_BLOCK_DIM 128
+#define CUMO_MAX_GRID_DIM 2147483647 // ref. http://docs.nvidia.com/cuda/cuda-c-programming-guide/index.html#compute-capabilities
 
 static inline size_t
-get_gridDim(size_t n)
+cumo_get_grid_dim(size_t n)
 {
-    size_t gridDim = (n / MAX_BLOCK_DIM) + 1;
-    if (gridDim > MAX_GRID_DIM) gridDim = MAX_GRID_DIM;
-    return gridDim;
+    size_t grid_dim = (n / CUMO_MAX_BLOCK_DIM) + 1;
+    if (grid_dim > CUMO_MAX_GRID_DIM) grid_dim = CUMO_MAX_GRID_DIM;
+    return grid_dim;
 }
 
 static inline size_t
-get_blockDim(size_t n)
+cumo_get_block_dim(size_t n)
 {
-    size_t blockDim = (n > MAX_BLOCK_DIM) ? MAX_BLOCK_DIM : n;
-    return blockDim;
+    size_t block_dim = (n > CUMO_MAX_BLOCK_DIM) ? CUMO_MAX_BLOCK_DIM : n;
+    return block_dim;
 }
 
 

data/ext/cumo/include/cumo/types/bit.h

@@ -1,5 +1,5 @@
-typedef BIT_DIGIT dtype;
-typedef BIT_DIGIT rtype;
+typedef CUMO_BIT_DIGIT dtype;
+typedef CUMO_BIT_DIGIT rtype;
 #define cT  cumo_cBit
 #define cRT cT
 
@@ -30,7 +30,7 @@ typedef BIT_DIGIT rtype;
 #define m_count_false(x) ((x)==0)
 #define m_count_false_cpu(x) ((x)==0)
 
-static inline BIT_DIGIT m_num_to_data(VALUE num) {
+static inline CUMO_BIT_DIGIT m_num_to_data(VALUE num) {
     if (RTEST(num)) {
         if (!RTEST(rb_equal(num,INT2FIX(0)))) {
             return 1;

data/ext/cumo/include/cumo/types/bit_kernel.h

@@ -1,8 +1,8 @@
 #ifndef CUMO_BIT_KERNEL_H
 #define CUMO_BIT_KERNEL_H
 
-typedef BIT_DIGIT dtype;
-typedef BIT_DIGIT rtype;
+typedef CUMO_BIT_DIGIT dtype;
+typedef CUMO_BIT_DIGIT rtype;
 
 #define m_zero 0
 #define m_one  1

data/ext/cumo/include/cumo/types/complex.h

@@ -1,165 +1,165 @@
 static inline dtype c_new(rtype r, rtype i) {
     dtype z;
-    REAL(z) = r;
-    IMAG(z) = i;
+    CUMO_REAL(z) = r;
+    CUMO_IMAG(z) = i;
     return z;
 }
 
 static inline dtype c_set_real(dtype x, rtype r) {
-    REAL(x)=r;
+    CUMO_REAL(x)=r;
     return x;
 }
 
 static inline dtype c_set_imag(dtype x, rtype i) {
-    IMAG(x)=i;
+    CUMO_IMAG(x)=i;
     return x;
 }
 
 static inline VALUE COMP2NUM(dtype x) {
     VALUE v;
     v = rb_funcall(rb_intern("Kernel"), rb_intern("Complex"), 2,
-                   rb_float_new(REAL(x)), rb_float_new(IMAG(x)));
+                   rb_float_new(CUMO_REAL(x)), rb_float_new(CUMO_IMAG(x)));
     return v;
 }
 
 static inline dtype NUM2COMP(VALUE v) {
     dtype z;
-    REAL(z) = NUM2DBL(rb_funcall(v,id_real,0));
-    IMAG(z) = NUM2DBL(rb_funcall(v,id_imag,0));
+    CUMO_REAL(z) = NUM2DBL(rb_funcall(v,cumo_id_real,0));
+    CUMO_IMAG(z) = NUM2DBL(rb_funcall(v,cumo_id_imag,0));
     return z;
 }
 
-#define c_is_zero(x) (REAL(x)==0 && IMAG(x)==0)
-#define c_eq(x,y) (REAL(x)==REAL(y) && IMAG(x)==IMAG(y))
-#define c_ne(x,y) (REAL(x)!=REAL(y) || IMAG(x)!=IMAG(y))
-#define c_isnan(x) (isnan(REAL(x)) || isnan(IMAG(x)))
-#define c_isinf(x) (isinf(REAL(x)) || isinf(IMAG(x)))
-#define c_isposinf(x) ((isinf(REAL(x)) && signbit(REAL(x))==0) || \
-                       (isinf(IMAG(x)) && signbit(IMAG(x))==0))
-#define c_isneginf(x) ((isinf(REAL(x)) && signbit(REAL(x))) || \
-                       (isinf(IMAG(x)) && signbit(IMAG(x))))
-#define c_isfinite(x) (isfinite(REAL(x)) && isfinite(IMAG(x)))
+#define c_is_zero(x) (CUMO_REAL(x)==0 && CUMO_IMAG(x)==0)
+#define c_eq(x,y) (CUMO_REAL(x)==CUMO_REAL(y) && CUMO_IMAG(x)==CUMO_IMAG(y))
+#define c_ne(x,y) (CUMO_REAL(x)!=CUMO_REAL(y) || CUMO_IMAG(x)!=CUMO_IMAG(y))
+#define c_isnan(x) (isnan(CUMO_REAL(x)) || isnan(CUMO_IMAG(x)))
+#define c_isinf(x) (isinf(CUMO_REAL(x)) || isinf(CUMO_IMAG(x)))
+#define c_isposinf(x) ((isinf(CUMO_REAL(x)) && signbit(CUMO_REAL(x))==0) || \
+                       (isinf(CUMO_IMAG(x)) && signbit(CUMO_IMAG(x))==0))
+#define c_isneginf(x) ((isinf(CUMO_REAL(x)) && signbit(CUMO_REAL(x))) || \
+                       (isinf(CUMO_IMAG(x)) && signbit(CUMO_IMAG(x))))
+#define c_isfinite(x) (isfinite(CUMO_REAL(x)) && isfinite(CUMO_IMAG(x)))
 
 static inline dtype c_zero() {
     dtype z;
-    REAL(z) = 0;
-    IMAG(z) = 0;
+    CUMO_REAL(z) = 0;
+    CUMO_IMAG(z) = 0;
     return z;
 }
 
 static inline dtype c_one() {
     dtype z;
-    REAL(z) = 1;
-    IMAG(z) = 0;
+    CUMO_REAL(z) = 1;
+    CUMO_IMAG(z) = 0;
     return z;
 }
 
 static inline dtype c_minus(dtype x) {
     dtype z;
-    REAL(z) = -REAL(x);
-    IMAG(z) = -IMAG(x);
+    CUMO_REAL(z) = -CUMO_REAL(x);
+    CUMO_IMAG(z) = -CUMO_IMAG(x);
     return z;
 }
 
 static inline dtype c_im(dtype x) {
     dtype z;
-    REAL(z) = -IMAG(x);
-    IMAG(z) = REAL(x);
+    CUMO_REAL(z) = -CUMO_IMAG(x);
+    CUMO_IMAG(z) = CUMO_REAL(x);
     return z;
 }
 
 static inline dtype c_add(dtype x, dtype y) {
     dtype z;
-    REAL(z) = REAL(x)+REAL(y);
-    IMAG(z) = IMAG(x)+IMAG(y);
+    CUMO_REAL(z) = CUMO_REAL(x)+CUMO_REAL(y);
+    CUMO_IMAG(z) = CUMO_IMAG(x)+CUMO_IMAG(y);
     return z;
 }
 
 static inline dtype c_sub(dtype x, dtype y) {
     dtype z;
-    REAL(z) = REAL(x)-REAL(y);
-    IMAG(z) = IMAG(x)-IMAG(y);
+    CUMO_REAL(z) = CUMO_REAL(x)-CUMO_REAL(y);
+    CUMO_IMAG(z) = CUMO_IMAG(x)-CUMO_IMAG(y);
     return z;
 }
 
 
 static inline dtype c_mul(dtype x, dtype y) {
     dtype z;
-    REAL(z) = REAL(x)*REAL(y)-IMAG(x)*IMAG(y);
-    IMAG(z) = REAL(x)*IMAG(y)+IMAG(x)*REAL(y);
+    CUMO_REAL(z) = CUMO_REAL(x)*CUMO_REAL(y)-CUMO_IMAG(x)*CUMO_IMAG(y);
+    CUMO_IMAG(z) = CUMO_REAL(x)*CUMO_IMAG(y)+CUMO_IMAG(x)*CUMO_REAL(y);
     return z;
 }
 
 static inline dtype c_mul_r(dtype x, rtype y) {
     dtype z;
-    REAL(z) = REAL(x)*y;
-    IMAG(z) = IMAG(x)*y;
+    CUMO_REAL(z) = CUMO_REAL(x)*y;
+    CUMO_IMAG(z) = CUMO_IMAG(x)*y;
     return z;
 }
 
 static inline dtype c_div(dtype x, dtype y) {
     dtype z;
     rtype s,yr,yi;
-    s  = r_hypot(REAL(y),IMAG(y));
-    yr = REAL(y)/s;
-    yi = IMAG(y)/s;
-    REAL(z) = (REAL(x)*yr+IMAG(x)*yi)/s;
-    IMAG(z) = (IMAG(x)*yr-REAL(x)*yi)/s;
+    s  = r_hypot(CUMO_REAL(y),CUMO_IMAG(y));
+    yr = CUMO_REAL(y)/s;
+    yi = CUMO_IMAG(y)/s;
+    CUMO_REAL(z) = (CUMO_REAL(x)*yr+CUMO_IMAG(x)*yi)/s;
+    CUMO_IMAG(z) = (CUMO_IMAG(x)*yr-CUMO_REAL(x)*yi)/s;
     return z;
 }
 
 static inline dtype c_div_r(dtype x, rtype y) {
     dtype z;
-    REAL(z) = REAL(x)/y;
-    IMAG(z) = IMAG(x)/y;
+    CUMO_REAL(z) = CUMO_REAL(x)/y;
+    CUMO_IMAG(z) = CUMO_IMAG(x)/y;
     return z;
 }
 
 static inline dtype c_reciprocal(dtype x) {
     dtype z;
-    if ( r_abs(REAL(x)) > r_abs(IMAG(x)) ) {
-        IMAG(z) = IMAG(x)/REAL(x);
-        REAL(z) = (1+IMAG(z)*IMAG(z))*REAL(x);
-        IMAG(z) /= -REAL(z);
-        REAL(z) = 1/REAL(z);
+    if ( r_abs(CUMO_REAL(x)) > r_abs(CUMO_IMAG(x)) ) {
+        CUMO_IMAG(z) = CUMO_IMAG(x)/CUMO_REAL(x);
+        CUMO_REAL(z) = (1+CUMO_IMAG(z)*CUMO_IMAG(z))*CUMO_REAL(x);
+        CUMO_IMAG(z) /= -CUMO_REAL(z);
+        CUMO_REAL(z) = 1/CUMO_REAL(z);
     } else {
-        REAL(z) = REAL(x)/IMAG(x);
-        IMAG(z) = (1+REAL(z)*REAL(z))*IMAG(x);
-        REAL(z) /= IMAG(z);
-        IMAG(z) = -1/IMAG(z);
+        CUMO_REAL(z) = CUMO_REAL(x)/CUMO_IMAG(x);
+        CUMO_IMAG(z) = (1+CUMO_REAL(z)*CUMO_REAL(z))*CUMO_IMAG(x);
+        CUMO_REAL(z) /= CUMO_IMAG(z);
+        CUMO_IMAG(z) = -1/CUMO_IMAG(z);
     }
     return z;
 }
 
 static inline dtype c_square(dtype x) {
     dtype z;
-    REAL(z) = REAL(x)*REAL(x)-IMAG(x)*IMAG(x);
-    IMAG(z) = 2*REAL(x)*IMAG(x);
+    CUMO_REAL(z) = CUMO_REAL(x)*CUMO_REAL(x)-CUMO_IMAG(x)*CUMO_IMAG(x);
+    CUMO_IMAG(z) = 2*CUMO_REAL(x)*CUMO_IMAG(x);
     return z;
 }
 
 static inline dtype c_sqrt(dtype x) {
     dtype z;
     rtype xr, xi, r;
-    xr = REAL(x)/2;
-    xi = IMAG(x)/2;
+    xr = CUMO_REAL(x)/2;
+    xi = CUMO_IMAG(x)/2;
     r  = r_hypot(xr,xi);
     if (xr>0) {
-        REAL(z) = sqrt(r+xr);
-        IMAG(z) = xi/REAL(z);
+        CUMO_REAL(z) = sqrt(r+xr);
+        CUMO_IMAG(z) = xi/CUMO_REAL(z);
     } else if ( (r-=xr)!=0 ) {
-        IMAG(z) = (xi>=0) ? sqrt(r):-sqrt(r);
-        REAL(z) = xi/IMAG(z);
+        CUMO_IMAG(z) = (xi>=0) ? sqrt(r):-sqrt(r);
+        CUMO_REAL(z) = xi/CUMO_IMAG(z);
     } else {
-        REAL(z) = IMAG(z) = 0;
+        CUMO_REAL(z) = CUMO_IMAG(z) = 0;
     }
     return z;
 }
 
 static inline dtype c_log(dtype x) {
     dtype z;
-    REAL(z) = r_log(r_hypot(REAL(x),IMAG(x)));
-    IMAG(z) = r_atan2(IMAG(x),REAL(x));
+    CUMO_REAL(z) = r_log(r_hypot(CUMO_REAL(x),CUMO_IMAG(x)));
+    CUMO_IMAG(z) = r_atan2(CUMO_IMAG(x),CUMO_REAL(x));
     return z;
 }
 
@@ -179,73 +179,73 @@ static inline dtype c_log10(dtype x) {
 
 static inline dtype c_exp(dtype x) {
     dtype z;
-    rtype a = r_exp(REAL(x));
-    REAL(z) = a*r_cos(IMAG(x));
-    IMAG(z) = a*r_sin(IMAG(x));
+    rtype a = r_exp(CUMO_REAL(x));
+    CUMO_REAL(z) = a*r_cos(CUMO_IMAG(x));
+    CUMO_IMAG(z) = a*r_sin(CUMO_IMAG(x));
     return z;
 }
 
 static inline dtype c_exp2(dtype x) {
     dtype z;
-    rtype a = r_exp(REAL(x)*M_LN2);
-    REAL(z) = a*r_cos(IMAG(x));
-    IMAG(z) = a*r_sin(IMAG(x));
+    rtype a = r_exp(CUMO_REAL(x)*M_LN2);
+    CUMO_REAL(z) = a*r_cos(CUMO_IMAG(x));
+    CUMO_IMAG(z) = a*r_sin(CUMO_IMAG(x));
     return z;
 }
 
 static inline dtype c_exp10(dtype x) {
     dtype z;
-    rtype a = r_exp(REAL(x)*M_LN10);
-    REAL(z) = a*r_cos(IMAG(x));
-    IMAG(z) = a*r_sin(IMAG(x));
+    rtype a = r_exp(CUMO_REAL(x)*M_LN10);
+    CUMO_REAL(z) = a*r_cos(CUMO_IMAG(x));
+    CUMO_IMAG(z) = a*r_sin(CUMO_IMAG(x));
     return z;
 }
 
 static inline dtype c_sin(dtype x) {
     dtype z;
-    REAL(z) = r_sin(REAL(x))*r_cosh(IMAG(x));
-    IMAG(z) = r_cos(REAL(x))*r_sinh(IMAG(x));
+    CUMO_REAL(z) = r_sin(CUMO_REAL(x))*r_cosh(CUMO_IMAG(x));
+    CUMO_IMAG(z) = r_cos(CUMO_REAL(x))*r_sinh(CUMO_IMAG(x));
     return z;
 }
 
 static inline dtype c_sinh(dtype x) {
     dtype z;
-    REAL(z) = r_sinh(REAL(x))*r_cos(IMAG(x));
-    IMAG(z) = r_cosh(REAL(x))*r_sin(IMAG(x));
+    CUMO_REAL(z) = r_sinh(CUMO_REAL(x))*r_cos(CUMO_IMAG(x));
+    CUMO_IMAG(z) = r_cosh(CUMO_REAL(x))*r_sin(CUMO_IMAG(x));
     return z;
 }
 
 static inline dtype c_cos(dtype x) {
     dtype z;
-    REAL(z) = r_cos(REAL(x))*r_cosh(IMAG(x));
-    IMAG(z) = -r_sin(REAL(x))*r_sinh(IMAG(x));
+    CUMO_REAL(z) = r_cos(CUMO_REAL(x))*r_cosh(CUMO_IMAG(x));
+    CUMO_IMAG(z) = -r_sin(CUMO_REAL(x))*r_sinh(CUMO_IMAG(x));
     return z;
 }
 
 static inline dtype c_cosh(dtype x) {
     dtype z;
-    REAL(z) = r_cosh(REAL(x))*r_cos(IMAG(x));
-    IMAG(z) = r_sinh(REAL(x))*r_sin(IMAG(x));
+    CUMO_REAL(z) = r_cosh(CUMO_REAL(x))*r_cos(CUMO_IMAG(x));
+    CUMO_IMAG(z) = r_sinh(CUMO_REAL(x))*r_sin(CUMO_IMAG(x));
     return z;
 }
 
 static inline dtype c_tan(dtype x) {
     dtype z;
     rtype c, d;
-    if (r_abs(IMAG(x))<1) {
-        c = r_cos(REAL(x));
-        d = r_sinh(IMAG(x));
+    if (r_abs(CUMO_IMAG(x))<1) {
+        c = r_cos(CUMO_REAL(x));
+        d = r_sinh(CUMO_IMAG(x));
         d = c*c + d*d;
-        REAL(z) = 0.5*r_sin(2*REAL(x))/d;
-        IMAG(z) = 0.5*r_sinh(2*IMAG(x))/d;
+        CUMO_REAL(z) = 0.5*r_sin(2*CUMO_REAL(x))/d;
+        CUMO_IMAG(z) = 0.5*r_sinh(2*CUMO_IMAG(x))/d;
     } else {
-        d = r_exp(-IMAG(x));
+        d = r_exp(-CUMO_IMAG(x));
         c = 2*d/(1-d*d);
         c = c*c;
-        d = r_cos(REAL(x));
+        d = r_cos(CUMO_REAL(x));
         d = 1.0 + d*d*c;
-        REAL(z) = 0.5*r_sin(2*REAL(x))*c/d;
-        IMAG(z) = 1/r_tanh(IMAG(x))/d;
+        CUMO_REAL(z) = 0.5*r_sin(2*CUMO_REAL(x))*c/d;
+        CUMO_IMAG(z) = 1/r_tanh(CUMO_IMAG(x))/d;
     }
     return z;
 }
@@ -253,17 +253,17 @@ static inline dtype c_tan(dtype x) {
 static inline dtype c_tanh(dtype x) {
     dtype z;
     rtype c, d, s;
-    c = r_cos(IMAG(x));
-    s = r_sinh(REAL(x));
+    c = r_cos(CUMO_IMAG(x));
+    s = r_sinh(CUMO_REAL(x));
     d = c*c + s*s;
-    if (r_abs(REAL(x))<1) {
-        REAL(z) = s*r_cosh(REAL(x))/d;
-        IMAG(z) = 0.5*r_sin(2*IMAG(x))/d;
+    if (r_abs(CUMO_REAL(x))<1) {
+        CUMO_REAL(z) = s*r_cosh(CUMO_REAL(x))/d;
+        CUMO_IMAG(z) = 0.5*r_sin(2*CUMO_IMAG(x))/d;
     } else {
         c = c / s;
         c = 1 + c*c;
-        REAL(z) = 1/(r_tanh(REAL(x))*c);
-        IMAG(z) = 0.5*r_sin(2*IMAG(x))/d;
+        CUMO_REAL(z) = 1/(r_tanh(CUMO_REAL(x))*c);
+        CUMO_IMAG(z) = 0.5*r_sin(2*CUMO_IMAG(x))/d;
     }
     return z;
 }
@@ -271,24 +271,24 @@ static inline dtype c_tanh(dtype x) {
 static inline dtype c_asin(dtype x) {
     dtype z, y;
     y = c_square(x);
-    REAL(y) = 1-REAL(y);
-    IMAG(y) = -IMAG(y);
+    CUMO_REAL(y) = 1-CUMO_REAL(y);
+    CUMO_IMAG(y) = -CUMO_IMAG(y);
     y = c_sqrt(y);
-    REAL(y) -= IMAG(x);
-    IMAG(y) += REAL(x);
+    CUMO_REAL(y) -= CUMO_IMAG(x);
+    CUMO_IMAG(y) += CUMO_REAL(x);
     y = c_log(y);
-    REAL(z) = IMAG(y);
-    IMAG(z) = -REAL(y);
+    CUMO_REAL(z) = CUMO_IMAG(y);
+    CUMO_IMAG(z) = -CUMO_REAL(y);
     return z;
 }
 
 static inline dtype c_asinh(dtype x) {
     dtype z, y;
     y = c_square(x);
-    REAL(y) += 1;
+    CUMO_REAL(y) += 1;
     y = c_sqrt(y);
-    REAL(y) += REAL(x);
-    IMAG(y) += IMAG(x);
+    CUMO_REAL(y) += CUMO_REAL(x);
+    CUMO_IMAG(y) += CUMO_IMAG(x);
     z = c_log(y);
     return z;
 }
@@ -296,51 +296,51 @@ static inline dtype c_asinh(dtype x) {
 static inline dtype c_acos(dtype x) {
     dtype z, y;
     y = c_square(x);
-    REAL(y) = 1-REAL(y);
-    IMAG(y) = -IMAG(y);
+    CUMO_REAL(y) = 1-CUMO_REAL(y);
+    CUMO_IMAG(y) = -CUMO_IMAG(y);
     y = c_sqrt(y);
-    REAL(z) = REAL(x)-IMAG(y);
-    IMAG(z) = IMAG(x)+REAL(y);
+    CUMO_REAL(z) = CUMO_REAL(x)-CUMO_IMAG(y);
+    CUMO_IMAG(z) = CUMO_IMAG(x)+CUMO_REAL(y);
     y = c_log(z);
-    REAL(z) = IMAG(y);
-    IMAG(z) = -REAL(y);
+    CUMO_REAL(z) = CUMO_IMAG(y);
+    CUMO_IMAG(z) = -CUMO_REAL(y);
     return z;
 }
 
 static inline dtype c_acosh(dtype x) {
     dtype z, y;
     y = c_square(x);
-    REAL(y) -= 1;
+    CUMO_REAL(y) -= 1;
     y = c_sqrt(y);
-    REAL(y) += REAL(x);
-    IMAG(y) += IMAG(x);
+    CUMO_REAL(y) += CUMO_REAL(x);
+    CUMO_IMAG(y) += CUMO_IMAG(x);
     z = c_log(y);
     return z;
 }
 
 static inline dtype c_atan(dtype x) {
     dtype z, y;
-    REAL(y) = -REAL(x);
-    IMAG(y) = 1-IMAG(x);
-    REAL(z) = REAL(x);
-    IMAG(z) = 1+IMAG(x);
+    CUMO_REAL(y) = -CUMO_REAL(x);
+    CUMO_IMAG(y) = 1-CUMO_IMAG(x);
+    CUMO_REAL(z) = CUMO_REAL(x);
+    CUMO_IMAG(z) = 1+CUMO_IMAG(x);
     y = c_div(z,y);
     y = c_log(y);
-    REAL(z) = -IMAG(y)/2;
-    IMAG(z) = REAL(y)/2;
+    CUMO_REAL(z) = -CUMO_IMAG(y)/2;
+    CUMO_IMAG(z) = CUMO_REAL(y)/2;
     return z;
 }
 
 static inline dtype c_atanh(dtype x) {
     dtype z, y;
-    REAL(y) = 1-REAL(x);
-    IMAG(y) = -IMAG(x);
-    REAL(z) = 1+REAL(x);
-    IMAG(z) = IMAG(x);
+    CUMO_REAL(y) = 1-CUMO_REAL(x);
+    CUMO_IMAG(y) = -CUMO_IMAG(x);
+    CUMO_REAL(z) = 1+CUMO_REAL(x);
+    CUMO_IMAG(z) = CUMO_IMAG(x);
     y = c_div(z,y);
     y = c_log(y);
-    REAL(z) = REAL(y)/2;
-    IMAG(z) = IMAG(y)/2;
+    CUMO_REAL(z) = CUMO_REAL(y)/2;
+    CUMO_IMAG(z) = CUMO_IMAG(y)/2;
     return z;
 }
 
@@ -349,7 +349,7 @@ static inline dtype c_pow(dtype x, dtype y)
     dtype z;
     if (c_is_zero(y)) {
         z = c_one();
-    } else if (c_is_zero(x) && REAL(y)>0 && IMAG(y)==0) {
+    } else if (c_is_zero(x) && CUMO_REAL(y)>0 && CUMO_IMAG(y)==0) {
         z = c_zero();
     } else {
         z = c_log(x);
@@ -386,11 +386,11 @@ static inline dtype c_cbrt(dtype x) {
 }
 
 static inline rtype c_abs(dtype x) {
-    return r_hypot(REAL(x),IMAG(x));
+    return r_hypot(CUMO_REAL(x),CUMO_IMAG(x));
 }
 
 static inline rtype c_abs_square(dtype x) {
-    return REAL(x)*REAL(x)+IMAG(x)*IMAG(x);
+    return CUMO_REAL(x)*CUMO_REAL(x)+CUMO_IMAG(x)*CUMO_IMAG(x);
 }