long-decimal 0.00.08 → 0.00.09

data/test/testlongdecimal.rb

@@ -2,8 +2,8 @@
 #
 # testlongdecimal.rb -- runit test for longdecimal.rb
 #
-# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testlongdecimal.rb,v 1.4 2006/03/02 20:20:12 bk1 Exp $
-# CVS-Label: $Name: PRE_ALPHA_0_08 $
+# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testlongdecimal.rb,v 1.5 2006/03/04 21:49:00 bk1 Exp $
+# CVS-Label: $Name: PRE_ALPHA_0_09 $
 # Author:    $Author: bk1 $ (Karl Brodowsky)
 #
 
@@ -18,7 +18,95 @@ load "lib/longdecimal.rb"
 #
 class TestLongDecimal_class < RUNIT::TestCase
 
-  @RCS_ID='-$Id: testlongdecimal.rb,v 1.4 2006/03/02 20:20:12 bk1 Exp $-'
+  @RCS_ID='-$Id: testlongdecimal.rb,v 1.5 2006/03/04 21:49:00 bk1 Exp $-'
+
+  def check_split_merge_words(x, l, wl)
+    w = LongMath.split_to_words(x, l)
+    y = LongMath.merge_from_words(w, l)
+    assert_equal(x, y, "#{x} splitted and merged should be equal but is #{y} l=#{l}")
+    assert_equal(wl, w.length, "#{x} splitted to l=#{l} should have length #{wl} but has #{w.length}")
+    w
+  end
+
+  #
+  # test split_to_words and merge_from_words
+  #
+  def test_split_merge_words
+    check_split_merge_words(0, 1, 1)
+    check_split_merge_words(0, 10, 1)
+    check_split_merge_words(0, 100, 1)
+    check_split_merge_words(0, 1000, 1)
+    check_split_merge_words(-1, 1, 1)
+    check_split_merge_words(-1, 10, 1)
+    check_split_merge_words(-1, 100, 1)
+    check_split_merge_words(-1, 1000, 1)
+    check_split_merge_words(1, 1, 1)
+    check_split_merge_words(1, 10, 1)
+    check_split_merge_words(1, 100, 1)
+    check_split_merge_words(1, 1000, 1)
+    check_split_merge_words(10, 1, 4)
+    check_split_merge_words(10, 10, 1)
+    check_split_merge_words(10, 100, 1)
+    check_split_merge_words(10, 1000, 1)
+    x = 10**12 # 10**12 = 1000**4 < 2**40
+    check_split_merge_words(x, 1, 40)
+    check_split_merge_words(x, 10, 4)
+    check_split_merge_words(x, 100, 1)
+    check_split_merge_words(x, 1000, 1)
+    x = 2**40
+    check_split_merge_words(x, 1, 41)
+    check_split_merge_words(x, 10, 5)
+    check_split_merge_words(x, 100, 1)
+    check_split_merge_words(x, 1000, 1)
+  end
+
+  #
+  # helper method for test_exp
+  #
+  def check_exp_floated(x, prec)
+    x0 = x
+    x = x.to_ld
+    y = LongMath.exp(x, prec)
+    z = Math.exp(x.to_f)
+    yf = y.to_f
+    assert((yf - z) / [yf.abs, z.abs, Float::MIN].max < 1e-9, "y=#{yf.to_s} and z=#{z.to_s} should be almost equal x=#{x}")
+    yd = LongMath.exp_internal(x, prec, nil, nil, nil, nil, LongDecimal::ROUND_DOWN)
+    yu = LongMath.exp_internal(x, prec, nil, nil, nil, nil, LongDecimal::ROUND_UP)
+    #    if (x >= 0)
+    assert_equal(yd, yu, "the result yd/yu should not depend on the internal rounding mode x0=#{x0} x=#{x} p=#{prec} d=#{(yd-yu).to_f.to_s}")
+    assert_equal(y,  yu, "the result y/yu  should not depend on the internal rounding mode x0=#{x0} x=#{x} p=#{prec} d=#{(y -yu).to_f.to_s}")
+    #   else
+    #      assert((yu - yd).abs <= y.unit, "the difference should not be too much x=#{x} d=#{(yd-yu).to_f.to_s}")
+    #      assert((yu - y ).abs <= y.unit, "the difference should not be too much x=#{x} d=#{(y -yu).to_f.to_s}")
+    #    end
+  end
+
+  #
+  # test the calculation of the exponential function
+  #
+  def test_exp
+    check_exp_floated(700, 10)
+    check_exp_floated(100, 10)
+    check_exp_floated(1, 10)
+    check_exp_floated(0.01, 10)
+    check_exp_floated(1e-10, 10)
+    check_exp_floated(1e-90, 10)
+    check_exp_floated(0, 10)
+    check_exp_floated(-1, 10)
+    check_exp_floated(-100, 10)
+    check_exp_floated(-700, 10)
+
+    check_exp_floated(700, 100)
+    check_exp_floated(100, 100)
+    check_exp_floated(1, 100)
+    check_exp_floated(0.01, 100)
+    check_exp_floated(1e-10, 100)
+    check_exp_floated(1e-90, 100)
+    check_exp_floated(0, 100)
+    check_exp_floated(-1, 100)
+    check_exp_floated(-100, 100)
+    check_exp_floated(-700, 100)
+  end
 
   #
   # helper method for test_sqrtb
@@ -45,22 +133,22 @@ class TestLongDecimal_class < RUNIT::TestCase
       y = check_sqrtb(x1, " i*i i=#{i}")
       assert_equal(i, y, "i=#{i} y=#{y}")
       if (i > 0) then
-	x2 = x1 + 1
-	y = check_sqrtb(x2, " i*i+1 i=#{i}")
-	assert_equal(i, y, "i=#{i} y=#{y}")
-	x0 = x1 - 1
-	y = check_sqrtb(x0, " i*i-1 i=#{i}")
-	assert_equal(i-1, y, "i=#{i} y=#{y}")
+        x2 = x1 + 1
+        y = check_sqrtb(x2, " i*i+1 i=#{i}")
+        assert_equal(i, y, "i=#{i} y=#{y}")
+        x0 = x1 - 1
+        y = check_sqrtb(x0, " i*i-1 i=#{i}")
+        assert_equal(i-1, y, "i=#{i} y=#{y}")
       end
-      
+
       x1 = 1 << i
       y = check_sqrtb(x1, " 2**i i=#{i}")
       if (i[0] == 0)
-	assert_equal(1 << (i>>1), y, "2^(i/2) i=#{i} y=#{y}")
+        assert_equal(1 << (i>>1), y, "2^(i/2) i=#{i} y=#{y}")
       end
       if (i > 0) then
-	check_sqrtb(x1-1, " 2**i-1 i=#{i}")
-	check_sqrtb(x1+1, " 2**i+1 i=#{i}")
+        check_sqrtb(x1-1, " 2**i-1 i=#{i}")
+        check_sqrtb(x1+1, " 2**i+1 i=#{i}")
       end
 
       x1 = 3 << i
@@ -95,22 +183,22 @@ class TestLongDecimal_class < RUNIT::TestCase
       y = check_sqrtw(x1, " i*i i=#{i}")
       assert_equal(i, y, "i=#{i} y=#{y}")
       if (i > 0) then
-	x2 = x1 + 1
-	y = check_sqrtw(x2, " i*i+1 i=#{i}")
-	assert_equal(i, y, "i=#{i} y=#{y}")
-	x0 = x1 - 1
-	y = check_sqrtw(x0, " i*i-1 i=#{i}")
-	assert_equal(i-1, y, "i=#{i} y=#{y}")
+        x2 = x1 + 1
+        y = check_sqrtw(x2, " i*i+1 i=#{i}")
+        assert_equal(i, y, "i=#{i} y=#{y}")
+        x0 = x1 - 1
+        y = check_sqrtw(x0, " i*i-1 i=#{i}")
+        assert_equal(i-1, y, "i=#{i} y=#{y}")
       end
-      
+
       x1 = 1 << i
       y = check_sqrtw(x1, " 2**i i=#{i}")
       if (i[0] == 0)
-	assert_equal(1 << (i>>1), y, "2^(i/2) i=#{i} y=#{y}")
+        assert_equal(1 << (i>>1), y, "2^(i/2) i=#{i} y=#{y}")
       end
       if (i > 0) then
-	check_sqrtw(x1-1, " 2**i-1 i=#{i}")
-	check_sqrtw(x1+1, " 2**i+1 i=#{i}")
+        check_sqrtw(x1-1, " 2**i-1 i=#{i}")
+        check_sqrtw(x1+1, " 2**i+1 i=#{i}")
       end
 
       x1 = 3 << i
@@ -154,7 +242,7 @@ class TestLongDecimal_class < RUNIT::TestCase
 
   #
   # test multiplicity_of_factor for rationals with numerator and
-  # denominator exceeding Float 
+  # denominator exceeding Float
   #
   def test_rat_long_multiplicity_of_factor
     n = Rational(224*(10**600+1), 225*(5**800))

metadata

@@ -3,8 +3,8 @@ rubygems_version: 0.8.11
 specification_version: 1
 name: long-decimal
 version: !ruby/object:Gem::Version 
-  version: 0.00.08
-date: 2006-03-02 00:00:00 +01:00
+  version: 0.00.09
+date: 2006-03-04 00:00:00 +01:00
 summary: LongDecimal for numbers with fixed point
 require_paths: 
 - lib
@@ -28,6 +28,7 @@ cert_chain:
 authors: 
 - Karl Brodowsky
 files: 
+- doc
 - lib
 - test
 - Rakefile
@@ -36,6 +37,152 @@ files:
 - VERSION
 - make_doc.rb
 - version.rb
+- doc/dot
+- doc/files
+- doc/index.html
+- doc/rdoc-style.css
+- doc/fr_method_index.html
+- doc/fr_class_index.html
+- doc/fr_file_index.html
+- doc/created.rid
+- doc/classes
+- doc/dot/m_0_0.dot
+- doc/dot/m_0_0.png
+- doc/dot/m_0_1.dot
+- doc/dot/m_0_1.png
+- doc/dot/f_0.dot
+- doc/dot/f_0.png
+- doc/files/lib
+- doc/files/lib/longdecimal_rb.src
+- doc/files/lib/longdecimal_rb.html
+- doc/files/lib/longdecimal_rb.src/M000001.html
+- doc/files/lib/longdecimal_rb.src/M000002.html
+- doc/classes/LongDecimal.src
+- doc/classes/LongMath.html
+- doc/classes/Numeric.src
+- doc/classes/LongDecimal.html
+- doc/classes/LongDecimalQuot.src
+- doc/classes/LongMath.src
+- doc/classes/LongDecimalQuot.html
+- doc/classes/LongDecimalRoundingMode.html
+- doc/classes/Numeric.html
+- doc/classes/LongDecimalRoundingMode
+- doc/classes/LongDecimal.src/M000076.html
+- doc/classes/LongDecimal.src/M000077.html
+- doc/classes/LongDecimal.src/M000078.html
+- doc/classes/LongDecimal.src/M000079.html
+- doc/classes/LongDecimal.src/M000080.html
+- doc/classes/LongDecimal.src/M000081.html
+- doc/classes/LongDecimal.src/M000082.html
+- doc/classes/LongDecimal.src/M000083.html
+- doc/classes/LongDecimal.src/M000084.html
+- doc/classes/LongDecimal.src/M000085.html
+- doc/classes/LongDecimal.src/M000086.html
+- doc/classes/LongDecimal.src/M000087.html
+- doc/classes/LongDecimal.src/M000088.html
+- doc/classes/LongDecimal.src/M000089.html
+- doc/classes/LongDecimal.src/M000090.html
+- doc/classes/LongDecimal.src/M000091.html
+- doc/classes/LongDecimal.src/M000092.html
+- doc/classes/LongDecimal.src/M000093.html
+- doc/classes/LongDecimal.src/M000094.html
+- doc/classes/LongDecimal.src/M000095.html
+- doc/classes/LongDecimal.src/M000096.html
+- doc/classes/LongDecimal.src/M000097.html
+- doc/classes/LongDecimal.src/M000098.html
+- doc/classes/LongDecimal.src/M000099.html
+- doc/classes/LongDecimal.src/M000040.html
+- doc/classes/LongDecimal.src/M000041.html
+- doc/classes/LongDecimal.src/M000100.html
+- doc/classes/LongDecimal.src/M000042.html
+- doc/classes/LongDecimal.src/M000101.html
+- doc/classes/LongDecimal.src/M000043.html
+- doc/classes/LongDecimal.src/M000102.html
+- doc/classes/LongDecimal.src/M000044.html
+- doc/classes/LongDecimal.src/M000103.html
+- doc/classes/LongDecimal.src/M000045.html
+- doc/classes/LongDecimal.src/M000104.html
+- doc/classes/LongDecimal.src/M000046.html
+- doc/classes/LongDecimal.src/M000047.html
+- doc/classes/LongDecimal.src/M000048.html
+- doc/classes/LongDecimal.src/M000107.html
+- doc/classes/LongDecimal.src/M000049.html
+- doc/classes/LongDecimal.src/M000108.html
+- doc/classes/LongDecimal.src/M000109.html
+- doc/classes/LongDecimal.src/M000050.html
+- doc/classes/LongDecimal.src/M000051.html
+- doc/classes/LongDecimal.src/M000110.html
+- doc/classes/LongDecimal.src/M000052.html
+- doc/classes/LongDecimal.src/M000053.html
+- doc/classes/LongDecimal.src/M000054.html
+- doc/classes/LongDecimal.src/M000055.html
+- doc/classes/LongDecimal.src/M000056.html
+- doc/classes/LongDecimal.src/M000057.html
+- doc/classes/LongDecimal.src/M000058.html
+- doc/classes/LongDecimal.src/M000059.html
+- doc/classes/LongDecimal.src/M000060.html
+- doc/classes/LongDecimal.src/M000062.html
+- doc/classes/LongDecimal.src/M000063.html
+- doc/classes/LongDecimal.src/M000064.html
+- doc/classes/LongDecimal.src/M000066.html
+- doc/classes/LongDecimal.src/M000067.html
+- doc/classes/LongDecimal.src/M000068.html
+- doc/classes/LongDecimal.src/M000069.html
+- doc/classes/LongDecimal.src/M000070.html
+- doc/classes/LongDecimal.src/M000071.html
+- doc/classes/LongDecimal.src/M000072.html
+- doc/classes/LongDecimal.src/M000073.html
+- doc/classes/LongDecimal.src/M000074.html
+- doc/classes/LongDecimal.src/M000075.html
+- doc/classes/Numeric.src/M000111.html
+- doc/classes/LongDecimalQuot.src/M000008.html
+- doc/classes/LongDecimalQuot.src/M000009.html
+- doc/classes/LongDecimalQuot.src/M000010.html
+- doc/classes/LongDecimalQuot.src/M000011.html
+- doc/classes/LongDecimalQuot.src/M000012.html
+- doc/classes/LongDecimalQuot.src/M000013.html
+- doc/classes/LongDecimalQuot.src/M000014.html
+- doc/classes/LongDecimalQuot.src/M000015.html
+- doc/classes/LongDecimalQuot.src/M000016.html
+- doc/classes/LongDecimalQuot.src/M000017.html
+- doc/classes/LongDecimalQuot.src/M000018.html
+- doc/classes/LongDecimalQuot.src/M000019.html
+- doc/classes/LongDecimalQuot.src/M000020.html
+- doc/classes/LongDecimalQuot.src/M000021.html
+- doc/classes/LongDecimalQuot.src/M000022.html
+- doc/classes/LongDecimalQuot.src/M000023.html
+- doc/classes/LongDecimalQuot.src/M000024.html
+- doc/classes/LongDecimalQuot.src/M000025.html
+- doc/classes/LongDecimalQuot.src/M000026.html
+- doc/classes/LongDecimalQuot.src/M000027.html
+- doc/classes/LongDecimalQuot.src/M000028.html
+- doc/classes/LongDecimalQuot.src/M000029.html
+- doc/classes/LongDecimalQuot.src/M000030.html
+- doc/classes/LongDecimalQuot.src/M000031.html
+- doc/classes/LongDecimalQuot.src/M000032.html
+- doc/classes/LongDecimalQuot.src/M000033.html
+- doc/classes/LongDecimalQuot.src/M000034.html
+- doc/classes/LongDecimalQuot.src/M000037.html
+- doc/classes/LongDecimalQuot.src/M000038.html
+- doc/classes/LongDecimalQuot.src/M000039.html
+- doc/classes/LongDecimalQuot.src/M000003.html
+- doc/classes/LongDecimalQuot.src/M000004.html
+- doc/classes/LongDecimalQuot.src/M000005.html
+- doc/classes/LongDecimalQuot.src/M000006.html
+- doc/classes/LongDecimalQuot.src/M000007.html
+- doc/classes/LongMath.src/M000112.html
+- doc/classes/LongMath.src/M000113.html
+- doc/classes/LongMath.src/M000114.html
+- doc/classes/LongMath.src/M000115.html
+- doc/classes/LongMath.src/M000116.html
+- doc/classes/LongMath.src/M000117.html
+- doc/classes/LongMath.src/M000118.html
+- doc/classes/LongMath.src/M000119.html
+- doc/classes/LongMath.src/M000120.html
+- doc/classes/LongMath.src/M000121.html
+- doc/classes/LongDecimalRoundingMode/RoundingModeClass.src
+- doc/classes/LongDecimalRoundingMode/RoundingModeClass.html
+- doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000122.html
 - lib/longdecimal.rb
 - test/testlongdecimal.rb
 test_files: