nmatrix 0.0.1

data/spec/nmatrix_yale_spec.rb

@@ -0,0 +1,239 @@
+# = NMatrix
+#
+# A linear algebra library for scientific computation in Ruby.
+# NMatrix is part of SciRuby.
+#
+# NMatrix was originally inspired by and derived from NArray, by
+# Masahiro Tanaka: http://narray.rubyforge.org
+#
+# == Copyright Information
+#
+# SciRuby is Copyright (c) 2010 - 2012, Ruby Science Foundation
+# NMatrix is Copyright (c) 2012, Ruby Science Foundation
+#
+# Please see LICENSE.txt for additional copyright notices.
+#
+# == Contributing
+#
+# By contributing source code to SciRuby, you agree to be bound by
+# our Contributor Agreement:
+#
+# * https://github.com/SciRuby/sciruby/wiki/Contributor-Agreement
+#
+# == nmatrix_yale_spec.rb
+#
+# Basic tests for NMatrix's Yale storage type.
+#
+require "./lib/nmatrix"
+
+describe NMatrix do
+  it "correctly compares two empty yale matrices" do
+    n = NMatrix.new(:yale, [4,4], :float64)
+    m = NMatrix.new(:yale, [4,4], :float64)
+    n.should.eql? m
+  end
+
+  it "correctly compares two yale matrices following basic assignments" do
+    n = NMatrix.new(:yale, [2,2], :float64)
+    m = NMatrix.new(:yale, [2,2], :float64)
+    m[0,0] = 1
+    m[0,1] = 1
+    n.should_not.eql? m
+    n[0,0] = 1
+    n.should_not.eql? m
+    n[0,1] = 1
+    n.should.eql? m
+  end
+
+  it "correctly compares two yale matrices following elementwise operations" do
+    n = NMatrix.new(:yale, [2,2], :float64)
+    n[0,1] = 1
+    m = NMatrix.new(:yale, [2,2], :float64)
+    m[0,1] = -1
+    r = NMatrix.new(:yale, [2,2], :float64)
+    r[0,1] = 0
+    (n+m).should.eql? r
+  end
+
+  it "correctly sets diagonal values in yale" do
+    n = NMatrix.new(:yale, [2,3], :float64)
+    n[1,1] = 0.1
+    n[0,0] = 0.2
+    n.__yale_d__.should == [0.2, 0.1]
+  end
+
+  it "does not resize yale until necessary" do
+    n = NMatrix.new(:yale, [2,3], :float64)
+    n.__yale_size__.should == 3
+    n.capacity.should == 5
+    n[0,0] = 0.1
+    n[0,1] = 0.2
+    n[1,0] = 0.3
+    n.__yale_size__.should == 5
+    n.capacity.should == 5
+  end
+
+
+  it "correctly sets when not resizing (yale)" do
+    n = NMatrix.new(:yale, [2,3], :float64)
+    n[0,0] = 0.1
+    n[0,1] = 0.2
+    n[1,0] = 0.3
+    n.__yale_a__ == [0.1, 0.0, 0.0, 0.2, 0.3]
+    n.__yale_ija__ == [3,4,5,1,0]
+  end
+
+  it "correctly sets when resizing (yale)" do
+    n = NMatrix.new(:yale, [2,3], :float64)
+    n[0,0] = 0.01
+    n[1,1] = 0.1
+    n[0,1] = 0.2
+    n[1,0] = 0.3
+    n[1,2] = 0.4
+    n.__yale_d__.should == [0.01, 0.1]
+    n.__yale_ia__.should == [3,4,6]
+    n.__yale_ja__.should == [1,0,2,nil]
+    n.__yale_lu__.should == [0.2, 0.3, 0.4, nil]
+  end
+
+  it "correctly sets values within rows of yale" do
+    n = NMatrix.new(:yale, [3,20], :float64)
+    n[2,1]   = 1.0
+    n[2,0]   = 1.5
+    n[2,15]  = 2.0
+    n.__yale_lu__.should == [1.5, 1.0, 2.0]
+    n.__yale_ja__.should == [0, 1, 15]
+  end
+
+  it "correctly gets values within rows of yale" do
+    n = NMatrix.new(:yale, [3,20], :float64)
+    n[2,1]   = 1.0
+    n[2,0]   = 1.5
+    n[2,15]  = 2.0
+    n[2,1].should == 1.0
+    n[2,0].should == 1.5
+    n[2,15].should == 2.0
+  end
+
+  it "correctly sets values within large rows of yale" do
+    n = NMatrix.new(:yale, [10,300], :float64)
+    n[5,1]   = 1.0
+    n[5,0]   = 1.5
+    n[5,15]  = 2.0
+    n[5,291] = 3.0
+    n[5,292] = 4.0
+    n[5,289] = 5.0
+    n[5,290] = 6.0
+    n[5,293] = 2.0
+    n[5,299] = 7.0
+    n[5,100] = 8.0
+    n.__yale_lu__.should == [1.5, 1.0, 2.0, 8.0, 5.0, 6.0, 3.0, 4.0, 2.0, 7.0]
+    n.__yale_ja__.should == [0,   1,   15,  100, 289, 290, 291, 292, 293, 299]
+  end
+
+  it "correctly gets values within large rows of yale" do
+    n = NMatrix.new(:yale, [10,300], :float64)
+    n[5,1]   = 1.0
+    n[5,0]   = 1.5
+    n[5,15]  = 2.0
+    n[5,291] = 3.0
+    n[5,292] = 4.0
+    n[5,289] = 5.0
+    n[5,290] = 6.0
+    n[5,293] = 2.0
+    n[5,299] = 7.0
+    n[5,100] = 8.0
+
+    n.__yale_ja__.each_index do |idx|
+      j = n.__yale_ja__[idx]
+      n[5,j].should == n.__yale_lu__[idx]
+    end
+  end
+
+  it "correctly dots two identical yale matrices" do
+    a = NMatrix.new(:yale, 4, :float64)
+    a[0,1] = 4.0
+    a[1,2] = 1.0
+    a[1,3] = 1.0
+    a[3,1] = 2.0
+
+    b = a.dup
+    c = a.dot b
+
+    c[0,0].should == 0.0
+    c[0,1].should == 0.0
+    c[0,2].should == 4.0
+    c[0,3].should == 4.0
+    c[1,0].should == 0.0
+    c[1,1].should == 2.0
+    c[1,2].should == 0.0
+    c[1,3].should == 0.0
+    c[2,0].should == 0.0
+    c[2,1].should == 0.0
+    c[2,2].should == 0.0
+    c[2,3].should == 0.0
+    c[3,0].should == 0.0
+    c[3,1].should == 0.0
+    c[3,2].should == 2.0
+    c[3,3].should == 2.0
+  end
+
+  it "correctly dots two identical yale matrices where a positive and negative partial sum cancel on the diagonal" do
+    a = NMatrix.new(:yale, 4, :float64)
+    a[0,0] = 1.0
+    a[0,1] = 4.0
+    a[1,2] = 2.0
+    a[1,3] = -4.0
+    a[3,1] = 4.0
+    a[3,3] = 4.0
+
+    b = a.dup
+    c = a.dot b
+
+    #c[0,0].should == 1.0
+    #c[0,1].should == 4.0
+    #c[0,2].should == 8.0
+    #c[0,3].should == -16.0
+    #c[1,0].should == 0.0
+    #c[1,1].should == -16.0
+    #c[1,2].should == 0.0
+    #c[1,3].should == -16.0
+    #c[2,0].should == 0.0
+    #c[2,1].should == 0.0
+    #c[2,2].should == 0.0
+    #c[2,3].should == 0.0
+    #c[3,0].should == 0.0
+    #c[3,1].should == 0.0
+    #c[3,2].should == 8.0
+    #c[3,3].should == 0.0 # this is the positive and negative partial sum cancel
+
+    c.__yale_ija__.reject { |i| i.nil? }.should == [5,8,9,9,11,1,2,3,3,1,2]
+    c.__yale_a__.reject { |i| i.nil? }.should == [1.0, -16.0, 0.0, 0.0, 0.0, 4.0, 8.0, -16.0, -16.0, 16.0, 8.0]
+
+  end
+
+  it "correctly transposes yale" do
+    a = NMatrix.new(:yale, 4, :float64)
+    a[0,0] = 1.0
+    a[0,1] = 4.0
+    a[1,2] = 2.0
+    a[1,3] = -4.0
+    a[3,1] = 5.0
+    a[3,3] = 6.0
+    b = a.transpose
+
+    b[0,0].should == 1.0
+    b[1,0].should == 4.0
+    b[2,0].should == 0.0
+    b[3,0].should == 0.0
+    b[0,1].should == 0.0
+    b[1,1].should == 0.0
+    b[2,1].should == 2.0
+    b[3,1].should == -4.0
+    b[0,3].should == 0.0
+    b[1,3].should == 5.0
+    b[2,3].should == 0.0
+    b[3,3].should == 6.0
+  end
+
+end

data/spec/nvector_spec.rb

@@ -0,0 +1,43 @@
+# = NMatrix
+#
+# A linear algebra library for scientific computation in Ruby.
+# NMatrix is part of SciRuby.
+#
+# NMatrix was originally inspired by and derived from NArray, by
+# Masahiro Tanaka: http://narray.rubyforge.org
+#
+# == Copyright Information
+#
+# SciRuby is Copyright (c) 2010 - 2012, Ruby Science Foundation
+# NMatrix is Copyright (c) 2012, Ruby Science Foundation
+#
+# Please see LICENSE.txt for additional copyright notices.
+#
+# == Contributing
+#
+# By contributing source code to SciRuby, you agree to be bound by
+# our Contributor Agreement:
+#
+# * https://github.com/SciRuby/sciruby/wiki/Contributor-Agreement
+#
+# == nvector_spec.rb
+#
+# Basic tests for NVector.
+#
+
+require "./lib/nmatrix"
+
+describe NVector do
+  it "correctly initializes" do
+    v = NVector.new 5, :float64
+    v.shape[0].should == 5
+    v.shape[1].should == 1
+  end
+
+  it "permits setting and getting contents" do
+    v = NVector.new 5, :float64
+    v[0] = 1.555
+    v[0].should == 1.555
+  end
+
+end

data/spec/syntax_tree_spec.rb

@@ -0,0 +1,46 @@
+# = NMatrix
+#
+# A linear algebra library for scientific computation in Ruby.
+# NMatrix is part of SciRuby.
+#
+# NMatrix was originally inspired by and derived from NArray, by
+# Masahiro Tanaka: http://narray.rubyforge.org
+#
+# == Copyright Information
+#
+# SciRuby is Copyright (c) 2010 - 2012, Ruby Science Foundation
+# NMatrix is Copyright (c) 2012, Ruby Science Foundation
+#
+# Please see LICENSE.txt for additional copyright notices.
+#
+# == Contributing
+#
+# By contributing source code to SciRuby, you agree to be bound by
+# our Contributor Agreement:
+#
+# * https://github.com/SciRuby/sciruby/wiki/Contributor-Agreement
+#
+# == syntax_tree_spec.rb
+#
+# Tests for the source code generator (syntax_tree.rb in the
+# ext/nmatrix/generator dir).
+#
+require "./ext/nmatrix/generator/syntax_tree.rb"
+
+describe SyntaxTree do
+  #it "correctly handles rational == 0" do
+  #  SyntaxTree.parse("x == 0").operate(:rational,:r128).should == "x.n == 0"
+  #end
+  #
+  #it "correctly handles rational = 0" do
+  #  SyntaxTree.parse("x == 0").operate(:rational,:r128).split("\n").should == ["x.n = 0", "x.d = 1;"]
+  #end
+  #
+  #it "correctly handles complex == 0" do
+  #  SyntaxTree.parse("x == 0").operate(:complex,:c128).should == "x.r == 0 && x.i == 0"
+  #end
+  #
+  #it "correctly handles complex = 0" do
+  #  SyntaxTree.parse("x = 0").operate(:complex,:c128).split("\n").should == ["x.r = 0", "x.i = 1;"]
+  #end
+end

metadata

@@ -0,0 +1,150 @@
+--- !ruby/object:Gem::Specification
+name: nmatrix
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+  prerelease: 
+platform: ruby
+authors:
+- John Woods
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2012-04-10 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: rdoc
+  requirement: &81217480 !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: '3.10'
+  type: :development
+  prerelease: false
+  version_requirements: *81217480
+- !ruby/object:Gem::Dependency
+  name: rake-compiler
+  requirement: &81217250 !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: '0.7'
+  type: :development
+  prerelease: false
+  version_requirements: *81217250
+- !ruby/object:Gem::Dependency
+  name: hoe
+  requirement: &81217030 !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: '3.0'
+  type: :development
+  prerelease: false
+  version_requirements: *81217030
+description: ! 'NMatrix is an experimental linear algebra library for Ruby, written
+  mostly in C. It can be used with or without SciRuby, but is part of the SciRuby
+  project.
+
+
+  NMatrix was inspired by and based heavily upon {NArray}[http://narray.rubyforge.org],
+  by Masahiro Tanaka.
+
+
+  {<img src=https://www.pledgie.com/campaigns/15783.png?skin_name=chrome>}[http://www.pledgie.com/campaigns/15783]'
+email:
+- john.o.woods@gmail.com
+executables: []
+extensions:
+- ext/nmatrix/extconf.rb
+extra_rdoc_files:
+- History.txt
+- Manifest.txt
+- README.rdoc
+- LICENSE.txt
+files:
+- .autotest
+- History.txt
+- Manifest.txt
+- README.rdoc
+- LICENSE.txt
+- Rakefile
+- Gemfile
+- spec/nmatrix_spec.rb
+- spec/nmatrix_yale_spec.rb
+- spec/nvector_spec.rb
+- spec/syntax_tree_spec.rb
+- lib/array.rb
+- lib/nmatrix.rb
+- lib/string.rb
+- ext/nmatrix/cblas.c
+- ext/nmatrix/dense.c
+- ext/nmatrix/depend
+- ext/nmatrix/extconf.rb
+- ext/nmatrix/generator.rb
+- ext/nmatrix/list.c
+- ext/nmatrix/nmatrix.c
+- ext/nmatrix/nmatrix.h
+- ext/nmatrix/rational.c
+- ext/nmatrix/yale.c
+- ext/nmatrix/generator/syntax_tree.rb
+- ext/nmatrix/dense/blas_header.template.c
+- ext/nmatrix/dense/elementwise.template.c
+- ext/nmatrix/dense/gemm.template.c
+- ext/nmatrix/dense/gemv.template.c
+- ext/nmatrix/dense/rationalmath.template.c
+- ext/nmatrix/yale/smmp1.template.c
+- ext/nmatrix/yale/smmp2.template.c
+- ext/nmatrix/yale/smmp1_header.template.c
+- ext/nmatrix/yale/smmp2_header.template.c
+- ext/nmatrix/yale/sort_columns.template.c
+- ext/nmatrix/yale/symbmm.template.c
+- ext/nmatrix/yale/numbmm.template.c
+- ext/nmatrix/yale/transp.template.c
+- ext/nmatrix/yale/complexmath.template.c
+- ext/nmatrix/yale/elementwise.template.c
+- ext/nmatrix/yale/elementwise_op.template.c
+- .gemtest
+homepage: http://sciruby.com
+licenses: []
+post_install_message: ! "***********************************************************\nWelcome
+  to SciRuby: Tools for Scientific Computing in Ruby!\n\n                     ***
+  WARNING ***\nPlease be aware that NMatrix is in ALPHA status. If you're\nthinking
+  of using NMatrix to write mission critical code,\nsuch as for driving a car or flying
+  a space shuttle, you\nmay wish to choose other software (for now).\n\nNMatrix requires
+  a C compiler, and has been tested only\nwith GCC 4.6.1. We are happy to accept contributions\nwhich
+  improve the portability of this project.\n\nAlso required is ATLAS. Most Linux distributions
+  and Mac\nversions include ATLAS, but you may wish to compile it\nyourself.\n\nMore
+  explicit instructions for NMatrix and SciRuby should\nbe available on the SciRuby
+  website, sciruby.com, or\nthrough our mailing list (which can be found on our web-\nsite).\n\nThanks
+  for trying out NMatrix! Happy coding!\n\n***********************************************************\n"
+rdoc_options:
+- --main
+- README.rdoc
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  none: false
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: '1.9'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  none: false
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+      segments:
+      - 0
+      hash: -102996319
+requirements: []
+rubyforge_project: nmatrix
+rubygems_version: 1.8.10
+signing_key: 
+specification_version: 3
+summary: NMatrix is an experimental linear algebra library for Ruby, written mostly
+  in C
+test_files: []