ruby-mpfr 0.0.2 → 0.0.4

data/spec/mpfr/set_value_spec.rb

@@ -0,0 +1,77 @@
+require File.dirname(__FILE__) + '/spec_helper.rb'
+
+MPFR.set_default_prec(300)
+
+describe MPFR, 'when setting abnormal number' do
+  it "should be set nan" do
+    a = MPFR.new(rand)
+    a.set_nan
+    a.nan?.should be_true
+    a.to_s.should == "nan"
+  end
+
+  it "should be set infinity" do
+    a = MPFR.new(rand)
+    a.set_inf(0)
+    a.infinite?.should be_true
+    a.to_s.should == "inf"
+    
+    a.set_inf(1)
+    a.infinite?.should be_true
+    a.to_s.should == "inf"
+    
+    a.set_inf(-1)
+    a.infinite?.should be_true
+    a.to_s.should == "-inf"
+  end  
+end
+
+describe MPFR, 'when setting value' do
+  it "should swap each other" do
+    num = 100
+    args = [GenerateNumber.mpfr_args(num), GenerateNumber.mpfr_args(num)]
+    args[0].each_index{ |i| check_swap(args[0][i], args[1][i]) }
+  end
+
+  it "should set other value" do
+    [0, 3, 2.2, -0.2, '1.0e-3', '-2e14', '2225'].each{ |val| check_set(val) }
+  end
+
+  it "should set value by set_fixnum_2exp" do
+    for i in 0..100
+      check_fixnum_2exp(rand(100), rand(100))
+    end
+  end
+end
+
+def check_swap(a, b)
+  tmp_a = a.dup
+  tmp_b = b.dup
+  tmp_a.swap(tmp_b)
+  a.should == tmp_b
+  b.should == tmp_a
+end
+
+def check_set(val)
+  tmp = MPFR.new
+  tmp.set(val)
+  tmp2 = tmp.dup
+  a = (tmp2 + 3) / Math::PI
+  tmp.nan?.should_not be_true
+  tmp.should_not == a
+end
+
+def check_fixnum_2exp(num, exp)
+  ret = MPFR.new(num)
+  if exp < 0
+    for i in 0...(exp.abs)
+      ret /= MPFR.new(2)
+    end
+  else
+    for i in 0...exp
+      ret *= MPFR.new(2)
+    end
+  end
+  MPFR.new.set_fixnum_2exp(num, exp).should == ret
+end
+

data/spec/mpfr/spec_helper.rb

@@ -0,0 +1,13 @@
+begin
+  require 'spec'
+rescue LoadError
+  require 'rubygems' unless ENV['NO_RUBYGEMS']
+  gem 'rspec'
+  require 'spec'
+end
+
+$:.unshift(File.dirname(__FILE__) + '/../../lib')
+$:.unshift(File.dirname(__FILE__) + '/../../ext')
+$:.unshift(File.dirname(__FILE__))
+require 'mpfr'
+require "generate_number_modulue"

data/spec/mpfr/string_spec.rb

@@ -0,0 +1,58 @@
+require File.dirname(__FILE__) + '/spec_helper.rb'
+
+describe MPFR, "when converted to string" do
+  before(:all) do
+    MPFR.set_default_prec(200)
+    @args = GenerateNumber.float(10).map{ |a| MPFR.new(a.to_s) }
+  end
+
+  it "should be to_s" do
+    @args.each do |a|
+      a.to_s.should match(/^(-|)[1-9]\.[0-9]+e[-+][0-9]+$/)
+    end
+  end
+
+  it "should be formatted string \"#<MPFR:%lx,'%.Re',%d>\" for inspect" do
+    @args.each do |a|
+      a.inspect.should match(/^#<MPFR:.+,'(-|)[1-9]\.[0-9]+e[-+][0-9]+',\d+>$/)
+    end
+  end
+
+  it "should be formatted string by to_strf" do
+    @args.each do |a|
+      a.to_strf("%.Re").should eql a.to_s
+      a.to_strf("%.RE").should match(/^(-|)[1-9]\.[0-9]+E[-+][0-9]+$/)
+      a.to_strf("%.Rf").should match(/^(-|)[0-9]*\.[0-9]*$/)
+    end
+
+    format = ["14.6", ".8", "10.10", "+10.4", "-10.4", "+012.2", "-012.2", "3.0"]
+    @args.each do |a|
+      MPFR.set_default_prec(100)
+      for i in 0...1000
+        a = rand(100) + rand
+        check_to_s_e(a, format)
+        check_to_s_e(-a, format)
+        # check_to_s_f(a, format)
+        # check_to_s_f(-a, format)
+      end
+    end
+  end
+
+end
+
+def check_to_s_e(float, format_ary)
+  a = MPFR.new(float)
+  format_ary.each do |format|
+    sprintf("%#{format}e", float).should == a.to_strf("%#{format}Re")
+    sprintf("%#{format}E", float).should == a.to_strf("%#{format}RE")
+  end
+end
+
+# Probably MPFR bug makes this test be failure. Need further investigations. (Tue Aug 18 10:43:51 2009)
+def check_to_s_f(float, format_ary)
+  a = MPFR.new(float)
+  format_ary.each do |format|
+    sprintf("%#{format}f", float).should == a.to_strf("%#{format}Rf")
+  end
+end
+

data/spec/mpfr_matrix/generate_matrix_arguments.rb

@@ -0,0 +1,55 @@
+module GenerateNumber
+
+  def self.float_matrix_arguments(number, row, column, max = 100)
+    ret = []
+    (0...number).each do
+      ret << Array.new(row){ |j| Array.new(column){ |i| rand(max) - rand } }
+      for i in 0...(ret[-1].size)
+        ret[-1][i].map!{ |a| (rand > 0.5 ? -a : a) }
+      end
+    end
+    ret
+  end
+
+  def self.float_arguments(size, number, max = 100)
+    ret = Array.new(number){ |j| Array.new(size){ |i| rand(max) - rand } }
+    for i in 0...(ret.size)
+      case i % 4
+      when 1
+        ret[i].map!{ |a| -a }
+      when 2, 3
+        ret[i].map!{ |a| (rand > 0.5 ? -a : a) }
+      end
+    end
+    ret
+  end
+
+  def self.float(num, max = 100, &block)
+    if block_given?
+      for i in 0...num
+        a = rand(max) - rand
+        a = -a if rand > 0.5
+        yield(a)
+      end
+    else
+      ary = Array.new(num){ |i| rand(max) - rand }
+      ary.map!{ |a| (rand > 0.5 ? -a : a) }
+    end
+  end
+
+  def self.string(number)
+    prec = MPFR.get_default_prec
+    max = 2 ** prec
+    Array.new(number) do |i|
+      sign = ((rand > 0.5 ? '-' : ''))
+      "#{sign}#{rand(max)}.#{rand(max)}"
+    end
+  end
+  
+  def self.mpfr_args(number)
+    ret = self.string(number)
+    ret.map!{ |a| MPFR.new(a) }
+  end
+  
+
+end

data/spec/mpfr_matrix/mpfr_matrix_alloc_spec.rb

@@ -0,0 +1,126 @@
+require File.dirname(__FILE__) + '/spec_helper.rb'
+
+describe "initialization of matrix" do
+  before(:all) do
+    MPFR.set_default_prec(128)
+    @sizes = [[1, 2], [2, 3], [8, 5], [7, 3]]
+    @matrixes = @sizes.map{ |a| MPFR::Matrix.new(*a) }
+  end
+
+  it "should has size which equals size of array." do
+    @matrixes.each_with_index do |a, i|
+      a.size.should == (@sizes[i][0] * @sizes[i][1])
+    end
+  end
+
+  it "should has size of row which equals size of array." do
+    @matrixes.each_with_index do |a, i|
+      a.row_size.should == @sizes[i][0]
+    end
+  end
+
+  it "should has size of column which equals size of array." do
+    @matrixes.each_with_index do |a, i|
+      a.column_size.should == @sizes[i][1]
+    end
+  end
+
+  it "should have MPFR instances as elements." do
+    @matrixes.each do |m|
+      m.each { |a| a.should be_an_instance_of MPFR }
+    end
+  end
+  
+end
+
+describe "argument error raises for invalid initialization argument" do
+  before(:all) do
+    MPFR.set_default_prec(256)
+    @data = [[[1, 2, 3], [0, 0]],
+             [[MPFR.new(0)], [MPFR.new(1), 3]]]
+  end
+  
+  it "should make error for invalid array argument" do
+    @data.each do |ary|
+      lambda { MPFR::Matrix.new(ary) }.should raise_error ArgumentError
+    end
+  end
+end
+
+describe "when methods of matrix 'dup' and 'transpose' are executed" do
+  before(:all) do
+    MPFR.set_default_prec(256)
+    @data = [[[1, 2], [3, 3]],
+             [[-1, 0.8], ['1.2342', 99]],
+             [[36.237582938, 7122], ['-1e7', 8237]],
+             [[-2.2738, 8.23, 9.237482], [11, '8237.2374', '1239']],
+             [[MPFR.new('23'), MPFR.new('123.23')], [MPFR.new(-3.2378), MPFR.new('1.2837')]]]
+    @matrixes = @data.map{ |a| MPFR::Matrix.new(a) }
+  end
+
+  it "should have MPFR instances as elements." do
+    @matrixes.each do |m|
+      m.each { |a| a.should be_an_instance_of MPFR }
+    end
+  end  
+
+  it "should create matrix which has same values as original and different object id." do
+    @matrixes.each do |m|
+      new = m.dup
+      m.object_id.should_not == new.object_id
+      m.should == new
+    end
+  end
+  
+  it "should has transposed values." do
+    @matrixes.each do |m|
+      id = m.object_id
+      old = m.dup
+      new = m.transpose
+      m.transpose!
+      m.should_not == old
+      m.should == new
+      m.object_id.should == id
+      m.object_id.should_not == new.object_id
+      m.column_size.should == old.row_size
+      m.row_size.should == old.column_size
+    end
+  end
+end
+
+describe "when methods of vector 'dup' and 'transpose' are executed" do
+  before(:all) do
+    MPFR.set_default_prec(768)
+    @data = [[1, 2, 3, 3],
+             [-1, 0.8, '1.2342', 99],
+             [36.237582938, 7122, '-1e7', 8237],
+             [-2.2738, 8.23, 9.237482, '8237.2374', '1239'],
+             [MPFR.new('23'), MPFR.new('123.23'), MPFR.new(-3.2378), MPFR.new('1.2837')]]
+    @column = @data.map{ |a| MPFR::ColumnVector.new(a) }
+    @row = @data.map{ |a| MPFR::RowVector.new(a) }
+  end
+
+  it "should create matrix which has same values as original and different object id." do
+    (@column + @row).each do |m|
+      new = m.dup
+      m.object_id.should_not == new.object_id
+      m.should == new
+    end
+  end
+
+  it "should change row vector to column vector and make the opposite changes" do
+    (@column + @row).each do |m|
+      id = m.object_id
+      old = m.dup
+      new = m.transpose
+      m.transpose!
+      m.should_not == old
+      m.should == new
+      m.object_id.should == id
+      m.object_id.should_not == new.object_id
+      m.column_size.should == old.row_size
+      m.row_size.should == old.column_size
+    end
+  end
+  
+end

data/spec/mpfr_matrix/mpfr_matrix_arithmetic_spec.rb

@@ -0,0 +1,93 @@
+require File.dirname(__FILE__) + '/spec_helper.rb'
+
+describe MPFR::Matrix, "when making matrix negative" do
+  before(:all) do
+    MPFR.set_default_prec(512)
+    @row = 8
+    @column = 3
+    @args = GenerateNumber.float_matrix_arguments(1000, @row, @column).map{ |a| MPFR::Matrix.new(a) }
+  end
+
+  it "should be negative" do
+    zero_matrix = MPFR::Matrix.new(@row, @column)
+    @args.each { |m| (zero_matrix - m).should == m.neg }
+  end
+  
+end
+
+describe MPFR::Matrix, "when arithmetic operator applies to 2x3 matrixies" do
+  before(:all) do
+    MPFR.set_default_prec(128)
+    @row = 2
+    @column = 3
+    @args = GenerateNumber.float_matrix_arguments(1000, @row, @column).map!{ |a| MPFR::Matrix.new(a) }
+  end
+
+  it "should be added" do
+    @args.each_index do |i|
+      if i > 0
+        res = @args[i-1] + @args[i]
+        (0...@row).each do |j|
+          (0...@column).each { |k| res[j, k].should eql(@args[i-1][j, k] + @args[i][j, k]) }
+        end
+      end
+    end
+  end
+
+  it "should be subtracted" do
+    @args.each_index do |i|
+      if i > 0
+        res = @args[i-1] + @args[i]
+        (0...@row).each do |j|
+          (0...@column).each { |k| res[j, k].should eql(@args[i-1][j, k] + @args[i][j, k]) }
+        end
+      end
+    end
+  end
+
+  it "should be multiplied" do
+    @args.each_index do |i|
+      if i > 0
+        trans = @args[i].transpose
+        res = @args[i-1] * trans
+        (0...@row).each do |j|
+          (0...@row).each do |k|
+            res2 = MPFR.new(0)
+            (0...@column).each { |l| res2 += @args[i-1][j, l] * trans[l, k] }
+            res[j, k].should eql res2
+          end
+        end
+        trans = @args[i-1].transpose
+        res = trans * @args[i]
+        (0...@column).each do |j|
+          (0...@column).each do |k|
+            res2 = MPFR.new(0)
+            (0...@row).each { |l| res2 += trans[j, l] * @args[i][l, k] }
+            res[j, k].should eql res2
+          end
+        end
+      end
+    end
+  end
+
+  it "should be multyplied by scalar" do
+    @args.each do |a|
+      scalar = MPFR.new(GenerateNumber.float(1)[0])
+      res = a.mul_scalar(scalar)
+      (0...@row).each do |j|
+        (0...@column).each { |k| res[j, k].should eql a[j, k] * scalar }
+      end
+    end
+  end
+  
+  it "should be divided by scalar" do
+    @args.each do |a|
+      scalar = MPFR.new(GenerateNumber.float(1)[0])
+      res = a.div_scalar(scalar)
+      (0...@row).each do |j|
+        (0...@column).each { |k| res[j, k].should eql a[j, k] / scalar }
+      end
+    end
+  end
+  
+end

data/spec/mpfr_matrix/mpfr_matrix_set_element_spec.rb

@@ -0,0 +1,55 @@
+require File.dirname(__FILE__) + '/spec_helper.rb'
+
+describe MPFR::Matrix, "when setting number to particular element" do
+  before(:all) do
+    MPFR.set_default_prec(128)
+    @row = 4
+    @column = 7
+    @matrix = MPFR::Matrix.new(@row, @column)
+    @source = GenerateNumber.float_matrix_arguments(1, @row, @column)[0]
+  end
+
+  it "should have elements which equals to source number." do
+    (0...(@matrix.row_size)).each do |i|
+      (0...(@matrix.column_size)).each do |j|
+        fr = MPFR.new(@source[i][j].to_s)
+        @matrix.set_element(i, j, fr)
+        @matrix[i, j].should eql fr
+      end
+    end
+  end
+
+end
+
+describe MPFR::Matrix, "when calculating row and column" do
+  before(:all) do
+    MPFR.set_default_prec(128)
+    @row = 6
+    @column = 4
+    @sources = GenerateNumber.float_matrix_arguments(100, @row, @column)
+    @matrixies = @sources.map { |a| MPFR::Matrix.new(a) }
+  end
+
+  it "should have row vector." do
+    @matrixies.each do |m|
+      (0...(m.column_size)).each do |i|
+        col = m.column(i)
+        (0...(m.row_size)).each do |j|
+          col[j].should eql m[j, i]
+        end
+      end
+    end
+  end
+
+  it "should have column vector." do
+    @matrixies.each do |m|
+      (0...(m.row_size)).each do |i|
+        row = m.row(i)
+        (0...(m.column_size)).each do |j|
+          row[j].should eql m[i, j]
+        end
+      end
+    end
+  end
+  
+end

data/spec/mpfr_matrix/mpfr_matrix_string_spec.rb

@@ -0,0 +1,31 @@
+require File.dirname(__FILE__) + '/spec_helper.rb'
+
+describe MPFR::Matrix, "when MPFR instance is converted to string by some methods." do
+  
+  before(:all) do
+    MPFR.set_default_prec(128)
+    @row = 2
+    @column = 3
+    @matrixes = ((0..100).map { |i| MPFR::Matrix.new(GenerateNumber.float_arguments(@column, @row)) })
+  end
+
+  it "should return array having strings." do
+    @matrixes.each do |m|
+      m.str_ary("%.Re").each_with_index do |s, i|
+        s.to_s.should eql m[i % @row, i / @row].to_s
+      end
+    end
+  end
+
+  it "should return inspect string." do
+    @matrixes.each do |m|
+      ins = m.inspect
+      ins.should match(/#{Regexp.escape(sprintf("%x", m.object_id))}/)
+      m.each do |a|
+        ins.should match(/#{Regexp.escape(a.to_s)}/)
+      end
+    end
+  end
+  
+end
+

data/spec/mpfr_matrix/mpfr_square_matrix_spec.rb

@@ -0,0 +1,75 @@
+require File.dirname(__FILE__) + '/spec_helper.rb'
+
+describe MPFR::SquareMatrix, "when calculating LU decomposition" do
+  before(:all) do
+    MPFR.set_default_prec(128)
+    @dim = 7
+    @matrixies = GenerateNumber.float_matrix_arguments(100, @dim, @dim).map{ |a| MPFR::SquareMatrix.new(a) }
+    @error = MPFR.new('1e-30')
+  end
+
+  it "should have lu decomposition" do
+    @matrixies.each do |mat|
+      l, u, indx = mat.lu_decomp
+      if l
+        mul = l * u
+        (0...(indx.size)).to_a.reverse.each do |i|
+          if i != indx[i]
+            row = mul.row(i)
+            (0...(mul.column_size)).each do |j|
+              mul.set_element(i, j, mul[indx[i], j])
+              mul.set_element(indx[i], j, row[j])
+            end
+          end
+        end
+        (0...(mat.row_size)).each do |i|
+          (0...(mat.column_size)).each do |j|
+            (mat[i, j] - mul[i, j]).abs.should < @error
+          end
+        end
+      end
+    end
+  end
+end
+
+describe MPFR::SquareMatrix, "when calculating determinant" do
+  before(:all) do
+    MPFR.set_default_prec(512)
+    @error = MPFR.new('1e-30')
+  end
+
+  it "should have the same values of determinant by formula and lu decomposition" do
+    [2, 3].each do |dim|
+      GenerateNumber.float_matrix_arguments(100, dim, dim).each do |src|
+        m = MPFR::SquareMatrix.new(src)
+        l, u, idx = m.lu_decomp
+        if l
+          det = MPFR.new(1)
+          (0...(u.row_size)).each { |i| det *= u[i, i] }
+          count = 0
+          idx.each_with_index { |ind, i| count += (ind != i ? 1 : 0) }
+          det = -det if count % 2 == 1
+          (det - m.determinant).should < @error
+        end
+      end
+    end
+  end
+
+  it "should have the same values of determinant by formula and lu decomposition" do
+    [4, 5, 6].each do |dim|
+      GenerateNumber.float_matrix_arguments(100, dim, dim).each do |src|
+        m = MPFR::SquareMatrix.new(src)
+        l, u, idx = m.lu_decomp
+        if l
+          det = MPFR.new(1)
+          (0...(u.row_size)).each { |i| det *= u[i, i] }
+          count = 0
+          idx.each_with_index { |ind, i| count += (ind != i ? 1 : 0) }
+          det = -det if count % 2 == 1
+          (det - m.determinant).should < @error
+        end
+      end
+    end
+  end
+
+end

data/spec/mpfr_matrix/spec_helper.rb

@@ -0,0 +1,16 @@
+begin
+  require 'spec'
+rescue LoadError
+  require 'rubygems' unless ENV['NO_RUBYGEMS']
+  gem 'rspec'
+  require 'spec'
+end
+
+$:.unshift(File.dirname(__FILE__) + '/../../lib')
+$:.unshift(*Dir.glob(File.dirname(__FILE__) + '/../../ext/*'))
+$:.unshift(File.dirname(__FILE__))
+
+require 'mpfr'
+require 'mpfr/matrix'
+require "generate_matrix_arguments"
+

data/tasks/extconf.rake

@@ -0,0 +1,36 @@
+begin
+  require 'spec'
+rescue LoadError
+  require 'rubygems' unless ENV['NO_RUBYGEMS']
+  require 'spec'
+end
+begin
+  require 'spec/rake/spectask'
+rescue LoadError
+  puts <<-EOS
+To use rspec for testing you must install rspec gem:
+    gem install rspec
+EOS
+  exit(0)
+end
+
+desc "Run 'make realclean' for extended libraries"
+task "ext:realclean" do
+  Dir.glob("ext/**/Makefile").each do |path|
+    system("cd #{File.dirname(path)}; make realclean")
+  end
+end
+
+desc "Run 'make clean' for extended libraries"
+task "ext:clean" do
+  Dir.glob("ext/**/Makefile").each do |path|
+    system("cd #{File.dirname(path)}; make clean")
+  end
+end
+
+desc "Run 'make realclean' for extended libraries"
+task 'ext:make' do
+  Dir.glob("ext/**/extconf.rb").each do |path|
+    system("cd #{File.dirname(path)}; ruby extconf.rb && make")
+  end
+end