ruby-decimal 0.1.0

data/tasks/test.rake

@@ -0,0 +1,40 @@
+
+if test(?e, PROJ.test.file) or not PROJ.test.files.to_a.empty?
+require 'rake/testtask'
+
+namespace :test do
+
+  Rake::TestTask.new(:run) do |t|
+    t.libs = PROJ.libs
+    t.test_files = if test(?f, PROJ.test.file) then [PROJ.test.file]
+                   else PROJ.test.files end
+    t.ruby_opts += PROJ.ruby_opts
+    t.ruby_opts += PROJ.test.opts
+  end
+
+  if HAVE_RCOV
+    desc 'Run rcov on the unit tests'
+    task :rcov => :clobber_rcov do
+      opts = PROJ.rcov.opts.dup << '-o' << PROJ.rcov.dir
+      opts = opts.join(' ')
+      files = if test(?f, PROJ.test.file) then [PROJ.test.file]
+              else PROJ.test.files end
+      files = files.join(' ')
+      sh "#{RCOV} #{files} #{opts}"
+    end
+
+    task :clobber_rcov do
+      rm_r 'coverage' rescue nil
+    end
+  end
+
+end  # namespace :test
+
+desc 'Alias to test:run'
+task :test => 'test:run'
+
+task :clobber => 'test:clobber_rcov' if HAVE_RCOV
+
+end
+
+# EOF

data/test/all_tests.rb

@@ -0,0 +1,11 @@
+require File.dirname(__FILE__) + '/test_flags.rb'
+require File.dirname(__FILE__) + '/test_basic.rb'
+require File.dirname(__FILE__) + '/test_dectest.rb'
+require File.dirname(__FILE__) + '/test_exact.rb'
+require File.dirname(__FILE__) + '/test_round.rb'
+require File.dirname(__FILE__) + '/test_multithreading.rb'
+require File.dirname(__FILE__) + '/test_comparisons.rb'
+require File.dirname(__FILE__) + '/test_coercion.rb'
+require File.dirname(__FILE__) + '/test_to_int.rb'
+require File.dirname(__FILE__) + '/test_to_rf.rb'
+require File.dirname(__FILE__) + '/test_define_conversions.rb'

data/test/helper.rb

@@ -0,0 +1,7 @@
+require 'test/unit'
+require File.dirname(__FILE__) + '/../lib/decimal'
+
+def initialize_context
+  Decimal.context = Decimal::ExtendedContext
+end
+

data/test/test_basic.rb

@@ -0,0 +1,398 @@
+require File.dirname(__FILE__) + '/helper.rb'
+
+def exp(x,c=nil)
+  i, lasts, s, fact, num = 0, 0, 1, 1, 1
+  Decimal.local_context(c) do
+    # result context
+    Decimal.local_context do |context|
+      # working context
+      context.precision += 2
+      context.rounding = Decimal::ROUND_HALF_EVEN
+      while s != lasts
+        lasts = s
+        i += 1
+        fact *= i
+        num *= x
+        s += num / fact
+      end
+    end
+    +s
+  end
+end
+
+def exp1(x, c=nil)
+  return Decimal(BigDecimal("NaN")) if x.infinite? || x.nan?
+  y = nil
+  ext = 2
+  Decimal.local_context(c) do |context|
+    n = (context.precision += ext)
+
+    one  = Decimal("1")
+    x1 = one
+    y  = one
+    d  = y
+    z  = one
+    i  = 0
+    while d.nonzero? && ((m = n - (y.fractional_exponent - d.fractional_exponent).abs) > 0)
+      m = ext if m < ext
+      x1  *= x
+      i += 1
+      z *= i
+
+      #d  = x1.divide(z,:precision=>m)
+      context.precision = m
+      d = x1/z
+      context.precision = n
+
+      y += d
+    end
+  end
+  return +y
+end
+
+
+
+class TestBasic < Test::Unit::TestCase
+
+
+  def setup
+    initialize_context
+  end
+
+
+  def test_basic
+
+    Decimal.context.precision = 4
+    assert_equal 4, Decimal.context.precision
+    assert_equal Decimal("0.3333"), Decimal(1)/Decimal(3)
+    Decimal.context.precision = 10
+    assert_equal 10, Decimal.context.precision
+    assert_equal Decimal("0.3333333333"), Decimal(1)/Decimal(3)
+    Decimal.local_context {
+      assert_equal 10, Decimal.context.precision
+      Decimal.context.precision = 3
+      assert_equal 3, Decimal.context.precision
+      assert_equal Decimal("0.333"), Decimal(1)/Decimal(3)
+    }
+    assert_equal 10, Decimal.context.precision
+    assert_equal "0.3333333333", (Decimal(1)/Decimal(3)).to_s
+    assert_equal Decimal("0.3333333333"), Decimal(1)/Decimal(3)
+
+    Decimal.local_context(:precision=>4) {
+      assert_equal 4, Decimal.context.precision
+    }
+    assert_equal 10, Decimal.context.precision
+
+    Decimal.context(:precision=>4) {
+      assert_equal 4, Decimal.context.precision
+    }
+    assert_equal 10, Decimal.context.precision
+
+    Decimal.local_context(Decimal::BasicContext) {
+      assert_equal :half_up, Decimal.context.rounding
+      assert_equal 9, Decimal.context.precision
+    }
+    assert_equal :half_even, Decimal.context.rounding
+    assert_equal 10, Decimal.context.precision
+
+    Decimal.context(Decimal::BasicContext) {
+      assert_equal :half_up, Decimal.context.rounding
+      assert_equal 9, Decimal.context.precision
+    }
+    assert_equal :half_even, Decimal.context.rounding
+    assert_equal 10, Decimal.context.precision
+
+    Decimal.local_context(Decimal::BasicContext, :precision=>4) {
+      assert_equal :half_up, Decimal.context.rounding
+      assert_equal 4, Decimal.context.precision
+    }
+    assert_equal :half_even, Decimal.context.rounding
+    assert_equal 10, Decimal.context.precision
+
+    Decimal.context(Decimal::BasicContext, :precision=>4) {
+      assert_equal :half_up, Decimal.context.rounding
+      assert_equal 4, Decimal.context.precision
+    }
+    assert_equal :half_even, Decimal.context.rounding
+    assert_equal 10, Decimal.context.precision
+
+
+    assert_equal Decimal("0."+"3"*100), Decimal(1)./(Decimal(3),Decimal.Context(:precision=>100))
+    assert_equal 10, Decimal.context.precision
+    assert_equal Decimal("0.33"), Decimal(1).divide(Decimal(3),Decimal.Context(:precision=>2))
+    assert_equal Decimal("0.3333333333"), Decimal(1)/Decimal(3)
+
+    assert_equal Decimal("11.0"), Decimal(11).abs
+    assert_equal Decimal("11.0"), Decimal(-11).abs
+    assert_equal Decimal("-11.0"), -Decimal(11)
+
+    assert_equal [-11,0], Decimal("-11").to_int_scale
+
+    assert_equal [-110,-1], Decimal("-11.0").to_int_scale
+    assert_equal [-11,-1], Decimal("-1.1").to_int_scale
+    assert_equal [-110,0], Decimal("-110").to_int_scale
+    assert_equal [110,-1], Decimal("11.0").to_int_scale
+    assert_equal [11,-1], Decimal("1.1").to_int_scale
+    assert_equal [110,0], Decimal("110").to_int_scale
+
+    assert_equal [-11,0], Decimal("-11.0").reduce.to_int_scale
+    assert_equal [-11,-1], Decimal("-1.1").reduce.to_int_scale
+    assert_equal [-11,1], Decimal("-110").reduce.to_int_scale
+    assert_equal [11,0], Decimal("11.0").reduce.to_int_scale
+    assert_equal [11,-1], Decimal("1.1").reduce.to_int_scale
+    assert_equal [11,1], Decimal("110").reduce.to_int_scale
+
+
+    assert_equal Decimal('2.1'), Decimal('2.1').remainder(Decimal('3'))
+    assert_equal Decimal('-2.1'), Decimal('-2.1').remainder(Decimal('3'))
+    assert_equal Decimal('2.1'), Decimal('2.1').remainder(Decimal('-3'))
+    assert_equal Decimal('-2.1'), Decimal('-2.1').remainder(Decimal('-3'))
+    assert_equal Decimal('1'), Decimal('10').remainder(Decimal('3'))
+    assert_equal Decimal('-1'), Decimal('-10').remainder(Decimal('3'))
+    assert_equal Decimal('1'), Decimal('10').remainder(Decimal('-3'))
+    assert_equal Decimal('-1'), Decimal('-10').remainder(Decimal('-3'))
+    assert_equal Decimal('0.2'), Decimal('10.2').remainder(Decimal('1'))
+    assert_equal Decimal('0.1'), Decimal('10').remainder(Decimal('0.3'))
+    assert_equal Decimal('1.0'), Decimal('3.6').remainder(Decimal('1.3'))
+    assert_equal Decimal('2'), Decimal('2').remainder(Decimal('3'))
+    assert_equal Decimal('1'), Decimal('10').remainder(Decimal('3'))
+    assert_equal Decimal('.1'), Decimal('1').remainder(Decimal('0.3'))
+
+    assert_equal Decimal('0'), Decimal('2.1').divide_int(Decimal('3'))
+    assert_equal Decimal('0'), Decimal('-2.1').divide_int(Decimal('3'))
+    assert_equal Decimal('0'), Decimal('2.1').divide_int(Decimal('-3'))
+    assert_equal Decimal('0'), Decimal('-2.1').divide_int(Decimal('-3'))
+    assert_equal Decimal('3'), Decimal('10').divide_int(Decimal('3'))
+    assert_equal Decimal('-3'), Decimal('-10').divide_int(Decimal('3'))
+    assert_equal Decimal('-3'), Decimal('10').divide_int(Decimal('-3'))
+    assert_equal Decimal('3'), Decimal('-10').divide_int(Decimal('-3'))
+    assert_equal Decimal('10'), Decimal('10.2').divide_int(Decimal('1'))
+    assert_equal Decimal('33'), Decimal('10').divide_int(Decimal('0.3'))
+    assert_equal Decimal('2'), Decimal('3.6').divide_int(Decimal('1.3'))
+    assert_equal Decimal('0'), Decimal('2').divide_int(Decimal('3'))
+    assert_equal Decimal('3'), Decimal('10').divide_int(Decimal('3'))
+    assert_equal Decimal('3'), Decimal('1').divide_int(Decimal('0.3'))
+
+    assert_equal Decimal('2.1'), Decimal('2.1').modulo(Decimal('3'))
+    assert_equal Decimal('0.9'), Decimal('-2.1').modulo(Decimal('3'))
+    assert_equal Decimal('-0.9'), Decimal('2.1').modulo(Decimal('-3'))
+    assert_equal Decimal('-2.1'), Decimal('-2.1').modulo(Decimal('-3'))
+    assert_equal Decimal('1'), Decimal('10').modulo(Decimal('3'))
+    assert_equal Decimal('2'), Decimal('-10').modulo(Decimal('3'))
+    assert_equal Decimal('-2'), Decimal('10').modulo(Decimal('-3'))
+    assert_equal Decimal('-1'), Decimal('-10').modulo(Decimal('-3'))
+    assert_equal Decimal('0.2'), Decimal('10.2').modulo(Decimal('1'))
+    assert_equal Decimal('0.1'), Decimal('10').modulo(Decimal('0.3'))
+    assert_equal Decimal('1.0'), Decimal('3.6').modulo(Decimal('1.3'))
+    assert_equal Decimal('2'), Decimal('2').modulo(Decimal('3'))
+    assert_equal Decimal('1'), Decimal('10').modulo(Decimal('3'))
+    assert_equal Decimal('.1'), Decimal('1').modulo(Decimal('0.3'))
+
+    assert_equal Decimal('0'), Decimal('2.1').div(Decimal('3'))
+    assert_equal Decimal('-1'), Decimal('-2.1').div(Decimal('3'))
+    assert_equal Decimal('-1'), Decimal('2.1').div(Decimal('-3'))
+    assert_equal Decimal('0'), Decimal('-2.1').div(Decimal('-3'))
+    assert_equal Decimal('3'), Decimal('10').div(Decimal('3'))
+    assert_equal Decimal('-4'), Decimal('-10').div(Decimal('3'))
+    assert_equal Decimal('-4'), Decimal('10').div(Decimal('-3'))
+    assert_equal Decimal('3'), Decimal('-10').div(Decimal('-3'))
+    assert_equal Decimal('10'), Decimal('10.2').div(Decimal('1'))
+    assert_equal Decimal('33'), Decimal('10').div(Decimal('0.3'))
+    assert_equal Decimal('2'), Decimal('3.6').div(Decimal('1.3'))
+    assert_equal Decimal('0'), Decimal('2').div(Decimal('3'))
+    assert_equal Decimal('3'), Decimal('10').div(Decimal('3'))
+    assert_equal Decimal('3'), Decimal('1').div(Decimal('0.3'))
+
+    assert_equal Decimal('-0.9'), Decimal('2.1').remainder_near(Decimal('3'))
+    assert_equal Decimal('-2'), Decimal('10').remainder_near(Decimal('6'))
+    assert_equal Decimal('1'), Decimal('10').remainder_near(Decimal('3'))
+    assert_equal Decimal('-1'), Decimal('-10').remainder_near(Decimal('3'))
+    assert_equal Decimal('0.2'), Decimal('10.2').remainder_near(Decimal('1'))
+    assert_equal Decimal('0.1'), Decimal('10').remainder_near(Decimal('0.3'))
+    assert_equal Decimal('-0.3'), Decimal('3.6').remainder_near(Decimal('1.3'))
+
+
+   assert_equal 2, Decimal('123.4567').adjusted_exponent
+   assert_equal 2, Decimal('123.45670').adjusted_exponent
+   assert_equal 2, Decimal.context.scaleb(Decimal('1.2345670'),2).adjusted_exponent
+   assert_equal 2, Decimal.context.scaleb(Decimal('1.2345670'),Decimal('2')).adjusted_exponent
+   assert_equal Decimal(2), Decimal.context.logb(Decimal('123.4567'))
+   assert_equal Decimal(2), Decimal.context.logb(Decimal('123.45670'))
+   assert_equal 3, Decimal('123.4567').fractional_exponent
+   assert_equal 3, Decimal('123.45670').fractional_exponent
+
+   assert_equal(-7, Decimal.context.normalized_integral_exponent(Decimal('123.4567')))
+   assert_equal(1234567000, Decimal.context.normalized_integral_significand(Decimal('123.4567')))
+
+   assert_equal 7, Decimal('123.4567').number_of_digits
+   assert_equal 8, Decimal('123.45670').number_of_digits
+   assert_equal 7, Decimal('123.45670').reduce.number_of_digits
+
+   assert_equal 1234567, Decimal('123.4567').integral_significand
+   assert_equal 12345670, Decimal('123.45670').integral_significand
+   assert_equal 1234567, Decimal('123.45670').reduce.integral_significand
+
+   assert_equal 2, Decimal('-123.4567').adjusted_exponent
+   assert_equal 2, Decimal('-123.45670').adjusted_exponent
+   assert_equal 2, Decimal.context.scaleb(Decimal('-1.2345670'),2).adjusted_exponent
+   assert_equal 2, Decimal.context.scaleb(Decimal('-1.2345670'),Decimal('2')).adjusted_exponent
+   assert_equal Decimal(2), Decimal.context.logb(Decimal('-123.4567'))
+   assert_equal Decimal(2), Decimal.context.logb(Decimal('-123.45670'))
+   assert_equal 3, Decimal('-123.4567').fractional_exponent
+   assert_equal 3, Decimal('-123.45670').fractional_exponent
+
+   assert_equal(-7, Decimal.context.normalized_integral_exponent(Decimal('-123.4567')))
+   assert_equal(1234567000, Decimal.context.normalized_integral_significand(Decimal('-123.4567')))
+
+   assert_equal 7, Decimal('-123.4567').number_of_digits
+   # assert_equal 9, Decimal('123.45670').number_of_digits # not with BigDecimal
+   assert_equal 7, Decimal('-123.45670').reduce.number_of_digits
+
+   assert_equal(1234567, Decimal('-123.4567').integral_significand)
+   #assert_equal(12345670, Decimal('-123.45670').integral_significand) # not with BigDecimal
+   assert_equal(1234567, Decimal('-123.45670').reduce.integral_significand)
+
+   x = Decimal('123.4567')
+   assert_equal x, Decimal(x.integral_significand)*10**x.integral_exponent
+   assert_equal x, Decimal(Decimal.context.normalized_integral_significand(x))*10**Decimal.context.normalized_integral_exponent(x)
+
+   assert_equal x, Decimal.context.scaleb(x.integral_significand, x.integral_exponent)
+   assert_equal x, Decimal.context.scaleb(Decimal.context.normalized_integral_significand(x),Decimal.context.normalized_integral_exponent(x))
+
+   x = Decimal('-123.4567')
+   assert_equal x.abs, Decimal(x.integral_significand)*10**x.integral_exponent
+   assert_equal x.abs, Decimal(Decimal.context.normalized_integral_significand(x))*10**Decimal.context.normalized_integral_exponent(x)
+
+   assert_equal x.abs, Decimal.context.scaleb(x.integral_significand, x.integral_exponent)
+   assert_equal x.abs, Decimal.context.scaleb(Decimal.context.normalized_integral_significand(x),Decimal.context.normalized_integral_exponent(x))
+
+   Decimal.context.precision = 3
+   Decimal.context.rounding = :half_up
+   assert_equal Decimal("100"), (+Decimal('100.4'))
+   assert_equal Decimal("101"), (+Decimal('101.4'))
+   assert_equal Decimal("101"), (+Decimal('100.5'))
+   assert_equal Decimal("102"), (+Decimal('101.5'))
+   assert_equal Decimal("-101"), (+Decimal('-100.5'))
+   assert_equal Decimal("-102"), (+Decimal('-101.5'))
+   assert_equal Decimal("-100"), (+Decimal('-100.4'))
+   assert_equal Decimal("-101"), (+Decimal('-101.4'))
+   Decimal.context.rounding = :half_even
+   assert_equal Decimal("100"), (+Decimal('100.5'))
+   assert_equal Decimal("101"), (+Decimal('100.51'))
+   assert_equal Decimal("101"), (+Decimal('100.6'))
+   assert_equal Decimal("102"), (+Decimal('101.5'))
+   assert_equal Decimal("101"), (+Decimal('101.4'))
+   assert_equal Decimal("-100"), (+Decimal('-100.5'))
+   assert_equal Decimal("-102"), (+Decimal('-101.5'))
+   assert_equal Decimal("-101"), (+Decimal('-101.4'))
+   Decimal.context.rounding = :half_down
+   assert_equal Decimal("100"), (+Decimal('100.5'))
+   assert_equal Decimal("101"), (+Decimal('101.5'))
+   assert_equal Decimal("-100"), (+Decimal('-100.5'))
+   assert_equal Decimal("-101"), (+Decimal('-101.5'))
+   Decimal.context.rounding = :down
+   assert_equal Decimal("100"), (+Decimal('100.9'))
+   assert_equal Decimal("101"), (+Decimal('101.9'))
+   assert_equal Decimal("-100"), (+Decimal('-100.9'))
+   assert_equal Decimal("-101"), (+Decimal('-101.9'))
+   Decimal.context.rounding = :up
+   assert_equal Decimal("101"), (+Decimal('100.1'))
+   assert_equal Decimal("102"), (+Decimal('101.1'))
+   assert_equal Decimal("-101"), (+Decimal('-100.1'))
+   assert_equal Decimal("-102"), (+Decimal('-101.1'))
+   Decimal.context.rounding = :floor
+   assert_equal Decimal("100"), (+Decimal('100.9'))
+   assert_equal Decimal("101"), (+Decimal('101.9'))
+   assert_equal Decimal("-101"), (+Decimal('-100.9'))
+   assert_equal Decimal("-102"), (+Decimal('-101.9'))
+   Decimal.context.rounding = :ceiling
+   assert_equal Decimal("101"), (+Decimal('100.1'))
+   assert_equal Decimal("102"), (+Decimal('101.1'))
+   assert_equal Decimal("-100"), (+Decimal('-100.1'))
+   assert_equal Decimal("-101"), (+Decimal('-101.1'))
+
+  end
+
+  def test_exp
+    Decimal.context.precision = 100
+    Decimal.context.rounding = :half_even
+    e_100 = "2.718281828459045235360287471352662497757247093699959574966967627724076630353547594571382178525166427"
+    assert_equal e_100, exp(Decimal(1)).to_s
+    assert_equal Decimal(e_100), exp(Decimal(1))
+    assert_equal e_100, exp1(Decimal(1)).to_s
+    assert_equal Decimal(e_100), exp1(Decimal(1))
+  end
+
+  def test_special
+    nan = Decimal(0)/Decimal(0)
+    inf_pos = Decimal(1)/Decimal(0)
+    inf_neg = Decimal(-1)/Decimal(0)
+    zero_pos = Decimal(0)
+    zero_neg = Decimal('-0')
+    pos = Decimal(1)
+    neg = Decimal(-1)
+
+    assert nan.nan?
+    assert nan.special?
+    assert !nan.zero?
+    assert !nan.finite?
+    assert !nan.infinite?
+    assert inf_pos.infinite?
+    assert !inf_pos.finite?
+    assert !inf_pos.nan?
+    assert inf_pos.special?
+    assert inf_neg.infinite?
+    assert !inf_neg.finite?
+    assert !inf_neg.nan?
+    assert inf_neg.special?
+    assert zero_pos.finite?
+    assert zero_neg.finite?
+    assert zero_pos.zero?
+    assert zero_neg.zero?
+    assert !pos.zero?
+    assert !neg.zero?
+    assert pos.nonzero?
+    assert neg.nonzero?
+
+    #assert nan.sign.nil?
+    assert_equal(+1, inf_pos.sign)
+    assert_equal(-1, inf_neg.sign)
+    assert_equal(+1, zero_pos.sign)
+    assert_equal(-1, zero_neg.sign)
+    assert_equal(+1, pos.sign)
+    assert_equal(-1, neg.sign)
+  end
+
+  def test_context_parameters
+    #Decimal.context = Decimal.Context
+    Decimal.context.precision = 3
+    Decimal.context.rounding = :half_even
+    x = Decimal(1)
+    y = Decimal(3)
+    assert_equal Decimal('0.333'), x.divide(y)
+    assert_equal Decimal('0.3333'), x.divide(y, Decimal.Context(:precision=>4))
+    assert_equal Decimal('0.333'), x.divide(y)
+    assert_equal Decimal('0.33333'), x.divide(y, :precision=>5)
+  end
+
+  def test_integer
+
+    %w{ 0 0E10 0E-10 12.0 12 120E-1 12.1E10 }.each do |x|
+      assert Decimal(x).integral?, "Decimal('#{x}').integral?"
+    end
+
+    [12, Rational(6,2)].each do |x|
+      assert Decimal(x).integral?, "Decimal(#{x}).integral?"
+    end
+
+    %w{ NaN Inf 12.1 121E-1 }.each do |x|
+      assert !Decimal(x).integral?, "!Decimal('#{x}').integral?"
+    end
+
+    [ Rational(121,10) ].each do |x|
+      assert !Decimal(x).integral?, "!Decimal(#{x}).integral?"
+    end
+
+  end
+
+end