numo-narray 0.9.0.1-x64-mingw32 → 0.9.0.3-x64-mingw32

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 6ba9387322df645be933895e871e8cfc22405fbb
-  data.tar.gz: 6aa7706cac725c3daec4dc417189fe364a37bb82
+  metadata.gz: 1214879d114107330abbe27eb03843915badfadc
+  data.tar.gz: '028c8513bd1dc3ac69c9da71400a497ee175fee2'
 SHA512:
-  metadata.gz: bdb2934fa9249928968f35cf6eb437995281f32180d482f36960422bd1b79d2e1cd55f9481c1d7cbeb7592ffd8ea407fa5cc3bd4386c3f276f5994a364d5038b
-  data.tar.gz: 956176ac05733f08f251cb5ff177b20f393ac159d129e7b59daf2e6c9ec2e12f3ee134037b834d6e6cfe539ffbf204e29c0149f9a97b8476b18314a17ad213a1
+  metadata.gz: a922d131ec23008c270adea1dee0054816d26399b9c1f3e96cf175aa0d3c0758922c51db5e835e5e1cb878214c1de1f03aa636605d396a544382cbc73c729ba8
+  data.tar.gz: c629c7a5328f381bf79d31571e7534e9bfb36ad01d949fb635464c9f1053ac54e8c56b814320d68abffa3a2efa14e693087f4cbab4d15c1afa45b35b4c06088b

data/README.md

@@ -1,13 +1,15 @@
 # Numo::NArray - New NArray class library for Ruby/Numo (NUmerical MOdule)
 
+[![Binder](http://mybinder.org/badge.svg)](http://mybinder.org/repo/ruby-numo/narray)
 [![Build Status](https://travis-ci.org/ruby-numo/narray.svg?branch=master)](https://travis-ci.org/ruby-numo/narray)
 
 under development
 
 ## Related Projects
-* [Numo::Linalg](https://github.com/ruby-numo/linalg): Linear Algebra library with [LAPACK](http://www.netlib.org/lapack/).
-* [Numo::GSL](https://github.com/ruby-numo/gsl): Ruby interface for [GSL (GNU Scientific Library)](http://www.gnu.org/software/gsl/).
-* [Numo::FFTE](https://github.com/ruby-numo/ffte): Ruby interface for [FFTE (A Fast Fourier Transform library with radix-2,3,5)](http://www.ffte.jp/).
+* [Numo::Linalg](https://github.com/ruby-numo/linalg) - Linear Algebra library with [LAPACK](http://www.netlib.org/lapack/).
+* [Numo::GSL](https://github.com/ruby-numo/gsl) - Ruby interface for [GSL (GNU Scientific Library)](http://www.gnu.org/software/gsl/).
+* [Numo::FFTE](https://github.com/ruby-numo/ffte) - Ruby interface for [FFTE (A Fast Fourier Transform library with radix-2,3,5)](http://www.ffte.jp/).
+* [Numo::Gnuplot](https://github.com/ruby-numo/gnuplot) - Simple and easy-to-use Gnuplot interface.
 
 ## Installation
 ### Ubuntu, Debian
@@ -16,7 +18,7 @@ apt install -y git ruby gcc ruby-dev rake make
 git clone git://github.com/ruby-numo/narray
 cd narray
 gem build numo-narray.gemspec
-gem install numo-narray-0.9.0.1.gem
+gem install numo-narray-0.9.0.2.gem
 ```
 
 ## Quick start
@@ -40,8 +42,10 @@ An example
 ```
 For more examples, check out this [narray version of 100 numpy exercises](https://github.com/ruby-numo/narray/wiki/100-narray-exercises).
 
-## numo-array status compared to numpy
+## Documentation
 
-https://github.com/ruby-numo/narray/wiki/Numo-vs-numpy
+All documents are primitive.
 
-## [NArray Tentative API Document](http://ruby-numo.github.io/narray/narray/frames.html)
+* [Numo::NArray API Doc](http://ruby-numo.github.io/narray/narray/frames.html)
+* [Numo::NArray概要](https://github.com/ruby-numo/narray/wiki/Numo::NArray%E6%A6%82%E8%A6%81) (in Japanese)
+* [Numo::NArray vs numpy](https://github.com/ruby-numo/narray/wiki/Numo-vs-numpy)

data/Rakefile

@@ -1,23 +1,27 @@
 require "bundler/gem_tasks"
+begin
 require "rake/extensiontask"
 require "rake_compiler_dock"
 require "shellwords"
 
 spec = Bundler::GemHelper.gemspec
 
+cross_platforms = ["x86-mingw32", "x64-mingw32"]
 Rake::ExtensionTask.new("numo/narray", spec) do |ext|
   ext.cross_compile = true
-  ext.cross_platform = ["x86-mingw32", "x64-mingw32"]
+  ext.cross_platform = cross_platforms
+end
+
+pkg_dir = "pkg"
+windows_gem_paths = cross_platforms.collect do |platform|
+  File.join(pkg_dir, "#{spec.full_name}-#{platform}.gem")
 end
 
 namespace :build do
-  pkg_dir = "pkg"
   directory pkg_dir
 
   desc "Build gems for Windows into the pkg directory"
   task :windows => pkg_dir do
-    ruby_versions = "2.1.6:2.2.2:2.3.0"
-
     build_dir = "tmp/windows"
     rm_rf build_dir
     mkdir_p build_dir
@@ -26,7 +30,7 @@ namespace :build do
       ["git", "clone", "file://#{Dir.pwd}/.git", build_dir],
       ["cd", build_dir],
       ["bundle"],
-      ["rake", "cross", "native", "gem", "RUBY_CC_VERSION=#{ruby_versions}"],
+      ["rake", "cross", "native", "gem"],
     ]
     raw_commands = commands.collect do |command|
       Shellwords.join(command)
@@ -39,3 +43,14 @@ namespace :build do
        "#{pkg_dir}/")
   end
 end
+
+namespace :release do
+  task :windows => "build:windows" do
+    windows_gem_paths.each do |path|
+      ruby("-S", "gem", "push", path)
+    end
+  end
+end
+
+rescue LoadError
+end

data/ext/numo/narray/array.c

@@ -1,7 +1,7 @@
 /*
   array.c
   Numerical Array Extension for Ruby
-    (C) Copyright 1999-2011 by Masahiro TANAKA
+    (C) Copyright 1999-2016 by Masahiro TANAKA
 
   This program is free software.
   You can distribute/modify this program
@@ -268,7 +268,7 @@ na_mdai_result(na_mdai_t *mdai, na_compose_t *nc)
         for (i=0; i<ndim; i++) {
             shape[i] = mdai->item[i].shape;
             //printf("shape[%d]=%d\n",i,shape[i]);
-            //rb_ary_push( shape, SIZE2NUM(mdai->item[i].shape) );
+            //rb_ary_push( shape, SIZET2NUM(mdai->item[i].shape) );
         }
 
 	// DataType
@@ -355,7 +355,7 @@ na_s_array_shape(VALUE mod, VALUE ary)
     Data_Get_Struct(vnc, na_compose_t, nc);
     shape = rb_ary_new2(nc->ndim);
     for (i=0; i<nc->ndim; i++) {
-        rb_ary_push( shape, SIZE2NUM(nc->shape[i]) );
+        rb_ary_push( shape, SIZET2NUM(nc->shape[i]) );
     }
     return shape;
 }

data/ext/numo/narray/data.c

@@ -1,7 +1,7 @@
 /*
   data.c
   Numerical Array Extension for Ruby
-    (C) Copyright 1999-2011,2013 by Masahiro TANAKA
+    (C) Copyright 1999-2016 by Masahiro TANAKA
 
   This program is free software.
   You can distribute/modify this program
@@ -76,7 +76,6 @@ na_copy(VALUE self)
     return v;
 }
 
-
 VALUE
 na_store(VALUE self, VALUE src)
 {
@@ -198,65 +197,72 @@ na_transpose_map(VALUE self, int *map)
 VALUE
 na_transpose(int argc, VALUE *argv, VALUE self)
 {
-    int ndim, *map, tmp;
-    int row_major;
-    int i, j, c, r;
-    size_t len;
-    ssize_t beg, step;
-    volatile VALUE v, view;
+    int ndim, *map, *permute;
+    int i, d;
+    bool is_positive, is_negative;
     narray_t *na1;
 
     GetNArray(self,na1);
     ndim = na1->ndim;
-    row_major = TEST_COLUMN_MAJOR( self );
-
+    if (ndim < 2) {
+        if (argc > 0) {
+            rb_raise(rb_eArgError, "unnecessary argument for 1-d array");
+        }
+        return na_make_view(self);
+    }
     map = ALLOCA_N(int,ndim);
-    for (i=0;i<ndim;i++) {
-	map[i] = i;
-    }
-    if (argc==0) {
-	SWAP(map[ndim-1], map[ndim-2], tmp);
-	goto new_object;
-    }
-    if (argc==2) {
-	if (TYPE(argv[0])==T_FIXNUM && TYPE(argv[1])==T_FIXNUM) {
-	    i = FIX2INT(argv[0]);
-	    j = FIX2INT(argv[1]);
-	    if (row_major) {
-		i = ndim-1-i;
-		j = ndim-1-j;
-	    }
-	    SWAP( map[i], map[j], tmp );
-	    goto new_object;
-	}
+    if (argc == 0) {
+        for (i=0; i < ndim; i++) {
+            map[i] = ndim-1-i;
+        }
+        return na_transpose_map(self,map);
     }
-    for (i=argc,c=ndim-1; i;) {
-	v = argv[--i];
-	if (TYPE(v)==T_FIXNUM) {
-	    beg = FIX2INT(v);
-	    len = 1;
-	    step = 0;
-	} else if (rb_obj_is_kind_of(v,rb_cRange) || rb_obj_is_kind_of(v,na_cStep)) {
-            // write me
-	    nary_step_array_index(v, ndim, &len, &beg, &step);
-            //printf("len=%d beg=%d step=%d\n",len,beg,step);
-	}
-	for (j=len; j; ) {
-	    r = beg + step*(--j);
-	    if (row_major) {
-		r = ndim-1-r;
+    // with argument
+    if (argc > ndim) {
+        rb_raise(rb_eArgError, "more arguments than ndim");
+    }
+    for (i=0; i < ndim; i++) {
+        map[i] = i;
+    }
+    permute = ALLOCA_N(int,argc);
+    for (i=0; i < argc; i++) {
+        permute[i] = 0;
+    }
+    is_positive = is_negative = 0;
+    for (i=0; i < argc; i++) {
+	if (TYPE(argv[i]) != T_FIXNUM) {
+            rb_raise(rb_eArgError, "invalid argument");
+        }
+        d = FIX2INT(argv[i]);
+        if (d >= 0) {
+            if (d >= argc) {
+                rb_raise(rb_eArgError, "out of dimension range");
             }
-	    if ( c < 0 ) {
-		rb_raise(rb_eArgError, "too many dims");
+            if (is_negative) {
+                rb_raise(rb_eArgError, "dimension must be non-negative only or negative only");
             }
-	    map[c--] = r;
-	    //printf("r=%d\n",r);
-	}
+            if (permute[d]) {
+                rb_raise(rb_eArgError, "not permutation");
+            }
+            map[i] = d;
+            permute[d] = 1;
+            is_positive = 1;
+        } else {
+            if (d < -argc) {
+                rb_raise(rb_eArgError, "out of dimension range");
+            }
+            if (is_positive) {
+                rb_raise(rb_eArgError, "dimension must be non-negative only or negative only");
+            }
+            if (permute[argc+d]) {
+                rb_raise(rb_eArgError, "not permutation");
+            }
+            map[ndim-argc+i] = ndim+d;
+            permute[argc+d] = 1;
+            is_negative = 1;
+        }
     }
-
- new_object:
-    view = na_transpose_map(self,map);
-    return view;
+    return na_transpose_map(self,map);
 }
 
 //----------------------------------------------------------------------
@@ -304,7 +310,7 @@ na_reshape(int argc, VALUE *argv, VALUE self)
         rb_raise(rb_eArgError, "Total size must be same");
     }
 
-    copy = na_copy(self);
+    copy = rb_funcall(self,rb_intern("copy"),0);
     GetNArray(copy,na);
     //shape_save = NA_SHAPE(na);
     na_setup_shape(na,argc,shape);
@@ -592,7 +598,7 @@ na_diagonal(int argc, VALUE *argv, VALUE self)
     }
 
     if (vofs) {
-        kofs = NUM2SSIZE(vofs);
+        kofs = NUM2SSIZET(vofs);
     } else {
         kofs = 0;
     }
@@ -627,15 +633,15 @@ na_diagonal(int argc, VALUE *argv, VALUE self)
         k0 = 0;
         k1 = kofs;
         if (k1 >= na->shape[ax[1]]) {
-            rb_raise(rb_eArgError,"invalid diagonal offset(%ld) for "
-                     "last dimension size(%ld)",kofs,na->shape[ax[1]]);
+            rb_raise(rb_eArgError,"invalid diagonal offset(%"SZF"d) for "
+                     "last dimension size(%"SZF"d)",kofs,na->shape[ax[1]]);
         }
     } else {
         k0 = -kofs;
         k1 = 0;
         if (k0 >= na->shape[ax[0]]) {
-            rb_raise(rb_eArgError,"invalid diagonal offset(=%ld) for "
-                     "last-1 dimension size(%ld)",kofs,na->shape[ax[0]]);
+            rb_raise(rb_eArgError,"invalid diagonal offset(=%"SZF"d) for "
+                     "last-1 dimension size(%"SZF"d)",kofs,na->shape[ax[0]]);
         }
     }
 
@@ -749,7 +755,7 @@ na_diagonal(int argc, VALUE *argv, VALUE self)
 static VALUE
 na_new_dimension_for_dot(VALUE self, int pos, int len, bool transpose)
 {
-    int i, k, nd;
+    int i, k, l, nd;
     size_t  j;
     size_t *idx1, *idx2;
     size_t *shape;
@@ -775,27 +781,25 @@ na_new_dimension_for_dot(VALUE self, int pos, int len, bool transpose)
     }
     nd += len;
     shape = ALLOCA_N(size_t,nd);
-    i = k = 0;
-    while (i < nd) {
-        if (i == pos) {
-            for (; len; len--) {
-                shape[i++] = 1;
-            }
-            pos = -1; // new axis done
-        } else {
-            shape[i++] = na->shape[k++];
-        }
-    }
-
-    na_setup_shape((narray_t*)na2, nd, shape);
     na2->stridx = ALLOC_N(stridx_t,nd);
 
     switch(na->type) {
     case NARRAY_DATA_T:
     case NARRAY_FILEMAP_T:
+        i = k = 0;
+        while (i < nd) {
+            if (i == pos && len > 0) {
+                for (l=0; l<len; l++) {
+                    shape[i++] = 1;
+                }
+            } else {
+                shape[i++] = na->shape[k++];
+            }
+        }
+        na_setup_shape((narray_t*)na2, nd, shape);
         stride = na_get_elmsz(self);
         for (i=nd; i--;) {
-            SDX_SET_STRIDE(na2->stridx[i],stride);
+            SDX_SET_STRIDE(na2->stridx[i], stride);
             stride *= shape[i];
         }
         na2->offset = 0;
@@ -803,18 +807,35 @@ na_new_dimension_for_dot(VALUE self, int pos, int len, bool transpose)
         break;
     case NARRAY_VIEW_T:
         GetNArrayView(self, na1);
-        for (i=0; i<nd; i++) {
-            if (SDX_IS_INDEX(na1->stridx[i])) {
-                idx1 = SDX_GET_INDEX(na1->stridx[i]);
-                idx2 = ALLOC_N(size_t,na1->base.shape[i]);
-                for (j=0; j<na1->base.shape[i]; j++) {
-                    idx2[j] = idx1[j];
+        i = k = 0;
+        while (i < nd) {
+            if (i == pos && len > 0) {
+                if (SDX_IS_INDEX(na1->stridx[k])) {
+                    stride = SDX_GET_INDEX(na1->stridx[k])[0];
+                } else {
+                    stride = SDX_GET_STRIDE(na1->stridx[k]);
+                }
+                for (l=0; l<len; l++) {
+                    shape[i] = 1;
+                    SDX_SET_STRIDE(na2->stridx[i], stride);
+                    i++;
                 }
-                SDX_SET_INDEX(na2->stridx[i],idx2);
             } else {
-                na2->stridx[i] = na1->stridx[i];
+                shape[i] = na1->base.shape[k];
+                if (SDX_IS_INDEX(na1->stridx[k])) {
+                    idx1 = SDX_GET_INDEX(na1->stridx[k]);
+                    idx2 = ALLOC_N(size_t,na1->base.shape[k]);
+                    for (j=0; j<na1->base.shape[k]; j++) {
+                        idx2[j] = idx1[j];
+                    }
+                    SDX_SET_INDEX(na2->stridx[i], idx2);
+                } else {
+                    na2->stridx[i] = na1->stridx[k];
+                }
+                i++; k++;
             }
         }
+        na_setup_shape((narray_t*)na2, nd, shape);
         na2->offset = na1->offset;
         na2->data = na1->data;
         break;

data/ext/numo/narray/depend.erb

@@ -17,17 +17,19 @@ 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)
+<%
+   type_c = ["types/bit.c"]
+   type_rb = Dir.glob("#{__dir__}/gen/def/*.rb")
+   type_rb.each do |s|
+     type_c << c = "types/"+File.basename(s,".rb")+".c"
+%>
+<%=c%>: <%=s%> $(DEPENDS)
 	$(MAKEDIRS) $(@D) types
-	ruby $(COGEN_DTYPE) <%=s%> > $@
+	ruby $(COGEN) -l -o $@ $(DTYPE) <%=s%>
 <% end %>
 
-types/bit.c: <%= __dir__ %>/gen/bit.erb.c $(DEPENDS)
-	$(MAKEDIRS) $(@D) types
-	ruby $(COGEN) -l <%= __dir__ %>/gen/bit.erb.c > $@
+src : <%= type_c.join(" ") %>
 
 CLEANOBJS = *.o */*.o *.bak types/*.c

data/ext/numo/narray/extconf.rb

@@ -1,9 +1,5 @@
 require 'rbconfig.rb'
-
-#RbConfig::MAKEFILE_CONFIG["optflags"] = "-g3 -gdwarf-2"
-
 require 'mkmf'
-
 require "erb"
 
 if RUBY_VERSION < "2.0.0"
@@ -11,12 +7,15 @@ if RUBY_VERSION < "2.0.0"
   exit(1)
 end
 
+rm_f 'numo/extconf.h'
+
 #$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)'] }
+$INSTALLFILES << ['numo/extconf.h','$(archdir)']
 if /cygwin|mingw/ =~ RUBY_PLATFORM
   $INSTALLFILES << ['libnarray.a', '$(archdir)']
 end
@@ -48,20 +47,6 @@ 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
@@ -91,18 +76,14 @@ 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")
+have_var("rb_cFixnum")
 
 $objs = srcs.collect{|i| i+".o"}
 
-create_header
+create_header('numo/extconf.h')
 
 depend_path = File.join(__dir__, "depend")
 File.open(depend_path, "w") do |depend|