numo-narray 0.9.0.1-x64-mingw32

data/ext/numo/narray/depend.erb

@@ -0,0 +1,33 @@
+TAGSRC = \
+ ../../ruby/include/ruby/*.h \
+ ../../ruby/*.c \
+ *.h \
+ types/*.h \
+ *.c \
+ types/*.c
+
+tags : TAGS
+TAGS : $(TAGSRC)
+	etags $(TAGSRC)
+
+doc :
+	yard doc *.c types/*.c
+
+C_TMPL = <%=Dir.glob("#{__dir__}/gen/tmpl/*.c").join(" ")%>
+
+COGEN = <%= __dir__ %>/gen/cogen.rb
+DTYPE = <%= __dir__ %>/gen/dtype.erb.c
+COGEN_DTYPE = $(COGEN) -l $(DTYPE)
+DEPENDS = $(C_TMPL) $(DTYPE) <%= __dir__ %>/gen/*.rb
+
+<% Dir.glob("#{__dir__}/gen/def/*.rb") do |s| %>
+types/<%=File.basename(s,".rb")%>.c: <%=s%> $(DEPENDS)
+	$(MAKEDIRS) $(@D) types
+	ruby $(COGEN_DTYPE) <%=s%> > $@
+<% end %>
+
+types/bit.c: <%= __dir__ %>/gen/bit.erb.c $(DEPENDS)
+	$(MAKEDIRS) $(@D) types
+	ruby $(COGEN) -l <%= __dir__ %>/gen/bit.erb.c > $@
+
+CLEANOBJS = *.o */*.o *.bak types/*.c

data/ext/numo/narray/extconf.rb

@@ -0,0 +1,117 @@
+require 'rbconfig.rb'
+
+#RbConfig::MAKEFILE_CONFIG["optflags"] = "-g3 -gdwarf-2"
+
+require 'mkmf'
+
+require "erb"
+
+if RUBY_VERSION < "2.0.0"
+  puts "Numo::NArray requires Ruby version 2.0 or later."
+  exit(1)
+end
+
+#$CFLAGS="-g3 -O0 -Wall"
+#$CFLAGS=" $(cflags) -O3 -m64 -msse2 -funroll-loops"
+#$CFLAGS=" $(cflags) -O3"
+$INCFLAGS = "-Itypes #$INCFLAGS"
+
+$INSTALLFILES = Dir.glob(%w[numo/*.h numo/types/*.h]).map{|x| [x,'$(archdir)'] }
+if /cygwin|mingw/ =~ RUBY_PLATFORM
+  $INSTALLFILES << ['libnarray.a', '$(archdir)']
+end
+
+srcs = %w(
+narray
+array
+step
+index
+ndloop
+data
+types/bit
+types/int8
+types/int16
+types/int32
+types/int64
+types/uint8
+types/uint16
+types/uint32
+types/uint64
+types/sfloat
+types/dfloat
+types/scomplex
+types/dcomplex
+types/robject
+math
+SFMT
+struct
+rand
+)
+
+=begin
+have_header("atlas/cblas.h")
+have_library("atlas")
+
+if have_library("blas")
+  if have_library("lapack")
+    srcs.push "linalg"
+    $defs.push "-DHAVE_LAPACK"
+  else
+    #$defs.delete "-DHAVE_LAPACK"
+  end
+end
+=end
+
+if have_header("stdbool.h")
+  stdbool = "stdbool.h"
+else
+  stdbool = nil
+end
+
+if have_header("stdint.h")
+  stdint = "stdint.h"
+elsif have_header("sys/types.h")
+  stdint = "sys/types.h"
+else
+  stdint = nil
+end
+
+have_type("bool", stdbool)
+unless have_type("u_int8_t", stdint)
+  have_type("uint8_t",stdint)
+end
+unless have_type("u_int16_t", stdint)
+  have_type("uint16_t",stdint)
+end
+have_type("int32_t", stdint)
+unless have_type("u_int32_t", stdint)
+  have_type("uint32_t",stdint)
+end
+have_type("int64_t", stdint)
+unless have_type("u_int64_t", stdint)
+  have_type("uint64_t", stdint)
+end
+#have_library("m")
+#have_func("sincos")
+#have_func("asinh")
+have_func("exp10")
+
+have_var("rb_cComplex")
+#have_func("rb_alloc_tmp_buffer", "ruby.h")
+#have_func("rb_free_tmp_buffer", "ruby.h")
+
+$objs = srcs.collect{|i| i+".o"}
+
+create_header
+
+depend_path = File.join(__dir__, "depend")
+File.open(depend_path, "w") do |depend|
+  depend_erb_path = File.join(__dir__, "depend.erb")
+  File.open(depend_erb_path, "r") do |depend_erb|
+    erb = ERB.new(depend_erb.read)
+    erb.filename = depend_erb_path
+    depend.print(erb.result)
+  end
+end
+
+create_makefile('numo/narray')

data/ext/numo/narray/gen/bit.erb.c

@@ -0,0 +1,811 @@
+/*
+  bit.c
+  Numerical Array Extension for Ruby
+    (C) Copyright 1999-2011,2013 by Masahiro TANAKA
+
+  This program is free software.
+  You can distribute/modify this program
+  under the same terms as Ruby itself.
+  NO WARRANTY.
+*/
+#include <ruby.h>
+#include "numo/narray.h"
+#include "numo/template.h"
+
+<%
+$embed = true
+class_name "Bit"
+%>
+
+#define cT <%=type_var%>
+
+typedef int dtype;
+
+#define m_load_data(ptr,pos) load_data(ptr,pos)
+#define m_sprintf(s,x)   sprintf(s,"%1d",(int)(x))
+
+#define m_copy(x)  (x)
+#define m_not(x)   (~(x))
+#define m_and(x,y) ((x)&(y))
+#define m_or(x,y)  ((x)|(y))
+#define m_xor(x,y) ((x)^(y))
+#define m_eq(x,y)  (~((x)^(y)))
+#define m_count_true(x)  (x!=0)
+#define m_count_false(x) (x==0)
+
+static inline dtype load_data(void *ptr, size_t pos) {
+    return (((BIT_DIGIT*)(ptr))[(pos)/NB]>>((pos)%NB)) & 1u;
+}
+
+VALUE <%=type_var%>;
+
+static VALUE
+numo_bit_allocate(VALUE self)
+{
+    narray_t *na;
+    char *ptr;
+
+    GetNArray(self,na);
+
+    switch(NA_TYPE(na)) {
+    case NARRAY_DATA_T:
+        ptr = NA_DATA_PTR(na);
+        if (na->size > 0 && ptr == NULL) {
+            ptr = xmalloc(((na->size-1)/sizeof(BIT_DIGIT)+1)*sizeof(BIT_DIGIT)/8);
+            NA_DATA_PTR(na) = ptr;
+        }
+        break;
+    case NARRAY_FILEMAP_T:
+        //ptr = ((narray_filemap_t*)na)->ptr;
+        // to be implemented
+        break;
+    case NARRAY_VIEW_T:
+        rb_funcall(NA_VIEW_DATA(na), rb_intern("allocate"), 0);
+        break;
+    default:
+        rb_raise(rb_eRuntimeError,"invalid narray type");
+    }
+    return self;
+}
+
+
+static VALUE numo_cast_array_to_bit(VALUE ary);
+
+static VALUE
+numo_bit_cast_numeric(VALUE val)
+{
+    VALUE v;
+    BIT_DIGIT *ptr, b=2;
+    size_t dig_ofs;
+    int    bit_ofs;
+
+    if (FIXNUM_P(val)) {
+        b = FIX2INT(val);
+    } else if (val==Qtrue) {
+        b = 1;
+    } else if (val==Qfalse) {
+        b = 0;
+    }
+    if (b!=0 && b!=1) {
+        rb_raise(rb_eArgError, "bit can be cast from 0 or 1 or true or false");
+    }
+
+    v = rb_narray_new(cT, 0, NULL);
+    //dig_ofs = na->offset / NB;
+    //bit_ofs = na->offset % NB;
+    dig_ofs = 0;
+    bit_ofs = 0;
+    ptr = (BIT_DIGIT*)na_get_pointer_for_write(v) + dig_ofs;
+    *ptr = (*ptr & ~(1u<<bit_ofs)) | (b<<bit_ofs);
+    na_release_lock(v);
+    return v;
+}
+
+
+static VALUE
+numo_bit_s_cast(VALUE type, VALUE obj)
+{
+    VALUE r;
+
+    if (CLASS_OF(obj)==cT) {
+        return obj;
+    } else if (TYPE(obj)==T_ARRAY) {
+        return numo_cast_array_to_bit(obj);
+    } else if (TYPE(obj)==T_FLOAT || FIXNUM_P(obj) || TYPE(obj)==T_BIGNUM) {
+        return numo_bit_cast_numeric(obj);
+    }
+
+    if (IsNArray(obj)) {
+        r = rb_funcall(obj, rb_intern("coerce_cast"), 1, cT);
+        if (RTEST(r)) {
+            return r;
+        }
+    }
+
+    rb_raise(nary_eCastError, "unknown conversion from %s to %s",
+             rb_class2name(CLASS_OF(obj)),
+             rb_class2name(type));
+    return Qnil;
+}
+
+static VALUE
+numo_bit_coerce_cast(VALUE value, VALUE type)
+{
+    return Qnil;
+}
+
+//----------------------------------------------------------------------
+
+static VALUE
+format_bit(VALUE fmt, int x)
+{
+    if (NIL_P(fmt)) {
+        char s[4];
+        int n;
+        n = sprintf(s,"%1d",x);
+        return rb_str_new(s,n);
+    }
+    return rb_funcall(fmt, '%', 1, INT2FIX(x));
+}
+
+static VALUE
+ bit_inspect_element(char *ptr, size_t pos, VALUE fmt)
+{
+    int x;
+    LOAD_BIT(ptr,pos,x);
+    return format_bit(fmt, x);
+}
+ VALUE
+ numo_bit_inspect(VALUE ary)
+{
+    return na_ndloop_inspect(ary, bit_inspect_element, Qnil);
+}
+
+
+static void
+iter_bit_format(na_loop_t *const lp)
+{
+    size_t  i;
+    BIT_DIGIT *a1, x=0;
+    size_t     p1;
+    char      *p2;
+    ssize_t    s1, s2;
+    size_t    *idx1;
+    VALUE y;
+    //VALUE fmt = *(VALUE*)(lp->opt_ptr);
+    VALUE fmt = lp->option;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR_BIT_IDX(lp, 0, a1, p1, s1, idx1);
+    INIT_PTR(lp, 1, p2, s2);
+
+    if (idx1) {
+        for (; i--;) {
+            LOAD_BIT(a1, p1 + *idx1, x); idx1++;
+            y = format_<%=tp%>(fmt, x);
+            SET_DATA_STRIDE(p2, s2, VALUE, y);
+        }
+    } else {
+        for (; i--;) {
+            LOAD_BIT(a1, p1, x);         p1+=s1;
+            y = format_<%=tp%>(fmt, x);
+            SET_DATA_STRIDE(p2, s2, VALUE, y);
+        }
+    }
+}
+
+static VALUE
+numo_bit_format(int argc, VALUE *argv, VALUE self)
+{
+    VALUE fmt=Qnil;
+    ndfunc_arg_in_t ain[2] = {{Qnil,0},{sym_option}};
+    ndfunc_arg_out_t aout[1] = {{numo_cRObject,0}};
+    ndfunc_t ndf = { iter_bit_format, FULL_LOOP, 2, 1, ain, aout };
+
+    rb_scan_args(argc, argv, "01", &fmt);
+    return na_ndloop(&ndf, 2, self, fmt);
+}
+
+
+static void
+iter_bit_format_to_a(na_loop_t *const lp)
+{
+    size_t  i;
+    BIT_DIGIT *a1, x=0;
+    size_t     p1;
+    ssize_t    s1;
+    size_t   *idx1;
+    VALUE y;
+    //VALUE fmt = *(VALUE*)(lp->opt_ptr);
+    VALUE fmt = lp->option;
+    volatile VALUE a;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR_BIT_IDX(lp, 0, a1, p1, s1, idx1);
+    a = rb_ary_new2(i);
+    rb_ary_push(lp->args[1].value, a);
+    if (idx1) {
+        for (; i--;) {
+            LOAD_BIT(a1, p1+*idx1, x); idx1++;
+            y = format_bit(fmt, x);
+            rb_ary_push(a,y);
+        }
+    } else {
+        for (; i--;) {
+            LOAD_BIT(a1, p1, x); p1+=s1;
+            y = format_bit(fmt, x);
+            rb_ary_push(a,y);
+        }
+    }
+}
+
+static VALUE
+numo_bit_format_to_a(int argc, VALUE *argv, VALUE self)
+{
+    volatile VALUE fmt=Qnil;
+    ndfunc_arg_in_t ain[3] = {{Qnil,0},{sym_loop_opt},{sym_option}};
+    ndfunc_arg_out_t aout[1] = {{rb_cArray,0}}; // dummy?
+    ndfunc_t ndf = { iter_bit_format_to_a, FULL_LOOP, 3, 1, ain, aout };
+
+    rb_scan_args(argc, argv, "01", &fmt);
+    return na_ndloop_cast_narray_to_rarray(&ndf, self, fmt);
+}
+
+
+static void
+iter_bit_fill(na_loop_t *const lp)
+{
+    size_t  n;
+    size_t  p3;
+    ssize_t s3;
+    size_t *idx3;
+    int     len;
+    BIT_DIGIT *a3;
+    BIT_DIGIT  y;
+    //VALUE x = *(VALUE*)(lp->opt_ptr);
+    VALUE x = lp->option;
+
+    if (x==INT2FIX(0) || x==Qfalse) {
+        y = 0;
+    } else
+    if (x==INT2FIX(1) || x==Qtrue) {
+        y = ~(BIT_DIGIT)0;
+    } else {
+        rb_raise(rb_eArgError, "invalid value for Bit");
+    }
+
+    INIT_COUNTER(lp, n);
+    INIT_PTR_BIT_IDX(lp, 0, a3, p3, s3, idx3);
+    if (idx3) {
+        y = y&1;
+        for (; n--;) {
+            STORE_BIT(a3, p3+*idx3, y); idx3++;
+        }
+    } else if (s3!=1) {
+        y = y&1;
+        for (; n--;) {
+            STORE_BIT(a3, p3, y); p3+=s3;
+        }
+    } else {
+        if (p3>0 || n<NB) {
+            len = NB - p3;
+            if ((int)n<len) len=n;
+            *a3 = (y & (SLB(len)<<p3)) | (*a3 & ~(SLB(len)<<p3));
+            a3++;
+            n -= len;
+        }
+        for (; n>=NB; n-=NB) {
+            *(a3++) = y;
+        }
+        if (n>0) {
+            *a3 = (y & SLB(n)) | (*a3 & BALL<<n);
+        }
+    }
+}
+
+
+static VALUE
+numo_bit_fill(VALUE self, VALUE val)
+{
+    ndfunc_arg_in_t ain[2] = {{OVERWRITE,0},{sym_option}};
+    ndfunc_t ndf = { iter_bit_fill, FULL_LOOP, 2, 0, ain, 0 };
+
+    na_ndloop(&ndf, 2, self, val);
+    return self;
+}
+
+
+void
+bit_cast_to_robj(na_loop_t *const lp)
+{
+    size_t     i;
+    ssize_t    s1;
+    size_t    *idx1;
+    BIT_DIGIT *a1;
+    size_t     p1;
+    BIT_DIGIT  x=0;
+    VALUE      y;
+    volatile VALUE a;
+
+    INIT_COUNTER(lp, i);
+    INIT_PTR_BIT_IDX(lp, 0, a1, p1, s1, idx1);
+    a = rb_ary_new2(i);
+    rb_ary_push(lp->args[1].value, a);
+    if (idx1) {
+        for (; i--;) {
+            LOAD_BIT(a1, p1+*idx1, x); idx1++;
+            y = INT2FIX(x);
+            rb_ary_push(a,y);
+        }
+    } else {
+        for (; i--;) {
+            LOAD_BIT(a1, p1, x); p1+=s1;
+            y = INT2FIX(x);
+            rb_ary_push(a,y);
+        }
+    }
+}
+
+static VALUE
+numo_bit_cast_to_rarray(VALUE self)
+{
+    ndfunc_arg_in_t ain[3] = {{Qnil,0},{sym_loop_opt},{sym_option}};
+    ndfunc_arg_out_t aout[1] = {{rb_cArray,0}}; // dummy?
+    ndfunc_t ndf = { bit_cast_to_robj, FULL_LOOP, 3, 1, ain, aout };
+    return na_ndloop_cast_narray_to_rarray(&ndf, self, Qnil);
+}
+
+
+static void
+iter_cast_rarray_to_bit(na_loop_t *const lp)
+{
+    size_t i, n, n1;
+    VALUE  v1, *ptr;
+    BIT_DIGIT *a2;
+    size_t p2;
+    size_t s2, *idx2;
+    VALUE  x;
+    BIT_DIGIT y;
+
+    INIT_COUNTER(lp, n);
+    INIT_PTR_BIT_IDX(lp, 1, a2, p2, s2, idx2);
+    v1 = lp->args[0].value;
+    ptr = &v1;
+
+    switch(TYPE(v1)) {
+    case T_ARRAY:
+        n1 = RARRAY_LEN(v1);
+        ptr = RARRAY_PTR(v1);
+        break;
+    case T_NIL:
+        n1 = 0;
+        break;
+    default:
+        n1 = 1;
+    }
+    if (idx2) {
+        <%
+        ["STORE_BIT(a2,p2+*idx2,y); idx2++;",
+         "STORE_BIT(a2,p2,y); p2+=s2;"].each_with_index do |s,i|
+        %>
+        for (i=0; i<n1 && i<n; i++) {
+            x = ptr[i];
+            y = 2;
+            if (FIXNUM_P(x)) {
+                y = FIX2INT(x);
+            } else if (x==Qtrue) {
+                y = 1;
+            } else if (x==Qfalse) {
+                y = 0;
+            } else if (x==Qnil) {
+                y = 0;
+            }
+            if (y!=0 && y!=1) {
+                rb_raise(rb_eArgError, "bit can be cast from 0 or 1 or true or false");
+            }
+            <%= s %>
+        }
+        y = 0;
+        for (; i<n; i++) {
+            <%= s %>
+        }
+        <% if i<1 %>
+    } else {
+        <% end; end %>
+    }
+}
+
+static VALUE
+numo_cast_array_to_bit(VALUE rary)
+{
+    volatile VALUE vnc, nary;
+    na_compose_t *nc;
+    ndfunc_arg_in_t ain[2] = {{Qnil,0},{rb_cArray,0}};
+    ndfunc_t ndf = { iter_cast_rarray_to_bit, FULL_LOOP, 2, 0, ain, 0 };
+
+    vnc = na_ary_composition(rary);
+    Data_Get_Struct(vnc, na_compose_t, nc);
+    nary = rb_narray_new(cT, nc->ndim, nc->shape);
+    numo_bit_allocate(nary);
+    na_ndloop_cast_rarray_to_narray(&ndf, rary, nary);
+    return nary;
+}
+
+
+static VALUE
+numo_bit_extract(VALUE self)
+{
+    BIT_DIGIT *ptr, val;
+    size_t pos;
+    narray_t *na;
+    GetNArray(self,na);
+
+    if (na->ndim==0) {
+        pos = na_get_offset(self);
+        ptr = (BIT_DIGIT*)na_get_pointer_for_read(self);
+        val = ((*((ptr)+(pos)/NB)) >> ((pos)%NB)) & 1u;
+        na_release_lock(self);
+        return INT2FIX(val);
+    }
+    return self;
+}
+
+<%
+bit_binary "and", "&"
+bit_binary "or" , "|"
+bit_binary "xor", "^"
+bit_binary "eq"
+bit_unary  "not", "~"
+bit_unary  "copy"
+bit_count  "count_true"
+bit_count  "count_false"
+Function.codes.each do |x|
+%>
+<%= x %>
+<% end %>
+
+
+/* !!! Bit#store: under construction !!! */
+VALUE
+numo_bit_store(VALUE dst, VALUE src)
+{
+    // check and fix me
+    ndfunc_arg_in_t ain[2] = {{INT2FIX(1),0},{Qnil,0}};
+    ndfunc_t ndf = { iter_bit_copy, FULL_LOOP, 2, 0, ain, 0 };
+
+    na_ndloop(&ndf, 2, src, dst);
+    return src;
+}
+
+//VALUE na_aref(int argc, VALUE *argv, VALUE self);
+
+/* method: []=(idx1,idx2,...,idxN,val) */
+static VALUE
+numo_bit_aset(int argc, VALUE *argv, VALUE self)
+{
+    VALUE a;
+    argc--;
+
+    if (argc==0)
+        numo_bit_store(self, argv[argc]);
+    else {
+        //a = na_aref(argc, argv, self);
+        a = rb_funcall2(self, rb_intern("[]"), argc, argv);
+        numo_bit_store(a, argv[argc]);
+    }
+    return argv[argc];
+}
+
+
+typedef struct {
+    size_t count;
+    char  *idx0;
+    char  *idx1;
+    size_t elmsz;
+} where_opt_t;
+
+#define STORE_INT(ptr, esz, x) memcpy(ptr,&(x),esz)
+
+static void
+iter_bit_where(na_loop_t *const lp)
+{
+    size_t  i;
+    BIT_DIGIT *a;
+    size_t  p;
+    ssize_t s;
+    size_t *idx;
+    BIT_DIGIT x=0;
+    char   *idx0, *idx1;
+    size_t  count;
+    size_t  e;
+    where_opt_t *g;
+
+    g = (where_opt_t*)(lp->opt_ptr);
+    count = g->count;
+    idx0  = g->idx0;
+    idx1  = g->idx1;
+    e     = g->elmsz;
+    INIT_COUNTER(lp, i);
+    INIT_PTR_BIT_IDX(lp, 0, a, p, s, idx);
+    if (idx) {
+        for (; i--;) {
+            LOAD_BIT(a, p+*idx, x);
+            idx++;
+            if (x==0) {
+                if (idx0) {
+                    STORE_INT(idx0,e,count);
+                    idx0 += e;
+                }
+            } else {
+                if (idx1) {
+                    STORE_INT(idx1,e,count);
+                    idx1 += e;
+                }
+            }
+            count++;
+        }
+    } else {
+        for (; i--;) {
+            LOAD_BIT(a, p, x);
+            p+=s;
+            if (x==0) {
+                if (idx0) {
+                    STORE_INT(idx0,e,count);
+                    idx0 += e;
+                }
+            } else {
+                if (idx1) {
+                    STORE_INT(idx1,e,count);
+                    idx1 += e;
+                }
+            }
+            count++;
+        }
+    }
+    g->count = count;
+    g->idx0  = idx0;
+    g->idx1  = idx1;
+}
+
+static VALUE
+ numo_bit_where(VALUE self)
+{
+    volatile VALUE idx_1;
+    size_t size, n_1;
+    where_opt_t *g;
+
+    ndfunc_arg_in_t ain[1] = {{cT,0}};
+    ndfunc_t ndf = { iter_bit_where, FULL_LOOP, 1, 0, ain, 0 };
+
+    size = RNARRAY_SIZE(self);
+    n_1 = NUM2SIZE(numo_bit_count_true(0, NULL, self));
+    g = ALLOCA_N(where_opt_t,1);
+    g->count = 0;
+    if (size>4294967295ul) {
+        idx_1 = rb_narray_new(numo_cInt64, 1, &n_1);
+        g->elmsz = 8;
+    } else {
+        idx_1 = rb_narray_new(numo_cInt32, 1, &n_1);
+        g->elmsz = 4;
+    }
+    g->idx1 = na_get_pointer_for_write(idx_1);
+    g->idx0 = NULL;
+    na_ndloop3(&ndf, g, 1, self);
+    na_release_lock(idx_1);
+    return idx_1;
+}
+
+static VALUE
+ numo_bit_where2(VALUE self)
+{
+    VALUE idx_1, idx_0;
+    size_t size, n_1;
+    where_opt_t *g;
+
+    ndfunc_arg_in_t ain[1] = {{cT,0}};
+    ndfunc_t ndf = { iter_bit_where, FULL_LOOP, 1, 0, ain, 0 };
+
+    size = RNARRAY_SIZE(self);
+    n_1 = NUM2SIZE(numo_bit_count_true(0, NULL, self));
+    g = ALLOCA_N(where_opt_t,1);
+    g->count = 0;
+    if (size>4294967295ul) {
+        idx_1 = rb_narray_new(numo_cInt64, 1, &n_1);
+        idx_0 = rb_narray_new(numo_cInt64, 1, &n_1);
+        g->elmsz = 8;
+    } else {
+        idx_1 = rb_narray_new(numo_cInt32, 1, &n_1);
+        idx_0 = rb_narray_new(numo_cInt32, 1, &n_1);
+        g->elmsz = 4;
+    }
+    g->idx1 = na_get_pointer_for_write(idx_1);
+    g->idx0 = na_get_pointer_for_write(idx_0);
+    na_ndloop3(&ndf, g, 1, self);
+    na_release_lock(idx_0);
+    na_release_lock(idx_1);
+    return rb_assoc_new(idx_1,idx_0);
+}
+
+
+static void
+iter_bit_pointer(na_loop_t *const lp)
+{
+    size_t  i;
+    BIT_DIGIT *a;
+    size_t  p1, p2;
+    ssize_t s1, s2;
+    size_t *idx1, *idx2, *pidx;
+    BIT_DIGIT x=0;
+    size_t  count;
+    where_opt_t *g;
+
+    g = (where_opt_t*)(lp->opt_ptr);
+    count = g->count;
+    pidx  = (size_t*)(g->idx1);
+    INIT_COUNTER(lp, i);
+    INIT_PTR_BIT_IDX(lp, 0, a, p1, s1, idx1);
+    //INIT_PTR_IDX(lp, 1, p2, s2, idx2);
+    p2 = lp->args[1].iter[0].pos;
+    s2 = lp->args[1].iter[0].step;
+    idx2 = lp->args[1].iter[0].idx;
+
+    if (idx1) {
+        if (idx2) {
+            for (; i--;) {
+                LOAD_BIT(a, p1+*idx1, x);
+                idx1++;
+                if (x) {
+                    *(pidx++) = p2+*idx2;
+                    count++;
+                }
+                idx2++;
+            }
+        } else {
+            for (; i--;) {
+                LOAD_BIT(a, p1+*idx1, x);
+                idx1++;
+                if (x) {
+                    *(pidx++) = p2;
+                    count++;
+                }
+                p2 += s2;
+            }
+        }
+    } else {
+        if (idx2) {
+            for (; i--;) {
+                LOAD_BIT(a, p1, x);
+                p1 += s1;
+                if (x) {
+                    *(pidx++) = p2+*idx2;
+                    count++;
+                }
+                idx2++;
+            }
+        } else {
+            for (; i--;) {
+                LOAD_BIT(a, p1, x);
+                p1 += s1;
+                if (x) {
+                    *(pidx++) = p2;
+                    count++;
+                }
+                p2 += s2;
+            }
+        }
+    }
+    g->count = count;
+    g->idx1  = (char*)pidx;
+}
+
+#if   SIZEOF_VOIDP == 8
+#define cIndex numo_cInt64
+#elif SIZEOF_VOIDP == 4
+#define cIndex numo_cInt32
+#endif
+
+static VALUE
+ numo_bit_mask(VALUE mask, VALUE val)
+{
+    volatile VALUE idx_1, view;
+    narray_data_t *nidx;
+    narray_view_t *nv;
+    stridx_t stridx0;
+    size_t n_1;
+    where_opt_t g;
+    ndfunc_arg_in_t ain[2] = {{cT,0},{Qnil,0}};
+    ndfunc_t ndf = {iter_bit_pointer, FULL_LOOP, 2, 0, ain, 0};
+
+    n_1 = NUM2SIZE(numo_bit_count_true(0, NULL, mask));
+    idx_1 = rb_narray_new(cIndex, 1, &n_1);
+    g.count = 0;
+    g.elmsz = SIZEOF_VOIDP;
+    g.idx1 = na_get_pointer_for_write(idx_1);
+    g.idx0 = NULL;
+    na_ndloop3(&ndf, &g, 2, mask, val);
+
+    view = na_s_allocate_view(CLASS_OF(val));
+    GetNArrayView(view, nv);
+    na_setup_shape((narray_t*)nv, 1, &n_1);
+
+    GetNArrayData(idx_1,nidx);
+    SDX_SET_INDEX(stridx0,(size_t*)nidx->ptr);
+    nidx->ptr = NULL;
+
+    nv->stridx = ALLOC_N(stridx_t,1);
+    nv->stridx[0] = stridx0;
+    nv->offset = 0;
+    nv->data = val;
+    return view;
+}
+
+
+VALUE
+numo_bit_all_p(VALUE self)
+{
+    return (rb_funcall(self, rb_intern("count_false"), 0)==INT2FIX(0)) ? Qtrue : Qfalse;
+}
+
+VALUE
+numo_bit_any_p(VALUE self)
+{
+    return (rb_funcall(self, rb_intern("count_true"), 0)!=INT2FIX(0)) ? Qtrue : Qfalse;
+}
+
+VALUE
+numo_bit_none_p(VALUE self)
+{
+    return (rb_funcall(self, rb_intern("count_true"), 0)==INT2FIX(0)) ? Qtrue : Qfalse;
+}
+
+
+void
+Init_nary_bit()
+{
+    volatile VALUE hCast;
+
+    cT = rb_define_class_under(mNumo, "Bit", cNArray);
+
+    rb_define_const(cT, "ELEMENT_BIT_SIZE",  INT2FIX(1));
+    rb_define_const(cT, "ELEMENT_BYTE_SIZE", rb_float_new(1.0/8));
+    rb_define_const(cT, "CONTIGUOUS_STRIDE", INT2FIX(1));
+
+    rb_define_method(cT, "allocate", numo_bit_allocate, 0);
+
+    rb_define_singleton_method(cT, "cast", numo_bit_s_cast, 1);
+    rb_define_singleton_method(cT, "[]", numo_bit_s_cast, -2);
+    rb_define_method(cT, "coerce_cast", numo_bit_coerce_cast, 1);
+
+    <% Function.definitions.each do |x| %>
+    <%= x %><% end %>
+
+    rb_define_alias (cT, "count_1","count_true");
+    rb_define_alias (cT, "count_0","count_false");
+    rb_define_method(cT, "where", numo_bit_where, 0);
+    rb_define_method(cT, "where2", numo_bit_where2, 0);
+    rb_define_method(cT, "mask", numo_bit_mask, 1);
+
+    rb_define_method(cT, "all?", numo_bit_all_p, 0);
+    rb_define_method(cT, "any?", numo_bit_any_p, 0);
+    rb_define_method(cT, "none?", numo_bit_none_p, 0);
+
+    rb_define_method(cT, "inspect", numo_bit_inspect, 0);
+    rb_define_method(cT, "format", numo_bit_format, -1);
+    rb_define_method(cT, "format_to_a", numo_bit_format_to_a, -1);
+
+    rb_define_method(cT, "fill", numo_bit_fill, 1);
+
+    rb_define_method(cT, "to_a", numo_bit_cast_to_rarray, 0);
+
+    rb_define_method(cT, "extract", numo_bit_extract, 0);
+
+    rb_define_method(cT, "copy",  numo_bit_copy, 0);
+    rb_define_method(cT, "store", numo_bit_store, 1);
+    rb_define_method(cT, "[]=",   numo_bit_aset, -1);
+
+    hCast = rb_hash_new();
+    rb_define_const(cT, "UPCAST", hCast);
+    rb_hash_aset(hCast, numo_cInt32, numo_cInt32);
+    rb_hash_aset(hCast, numo_cInt16, numo_cInt16);
+    rb_hash_aset(hCast, numo_cInt8,  numo_cInt8);
+}