long-decimal 0.01.01 → 0.01.02

data/test/testlongdeclib.rb

@@ -2,8 +2,8 @@
 #
 # library for testlongdecimal.rb
 #
-# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testlongdeclib.rb,v 1.15 2006/05/01 12:22:12 bk1 Exp $
-# CVS-Label: $Name: ALPHA_01_01 $
+# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testlongdeclib.rb,v 1.18 2007/06/19 01:36:23 bk1 Exp $
+# CVS-Label: $Name: ALPHA_01_02 $
 # Author:    $Author: bk1 $ (Karl Brodowsky)
 #
 
@@ -14,7 +14,7 @@ load "lib/long-decimal.rb"
 #
 module TestLongDecHelper
 
-  @RCS_ID='-$Id: testlongdeclib.rb,v 1.15 2006/05/01 12:22:12 bk1 Exp $-'
+  @RCS_ID='-$Id: testlongdeclib.rb,v 1.18 2007/06/19 01:36:23 bk1 Exp $-'
 
 
   #
@@ -146,6 +146,7 @@ module TestLongDecHelper
     # 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}")
     assert_small_interval(yd, yu, y, "the result y/yu  should not depend on the internal rounding mode x0=#{x0} x=#{x} p=#{prec} d=#{(yd-yu).to_f.to_s}")
+    y
   end
 
   #
@@ -216,6 +217,7 @@ module TestLongDecHelper
       delta = z.unit * df
       assert((x - z).abs <= delta, "x=#{x.to_s} and z=#{z.to_s} should be almost equal (#{(x-z).abs.inspect} < d=#{delta.inspect} y=#{y.to_s} lprec=#{lprec} prec=#{prec})")
     end
+    y
 
   end
 
@@ -264,6 +266,7 @@ module TestLongDecHelper
       delta = z.unit * df
       assert((x - z).abs <= delta, "x=#{x.to_s} and z=#{z.to_s} should be almost equal (#{(x-z).abs.inspect} < d=#{delta.inspect} y=#{y.to_s} lprec=#{lprec} prec=#{prec})")
     end
+    y
 
   end
 
@@ -313,6 +316,31 @@ module TestLongDecHelper
     assert((zf - zz).abs < [ z.unit.to_f, zf.abs / 1e9 ].max, "zf=#{zf.to_s} and zz=#{zz.to_s} should be almost equal x=#{x} y=#{y} z=#{z} zz=#{zz}")
   end
 
+  #
+  # helper method for test_lm_power_yint
+  # tests if exp(x) with precision prec is calculated correctly
+  #
+  def check_power_yhalfint(x, y2, prec)
+
+    print "."
+    $stdout.flush
+
+    y2i = y2.to_i
+    x = x.to_ld
+    y = LongDecimal("0.5") * y2i
+
+    z  = LongMath.power(x, y, prec, LongMath::ROUND_HALF_UP)
+    zz = LongMath.sqrt(x ** y2i, prec, LongMath::ROUND_HALF_UP)
+    assert_equal(z, zz, "power with ** or power-method x=#{x} y=#{y} z=#{z} zz=#{zz}")
+    zz = LongMath.power_internal(x, y, prec)
+    assert((zz - z).abs <= z.unit, "power with and without optimizations x=#{x} y=#{y} z=#{z} zz=#{zz}")
+
+    zz = (x.to_f) ** (y.to_f)
+    zf = z.to_f
+    # assert((zf - zz).abs / [zf.abs, zz.abs, Float::MIN].max < 1e-9, "z=#{zf.to_s} and zz=#{zz.to_s} should be almost equal x=#{x} y=#{y} z=#{z} zz=#{zz}")
+    assert((zf - zz).abs < [ z.unit.to_f, zf.abs / 1e9 ].max, "zf=#{zf.to_s} and zz=#{zz.to_s} should be almost equal x=#{x} y=#{y} z=#{z} zz=#{zz}")
+  end
+
   #
   # helper method for test_log
   # tests if log(x) with precision prec is calculated correctly
@@ -399,8 +427,8 @@ module TestLongDecHelper
   def check_power_floated(x, y, prec)
 
     print "."
-    # print("\nstart: check_power_floated: x=#{x} y=#{y} prec=#{prec}")
-    # t0 = Time.new
+    t0 = Time.new
+    # print("\nstart: check_power_floated: x=#{x} y=#{y} prec=#{prec} t0=#{t0} ")
     $stdout.flush
 
     # make sure x and y are LongDecimal
@@ -410,13 +438,26 @@ module TestLongDecHelper
     y = y.to_ld
     # calculate z = x**y
     z = LongMath.power(x, y, prec)
+    # puts("z=#{z}")
 
-    if (z.abs < LongMath::MAX_FLOATABLE)
+    if (z.abs < LongMath::MAX_FLOATABLE && (x - 1).abs*1000000000 > 1)
       # compare y against w = x**y calculated using regular floating point arithmetic
-      w = (x.to_f) ** (y.to_f)
+      xf = x.to_f
+      yf = y.to_f
+      wf = (xf) ** (yf)
       zf = z.to_f
-      # assert((zf - w).abs / [zf.abs, w.abs, Float::MIN].max < 1e-9, "z=#{zf.to_s} and w=#{w.to_s} should be almost equal x=#{x} y=#{y}")
-      assert((zf - w).abs <= [ z.unit, zf.abs / 1e9 ].max, "z=#{zf.to_s} and w=#{w.to_s} should be almost equal x=#{x} y=#{y}")
+      # puts("zf=#{zf} w=#{w}")
+      # assert((zf - w).abs / [zf.abs, w.abs, Float::MIN].max < 1e-9, "z=#{zf.to_s} and w=#{w.to_s} should be almost equal x=#{x} y=#{y} prec=#{prec}")
+      qf = 1e9
+      delta = [ z.unit, zf.abs / qf ].max
+      if (yf.abs > 1) 
+	l = Math.log(yf.abs)
+	if (l > 1)
+	  delta *= l
+        end
+      end
+      diff  = (zf - wf).abs
+      assert(diff <= delta, "zf=#{zf.to_s} and wf=#{wf.to_s} should be almost equal x=#{x} y=#{y} delta=#{delta} l=#{l} diff=#{diff} prec=#{prec}")
     end
 
     # check by taking log(z) = y * log(x)
@@ -445,7 +486,6 @@ module TestLongDecHelper
       yv = (y*v).round_to_scale(lprec, LongDecimal::ROUND_HALF_DOWN)
       assert((u - yv).abs <= unit, "u=#{u} and yv=y*v=#{yv} should be almost equal (unit=#{unit} x=#{x.to_s} y=#{y.to_s} z=#{z.to_s} u=#{u.to_s} v=#{v.to_s} lprec=#{lprec} prec=#{prec})")
     end
-    # puts("ok check_power_floated: x=#{x} y=#{y} prec=#{prec} t=#{Time.new - t0}\n")
 
   end
 
@@ -510,7 +550,7 @@ module TestLongDecHelper
       assert((x - z).abs <= delta, "x=#{x.to_s} and z=#{z.to_s} should be almost equal (y=#{y.to_s} eprec=#{eprec} prec=#{prec} delta=#{delta})")
       # assert((x - z).abs <= z.unit, "x=#{x.to_s} and z=#{z.to_s} should be almost equal (y=#{y.to_s} eprec=#{eprec} prec=#{prec})")
     end
-    
+
     return y
 
   end

data/test/testrandlib.rb

@@ -2,8 +2,8 @@
 #
 # testrandlib.rb -- runit test for long-decimal.rb
 #
-# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandlib.rb,v 1.3 2006/04/30 20:49:54 bk1 Exp $
-# CVS-Label: $Name: ALPHA_01_01 $
+# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandlib.rb,v 1.7 2006/05/07 20:00:24 bk1 Exp $
+# CVS-Label: $Name: ALPHA_01_02 $
 # Author:    $Author: bk1 $ (Karl Brodowsky)
 #
 
@@ -14,7 +14,7 @@ require "crypt/ISAAC"
 #
 module TestRandomHelper
 
-  @RCS_ID='-$Id: testrandlib.rb,v 1.3 2006/04/30 20:49:54 bk1 Exp $-'
+  @RCS_ID='-$Id: testrandlib.rb,v 1.7 2006/05/07 20:00:24 bk1 Exp $-'
 
   @@r1 = Crypt::ISAAC.new
   @@r2 = Crypt::ISAAC.new
@@ -24,9 +24,10 @@ module TestRandomHelper
   def arr_with_sq(arr)
     result = arr.map do |x|
       if (x >= 0) then
-	[ LongMath.sqrt(x, x.scale, LongMath::ROUND_HALF_UP), x, x.square ]
+        r = LongMath.sqrt(x, x.scale, LongMath::ROUND_HALF_UP)
+        [ r, x, r*x, x.square ]
       else
-	x
+        x
       end
     end
     result.flatten
@@ -72,13 +73,28 @@ module TestRandomHelper
     x5 = @@r4.rand(1000)
     x6 = @@r4.rand(100)+5
     xs = LongDecimal(x1, x2) + LongDecimal(x3, x4) + LongDecimal(x5, x6)
-    xm = 1 + xs
+    xt = (xs / LongMath.pi(xs.scale)).to_ld(xs.scale*2)
+    xu = 1 + xs
+    xv = 1 + xt
+    xw = 2 + xs
+    xe = LongMath.exp(1, xs.scale) + xs
+    xp = LongMath.pi(xs.scale) + xs
+    xn = 10 + xs
     xb = x0 + xs
+    xg = LongMath::MAX_EXP_ABLE - xs
+    xh = LongMath::MAX_EXP_ABLE + xs
     eprec = @@r1.rand(120)
     lprec = eprec+1
     sprec = eprec+((xb.scale+1)>>1)
     pprec = lprec >> 1
-    arr = arr_with_neg(arr_with_inv(arr_with_sq([xb, xm, xs]))).sort.uniq
+    arr1 = [ xt, xs, xu, xv, xw, xe, xp, xn, xb, xg, xh ]
+#     f    = LongDecimal("0.25")
+#     arr2 = (0..8).map do |j|
+#       ff = f*j
+#       [ ff + xs, ff**5 + xs ]
+#     end
+#     arr3 = [ arr1, arr2 ].flatten
+    arr  = arr_with_neg(arr_with_inv(arr_with_sq(arr1))).sort.uniq
     return [ arr, eprec, lprec, sprec, pprec, xb.scale ]
   end
 

data/test/testrandom.rb

@@ -2,8 +2,8 @@
 #
 # testrandom.rb -- runit test for long-decimal.rb
 #
-# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandom.rb,v 1.10 2006/04/30 20:49:54 bk1 Exp $
-# CVS-Label: $Name: ALPHA_01_01 $
+# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandom.rb,v 1.12 2006/05/05 22:15:36 bk1 Exp $
+# CVS-Label: $Name: ALPHA_01_02 $
 # Author:    $Author: bk1 $ (Karl Brodowsky)
 #
 
@@ -22,7 +22,7 @@ class TestRandom_class < RUNIT::TestCase
   include TestLongDecHelper
   include TestRandomHelper
 
-  @RCS_ID='-$Id: testrandom.rb,v 1.10 2006/04/30 20:49:54 bk1 Exp $-'
+  @RCS_ID='-$Id: testrandom.rb,v 1.12 2006/05/05 22:15:36 bk1 Exp $-'
 
   # for how many seconds should this test run? change to different
   # value on demand
@@ -32,7 +32,7 @@ class TestRandom_class < RUNIT::TestCase
     arr.each do |x|
       @scnt += 1
       puts("\ncnt=#{cnt} scnt=#{@scnt} x=#{x} ep=#{eprec} lp=#{lprec} sp=#{sprec} pp=#{pprec}\n")
-      if (x.abs <= LongMath::MAX_EXP_ABLE)
+      if (x <= LongMath::MAX_EXP_ABLE) then
 	check_exp_floated(x, eprec)
 	check_exp2_floated(x, pprec)
 	check_exp10_floated(x, pprec)

data/test/testrandpower.rb

@@ -2,8 +2,8 @@
 #
 # testrandpower.rb -- runit test for long-decimal.rb
 #
-# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandpower.rb,v 1.7 2006/05/01 12:22:12 bk1 Exp $
-# CVS-Label: $Name: ALPHA_01_01 $
+# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandpower.rb,v 1.9 2006/05/11 20:57:22 bk1 Exp $
+# CVS-Label: $Name: ALPHA_01_02 $
 # Author:    $Author: bk1 $ (Karl Brodowsky)
 #
 
@@ -22,7 +22,7 @@ class TestRandomPower_class < RUNIT::TestCase
   include TestLongDecHelper
   include TestRandomHelper
 
-  @RCS_ID='-$Id: testrandpower.rb,v 1.7 2006/05/01 12:22:12 bk1 Exp $-'
+  @RCS_ID='-$Id: testrandpower.rb,v 1.9 2006/05/11 20:57:22 bk1 Exp $-'
 
   # for how many seconds should this test run? change to different
   # value on demand
@@ -41,6 +41,7 @@ class TestRandomPower_class < RUNIT::TestCase
       break if d >= @@duration
       xarr, eprec1, lprec1, sprec1, pprec1, xs = random_arr
       yarr, eprec2, lprec2, sprec2, pprec2, ys = random_arr
+
       prec = (pprec1+pprec2) >> 1
       xarr.each do |x|
         puts(" t=#{Time.new-t1} x=#{x}: ")
@@ -48,15 +49,16 @@ class TestRandomPower_class < RUNIT::TestCase
         # next if (x < 1.1)
         # next if (x < 1)
         # next if (x < 0.99)
+	next if x.abs > LongMath::MAX_FLOATABLE || (1/x).abs > LongMath::MAX_FLOATABLE
         yarr.each do |y|
-          # next if x < 0.001 && y < 0
-          next if y.abs > LongMath::MAX_EXP_ABLE * 2
-          next if ((1+x.to_i) ** y.to_i.abs).abs > LongMath::MAX_FLOATABLE
-          next if ((1+(1/x).to_i) ** y.to_i.abs).abs > LongMath::MAX_FLOATABLE
-          next if (Math.log(x) * y.to_f).abs > LongMath::MAX_EXP_ABLE * 2
           puts("\ncnt=#{cnt} scnt=#{scnt} x=#{x} y=#{y} scx=#{x.scale} scy=#{y.scale} prec=#{prec}")
+          # next if x < 0.001 && y < 0
+          next if x > 1 && y > LongMath::MAX_EXP_ABLE * 2
+          next if x < 1 && y < -LongMath::MAX_EXP_ABLE * 2
+          next if Math.log(1+x.to_i) * y.to_i > LongMath::MAX_EXP_ABLE
+          next if -Math.log(1+(1/x).to_i) * y.to_i > LongMath::MAX_EXP_ABLE
+          next if (Math.log(x.to_f) * y.to_f).abs > LongMath::MAX_EXP_ABLE * 2
           scnt += 1
-          # puts(" t=#{Time.new-t1} x=#{x} y=#{y}: ")
           check_power_floated(x, y, prec)
         end
       end

metadata

@@ -1,10 +1,10 @@
 --- !ruby/object:Gem::Specification 
-rubygems_version: 0.8.11
+rubygems_version: 0.9.4
 specification_version: 1
 name: long-decimal
 version: !ruby/object:Gem::Version 
-  version: 0.01.01
-date: 2006-05-01 00:00:00 +02:00
+  version: 0.01.02
+date: 2007-07-13 00:00:00 +02:00
 summary: LongDecimal for numbers with fixed point
 require_paths: 
 - lib
@@ -25,6 +25,7 @@ required_ruby_version: !ruby/object:Gem::Version::Requirement
 platform: ruby
 signing_key: 
 cert_chain: 
+post_install_message: 
 authors: 
 - Karl Brodowsky
 files: 
@@ -58,6 +59,8 @@ files:
 - doc/files/lib/long-decimal_rb.src/M000001.html
 - doc/files/lib/long-decimal_rb.src/M000002.html
 - doc/classes/LongDecimal.src
+- doc/classes/Integer.src
+- doc/classes/Integer.html
 - doc/classes/LongMath.html
 - doc/classes/Numeric.src
 - doc/classes/Rational.src
@@ -112,10 +115,16 @@ files:
 - 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/M000105.html
 - doc/classes/LongDecimal.src/M000047.html
+- doc/classes/LongDecimal.src/M000106.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/M000052.html
@@ -142,15 +151,22 @@ files:
 - doc/classes/LongDecimal.src/M000073.html
 - doc/classes/LongDecimal.src/M000074.html
 - doc/classes/LongDecimal.src/M000075.html
+- doc/classes/Integer.src/M000036.html
+- doc/classes/Integer.src/M000039.html
 - doc/classes/Numeric.src/M000100.html
 - doc/classes/Numeric.src/M000102.html
 - doc/classes/Numeric.src/M000103.html
 - doc/classes/Numeric.src/M000104.html
 - doc/classes/Numeric.src/M000105.html
+- doc/classes/Numeric.src/M000106.html
+- doc/classes/Numeric.src/M000110.html
+- doc/classes/Numeric.src/M000111.html
 - doc/classes/Rational.src/M000101.html
 - doc/classes/Rational.src/M000103.html
 - doc/classes/Rational.src/M000104.html
 - doc/classes/Rational.src/M000106.html
+- doc/classes/Rational.src/M000107.html
+- doc/classes/Rational.src/M000112.html
 - doc/classes/LongDecimalQuot.src/M000008.html
 - doc/classes/LongDecimalQuot.src/M000009.html
 - doc/classes/LongDecimalQuot.src/M000010.html
@@ -201,12 +217,18 @@ files:
 - doc/classes/LongDecimalBase.src/M000115.html
 - doc/classes/LongDecimalBase.src/M000116.html
 - doc/classes/LongDecimalBase.src/M000117.html
+- doc/classes/LongDecimalBase.src/M000118.html
+- doc/classes/LongDecimalBase.src/M000121.html
+- doc/classes/LongDecimalBase.src/M000122.html
+- doc/classes/LongDecimalBase.src/M000123.html
 - doc/classes/LongMath/CacheKey.html
 - doc/classes/LongMath/CacheKey.src
 - doc/classes/LongMath/CacheKey.src/M000140.html
 - doc/classes/LongMath/CacheKey.src/M000145.html
 - doc/classes/LongMath/CacheKey.src/M000147.html
 - doc/classes/LongMath/CacheKey.src/M000149.html
+- doc/classes/LongMath/CacheKey.src/M000155.html
+- doc/classes/LongMath/CacheKey.src/M000160.html
 - doc/classes/LongMath.src/M000135.html
 - doc/classes/LongMath.src/M000136.html
 - doc/classes/LongMath.src/M000137.html
@@ -221,6 +243,17 @@ files:
 - doc/classes/LongMath.src/M000146.html
 - doc/classes/LongMath.src/M000147.html
 - doc/classes/LongMath.src/M000148.html
+- doc/classes/LongMath.src/M000149.html
+- doc/classes/LongMath.src/M000150.html
+- doc/classes/LongMath.src/M000151.html
+- doc/classes/LongMath.src/M000152.html
+- doc/classes/LongMath.src/M000153.html
+- doc/classes/LongMath.src/M000154.html
+- doc/classes/LongMath.src/M000155.html
+- doc/classes/LongMath.src/M000156.html
+- doc/classes/LongMath.src/M000157.html
+- doc/classes/LongMath.src/M000158.html
+- doc/classes/LongMath.src/M000159.html
 - doc/classes/LongMath.src/M000112.html
 - doc/classes/LongMath.src/M000113.html
 - doc/classes/LongMath.src/M000114.html
@@ -245,6 +278,8 @@ files:
 - doc/classes/LongMath.src/M000133.html
 - doc/classes/LongMath.src/M000134.html
 - doc/classes/LongDecimalRoundingMode/RoundingModeClass.src
+- doc/classes/LongDecimalRoundingMode/ZeroRoundingModeClass.src
+- doc/classes/LongDecimalRoundingMode/ZeroRoundingModeClass.html
 - doc/classes/LongDecimalRoundingMode/RoundingModeClass.html
 - doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000135.html
 - doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000137.html
@@ -258,6 +293,16 @@ files:
 - doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000150.html
 - doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000151.html
 - doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000152.html
+- doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000156.html
+- doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000157.html
+- doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000158.html
+- doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000159.html
+- doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000161.html
+- doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000162.html
+- doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000163.html
+- doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000164.html
+- doc/classes/LongDecimalRoundingMode/ZeroRoundingModeClass.src/M000165.html
+- doc/classes/LongDecimalRoundingMode/ZeroRoundingModeClass.src/M000166.html
 - lib/long-decimal.rb
 - test/testlongdecimal.rb
 - test/testrandpower.rb