numo-narray 0.9.1.4 → 0.9.1.5

data/ext/numo/narray/math.c

@@ -1,7 +1,7 @@
 /*
   math.c
-  Numerical Array Extension for Ruby
-    (C) Copyright 1999-2017 by Masahiro TANAKA
+  Ruby/Numo::NArray - Numerical Array class for Ruby
+    Copyright (C) 1999-2019 Masahiro TANAKA
 */
 #include <ruby.h>
 #include "numo/narray.h"

data/ext/numo/narray/narray.c

@@ -1,7 +1,7 @@
 /*
   narray.c
-  Numerical Array Extension for Ruby
-    (C) Copyright 1999-2017 by Masahiro TANAKA
+  Ruby/Numo::NArray - Numerical Array class for Ruby
+    Copyright (C) 1999-2019 Masahiro TANAKA
 */
 #define NARRAY_C
 #include <ruby.h>
@@ -42,17 +42,21 @@ VALUE sym_option;
 VALUE sym_loop_opt;
 VALUE sym_init;
 
-VALUE na_cStep;
 #ifndef HAVE_RB_CCOMPLEX
 VALUE rb_cComplex;
 #endif
+#ifdef HAVE_RB_ARITHMETIC_SEQUENCE_EXTRACT
+VALUE rb_cArithSeq;
+#endif
 
 int numo_na_inspect_rows=20;
 int numo_na_inspect_cols=80;
 
 void Init_nary_data();
 void Init_nary_ndloop();
+#ifndef HAVE_RB_ARITHMETIC_SEQUENCE_EXTRACT
 void Init_nary_step();
+#endif
 void Init_nary_index();
 void Init_numo_bit();
 void Init_numo_int8();
@@ -348,19 +352,19 @@ na_setup(VALUE self, int ndim, size_t *shape)
 
   @example
     i = Numo::Int64.new([2,4,3])
-    #=> Numo::Int64#shape=[2,4,3](empty)
+    # => Numo::Int64#shape=[2,4,3](empty)
 
     f = Numo::DFloat.new(3,4)
-    #=> Numo::DFloat#shape=[3,4](empty)
+    # => Numo::DFloat#shape=[3,4](empty)
 
     f.fill(2)
-    #=> Numo::DFloat#shape=[3,4]
+    # => Numo::DFloat#shape=[3,4]
     # [[2, 2, 2, 2],
     #  [2, 2, 2, 2],
     #  [2, 2, 2, 2]]
 
     x = Numo::NArray.new(5)
-    #=> in `new': allocator undefined for Numo::NArray (TypeError)
+    # => in `new': allocator undefined for Numo::NArray (TypeError)
     #   	from t.rb:9:in `<main>'
 
 */
@@ -444,10 +448,10 @@ na_initialize_copy(VALUE self, VALUE orig)
  *  but for typed NArray subclasses, e.g., DFloat, Int64.
  *  @example
  *    a = Numo::DFloat.zeros(3,5)
- *    => Numo::DFloat#shape=[3,5]
- *    [[0, 0, 0, 0, 0],
- *     [0, 0, 0, 0, 0],
- *     [0, 0, 0, 0, 0]]
+ *    # => Numo::DFloat#shape=[3,5]
+ *    # [[0, 0, 0, 0, 0],
+ *    #  [0, 0, 0, 0, 0],
+ *    #  [0, 0, 0, 0, 0]]
  */
 static VALUE
 na_s_zeros(int argc, VALUE *argv, VALUE klass)
@@ -468,10 +472,10 @@ na_s_zeros(int argc, VALUE *argv, VALUE klass)
  *  but for typed NArray subclasses, e.g., DFloat, Int64.
  *  @example
  *    a = Numo::DFloat.ones(3,5)
- *    => Numo::DFloat#shape=[3,5]
- *    [[1, 1, 1, 1, 1],
- *     [1, 1, 1, 1, 1],
- *     [1, 1, 1, 1, 1]]
+ *    # => Numo::DFloat#shape=[3,5]
+ *    # [[1, 1, 1, 1, 1],
+ *    #  [1, 1, 1, 1, 1],
+ *    #  [1, 1, 1, 1, 1]]
  */
 static VALUE
 na_s_ones(int argc, VALUE *argv, VALUE klass)
@@ -495,8 +499,8 @@ na_s_ones(int argc, VALUE *argv, VALUE klass)
 
   @example
     a = Numo::DFloat.linspace(-5,5,7)
-    => Numo::DFloat#shape=[7]
-    [-5, -3.33333, -1.66667, 0, 1.66667, 3.33333, 5]
+    # => Numo::DFloat#shape=[7]
+    # [-5, -3.33333, -1.66667, 0, 1.66667, 3.33333, 5]
  */
 static VALUE
 na_s_linspace(int argc, VALUE *argv, VALUE klass)
@@ -534,11 +538,12 @@ na_s_linspace(int argc, VALUE *argv, VALUE klass)
 
   @example
     Numo::DFloat.logspace(4,0,5,2)
-    => Numo::DFloat#shape=[5]
-       [16, 8, 4, 2, 1]
+    # => Numo::DFloat#shape=[5]
+    # [16, 8, 4, 2, 1]
+
     Numo::DComplex.logspace(0,1i*Math::PI,5,Math::E)
-    => Numo::DComplex#shape=[5]
-       [1+4.44659e-323i, 0.707107+0.707107i, 6.12323e-17+1i, -0.707107+0.707107i, ...]
+    # => Numo::DComplex#shape=[5]
+    # [1+4.44659e-323i, 0.707107+0.707107i, 6.12323e-17+1i, -0.707107+0.707107i, ...]
  */
 static VALUE
 na_s_logspace(int argc, VALUE *argv, VALUE klass)
@@ -572,10 +577,10 @@ na_s_logspace(int argc, VALUE *argv, VALUE klass)
   @return [Numo::NArray]  created NArray.
   @example
     a = Numo::DFloat.eye(3)
-    => Numo::DFloat#shape=[3,3]
-    [[1, 0, 0],
-     [0, 1, 0],
-     [0, 0, 1]]
+    # => Numo::DFloat#shape=[3,3]
+    # [[1, 0, 0],
+    #  [0, 1, 0],
+    #  [0, 0, 1]]
 */
 static VALUE
 na_s_eye(int argc, VALUE *argv, VALUE klass)
@@ -1583,8 +1588,13 @@ na_get_reduce_flag_from_axes(VALUE na_obj, VALUE axes)
             len = 1;
             step = 0;
             //printf("beg=%d step=%d len=%d\n",beg,step,len);
-        } else if (rb_obj_is_kind_of(v,rb_cRange) ||
-                   rb_obj_is_kind_of(v,na_cStep)) {
+        } else if (rb_obj_is_kind_of(v,rb_cRange)
+#ifdef HAVE_RB_ARITHMETIC_SEQUENCE_EXTRACT
+                   || rb_obj_is_kind_of(v,rb_cArithSeq)
+#else
+                   || rb_obj_is_kind_of(v,rb_cEnumerator)
+#endif
+                   ) {
             nary_step_array_index( v, ndim, &len, &beg, &step );
         } else {
             rb_raise(nary_eDimensionError, "invalid dimension argument %s",
@@ -1870,6 +1880,9 @@ na_equal(VALUE self, volatile VALUE other)
             return Qfalse;
         }
     }
+    if (na1->size == 0) {
+      return Qtrue;
+    }
     vbool = rb_funcall(self, id_eq, 1, other);
     return (rb_funcall(vbool, id_count_false, 0)==INT2FIX(0)) ? Qtrue : Qfalse;
 }
@@ -1896,6 +1909,9 @@ Init_narray()
     rb_require("complex");
     rb_cComplex = rb_const_get(rb_cObject, rb_intern("Complex"));
 #endif
+#ifdef HAVE_RB_ARITHMETIC_SEQUENCE_EXTRACT
+    rb_cArithSeq = rb_path2class("Enumerator::ArithmeticSequence");
+#endif
 
     rb_define_const(cNArray, "VERSION", rb_str_new2(NARRAY_VERSION));
 
@@ -1999,7 +2015,9 @@ Init_narray()
     sym_loop_opt = ID2SYM(rb_intern("loop_opt"));
     sym_init     = ID2SYM(rb_intern("init"));
 
+#ifndef HAVE_RB_ARITHMETIC_SEQUENCE_EXTRACT
     Init_nary_step();
+#endif
     Init_nary_index();
 
     Init_nary_data();

data/ext/numo/narray/ndloop.c

@@ -1,7 +1,7 @@
 /*
   ndloop.c
-  Numerical Array Extension for Ruby
-    (C) Copyright 1999-2017 by Masahiro TANAKA
+  Ruby/Numo::NArray - Numerical Array class for Ruby
+    Copyright (C) 1999-2019 Masahiro TANAKA
 */
 
 #include <ruby.h>

data/ext/numo/narray/numo/intern.h

@@ -1,7 +1,7 @@
 /*
   intern.h
-  Numerical Array Extension for Ruby
-    (C) Copyright 1999-2017 by Masahiro TANAKA
+  Ruby/Numo::NArray - Numerical Array class for Ruby
+    Copyright (C) 1999-2019 Masahiro TANAKA
 */
 #ifndef INTERN_H
 #define INTERN_H
@@ -99,6 +99,7 @@ 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);
+void na_parse_enumerator_step(VALUE enum_obj, VALUE *pstep );
 
 // used in aref, aset
 #define na_get_result_dimension nary_get_result_dimension

data/ext/numo/narray/numo/narray.h

@@ -1,7 +1,7 @@
 /*
   narray.h
-  Numerical Array Extension for Ruby
-    (C) Copyright 1999-2017 by Masahiro TANAKA
+  Ruby/Numo::NArray - Numerical Array class for Ruby
+    Copyright (C) 1999-2019 Masahiro TANAKA
 */
 #ifndef NARRAY_H
 #define NARRAY_H
@@ -13,8 +13,8 @@ extern "C" {
 #endif
 #endif
 
-#define NARRAY_VERSION "0.9.1.4"
-#define NARRAY_VERSION_CODE 914
+#define NARRAY_VERSION "0.9.1.5"
+#define NARRAY_VERSION_CODE 915
 
 #include <math.h>
 #include "numo/compat.h"
@@ -165,10 +165,12 @@ 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
+#ifdef HAVE_RB_ARITHMETIC_SEQUENCE_EXTRACT
+extern VALUE rb_cArithSeq;
+#endif
 
 extern VALUE sym_reduce;
 extern VALUE sym_option;
@@ -234,6 +236,23 @@ typedef struct {
     unsigned int element_stride;
 } narray_type_info_t;
 
+// from ruby/enumerator.c
+struct enumerator {
+    VALUE obj;
+    ID    meth;
+    VALUE args;
+    // use only above in this source
+    VALUE fib;
+    VALUE dst;
+    VALUE lookahead;
+    VALUE feedvalue;
+    VALUE stop_exc;
+    VALUE size;
+    // incompatible below depending on ruby version
+    //VALUE procs;                      // ruby 2.4
+    //rb_enumerator_size_func *size_fn; // ruby 2.1-2.4
+    //VALUE (*size_fn)(ANYARGS);        // ruby 2.0
+};
 
 static inline narray_t *
 na_get_narray_t(VALUE obj)

data/ext/numo/narray/numo/ndloop.h

@@ -1,7 +1,7 @@
 /*
   ndloop.h
-  Numerical Array Extension for Ruby
-    (C) Copyright 1999-2017 by Masahiro TANAKA
+  Ruby/Numo::NArray - Numerical Array class for Ruby
+    Copyright (C) 1999-2019 Masahiro TANAKA
 */
 #ifndef NDLOOP_H
 #define NDLOOP_H

data/ext/numo/narray/numo/template.h

@@ -1,7 +1,7 @@
 /*
   template.h
-  Numerical Array Extension for Ruby
-    (C) Copyright 1999-2017 by Masahiro TANAKA
+  Ruby/Numo::NArray - Numerical Array class for Ruby
+    Copyright (C) 1999-2019 Masahiro TANAKA
 */
 #ifndef TEMPLATE_H
 #define TEMPLATE_H
@@ -112,9 +112,8 @@
         size_t dig = (pos) / NB;                \
         int    bit = (pos) % NB;                \
         ((BIT_DIGIT*)(adr))[dig] =              \
-            (((BIT_DIGIT*)(adr))[dig] & ~(1u<<(bit))) | ((val)<<(bit)); \
+            (((BIT_DIGIT*)(adr))[dig] & ~(1u<<(bit))) | (((val)&1u)<<(bit)); \
     }
-// val -> val&1 ??
 
 #define STORE_BIT_STEP( adr, pos, step, idx, val )\
     {                                           \
@@ -129,9 +128,8 @@
             pos += step;                        \
         }                                       \
         ((BIT_DIGIT*)(adr))[dig] =              \
-            (((BIT_DIGIT*)(adr))[dig] & ~(1u<<(bit))) | ((val)<<(bit)); \
+            (((BIT_DIGIT*)(adr))[dig] & ~(1u<<(bit))) | (((val)&1u)<<(bit)); \
     }
-// val -> val&1 ??
 
 static inline int
 is_aligned(const void *ptr, const size_t alignment)

data/ext/numo/narray/numo/types/complex.h

@@ -1,7 +1,7 @@
 /*
   complex.h
-  Numerical Array Extension for Ruby
-    (C) Copyright 1999-2017 by Masahiro TANAKA
+  Ruby/Numo::NArray - Numerical Array class for Ruby
+    Copyright (C) 1999-2019 Masahiro TANAKA
 */
 
 

data/ext/numo/narray/step.c

@@ -1,7 +1,7 @@
 /*
   step.c
-  Numerical Array Extension for Ruby
-    (C) Copyright 2007,2013 by Masahiro TANAKA
+  Ruby/Numo::NArray - Numerical Array class for Ruby
+    Copyright (C) 2007-2019 Masahiro TANAKA
 */
 #include <ruby.h>
 #include <math.h>
@@ -24,164 +24,10 @@
 #define DBL_EPSILON 2.2204460492503131e-16
 #endif
 
-static ID id_beg, id_end, id_len, id_step, id_excl;
+static ID id_beg, id_end, id_len, id_step;
 
-//#define EXCL(r) RTEST(rb_ivar_get((r), id_excl))
-#define EXCL(r) RTEST(rb_funcall((r), rb_intern("exclude_end?"), 0))
-
-#define SET_EXCL(r,v) rb_ivar_set((r), id_excl, (v) ? Qtrue : Qfalse)
-
-static void
-step_init(
-  VALUE self,
-  VALUE beg,
-  VALUE end,
-  VALUE step,
-  VALUE len,
-  VALUE excl
-)
-{
-    if (RTEST(len)) {
-        if (!(FIXNUM_P(len) || TYPE(len)==T_BIGNUM)) {
-            rb_raise(rb_eArgError, "length must be Integer");
-        }
-        if (RTEST(rb_funcall(len,rb_intern("<"),1,INT2FIX(0)))) {
-            rb_raise(rb_eRangeError,"length must be non negative");
-        }
-    }
-    rb_ivar_set(self, id_beg, beg);
-    rb_ivar_set(self, id_end, end);
-    rb_ivar_set(self, id_len, len);
-    rb_ivar_set(self, id_step, step);
-    SET_EXCL(self, excl);
-}
-
-static VALUE
-nary_step_new2(
-  VALUE range,
-  VALUE step,
-  VALUE len
-)
-{
-    VALUE beg, end, excl;
-    VALUE self = rb_obj_alloc(na_cStep);
-
-    //beg = rb_ivar_get(range, id_beg);
-    beg = rb_funcall(range, id_beg, 0);
-    //end = rb_ivar_get(range, id_end);
-    end = rb_funcall(range, id_end, 0);
-    excl = rb_funcall(range, rb_intern("exclude_end?"), 0);
-
-    step_init(self, beg, end, step, len, excl);
-    return self;
-}
-
-
-/*
- *  call-seq:
- *     Step.new(start, end, step=nil, length=nil)    => step
- *     Step.new(range, step=nil, length=nil)         => step
- *
- *  Constructs a step using three parameters among <i>start</i>,
- *  <i>end</i>, <i>step</i> and <i>length</i>.  <i>start</i>,
- *  <i>end</i> parameters can be replaced with <i>range</i>.  If the
- *  <i>step</i> is omitted (or supplied with nil), then calculated
- *  from <i>length</i> or definded as 1.
- */
-
-static VALUE
-step_initialize( int argc, VALUE *argv, VALUE self )
-{
-    VALUE a, b=Qnil, c=Qnil, d=Qnil, e=Qnil;
-
-    rb_scan_args(argc, argv, "13", &a, &b, &c, &d);
-    /* Selfs are immutable, so that they should be initialized only once. */
-    if (rb_ivar_defined(self, id_beg)) {
-        rb_name_error(rb_intern("initialize"), "`initialize' called twice");
-    }
-    if (rb_obj_is_kind_of(a,rb_cRange)) {
-        if (argc>3) {
-            rb_raise(rb_eArgError, "extra argument");
-        }
-        d = c;
-        c = b;
-        e = rb_funcall(a, rb_intern("exclude_end?"), 0);
-        //b = rb_ivar_get(a, id_end);
-        b = rb_funcall(a, id_end, 0);
-        //a = rb_ivar_get(a, id_beg);
-        a = rb_funcall(a, id_beg, 0);
-    }
-    step_init(self, a, b, c, d, e);
-    return Qnil;
-}
-
-/*
- *  call-seq:
- *     step.begin  => obj
- *     step.first  => obj
- *
- *  Returns the start of <i>step</i>.
- */
-
-static VALUE
-step_first( VALUE self )
-{
-    return rb_ivar_get(self, id_beg);
-}
-
-/*
- *  call-seq:
- *     step.end    => obj
- *     step.last   => obj
- *
- *  Returns the object that defines the end of <i>step</i>.
- */
-
-static VALUE
-step_last( VALUE self )
-{
-    return rb_ivar_get(self, id_end);
-}
-
-/*
- *  call-seq:
- *     step.length  => obj
- *     step.size    => obj
- *
- *  Returns the length of <i>step</i>.
- */
-
-static VALUE
-step_length( VALUE self )
-{
-    return rb_ivar_get(self, id_len);
-}
-
-/*
- *  call-seq:
- *     step.step    => obj
- *
- *  Returns the step of <i>step</i>.
- */
-
-static VALUE
-step_step( VALUE self )
-{
-    return rb_ivar_get(self, id_step);
-}
-
-/*
- *  call-seq:
- *     step.exclude_end?    => true or false
- *
- *  Returns <code>true</code> if <i>step</i> excludes its end value.
- */
-static VALUE
-step_exclude_end_p(VALUE self)
-{
-    return RTEST(rb_ivar_get(self, id_excl)) ? Qtrue : Qfalse;
-}
 
+#define EXCL(r) RTEST(rb_funcall((r), rb_intern("exclude_end?"), 0))
 
 /*
  *  call-seq:
@@ -192,7 +38,7 @@ step_exclude_end_p(VALUE self)
  */
 
 void
-nary_step_array_index(VALUE self, size_t ary_size,
+nary_step_array_index(VALUE obj, size_t ary_size,
                       size_t *plen, ssize_t *pbeg, ssize_t *pstep)
 {
     size_t len;
@@ -200,14 +46,28 @@ nary_step_array_index(VALUE self, size_t ary_size,
     VALUE vbeg, vend, vstep, vlen;
     ssize_t end=ary_size;
 
-    //vbeg = rb_ivar_get(self, id_beg);
-    //vend = rb_ivar_get(self, id_end);
-    vlen = rb_ivar_get(self, id_len);
-    vstep = rb_ivar_get(self, id_step);
-    vbeg = rb_funcall(self, id_beg, 0);
-    vend = rb_funcall(self, id_end, 0);
-    //vlen = rb_funcall(self, id_len, 0);
-    //vstep = rb_funcall(self, id_step, 0);
+#ifdef HAVE_RB_ARITHMETIC_SEQUENCE_EXTRACT
+    rb_arithmetic_sequence_components_t x;
+    rb_arithmetic_sequence_extract(obj, &x);
+
+    vstep = x.step;
+    vbeg = x.begin;
+    vend = x.end;
+#else
+    struct enumerator *e;
+
+    if (rb_obj_is_kind_of(obj, rb_cRange)) {
+        vstep = rb_ivar_get(obj, id_step);
+    } else { // Enumerator
+        na_parse_enumerator_step(obj, &vstep);
+        e = (struct enumerator *)DATA_PTR(obj);
+        obj =  e->obj; // Range
+    }
+
+    vbeg = rb_funcall(obj, id_beg, 0);
+    vend = rb_funcall(obj, id_end, 0);
+#endif
+    vlen = rb_ivar_get(obj, id_len);
 
     if (RTEST(vbeg)) {
         beg = NUM2SSIZET(vbeg);
@@ -237,7 +97,7 @@ nary_step_array_index(VALUE self, size_t ary_size,
                     }
                 } else {
                     if (RTEST(vend)) {
-                        if (EXCL(self)) {
+                        if (EXCL(obj)) {
                             if (step>0) end--;
                             if (step<0) end++;
                         }
@@ -251,7 +111,7 @@ nary_step_array_index(VALUE self, size_t ary_size,
                 step = 1;
                 if (RTEST(vbeg)) {
                     if (RTEST(vend)) {
-                        if (EXCL(self)) {
+                        if (EXCL(obj)) {
                             if (beg<end) end--;
                             if (beg>end) end++;
                         }
@@ -262,7 +122,7 @@ nary_step_array_index(VALUE self, size_t ary_size,
                     }
                 } else {
                     if (RTEST(vend)) {
-                        if (EXCL(self)) {
+                        if (EXCL(obj)) {
                             end--;
                         }
                         beg = end - (len-1);
@@ -286,7 +146,7 @@ nary_step_array_index(VALUE self, size_t ary_size,
             if (!RTEST(vend)) {
                 end = ary_size-1;
             }
-            else if (EXCL(self)) {
+            else if (EXCL(obj)) {
                 end--;
             }
             if (beg<=end) {
@@ -301,7 +161,7 @@ nary_step_array_index(VALUE self, size_t ary_size,
             if (!RTEST(vend)) {
                 end = 0;
             }
-            else if (EXCL(self)) {
+            else if (EXCL(obj)) {
                 end++;
             }
             if (beg>=end) {
@@ -327,25 +187,35 @@ nary_step_array_index(VALUE self, size_t ary_size,
     if (pstep) *pstep = step;
 }
 
-
 void
-nary_step_sequence( VALUE self, size_t *plen, double *pbeg, double *pstep )
+nary_step_sequence( VALUE obj, size_t *plen, double *pbeg, double *pstep )
 {
-    VALUE vbeg, vend, vstep, vlen;
+    VALUE vend, vstep, vlen;
     double dbeg, dend, dstep=1, dsize, err;
     size_t size, n;
 
-    //vbeg = rb_ivar_get(self, id_beg);
-    vbeg = rb_funcall(self, id_beg, 0);
-    dbeg = NUM2DBL(vbeg);
-
-    //vend = rb_ivar_get(self, id_end);
-    vend = rb_funcall(self, id_end, 0);
+#ifdef HAVE_RB_ARITHMETIC_SEQUENCE_EXTRACT
+    rb_arithmetic_sequence_components_t x;
+    rb_arithmetic_sequence_extract(obj, &x);
+
+    vstep = x.step;
+    dbeg = NUM2DBL(x.begin);
+    vend = x.end;
+#else
+    struct enumerator *e;
+
+    if (rb_obj_is_kind_of(obj, rb_cRange)) {
+        vstep = rb_ivar_get(obj, id_step);
+    } else { // Enumerator
+        na_parse_enumerator_step(obj, &vstep);
+        e = (struct enumerator *)DATA_PTR(obj);
+        obj =  e->obj; // Range
+    }
 
-    vlen = rb_ivar_get(self, id_len);
-    vstep = rb_ivar_get(self, id_step);
-    //vlen  = rb_funcall(self, id_len ,0);
-    //vstep = rb_funcall(self, id_step,0);
+    dbeg = NUM2DBL(rb_funcall(obj, id_beg, 0));
+    vend = rb_funcall(obj, id_end, 0);
+#endif
+    vlen = rb_ivar_get(obj, id_len);
 
     if (RTEST(vlen)) {
         size = NUM2SIZET(vlen);
@@ -353,7 +223,7 @@ nary_step_sequence( VALUE self, size_t *plen, double *pbeg, double *pstep )
         if (!RTEST(vstep)) {
             if (RTEST(vend)) {
                 dend = NUM2DBL(vend);
-                if (EXCL(self)) {
+                if (EXCL(obj)) {
                     n = size;
                 } else {
                     n = size-1;
@@ -378,7 +248,7 @@ nary_step_sequence( VALUE self, size_t *plen, double *pbeg, double *pstep )
             err = (fabs(dbeg)+fabs(dend)+fabs(dend-dbeg))/fabs(dstep)*DBL_EPSILON;
             if (err>0.5) err=0.5;
             dsize = (dend-dbeg)/dstep;
-            if (EXCL(self))
+            if (EXCL(obj))
                 dsize -= err;
             else
                 dsize += err;
@@ -398,77 +268,13 @@ nary_step_sequence( VALUE self, size_t *plen, double *pbeg, double *pstep )
     if (pstep) *pstep = dstep;
 }
 
-/*
-static VALUE
-step_each( VALUE self )
-{
-    VALUE  a;
-    double beg, step;
-    size_t i, size;
-
-    a = nary_step_parameters( self, Qnil );
-    beg  = NUM2DBL(RARRAY_PTR(a)[0]);
-    step = NUM2DBL(RARRAY_PTR(a)[1]);
-    size = NUM2SIZET(RARRAY_PTR(a)[2]);
-
-    for (i=0; i<size; i++) {
-        rb_yield(rb_float_new(beg+i*step));
-    }
-    return self;
-}
-*/
-
-static VALUE
-range_with_step( VALUE range, VALUE step )
-{
-    return nary_step_new2( range, step, Qnil );
-}
-
-static VALUE
-range_with_length( VALUE range, VALUE len )
-{
-    return nary_step_new2( range, Qnil, len );
-}
-
-
-static VALUE
-nary_s_step( int argc, VALUE *argv, VALUE mod )
-{
-    VALUE self = rb_obj_alloc(na_cStep);
-    step_initialize(argc, argv, self);
-    return self;
-}
-
-
 void
 Init_nary_step()
 {
-    na_cStep = rb_define_class_under(cNArray, "Step", rb_cObject);
-    rb_include_module(na_cStep, rb_mEnumerable);
-    rb_define_method(na_cStep, "initialize", step_initialize, -1);
-
-    //rb_define_method(na_cStep, "each", step_each, 0);
-
-    rb_define_method(na_cStep, "first", step_first, 0);
-    rb_define_method(na_cStep, "last", step_last, 0);
-    rb_define_method(na_cStep, "begin", step_first, 0);
-    rb_define_method(na_cStep, "end", step_last, 0);
-    rb_define_method(na_cStep, "step", step_step, 0);
-    rb_define_method(na_cStep, "length", step_length, 0);
-    rb_define_method(na_cStep, "size", step_length, 0);
-    rb_define_method(na_cStep, "exclude_end?", step_exclude_end_p, 0);
-    //rb_define_method(na_cStep, "to_s", step_to_s, 0);
-    //rb_define_method(na_cStep, "inspect", step_inspect, 0);
-    //rb_define_method(na_cStep, "parameters", nary_step_parameters, 1);
-
-    rb_define_method(rb_cRange, "%", range_with_step, 1);
-    rb_define_method(rb_cRange, "*", range_with_length, 1);
-
-    rb_define_singleton_method(cNArray, "step", nary_s_step, -1);
+    rb_define_alias(rb_cRange, "%", "step");
 
     id_beg  = rb_intern("begin");
     id_end  = rb_intern("end");
     id_len  = rb_intern("length");
     id_step = rb_intern("step");
-    id_excl = rb_intern("excl");
 }