narray 0.5.9.4 → 0.5.9.5

data/src/speed/add_int.py

@@ -1,12 +0,0 @@
-from mybench import *
-
-a = arrayrange(ARRSZ)
-b = arrayrange(ARRSZ)
-
-print "a.typecode:",a.typecode(),",  a.shape:",a.shape
-print "b.typecode:",b.typecode(),",  b.shape:",b.shape
-print "calculating c = a+b ..."
-
-def bench_body(a=a,b=b): c=a+b
-
-bench_time(bench_body)

data/src/speed/add_int.rb

@@ -1,9 +0,0 @@
-require 'mybench'
-
-a, b = bench_int
-
-print "a = "; p a;
-print "b = "; p b;
-
-print "calculating c = a+b ...\n"
-bench_time{ c=a+b }

data/src/speed/lu.m

@@ -1,14 +0,0 @@
-# mesuring performance of Octave
-
-n = 1000;
-m = linspace(0,n*n-1,n*n);
-m = rem(m, n+1) + 1;
-m = reshape(m,n,n);
-
-printf ("executing %ix%i LU...\n",n,n);
-[t1, u1, s1] = cputime ();
-[l,u,p] = lu(m);
-[t2, u2, s2] = cputime ();
-printf (" Time: %5.2f sec\n", u2 - u1);
-# Time: 20.43 sec, GNU Octave, version 2.0.16 (sparc-sun-solaris2.7).
-exit

data/src/speed/lu.rb

@@ -1,22 +0,0 @@
-require 'mybench'
-
-n = 1000
-m = NArray.float(n,n).indgen!
-m = m % (n+1) + 1
-m = NMatrix.ref(m).transpose
-
-puts 'm ='
-p m
-
-printf "executing %ix%i LU...\n",n,n
-bench_time(1){lu = m.lu}
-exit
-
-puts 'lu='
-p lu
-puts 'y='
-p y
-puts 'm*y='
-p m*y
-
-#system "ps -v"

data/src/speed/mat.m

@@ -1,23 +0,0 @@
-# mesuring performance of Octave
-
-n = 500;
-
-a = linspace(0,n*n-1,n*n);
-a = rem(a, n+1) + 1;
-a = reshape(a,n,n);
-
-b = linspace(0,n*n-1,n*n);
-b = rem(b, n-1) + 1;
-b = reshape(b,n,n);
-
-printf ("executing product of %ix%i Matrix...\n",n,n);
-printf ("c=a*b\n");
-
-[t1, u1, s1] = cputime ();
-c = a*b;
-[t2, u2, s2] = cputime ();
-
-printf (" Time: %f\n", u2 - u1);
-# Time: 7.280000
-# GNU Octave, version 2.0.16 (sparc-sun-solaris2.7).
-exit

data/src/speed/mat.rb

@@ -1,28 +0,0 @@
-require 'mybench'
-
-n = 500
-
-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
-
-c = 0
-
-puts 'a='
-p a
-puts 'b='
-p b
-
-printf "executing %ix%i Matrix product...\n",n,n
-puts 'c=a*b'
-
-bench_time(1) { c=a*b }
-
-# time: 4.21 sec
-
-puts 'c='
-p c

data/src/speed/mul.py

@@ -1,12 +0,0 @@
-from mybench import *
-
-a = bench_array()
-b = bench_array()
-
-print "a.typecode:",a.typecode(),",  a.shape:",a.shape
-print "b.typecode:",b.typecode(),",  b.shape:",b.shape
-print "calculating c = a*b ..."
-
-def bench_body(a=a,b=b): c=a*b
-
-bench_time(bench_body)

data/src/speed/mul.rb

@@ -1,9 +0,0 @@
-require 'mybench'
-
-a, b = bench_float
-
-print "a = "; p a;
-print "b = "; p b;
-print "calculating c = a*b ...\n"
-
-bench_time{ c = a*b }

data/src/speed/mul2.py

@@ -1,15 +0,0 @@
-from mybench import *
-
-a = reshape(arrayrange(1000),(1000,1)).astype(Float64)
-b = reshape(arrayrange(1000),(1,1000)).astype(Float64)
-
-print "a.typecode:",a.typecode(),",  a.shape:",a.shape
-print "b.typecode:",b.typecode(),",  b.shape:",b.shape
-print "calculating  c = a*b ..."
-
-def bench_body(a=a,b=b): c=a*b
-
-bench_time(bench_body)
-
-c = a*b
-print "c.shape:", c.shape,"\n"

data/src/speed/mul2.rb

@@ -1,13 +0,0 @@
-require 'mybench'
-
-a = NArray.float(1000,1).indgen!
-b = NArray.float(1,1000).indgen!
-c = 0
-
-print "a = "; p a;
-print "b = "; p b;
-print "calculating c = a*b ...\n"
-
-bench_time{ c = a*b }
-
-print "c = "; p c;

data/src/speed/mul_comp.py

@@ -1,12 +0,0 @@
-from mybench import *
-
-a = bench_array(Complex64)
-b = bench_array(Complex64)
-
-print "a.typecode:",a.typecode(),",  a.shape:",a.shape
-print "b.typecode:",b.typecode(),",  b.shape:",b.shape
-print "calculating c = a*b ..."
-
-def bench_body(a=a,b=b): c=a*b
-
-bench_time(bench_body)

data/src/speed/mul_comp.rb

@@ -1,9 +0,0 @@
-require 'mybench'
-
-a, b = bench_complex
-
-print "a = "; p a;
-print "b = "; p b;
-print "calculating c = a*b ...\n"
-
-bench_time{ c = a*b }

data/src/speed/mul_int.py

@@ -1,12 +0,0 @@
-from mybench import *
-
-a = arrayrange(ARRSZ) # % 1000
-b = arrayrange(ARRSZ) # / 1000
-
-print "a.typecode:",a.typecode(),",  a.shape:",a.shape
-print "b.typecode:",b.typecode(),",  b.shape:",b.shape
-print "calculating c = a*b ..."
-
-def bench_body(a=a,b=b): c=a*b
-
-bench_time(bench_body)

data/src/speed/mul_int.rb

@@ -1,9 +0,0 @@
-require 'mybench'
-
-a, b = bench_int
-
-print "a = "; p a;
-print "b = "; p b;
-
-print "calculating c = a*b ...\n"
-bench_time{ c=a*b }

data/src/speed/mybench.py

@@ -1,15 +0,0 @@
-from Numeric import *
-import time
-
-REPEAT = 100
-ARRSZ = 1000000
-
-def bench_array(type=Float64):
-  return arrayrange(ARRSZ).astype(type)
-
-def bench_time(func,repeat=REPEAT):
-  start = time.clock()
-  for i in range(repeat):
-    func()
-  stop = time.clock()
-  print "  Time: %7.3f sec" % (stop-start)

data/src/speed/mybench.rb

@@ -1,31 +0,0 @@
-require 'narray'
-
-REPEAT = 100
-ARRSZ = 1_000_0000
-T = (RUBY_VERSION<"1.8.0") ? Time : Process
-
-def bench_array(type=Float)
-  [ NArray.new(type,ARRSZ).indgen!,
-    NArray.new(type,ARRSZ).indgen!  ]
-end
-
-def bench_float
-  bench_array(Float)
-end
-
-def bench_int
-  bench_array(Integer)
-end
-
-def bench_complex
-  bench_array(Complex)
-end
-
-def bench_time(n=REPEAT)
-  t1 = T.times.utime
-   for i in 1..n
-     yield
-   end
-  t2 = T.times.utime 
-  puts " Time: %.2f sec\n\n" % [t2-t1]
-end

data/src/speed/solve.m

@@ -1,18 +0,0 @@
-# mesuring performance of Octave
-n = 500;
-m = linspace(0,n*n-1,n*n);
-m = rem(m, n+1) + 1;
-m = reshape(m,n,n);
-x = reshape(linspace(1,n*n,n*n),n,n);
-
-printf ("solving %ix%i matrix...\n",n,n);
-[t1, u1, s1] = cputime ();
-
-y = m \ x;
-
-[t2, u2, s2] = cputime ();
-printf (" Time: %5.2f sec\n", u2 - u1);
-
-#  Time:  6.87 sec, GNU Octave, version 2.0.16 (sparc-sun-solaris2.7).
-#y
-exit

data/src/speed/solve.py

@@ -1,16 +0,0 @@
-from mybench import *
-from LinearAlgebra import *
-
-n = 500
-
-a = reshape(arrayrange(n*n)%(n+1)+1,(n,n)).astype(Float64)
-
-b = reshape(arrayrange(n*n)+1,(n,n)).astype(Float64)
-#b = arrayrange(n).astype(Float64)+1
-
-print "LU factorize & solve %ix%i matrix ..."%(n,n)
-
-def bench_body(a=a,b=b): c=solve_linear_equations(a,b)
-bench_time(bench_body,1)
-
-# Time:   8.120 sec

data/src/speed/solve.rb

@@ -1,21 +0,0 @@
-require 'mybench'
-
-n = 500
-m = NArray.float(n,n).indgen!
-m = m % (n+1) + 1
-m = NMatrix.ref(m).transpose
-x = NMatrix.float(n,n).indgen!(1).transpose
-#x = NVector.float(n).indgen!(1)
-y = 0
-
-puts 'm ='
-p m
-puts 'x ='
-p x
-
-printf "solving  y=x/m ...\n",n,n
-bench_time(1){ y = x/m }
-#  Time: 5.87 sec
-
-puts 'y='
-p y