numo-narray 0.9.0.1 → 0.9.0.3

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: numo-narray
 version: !ruby/object:Gem::Version
-  version: 0.9.0.1
+  version: 0.9.0.3
 platform: ruby
 authors:
 - Masahiro TANAKA
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2016-06-21 00:00:00.000000000 Z
+date: 2016-08-13 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -99,8 +99,8 @@ files:
 - ext/numo/narray/data.c
 - ext/numo/narray/depend.erb
 - ext/numo/narray/extconf.rb
-- ext/numo/narray/gen/bit.erb.c
 - ext/numo/narray/gen/cogen.rb
+- ext/numo/narray/gen/def/bit.rb
 - ext/numo/narray/gen/def/dcomplex.rb
 - ext/numo/narray/gen/def/dfloat.rb
 - ext/numo/narray/gen/def/int16.rb
@@ -115,6 +115,7 @@ files:
 - ext/numo/narray/gen/def/uint64.rb
 - ext/numo/narray/gen/def/uint8.rb
 - ext/numo/narray/gen/dtype.erb.c
+- ext/numo/narray/gen/spec.rb
 - ext/numo/narray/gen/tmpl/accum.c
 - ext/numo/narray/gen/tmpl/accum_binary.c
 - ext/numo/narray/gen/tmpl/accum_index.c
@@ -124,9 +125,6 @@ files:
 - ext/numo/narray/gen/tmpl/binary.c
 - ext/numo/narray/gen/tmpl/binary2.c
 - ext/numo/narray/gen/tmpl/binary_s.c
-- ext/numo/narray/gen/tmpl/bit_binary.c
-- ext/numo/narray/gen/tmpl/bit_count.c
-- ext/numo/narray/gen/tmpl/bit_unary.c
 - ext/numo/narray/gen/tmpl/cast.c
 - ext/numo/narray/gen/tmpl/cast_array.c
 - ext/numo/narray/gen/tmpl/cast_numeric.c
@@ -143,6 +141,7 @@ files:
 - ext/numo/narray/gen/tmpl/format_to_a.c
 - ext/numo/narray/gen/tmpl/head.c
 - ext/numo/narray/gen/tmpl/inspect.c
+- ext/numo/narray/gen/tmpl/logseq.c
 - ext/numo/narray/gen/tmpl/map_with_index.c
 - ext/numo/narray/gen/tmpl/median.c
 - ext/numo/narray/gen/tmpl/minmax.c
@@ -159,12 +158,42 @@ files:
 - ext/numo/narray/gen/tmpl/sort_index.c
 - ext/numo/narray/gen/tmpl/store.c
 - ext/numo/narray/gen/tmpl/store_array.c
+- ext/numo/narray/gen/tmpl/store_bit.c
 - ext/numo/narray/gen/tmpl/store_from.c
 - ext/numo/narray/gen/tmpl/store_numeric.c
 - ext/numo/narray/gen/tmpl/to_a.c
 - ext/numo/narray/gen/tmpl/unary.c
 - ext/numo/narray/gen/tmpl/unary2.c
 - ext/numo/narray/gen/tmpl/unary_s.c
+- ext/numo/narray/gen/tmpl_bit/allocate.c
+- ext/numo/narray/gen/tmpl_bit/aref.c
+- ext/numo/narray/gen/tmpl_bit/aset.c
+- ext/numo/narray/gen/tmpl_bit/binary.c
+- ext/numo/narray/gen/tmpl_bit/bit_count.c
+- ext/numo/narray/gen/tmpl_bit/bit_reduce.c
+- ext/numo/narray/gen/tmpl_bit/cast.c
+- ext/numo/narray/gen/tmpl_bit/cast_array.c
+- ext/numo/narray/gen/tmpl_bit/cast_numeric.c
+- ext/numo/narray/gen/tmpl_bit/coerce_cast.c
+- ext/numo/narray/gen/tmpl_bit/each.c
+- ext/numo/narray/gen/tmpl_bit/each_with_index.c
+- ext/numo/narray/gen/tmpl_bit/extract.c
+- ext/numo/narray/gen/tmpl_bit/fill.c
+- ext/numo/narray/gen/tmpl_bit/format.c
+- ext/numo/narray/gen/tmpl_bit/format_to_a.c
+- ext/numo/narray/gen/tmpl_bit/inspect.c
+- ext/numo/narray/gen/tmpl_bit/map_with_index.c
+- ext/numo/narray/gen/tmpl_bit/mask.c
+- ext/numo/narray/gen/tmpl_bit/none_p.c
+- ext/numo/narray/gen/tmpl_bit/store.c
+- ext/numo/narray/gen/tmpl_bit/store_array.c
+- ext/numo/narray/gen/tmpl_bit/store_bit.c
+- ext/numo/narray/gen/tmpl_bit/store_from.c
+- ext/numo/narray/gen/tmpl_bit/store_numeric.c
+- ext/numo/narray/gen/tmpl_bit/to_a.c
+- ext/numo/narray/gen/tmpl_bit/unary.c
+- ext/numo/narray/gen/tmpl_bit/where.c
+- ext/numo/narray/gen/tmpl_bit/where2.c
 - ext/numo/narray/index.c
 - ext/numo/narray/kwarg.c
 - ext/numo/narray/math.c

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

@@ -1,811 +0,0 @@
-/*
-  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);
-}