narray 0.5.9.9 → 0.6.0.0

data/src/bench/bench.rb

@@ -1,59 +0,0 @@
-require 'narray'
-T = (RUBY_VERSION<"1.8.0") ? Time : Process
-
-TYPE   = ARGV[0]
-OP     = ARGV[1]
-ARRSZ  = Integer(ARGV[2])
-REPEAT = Integer(ARGV[3])
-
-def bench_array(type=Float)
-  [ NArray.new(type,ARRSZ).indgen!,
-    NArray.new(type,ARRSZ).indgen!  ]
-end
-
-def bench_time(n=REPEAT)
-  t1 = T.times.utime
-   for i in 1..n
-     yield
-   end
-  t = T.times.utime - t1
-  printf "Ruby NArray type=%s size=%d op=%s repeat=%d  Time: %.2f sec\n",
-    TYPE,ARRSZ,OP,REPEAT,t
-end
-
-n = ARRSZ
-
-case TYPE
-when "float"
-  a,b = bench_array(Float)
-when "int"
-  a,b = bench_array(Integer)
-when "complex"
-  a,b = bench_array(Complex)
-when "float_cross"
-  a = NArray.float(n,1).indgen!
-  b = NArray.float(1,n).indgen!
-when "float_matrix"
-  a = NArray.float(n,n).indgen!
-  a = a % (n+1) + 1
-  a = NMatrix.ref(a)#.transpose
-  b = NArray.float(n,n).indgen!
-  b = b % (n-1) + 1
-  b = NMatrix.ref(b)#.transpose
-when "float_solve"
-  a = NMatrix.float(n,n).indgen!(1).transpose
-  b = NArray.float(n,n).indgen!
-  b = b % (n+1) + 1
-  b = NMatrix.ref(b).transpose
-end
-
-c = 0
-
-case OP
-when "add"
-  bench_time{ c = a+b }
-when "mul"
-  bench_time{ c = a*b }
-when "solve"
-  bench_time{ c = a/b }
-end

data/src/bench/dummy.py

@@ -1,13 +0,0 @@
-import sys
-
-MODULE = sys.argv[1]
-
-if MODULE=="numeric":
-    from Numeric import *
-    from LinearAlgebra import *
-elif MODULE=="numarray":
-    from numarray import *
-    from LinearAlgebra import *
-elif MODULE=="numpy":
-    from numpy import *
-    from numpy.linalg import solve

data/src/test/statistics.rb

@@ -1,22 +0,0 @@
-require "narray"
-include NMath
-
-x = NArray[65, 63, 67, 64, 68, 62, 70, 66, 68, 67, 69, 71]
-y = NArray[68, 66, 68, 65, 69, 66, 68, 65, 71, 67, 68, 70]
-
-def test str, x, y=nil
-  print str," #=> "
-  p eval(str)
-end
-
-test "x",x,y
-test "y",x,y
-test "covariance(x,y)",x,y
-
-a = covariance(x,y)
-
-test "x.stddev",x
-
-test "x.sort",x
-test "x.median",x
-test "(x+y.newrank!(0)).median(0)",x,y

data/src/test/testarray.rb

@@ -1,20 +0,0 @@
-require 'narray'
-
-a = NArray.float(3,3).indgen
-
-p NArray[a,[100,101]]
-
-b = NArray[ [ [ 0.0, 1.0, 2.0 ], 
-              [ 3.0, 4.0, 5.0 ], 
-              [ 6.0, 7.0, 8.0 ] ],
-            [100,101] ]
-p b
-
-a = NArray.float(2,2).indgen
-
-b = NArray[ a,[a] ]
-p b
-
-b = NArray[ [ 0.0, 1.0 ], 
-            [ [ 0.0, 1.0 ] ] ]
-p b

data/src/test/testbit.rb

@@ -1,27 +0,0 @@
-require "narray"
-
-a = NArray.byte(10).indgen!
-
-def test a, str
-  print str," #=>\n"
-  p eval(str)
-end
-
-test a, "a"
-test a, "a & 1"
-test a, "a & 2"
-test a, "a & -1"
-
-test a, "a | 1"
-test a, "a | 2"
-test a, "a | -1"
-
-test a, "a ^ 1"
-test a, "a ^ 2"
-test a, "a ^ -1"
-
-test a, "~a"
-
-a = NArray.int(10).indgen!
-test a, "a"
-test a, "~a"

data/src/test/testcast.rb

@@ -1,14 +0,0 @@
-require 'narray'
-
-def test a,b
-  print a," #=> "
-  p b
-end
-
-a = NArray.int(3,3).indgen!
-
-test "a",a
-test "a+1.5",a+1.5
-test "1.5+a",1.5+a
-test "a+NArray[1.2,3.4,5.6]",a+NArray[1.2,3.4,5.6]
-test "a+NArray[Complex(0.5,1.5)]",a+NArray[Complex(0.5,1.5)]

data/src/test/testcomplex.rb

@@ -1,35 +0,0 @@
-require 'narray'
-
-def testop(a,b)
-  print "a = "; p a
-  print "b = "; p b
-  print "a+b = "; p a+b
-  print "a-b = "; p a-b
-  print "a*b = "; p a*b
-  print "a/b = "; p a/b
-  print "a**b = "; p a**b
-end
-
-a = NArray.complex(4,1).indgen!.sbt!(-1) + 0.25.im
-b = NArray.complex(1,3).indgen!.add!(-0.5).mul!(0.5.im)
-testop(a,b)
-
-# compare a/b with real-number operation
-# a = NArray(4,1).indgen!.sbt!(-1)
-# b = NArray(1,3).fill!(1)
-# c = b*b + b*b
-# p ( (a*b + 0*b)/c )
-# p ( (0*b - a*b)/c )
-
-def testimag(a)
-  print "a.real = "
-  p a.real
-  print "a.imag = "
-  p a.imag
-  print "a.angle = "
-  p a.angle
-  print "a.conj = "
-  p a.conj
-end
-
-testimag a

data/src/test/testfftw.rb

@@ -1,16 +0,0 @@
-require 'narray'
-
-x = NArray.complex(1024,1024).indgen!
-
-p x
-puts 'executing fftw ...'
-t1 = Time.times.utime
-y = FFTW.fftw( x, 1 )
-t2 = Time.times.utime
-print "time: ",t2 - t1,"\n"
-p y
-
-exit
-
-x = NArray.complex(128,128).indgen!
-10000.times{ x = FFTW.fftw( x, 1 ) }

data/src/test/testindex.rb

@@ -1,11 +0,0 @@
-require 'narray'
-
-a = NArray.float(3,3).indgen!
-
-print "a #=> "; p a
-print "a[1, -1] #=> "; p a[1, -1]
-print "a[1, 2..0] #=> "; p a[1, 2..0]
-print "a[1...2, 0..1] #=> "; p a[1...2, 0..1]  # without rank-reduce
-print "a[true, 0..1] #=> "; p a[true, 0..1]
-print "a[0..5] #=> "; p a[0..5]
-

data/src/test/testindexary.rb

@@ -1,26 +0,0 @@
-require 'narray'
-
-a = NArray.new(NArray::DFLOAT,5,5).indgen!
-print "a #=> "; p a
-
-idx = [2,0,1]
-print "\nidx #=> "
- p idx
-
-print "\na[idx,idx]  # Index Array\n #=> "
- p a[idx,idx]
-
-idx = NArray.to_na(idx)
-idx += 10
-print "\nidx #=> "
- p idx
-
-print "\na[idx] #=> "
- p a[idx]
-
-print "\na[1,[]]  # Empty Array\n #=> "
- p a[1,[]]
-
-print "\nFollowing will fail...\n\n"
-print "a[idx,0] #=> \n"
-p a[idx,0]

data/src/test/testindexset.rb

@@ -1,55 +0,0 @@
-require 'narray'
-
-a = NArray.float(3,3).indgen!
-
-print "a #=> "; p a
-
-print "\na[[]] = []  # Do nothing \na #=> ";
-a[[]] = []; p a
-
-print "\na[1,1] = 0  # 1-element replace\na #=> ";
-b=a.dup; b[1,1] = 0
-p b
-
-print "\na[0..-1, 1] = 0  # Range replace\na #=> ";
-b=a.dup; b[0..-1, 1] = 0
-p b
-
-print "\na[1,2..0] = [100,101,102]  # Array replace\na #=> ";
-b=a.dup; b[1,2..0] = [100,101,102]
-p b
-
-print "\na[0,1] = [100,101,102]  # Specifing Starting point \na #=> ";
-b=a.dup; b[0,1] = [100,101,102]
-p b
-
-print "\na[1,0] = [[100],[101],[102]]  # Specifing Starting point \na #=> ";
-b=a.dup; b[1,0] = [[100],[101],[102]]
-p b
-
-print "\na[true,1] = [100,101,102]  # `true' means entire range\na #=> ";
-b=a.dup; b[true,1] = [100,101,102]
-p b
-
-print "\na[true,true] = [[100,101,102]] \na #=> ";
-b=a.dup; b[true,true] = [[100,101,102]]
-p b
-
-
-print "\nFollowing will fail ...\n"
-
-print "\na[true,1] = [[100,101,102]] \na #=> ";
-b=a.dup; b[true,1] = [[100,101,102]]
-p b
-
-print "\na[true,1] = [[100],[101],[102]] \na #=> ";
-b=a.dup; b[true,1] = [[100],[101],[102]]
-p b
-
-print "\na[true,true] = [100,101,102] \na #=> ";
-b=a.dup; b[true,true] = [100,101,102]
-p b
-
-print "\na[1,0] = [100,101,102] \na #=> ";
-b=a.dup; b[1,0] = [100,101,102]
-p b

data/src/test/testmask.rb

@@ -1,40 +0,0 @@
-require "narray"
-
-a = NArray.byte(10)
-a[2..4] = 1
-p a.class, a.count_true, a.count_false
-
-begin
-   a = NArray.float(10)
-   a[2..4] = 1
-   p a.class, a.count_true, a.count_false
-rescue
-   print a.class," -- Exception raised as expected. The message was: ", $!,"\n"
-end
-
-#-------------------
-print "\n--- test masking (float) ---\n"
-
-a = NArray.float(5,3).indgen!
-b = (a-2)*2
-c = a.lt(b)
-p c, c.typecode
-
-p a, b, a.mask( c ), a[c]
-
-#a[c] = ( NArray.int(c.length).indgen!+100 )
-#p a
-
-a[c] = 10000
-p a
-
-#-------------------
-print "\n--- test masking (complex) ---\n"
-p a = NArray.complex(5).indgen! + Complex::I
-
-m = NArray.byte(5)
-m[true] = [0,0,1,1,0]
-
-p a[m]
-a[m] = 100.0
-p a

data/src/test/testmath.rb

@@ -1,48 +0,0 @@
-require 'narray'
-include NMath
-
-def pr x
-  x.each{|i|
-    if i.kind_of?(Complex)
-      printf("%.3f%+.3fi ",i.real,i.image)
-    else
-      printf("%.3f ",i)
-    end      
-  }
-  print "\n"
-end
-
-def testmath(x)
-  print "x = "
-  pr x
-  print "sqrt(x) = "
-  pr sqrt(x)
-  print "sin(x) = "
-  pr sin(x)
-  print "cos(x) = "
-  pr cos(x)
-  print "tan(x) = "
-  pr tan(x)
-  print "sinh(x) = "
-  pr sinh(x)
-  print "cosh(x) = "
-  pr cosh(x)
-  print "tanh(x) = "
-  pr tanh(x)
-  print "exp(x) = "
-  pr exp(x)
-  print "log(x) = "
-  pr log(x)
-  print "log10(x) = "
-  pr log10(x)
-  print "atan(x) = "
-  pr atan(x)
-  print "atan(tan(x)) = "
-  pr atan(tan(x))
-end
-
-testmath NArray.sfloat(6).indgen.div!(2)
-testmath NArray.float(6).indgen.div!(2)
-
-testmath NArray.scomplex(6).indgen.div!(2)-2 - 1.im
-testmath NArray.complex(6).indgen!/5-0.5# - 0.3.im

data/src/test/testmath2.rb

@@ -1,46 +0,0 @@
-require "narray"
-include NMath
-
-def testm x,name
-  f = "a#{name}(#{name}(x))"
-  print "\n### #{f} ###\n"
-  y = eval(f)
-  p y
-  #d = x-y
-  #p d.abs
-end
-
-PI=Math::PI
-
-i = NArray.complex(1)
-i.imag=1
-
-n=7
-x = NArray.complex(n).random!(1)
-x.imag= NArray.float(n).random!(1)
-
-p x
-testm x,"sin"
-testm x,"cos"
-testm x,"tan"
-testm x,"sinh"
-testm x,"cosh"
-testm x,"tanh"
-testm x,"sec"
-testm x,"sech"
-
-p cot(1)
-
-exit
-
-a= -i * log( sqrt(1-x**2)*i + x )
-b= -i * log( sqrt(x**2-1) + x )
-p a
-p b
-p a-b
-
-a= -i * log( sqrt(1-x**2) + x*i )
-b= -i * log( (-sqrt(x**2-1) + x)*i )
-p a
-p b
-p a-b