ruby-dnn 0.1.6 → 0.1.7

data/lib/dnn/ext/cifar10/cifar10_ext.c

@@ -1,52 +0,0 @@
-#include <ruby.h>
-#include <stdlib.h>
-#include "numo/narray.h"
-
-#define CIFAR10_WIDTH 32
-#define CIFAR10_HEIGHT 32
-#define CIFAR10_CHANNEL 3
-#define CIFAR10_CLASSES 10
-
-VALUE cifar10_load(VALUE self, VALUE rb_bin, VALUE rb_num_datas) {
-  char* bin = StringValuePtr(rb_bin);
-  int num_datas = FIX2INT(rb_num_datas);
-  char script[64];
-  VALUE rb_na_x;
-  VALUE rb_na_y;
-  VALUE rb_xy;
-  narray_data_t* na_data_x;
-  narray_data_t* na_data_y;
-  int i;
-  int j = 0;
-  int k = 0;
-  int size = CIFAR10_WIDTH * CIFAR10_HEIGHT * CIFAR10_CHANNEL;
-
-  sprintf(script, "Numo::UInt8.zeros(%d, %d, %d, %d)", num_datas, CIFAR10_WIDTH, CIFAR10_HEIGHT, CIFAR10_CHANNEL);
-  rb_na_x = rb_eval_string(&script[0]);
-  na_data_x = RNARRAY_DATA(rb_na_x);
-  for(i = 0; i < 64; i++) {
-    script[i] = 0;
-  }
-  sprintf(script, "Numo::UInt8.zeros(%d, %d)", num_datas, CIFAR10_CLASSES);
-  rb_na_y = rb_eval_string(&script[0]);
-  na_data_y = RNARRAY_DATA(rb_na_y);
-
-  for (i = 0; i < num_datas; i++) {
-    na_data_y->ptr[i] = bin[j];
-    j++;
-    memcpy(&na_data_x->ptr[k], &bin[j], size);
-    j += size;
-    k += size;
-  }
-
-  rb_xy = rb_ary_new3(2, rb_na_x, rb_na_y);
-  return rb_xy;
-}
-
-void Init_cifar10_ext() {
-  VALUE rb_dnn;
-  VALUE rb_cifar10;
-  rb_dnn = rb_define_module("DNN");
-  rb_cifar10 = rb_define_module_under(rb_dnn, "CIFAR10");
-  rb_define_singleton_method(rb_cifar10, "_cifar10_load", cifar10_load, 2);
-}

data/lib/dnn/ext/cifar10/extconf.rb

@@ -1,3 +0,0 @@
-require "mkmf"
-
-create_makefile("cifar10_ext")

data/lib/dnn/ext/mnist/extconf.rb

@@ -1,3 +0,0 @@
-require "mkmf"
-
-create_makefile("mnist_ext")

data/lib/dnn/ext/mnist/mnist_ext.c

@@ -1,49 +0,0 @@
-#include <ruby.h>
-#include <stdlib.h>
-#include "numo/narray.h"
-
-VALUE mnist_load_images(VALUE self, VALUE rb_bin, VALUE rb_num_images, VALUE rb_cols, VALUE rb_rows) {
-  char* bin = StringValuePtr(rb_bin);
-  int num_images = FIX2INT(rb_num_images);
-  int cols = FIX2INT(rb_cols);
-  int rows = FIX2INT(rb_rows);
-  int i;
-  int j;
-  char script[64];
-  VALUE rb_na;
-  narray_data_t* na_data;
-
-  sprintf(script, "Numo::UInt8.zeros(%d, %d, %d)", num_images, cols, rows);
-  rb_na = rb_eval_string(&script[0]);
-  na_data = RNARRAY_DATA(rb_na);
-
-  for (i = 0; i < num_images; i++) {
-    j = i * cols * rows;
-    memcpy(&na_data->ptr[j], &bin[j], cols * rows);
-  }
-  return rb_na;
-}
-
-VALUE mnist_load_labels(VALUE self, VALUE rb_bin, VALUE rb_num_labels) {
-  char* bin = StringValuePtr(rb_bin);
-  int num_labels = FIX2INT(rb_num_labels);
-  char script[64];
-  VALUE rb_na;
-  narray_data_t* na_data;
-
-  sprintf(script, "Numo::UInt8.zeros(%d)", num_labels);
-  rb_na = rb_eval_string(&script[0]);
-  na_data = RNARRAY_DATA(rb_na);
-
-  memcpy(na_data->ptr, bin, num_labels);
-  return rb_na;
-}
-
-void Init_mnist_ext() {
-  VALUE rb_dnn;
-  VALUE rb_mnist;
-  rb_dnn = rb_define_module("DNN");
-  rb_mnist = rb_define_module_under(rb_dnn, "MNIST");
-  rb_define_singleton_method(rb_mnist, "_mnist_load_images", mnist_load_images, 4);
-  rb_define_singleton_method(rb_mnist, "_mnist_load_labels", mnist_load_labels, 2);
-}

data/lib/dnn/ext/mnist/numo/compat.h

@@ -1,23 +0,0 @@
-#ifndef COMPAT_H
-#define COMPAT_H
-
-#if !defined RSTRING_LEN
-#define RSTRING_LEN(a) RSTRING(a)->len
-#endif
-#if !defined RSTRING_PTR
-#define RSTRING_PTR(a) RSTRING(a)->ptr
-#endif
-#if !defined RARRAY_LEN
-#define RARRAY_LEN(a) RARRAY(a)->len
-#endif
-#if !defined RARRAY_PTR
-#define RARRAY_PTR(a) RARRAY(a)->ptr
-#endif
-#if !defined RARRAY_AREF
-#define RARRAY_AREF(a,i) RARRAY_PTR(a)[i]
-#endif
-#if !defined RARRAY_ASET
-#define RARRAY_ASET(a,i,v) (RARRAY_PTR(a)[i] = v)
-#endif
-
-#endif /* ifndef COMPAT_H */

data/lib/dnn/ext/mnist/numo/extconf.h

@@ -1,13 +0,0 @@
-#ifndef NUMO_EXTCONF_H
-#define NUMO_EXTCONF_H
-#define HAVE_STDBOOL_H 1
-#define HAVE_STDINT_H 1
-#define HAVE_TYPE_BOOL 1
-#define HAVE_TYPE_UINT8_T 1
-#define HAVE_TYPE_UINT16_T 1
-#define HAVE_TYPE_INT32_T 1
-#define HAVE_TYPE_UINT32_T 1
-#define HAVE_TYPE_INT64_T 1
-#define HAVE_TYPE_UINT64_T 1
-#define HAVE_RB_CCOMPLEX 1
-#endif

data/lib/dnn/ext/mnist/numo/intern.h

@@ -1,117 +0,0 @@
-/*
-  intern.h
-  Numerical Array Extension for Ruby
-    (C) Copyright 1999-2017 by Masahiro TANAKA
-*/
-#ifndef INTERN_H
-#define INTERN_H
-
-#define rb_narray_new nary_new
-VALUE nary_new(VALUE elem, int ndim, size_t *shape);
-#define rb_narray_view_new nary_view_new
-VALUE nary_view_new(VALUE elem, int ndim, size_t *shape);
-#define rb_narray_debug_info nary_debug_info
-VALUE nary_debug_info(VALUE);
-
-#define na_make_view nary_make_view
-VALUE nary_make_view(VALUE self);
-
-#define na_s_allocate nary_s_allocate
-VALUE nary_s_allocate(VALUE klass);
-#define na_s_allocate_view nary_s_allocate_view
-VALUE nary_s_allocate_view(VALUE klass);
-#define na_s_new_like nary_s_new_like
-VALUE nary_s_new_like(VALUE type, VALUE obj);
-
-void na_alloc_shape(narray_t *na, int ndim);
-void na_array_to_internal_shape(VALUE self, VALUE ary, size_t *shape);
-void na_index_arg_to_internal_order(int argc, VALUE *argv, VALUE self);
-void na_setup_shape(narray_t *na, int ndim, size_t *shape);
-
-#define na_get_elmsz nary_element_stride
-//#define na_element_stride nary_element_stride
-unsigned int nary_element_stride(VALUE nary);
-#define na_dtype_elmsz nary_dtype_element_stride
-size_t nary_dtype_element_stride(VALUE klass);
-
-#define na_get_pointer nary_get_pointer
-char *nary_get_pointer(VALUE);
-#define na_get_pointer_for_write nary_get_pointer_for_write
-char *nary_get_pointer_for_write(VALUE);
-#define na_get_pointer_for_read nary_get_pointer_for_read
-char *nary_get_pointer_for_read(VALUE);
-#define na_get_pointer_for_read_write nary_get_pointer_for_read_write
-char *nary_get_pointer_for_read_write(VALUE);
-#define na_get_offset nary_get_offset
-size_t nary_get_offset(VALUE self);
-
-#define na_copy_flags nary_copy_flags
-void nary_copy_flags(VALUE src, VALUE dst);
-
-#define na_check_ladder nary_check_ladder
-VALUE nary_check_ladder(VALUE self, int start_dim);
-#define na_check_contiguous nary_check_contiguous
-VALUE nary_check_contiguous(VALUE self);
-
-#define na_flatten_dim nary_flatten_dim
-VALUE nary_flatten_dim(VALUE self, int sd);
-
-#define na_flatten nary_flatten
-VALUE nary_flatten(VALUE);
-
-#define na_copy nary_dup
-VALUE nary_dup(VALUE);
-
-#define na_store nary_store
-VALUE nary_store(VALUE self, VALUE src);
-
-#define na_upcast numo_na_upcast
-VALUE numo_na_upcast(VALUE type1, VALUE type2);
-
-void na_release_lock(VALUE); // currently do nothing
-
-// used in reduce methods
-#define na_reduce_dimension nary_reduce_dimension
-VALUE nary_reduce_dimension(int argc, VALUE *argv, int naryc, VALUE *naryv,
-                            ndfunc_t *ndf, na_iter_func_t nan_iter);
-
-#define na_reduce_options nary_reduce_options
-VALUE nary_reduce_options(VALUE axes, VALUE *opts, int naryc, VALUE *naryv,
-                          ndfunc_t *ndf);
-
-// ndloop
-VALUE na_ndloop(ndfunc_t *nf, int argc, ...);
-VALUE na_ndloop2(ndfunc_t *nf, VALUE args);
-VALUE na_ndloop3(ndfunc_t *nf, void *ptr, int argc, ...);
-VALUE na_ndloop4(ndfunc_t *nf, void *ptr, VALUE args);
-
-VALUE na_ndloop_cast_narray_to_rarray(ndfunc_t *nf, VALUE nary, VALUE fmt);
-VALUE na_ndloop_store_rarray(ndfunc_t *nf, VALUE nary, VALUE rary);
-VALUE na_ndloop_store_rarray2(ndfunc_t *nf, VALUE nary, VALUE rary, VALUE opt);
-VALUE na_ndloop_inspect(VALUE nary, na_text_func_t func, VALUE opt);
-VALUE na_ndloop_with_index(ndfunc_t *nf, int argc, ...);
-
-#define na_info_str nary_info_str
-VALUE nary_info_str(VALUE);
-
-#define na_test_reduce nary_test_reduce
-bool nary_test_reduce(VALUE reduce, int dim);
-
-void nary_step_array_index(VALUE self, size_t ary_size, size_t *plen, ssize_t *pbeg, ssize_t *pstep);
-void nary_step_sequence(VALUE self, size_t *plen, double *pbeg, double *pstep);
-
-// used in aref, aset
-#define na_get_result_dimension nary_get_result_dimension
-int nary_get_result_dimension(VALUE self, int argc, VALUE *argv, ssize_t stride, size_t *pos_idx);
-#define na_aref_main nary_aref_main
-VALUE nary_aref_main(int nidx, VALUE *idx, VALUE self, int keep_dim, int nd);
-
-#include "ruby/version.h"
-
-#if RUBY_API_VERSION_CODE == 20100 // 2.1.0
-int rb_get_kwargs(VALUE keyword_hash, const ID *table, int required, int optional, VALUE *);
-VALUE rb_extract_keywords(VALUE *orighash);
-#endif
-
-
-#endif /* ifndef INTERN_H */

data/lib/dnn/ext/mnist/numo/narray.h

@@ -1,430 +0,0 @@
-/*
-  narray.h
-  Numerical Array Extension for Ruby
-    (C) Copyright 1999-2017 by Masahiro TANAKA
-*/
-#ifndef NARRAY_H
-#define NARRAY_H
-
-#if defined(__cplusplus)
-extern "C" {
-#if 0
-} /* satisfy cc-mode */
-#endif
-#endif
-
-#define NARRAY_VERSION "0.9.1.1"
-#define NARRAY_VERSION_CODE 911
-
-#include <math.h>
-#include "numo/compat.h"
-#include "numo/template.h"
-#include "numo/extconf.h"
-
-#ifdef HAVE_STDBOOL_H
-# include <stdbool.h>
-#endif
-
-#ifdef HAVE_STDINT_H
-# include <stdint.h>
-#endif
-
-#ifdef HAVE_SYS_TYPES_H
-# include <sys/types.h>
-#endif
-
-#ifndef HAVE_U_INT8_T
-# ifdef HAVE_UINT8_T
-    typedef uint8_t u_int8_t;
-# endif
-#endif
-
-#ifndef HAVE_U_INT16_T
-# ifdef HAVE_UINT16_T
-    typedef uint16_t u_int16_t;
-# endif
-#endif
-
-#ifndef HAVE_U_INT32_T
-# ifdef HAVE_UINT32_T
-    typedef uint32_t u_int32_t;
-# endif
-#endif
-
-#ifndef HAVE_U_INT64_T
-# ifdef HAVE_UINT64_T
-    typedef uint64_t u_int64_t;
-# endif
-#endif
-
-#define SZF PRI_SIZE_PREFIX // defined in ruby.h
-
-#if   SIZEOF_LONG==8
-# define NUM2INT64(x) NUM2LONG(x)
-# define INT642NUM(x) LONG2NUM(x)
-# define NUM2UINT64(x) NUM2ULONG(x)
-# define UINT642NUM(x) ULONG2NUM(x)
-# ifndef PRId64
-#  define PRId64 "ld"
-# endif
-# ifndef PRIu64
-#  define PRIu64 "lu"
-# endif
-#elif SIZEOF_LONG_LONG==8
-# define NUM2INT64(x) NUM2LL(x)
-# define INT642NUM(x) LL2NUM(x)
-# define NUM2UINT64(x) NUM2ULL(x)
-# define UINT642NUM(x) ULL2NUM(x)
-# ifndef PRId64
-#  define PRId64 "lld"
-# endif
-# ifndef PRIu64
-#  define PRIu64 "llu"
-# endif
-#endif
-
-#if   SIZEOF_LONG==4
-# define NUM2INT32(x) NUM2LONG(x)
-# define INT322NUM(x) LONG2NUM(x)
-# define NUM2UINT32(x) NUM2ULONG(x)
-# define UINT322NUM(x) ULONG2NUM(x)
-# ifndef PRId32
-#  define PRId32 "ld"
-# endif
-# ifndef PRIu32
-#  define PRIu32 "lu"
-# endif
-#elif SIZEOF_INT==4
-# define NUM2INT32(x) NUM2INT(x)
-# define INT322NUM(x) INT2NUM(x)
-# define NUM2UINT32(x) NUM2UINT(x)
-# define UINT322NUM(x) UINT2NUM(x)
-# ifndef PRId32
-#  define PRId32 "d"
-# endif
-# ifndef PRIu32
-#  define PRIu32 "u"
-# endif
-#endif
-
-#ifndef HAVE_TYPE_BOOL
-  typedef int bool;
-#endif
-#ifndef FALSE                   /* in case these macros already exist */
-# define FALSE   0              /* values of bool */
-#endif
-#ifndef TRUE
-# define TRUE    1
-#endif
-
-typedef struct { float dat[2]; }  scomplex;
-typedef struct { double dat[2]; } dcomplex;
-typedef int fortran_integer;
-
-#define REAL(x) ((x).dat[0])
-#define IMAG(x) ((x).dat[1])
-
-extern int na_debug_flag;
-
-#ifndef NARRAY_C
-extern VALUE numo_cNArray;
-extern VALUE rb_mNumo;
-extern VALUE nary_eCastError;
-extern VALUE nary_eShapeError;
-extern VALUE nary_eOperationError;
-extern VALUE nary_eDimensionError;
-extern VALUE nary_eValueError;
-extern const rb_data_type_t na_data_type;
-
-//EXTERN const int na_sizeof[NA_NTYPES+1];
-#endif
-
-#define cNArray numo_cNArray
-#define mNumo rb_mNumo
-//#define na_upcast(x,y) numo_na_upcast(x,y)
-
-/* global variables within this module */
-extern VALUE numo_cBit;
-extern VALUE numo_cDFloat;
-extern VALUE numo_cSFloat;
-extern VALUE numo_cDComplex;
-extern VALUE numo_cSComplex;
-extern VALUE numo_cInt64;
-extern VALUE numo_cInt32;
-extern VALUE numo_cInt16;
-extern VALUE numo_cInt8;
-extern VALUE numo_cUInt64;
-extern VALUE numo_cUInt32;
-extern VALUE numo_cUInt16;
-extern VALUE numo_cUInt8;
-extern VALUE numo_cRObject;
-extern VALUE na_cStep;
-#ifndef HAVE_RB_CCOMPLEX
-extern VALUE rb_cComplex;
-#endif
-
-extern VALUE sym_reduce;
-extern VALUE sym_option;
-extern VALUE sym_loop_opt;
-extern VALUE sym_init;
-
-#define NARRAY_DATA_T     0x1
-#define NARRAY_VIEW_T     0x2
-#define NARRAY_FILEMAP_T  0x3
-
-typedef struct RNArray {
-    unsigned char ndim;     // # of dimensions
-    unsigned char type;
-    unsigned char flag[2];  // flags
-    unsigned short elmsz;    // element size
-    size_t   size;          // # of total elements
-    size_t  *shape;         // # of elements for each dimension
-    VALUE    reduce;
-} narray_t;
-
-
-typedef struct RNArrayData {
-    narray_t base;
-    char    *ptr;
-} narray_data_t;
-
-
-typedef union {
-    ssize_t stride;
-    size_t *index;
-} stridx_t;
-
-typedef struct RNArrayView {
-    narray_t base;
-    VALUE    data;       // data object
-    size_t   offset;     // offset of start point from data pointer
-                         // :in units of elm.unit_bits
-                         // address_unit  pointer_unit access_unit data_unit
-                         // elm.step_unit = elm.bit_size / elm.access_unit
-                         // elm.step_unit = elm.size_bits / elm.unit_bits
-    stridx_t *stridx;    // stride or indices of data pointer for each dimension
-} narray_view_t;
-
-
-// filemap is unimplemented
-typedef struct RNArrayFileMap {
-    narray_t base;
-    char    *ptr;
-#ifdef WIN32
-    HANDLE hFile;
-    HANDLE hMap;
-#else // POSIX mmap
-    int prot;
-    int flag;
-#endif
-} narray_filemap_t;
-
-
-// this will be revised in future.
-typedef struct {
-    unsigned int element_bits;
-    unsigned int element_bytes;
-    unsigned int element_stride;
-} narray_type_info_t;
-
-
-static inline narray_t *
-na_get_narray_t(VALUE obj)
-{
-    narray_t *na;
-
-    Check_TypedStruct(obj,&na_data_type);
-    na = (narray_t*)DATA_PTR(obj);
-    return na;
-}
-
-static inline narray_t *
-_na_get_narray_t(VALUE obj, unsigned char na_type)
-{
-    narray_t *na;
-
-    Check_TypedStruct(obj,&na_data_type);
-    na = (narray_t*)DATA_PTR(obj);
-    if (na->type != na_type) {
-        rb_bug("unknown type 0x%x (0x%x given)", na_type, na->type);
-    }
-    return na;
-}
-
-#define na_get_narray_data_t(obj) (narray_data_t*)_na_get_narray_t(obj,NARRAY_DATA_T)
-#define na_get_narray_view_t(obj) (narray_view_t*)_na_get_narray_t(obj,NARRAY_VIEW_T)
-#define na_get_narray_filemap_t(obj) (narray_filemap_t*)_na_get_narray_t(obj,NARRAY_FILEMAP_T)
-
-#define GetNArray(obj,var)      TypedData_Get_Struct(obj, narray_t, &na_data_type, var)
-#define GetNArrayView(obj,var)  TypedData_Get_Struct(obj, narray_view_t, &na_data_type, var)
-#define GetNArrayData(obj,var)  TypedData_Get_Struct(obj, narray_data_t, &na_data_type, var)
-
-#define SDX_IS_STRIDE(x) ((x).stride&0x1)
-#define SDX_IS_INDEX(x)  (!SDX_IS_STRIDE(x))
-#define SDX_GET_STRIDE(x) ((x).stride>>1)
-#define SDX_GET_INDEX(x)  ((x).index)
-
-#define SDX_SET_STRIDE(x,s) ((x).stride=((s)<<1)|0x1)
-#define SDX_SET_INDEX(x,idx) ((x).index=idx)
-
-#define RNARRAY(val)            ((narray_t*)DATA_PTR(val))
-#define RNARRAY_DATA(val)       ((narray_data_t*)DATA_PTR(val))
-#define RNARRAY_VIEW(val)       ((narray_view_t*)DATA_PTR(val))
-#define RNARRAY_FILEMAP(val)    ((narray_filemap_t*)DATA_PTR(val))
-
-#define RNARRAY_NDIM(val)       (RNARRAY(val)->ndim)
-#define RNARRAY_TYPE(val)       (RNARRAY(val)->type)
-#define RNARRAY_FLAG(val)       (RNARRAY(val)->flag)
-#define RNARRAY_SIZE(val)       (RNARRAY(val)->size)
-#define RNARRAY_SHAPE(val)      (RNARRAY(val)->shape)
-#define RNARRAY_REDUCE(val)     (RNARRAY(val)->reduce)
-
-#define RNARRAY_DATA_PTR(val)    (RNARRAY_DATA(val)->ptr)
-#define RNARRAY_VIEW_DATA(val)   (RNARRAY_VIEW(val)->data)
-#define RNARRAY_VIEW_OFFSET(val) (RNARRAY_VIEW(val)->offset)
-#define RNARRAY_VIEW_STRIDX(val) (RNARRAY_VIEW(val)->stridx)
-
-#define NA_NDIM(na)     (((narray_t*)na)->ndim)
-#define NA_TYPE(na)     (((narray_t*)na)->type)
-#define NA_SIZE(na)     (((narray_t*)na)->size)
-#define NA_SHAPE(na)    (((narray_t*)na)->shape)
-#define NA_REDUCE(na)   (((narray_t*)na)->reduce)
-
-#define NA_FLAG(obj)    (na_get_narray_t(obj)->flag)
-#define NA_FLAG0(obj)   (NA_FLAG(obj)[0])
-#define NA_FLAG1(obj)   (NA_FLAG(obj)[1])
-
-#define NA_DATA(na)             ((narray_data_t*)(na))
-#define NA_VIEW(na)             ((narray_view_t*)(na))
-#define NA_DATA_PTR(na)         (NA_DATA(na)->ptr)
-#define NA_VIEW_DATA(na)        (NA_VIEW(na)->data)
-#define NA_VIEW_OFFSET(na)      (NA_VIEW(na)->offset)
-#define NA_VIEW_STRIDX(na)      (NA_VIEW(na)->stridx)
-
-#define NA_IS_INDEX_AT(na,i)    (SDX_IS_INDEX(NA_VIEW_STRIDX(na)[i]))
-#define NA_IS_STRIDE_AT(na,i)   (SDX_IS_STRIDE(NA_VIEW_STRIDX(na)[i]))
-#define NA_INDEX_AT(na,i)       (SDX_GET_INDEX(NA_VIEW_STRIDX(na)[i]))
-#define NA_STRIDE_AT(na,i)      (SDX_GET_STRIDE(NA_VIEW_STRIDX(na)[i]))
-
-#define NA_FILEMAP_PTR(na)      (((narray_filemap_t*)na)->ptr)
-
-
-#define NA_FL0_TEST(x,f) (NA_FLAG0(x)&(f))
-#define NA_FL1_TEST(x,f) (NA_FLAG1(x)&(f))
-
-#define NA_FL0_SET(x,f) do {NA_FLAG0(x) |= (f);} while(0)
-#define NA_FL1_SET(x,f) do {NA_FLAG1(x) |= (f);} while(0)
-
-#define NA_FL0_UNSET(x,f) do {NA_FLAG0(x) &= ~(f);} while(0)
-#define NA_FL1_UNSET(x,f) do {NA_FLAG1(x) &= ~(f);} while(0)
-
-#define NA_FL0_REVERSE(x,f) do {NA_FLAG0(x) ^= (f);} while(0)
-#define NA_FL1_REVERSE(x,f) do {NA_FLAG1(x) ^= (f);} while(0)
-
-
-/* FLAGS
-   - row-major / column-major
-   - Overwrite or not
-   - byteswapp
-   - Extensible?
-   - matrix or not
-*/
-
-#define NA_FL0_BIG_ENDIAN     (0x1<<0)
-#define NA_FL0_COLUMN_MAJOR   (0x1<<1)
-#define NA_FL1_LOCK           (0x1<<0)
-#define NA_FL1_INPLACE        (0x1<<1)
-
-#define TEST_COLUMN_MAJOR(x)   NA_FL0_TEST(x,NA_FL0_COLUMN_MAJOR)
-#define SET_COLUMN_MAJOR(x)    NA_FL0_SET(x,NA_FL0_COLUMN_MAJOR)
-#define UNSET_COLUMN_MAJOR(x)  NA_FL0_UNSET(x,NA_FL0_COLUMN_MAJOR)
-
-#define TEST_ROW_MAJOR(x)      (!TEST_COLUMN_MAJOR(x))
-#define SET_ROW_MAJOR(x)       UNSET_COLUMN_MAJOR(x)
-#define UNSET_ROW_MAJOR(x)     SET_COLUMN_MAJOR(x)
-
-#define TEST_BIG_ENDIAN(x)     NA_FL0_TEST(x,NA_FL0_BIG_ENDIAN)
-#define SET_BIG_ENDIAN(x)      NA_FL0_SET(x,NA_FL0_BIG_ENDIAN)
-#define UNSET_BIG_ENDIAN(x)    NA_FL0_UNSET(x,NA_FL0_BIG_ENDIAN)
-
-#define TEST_LITTLE_ENDIAN(x)  (!TEST_BIG_ENDIAN(x))
-#define SET_LITTLE_ENDIAN(x)   UNSET_BIG_ENDIAN(x)
-#define UNSET_LITTLE_ENDIAN(x) SET_BIG_ENDIAN(x)
-
-#define REVERSE_ENDIAN(x)      NA_FL0_REVERSE((x),NA_FL0_BIG_ENDIAN)
-
-#define TEST_LOCK(x)           NA_FL1_TEST(x,NA_FL1_LOCK)
-#define SET_LOCK(x)            NA_FL1_SET(x,NA_FL1_LOCK)
-#define UNSET_LOCK(x)          NA_FL1_UNSET(x,NA_FL1_LOCK)
-
-#define TEST_INPLACE(x)        NA_FL1_TEST(x,NA_FL1_INPLACE)
-#define SET_INPLACE(x)         NA_FL1_SET(x,NA_FL1_INPLACE)
-#define UNSET_INPLACE(x)       NA_FL1_UNSET(x,NA_FL1_INPLACE)
-
-#ifdef DYNAMIC_ENDIAN
-// not supported
-#else
-#ifdef WORDS_BIGENDIAN
-#define TEST_HOST_ORDER(x)     TEST_BIG_ENDIAN(x)
-#define SET_HOST_ORDER(x)      SET_BIG_ENDIAN(x)
-#define UNSET_HOST_ORDER(x)    UNSET_BIG_ENDIAN(x)
-#define TEST_BYTE_SWAPPED(x)   TEST_LITTLE_ENDIAN(x)
-#define SET_BYTE_SWAPPED(x)    SET_LITTLE_ENDIAN(x)
-#define UNSET_BYTE_SWAPPED(x)  UNSET_LITTLE_ENDIAN(x)
-#define NA_FL0_INIT            NA_FL0_BIG_ENDIAN
-#else // LITTLE ENDIAN
-#define TEST_HOST_ORDER(x)     TEST_LITTLE_ENDIAN(x)
-#define SET_HOST_ORDER(x)      SET_LITTLE_ENDIAN(x)
-#define UNSET_HOST_ORDER(x)    UNSET_LITTLE_ENDIAN(x)
-#define TEST_BYTE_SWAPPED(x)   TEST_BIG_ENDIAN(x)
-#define SET_BYTE_SWAPPED(x)    SET_BIG_ENDIAN(x)
-#define UNSET_BYTE_SWAPPED(x)  UNSET_BIG_ENDIAN(x)
-#define NA_FL0_INIT            0
-#endif
-#endif
-#define NA_FL1_INIT            0
-
-
-#define IsNArray(obj) (rb_obj_is_kind_of(obj,cNArray)==Qtrue)
-
-#define DEBUG_PRINT(v) puts(StringValueCStr(rb_funcall(v,rb_intern("inspect"),0)))
-
-#define NA_IsNArray(obj) \
-  (rb_obj_is_kind_of(obj,cNArray)==Qtrue)
-#define NA_IsArray(obj) \
-  (TYPE(obj)==T_ARRAY || rb_obj_is_kind_of(obj,cNArray)==Qtrue)
-
-#define NUM2REAL(v)  NUM2DBL( rb_funcall((v),na_id_real,0) )
-#define NUM2IMAG(v)  NUM2DBL( rb_funcall((v),na_id_imag,0) )
-
-#define NA_MAX_DIMENSION (int)(sizeof(VALUE)*8-2)
-#define NA_MAX_ELMSZ     65535
-
-typedef unsigned int BIT_DIGIT;
-//#define BYTE_BIT_DIGIT sizeof(BIT_DIGIT)
-#define NB     (sizeof(BIT_DIGIT)*8)
-#define BALL   (~(BIT_DIGIT)0)
-#define SLB(n) (((n)==NB)?~(BIT_DIGIT)0:(~(~(BIT_DIGIT)0<<(n))))
-
-#define ELEMENT_BIT_SIZE  "ELEMENT_BIT_SIZE"
-#define ELEMENT_BYTE_SIZE "ELEMENT_BYTE_SIZE"
-#define CONTIGUOUS_STRIDE "CONTIGUOUS_STRIDE"
-
-
-#ifdef RUBY_INTEGER_UNIFICATION
-#define IS_INTEGER_CLASS(c) ((c)==rb_cInteger)
-#else
-#define IS_INTEGER_CLASS(c) ((c)==rb_cFixnum||(c)==rb_cBignum)
-#endif
-
-#include "numo/ndloop.h"
-#include "numo/intern.h"
-
-#if defined(__cplusplus)
-#if 0
-{ /* satisfy cc-mode */
-#endif
-}  /* extern "C" { */
-#endif
-
-#endif /* ifndef NARRAY_H */