long-decimal 0.01.03 → 0.02.01

data/test/testlongdeclib.rb

@@ -1,21 +1,33 @@
 #!/usr/bin/env ruby
 #
-# library for testlongdecimal.rb
+# library for tests for long-decimal.rb and long-decimal-extra.rb
 #
-# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testlongdeclib.rb,v 1.19 2007/08/19 19:25:59 bk1 Exp $
-# CVS-Label: $Name: ALPHA_01_03 $
+# (C) Karl Brodowsky (IT Sky Consulting GmbH) 2006-2009
+#
+# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testlongdeclib.rb,v 1.29 2009/04/21 16:56:49 bk1 Exp $
+# CVS-Label: $Name: BETA_02_01 $
 # Author:    $Author: bk1 $ (Karl Brodowsky)
 #
 
 # load "lib/long-decimal.rb"
 
+class Integer
+  def cube
+    self*self*self
+  end
+
+  def square
+    self*self
+  end
+end
+
+
 #
 # test class for LongDecimal and LongDecimalQuot
 #
 module TestLongDecHelper
 
-  @RCS_ID='-$Id: testlongdeclib.rb,v 1.19 2007/08/19 19:25:59 bk1 Exp $-'
-
+  @RCS_ID='-$Id: testlongdeclib.rb,v 1.29 2009/04/21 16:56:49 bk1 Exp $-'
 
   #
   # convenience method for comparing two numbers. true if and only if
@@ -405,7 +417,6 @@ module TestLongDecHelper
       v = y.unit
       # delta = (u + u.move_point_left(1)) * df
       delta = [ v * z, u * df ].max
-      # puts("x=#{x.to_s} and z=#{z.to_s} should be almost equal (#{(x-z).abs} < d=#{delta} y=#{y.to_s} eprec=#{eprec} prec=#{prec} u=#{u.inspect} v=#{v.inspect} df=#{df})")
       assert((x - z).abs <= delta, "x=#{x.to_s} and z=#{z.to_s} should be almost equal (#{(x-z).abs} < d=#{delta} y=#{y.to_s} eprec=#{eprec} prec=#{prec})")
     end
 
@@ -428,7 +439,6 @@ module TestLongDecHelper
 
     print "."
     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
@@ -438,24 +448,31 @@ 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 && (x - 1).abs*1000000000 > 1)
+    corr2 = (x - 1).abs*1000000000
+    if (z.abs < LongMath::MAX_FLOATABLE && corr2 > 1)
       # compare y against w = x**y calculated using regular floating point arithmetic
       xf = x.to_f
       yf = y.to_f
       wf = (xf) ** (yf)
       zf = z.to_f
-      # 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
+      delta = [ z.unit.to_f, zf.abs / qf ].max
+      # puts "delta=#{delta} z=#{z} zu=#{z.unit} zuf=#{z.unit.to_f} zf=#{zf} |zf/qf|=#{zf.abs/qf}"
+      if (yf.abs > 1)
+        l = Math.log(yf.abs)
+        if (l > 1)
+          delta *= l
         end
+        # puts "delta=#{delta} l=#{l}"
       end
+      corr = corr2 * 0.5
+      if corr > 1
+        corr_f = [ corr.to_f, 2 ].min.to_f
+        delta *= corr_f
+      end
+      # puts "delta=#{delta} corr_f=#{corr_f} corr=#{corr}"
+
       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
@@ -504,7 +521,6 @@ module TestLongDecHelper
     # calculate y = log10(x)
     y  = LongMath.log10(x, prec)
     yy = LongMath.log10(x, prec+10)
-    # assert_equal(yy.round_to_scale(y.scale, LongDecimal::ROUND_HALF_DOWN), y, "x=#{x} y=#{y} yy=#{yy}")
     assert_equal_rounded(yy, y, "x=#{x} y=#{y} yy=#{yy}")
 
     # compare y against z = log10(x) calculated using regular floating
@@ -519,7 +535,6 @@ module TestLongDecHelper
         df = [ 1e-13, z.abs / 1e10 ].max
         dl = y.unit
         delta = [ df, dl ].max
-        # assert((yf - z).abs / [yf.abs, z.abs, Float::MIN].max < 1e-9, "y=#{yf.to_s} and z=#{z.to_s} should be almost equal x=#{x}")
         assert((yf - z).abs < delta, "y=#{yf.to_s} and z=#{z.to_s} should be almost equal (x=#{x} delta=#{delta}")
       end
     end
@@ -548,7 +563,6 @@ module TestLongDecHelper
       assert((zz - z).abs <= u, "zz=#{zz.to_s} and z=#{z.to_s} should be almost equal (y=#{y.to_s} eprec=#{eprec} prec=#{prec})")
       delta = [ v*z*LongMath::LOG10*1.2, u * 1.1 ].max
       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
@@ -600,7 +614,6 @@ module TestLongDecHelper
         df = [ 1e-13, z.abs / 1e10 ].max
         dl = y.unit.abs
         delta = [ df, dl ].max
-        # assert((yf - z).abs / [yf.abs, z.abs, Float::MIN].max < 1e-9, "y=#{yf.to_s} and z=#{z.to_s} should be almost equal x=#{x}")
         assert((yf - z).abs < delta.to_f, "y=#{yf.to_s} and z=#{z.to_s} should be almost equal (x=#{x} delta=#{delta}")
       end
     end
@@ -654,8 +667,8 @@ module TestLongDecHelper
   #
   def check_sqrtb(x, s)
     y = LongMath.sqrtb(x)
-    z = y * y
-    zz = (y+1)*(y+1)
+    z = y.square
+    zz = (y+1).square
     assert(0 <= y, "sqrt must be >= 0" + s)
     assert(z <= x && x < zz, "y=#{y}=sqrt(#{x}) and x in [#{z}, #{zz})" + s)
     y
@@ -730,6 +743,62 @@ module TestLongDecHelper
     r
   end
 
+  #
+  # helper method for test_int_cbrtb
+  #
+  def check_cbrtb(x, s)
+    y = LongMath.cbrtb(x)
+    z = y.cube
+    zz = (y+1).cube
+    assert(0 <= y, "cbrt must be >= 0" + s)
+    assert(z <= x && x < zz, "y=#{y}=cbrt(#{x}) and x in [#{z}, #{zz})" + s)
+    y
+  end
+
+  #
+  # helper method for test_int_cbrtb_with_remainder
+  #
+  def check_cbrtb_with_remainder(x, s)
+    y, r = LongMath.cbrtb_with_remainder(x)
+    z0 = y.cube
+    # puts "x=#{x} y=#{y} z0=#{z0} r=#{r}"
+    z1 = z0 + r
+    z2 = (y+1).cube
+    assert(0 <= y, "cbrt _with_remainder must be >= 0" + s)
+    assert_equal(z1, x, "x=#{x} y=#{y} r=#{r} z0=#{z0} z1=#{z1} z2=#{z2}" + s)
+    assert(z0 <= x && x < z2, "y=#{y}=cbrt(_with_remainder#{x}) and x in [#{z0}, #{z2}) r=#{r} z0=#{z0} z1=#{z1} z2=#{z2}" + s)
+    y
+  end
+
+  #
+  # helper method of test_cbrt
+  #
+  def check_cbrt(x, scale, mode, su0, su1, str)
+    y  = x.cbrt(scale, mode)
+    if (mode == LongMath::ROUND_HALF_UP || mode == LongMath::ROUND_HALF_DOWN || mode == LongMath::ROUND_HALF_EVEN)
+      yy = x.cbrt(scale+10, mode)
+      assert_equal(yy.round_to_scale(y.scale, mode), y, "x=#{x} y=#{y} yy=#{yy}")
+    end
+    z0 = (y+su0*y.unit).cube
+    z1 = (y+su1*y.unit).cube
+    assert(0 <= y.sign, "cbrt must be >= 0" + str)
+    assert(z0 <= x && x <= z1, "y=#{y}=cbrt(#{x}) and x in [#{z0}, #{z1}) su0=#{su0} su1=#{su1}" + str)
+    y
+  end
+
+  #
+  # helper method of test_cbrt_with_remainder
+  #
+  def check_cbrt_with_remainder(x, scale, str)
+    y, r = x.cbrt_with_remainder(scale)
+    z0 = y.cube
+    z1 = y.succ.cube
+    assert(0 <= y.sign, "cbrt must be >= 0" + str)
+    assert(z0 <= x && x < z1, "y=#{y}=cbrt(#{x}) and x in [#{z0}, #{z1}) " + str)
+    assert((x - z0 - r).zero?, "x=y**3+r x=#{x} z0=#{z0} z1=#{z1} y=#{y} r=#{r} total=#{x - z0 - r} " + str)
+    r
+  end
+
 end
 
 # end of file testlongdecimal.rb

data/test/testrandlib.rb

@@ -1,12 +1,15 @@
 #!/usr/bin/env ruby
 #
-# testrandlib.rb -- runit test for long-decimal.rb
+# testrandlib.rb -- library for random tests for long-decimal.rb and long-decimal-extra.rb
 #
-# 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_03 $
+# (C) Karl Brodowsky (IT Sky Consulting GmbH) 2006-2009
+#
+# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandlib.rb,v 1.9 2009/04/15 19:29:38 bk1 Exp $
+# CVS-Label: $Name: BETA_02_01 $
 # Author:    $Author: bk1 $ (Karl Brodowsky)
 #
 
+#require "rubygems"
 require "crypt/ISAAC"
 
 #
@@ -14,7 +17,7 @@ require "crypt/ISAAC"
 #
 module TestRandomHelper
 
-  @RCS_ID='-$Id: testrandlib.rb,v 1.7 2006/05/07 20:00:24 bk1 Exp $-'
+  @RCS_ID='-$Id: testrandlib.rb,v 1.9 2009/04/15 19:29:38 bk1 Exp $-'
 
   @@r1 = Crypt::ISAAC.new
   @@r2 = Crypt::ISAAC.new

data/test/testrandom-extra.rb

@@ -1,9 +1,11 @@
 #!/usr/bin/env ruby
 #
-# testrandom.rb -- runit test for long-decimal.rb
+# testrandom-extra.rb -- randem tests for long-decimal-extra.rb
 #
-# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandom-extra.rb,v 1.1 2007/08/19 19:25:59 bk1 Exp $
-# CVS-Label: $Name: ALPHA_01_03 $
+# (C) Karl Brodowsky (IT Sky Consulting GmbH) 2006-2009
+#
+# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandom-extra.rb,v 1.2 2009/04/15 19:29:38 bk1 Exp $
+# CVS-Label: $Name: BETA_02_01 $
 # Author:    $Author: bk1 $ (Karl Brodowsky)
 #
 
@@ -25,7 +27,7 @@ class TestRandom_class < RUNIT::TestCase
   include TestLongDecHelper
   include TestRandomHelper
 
-  @RCS_ID='-$Id: testrandom-extra.rb,v 1.1 2007/08/19 19:25:59 bk1 Exp $-'
+  @RCS_ID='-$Id: testrandom-extra.rb,v 1.2 2009/04/15 19:29:38 bk1 Exp $-'
 
   # for how many seconds should this test run? change to different
   # value on demand

data/test/testrandom.rb

@@ -1,15 +1,19 @@
 #!/usr/bin/env ruby
 #
-# testrandom.rb -- runit test for long-decimal.rb
+# testrandom.rb -- random tests for long-decimal.rb
 #
-# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandom.rb,v 1.13 2007/08/19 19:25:59 bk1 Exp $
-# CVS-Label: $Name: ALPHA_01_03 $
+# (C) Karl Brodowsky (IT Sky Consulting GmbH) 2006-2009
+#
+# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandom.rb,v 1.15 2009/04/15 19:29:38 bk1 Exp $
+# CVS-Label: $Name: BETA_02_01 $
 # Author:    $Author: bk1 $ (Karl Brodowsky)
 #
 
 require "runit/testcase"
 require "runit/cui/testrunner"
 require "runit/testsuite"
+
+#require "rubygems"
 require "crypt/ISAAC"
 
 load "lib/long-decimal.rb"
@@ -23,7 +27,7 @@ class TestRandom_class < RUNIT::TestCase
   include TestLongDecHelper
   include TestRandomHelper
 
-  @RCS_ID='-$Id: testrandom.rb,v 1.13 2007/08/19 19:25:59 bk1 Exp $-'
+  @RCS_ID='-$Id: testrandom.rb,v 1.15 2009/04/15 19:29:38 bk1 Exp $-'
 
   # for how many seconds should this test run? change to different
   # value on demand

data/test/testrandpower.rb

@@ -1,9 +1,11 @@
 #!/usr/bin/env ruby
 #
-# testrandpower.rb -- runit test for long-decimal.rb
+# testrandpower.rb -- random tests for power-method from long-decimal-extra.rb
 #
-# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandpower.rb,v 1.10 2007/08/19 19:25:59 bk1 Exp $
-# CVS-Label: $Name: ALPHA_01_03 $
+# (C) Karl Brodowsky (IT Sky Consulting GmbH) 2006-2009
+#
+# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/test/testrandpower.rb,v 1.12 2009/04/18 05:51:14 bk1 Exp $
+# CVS-Label: $Name: BETA_02_01 $
 # Author:    $Author: bk1 $ (Karl Brodowsky)
 #
 
@@ -25,7 +27,7 @@ class TestRandomPower_class < RUNIT::TestCase
   include TestLongDecHelper
   include TestRandomHelper
 
-  @RCS_ID='-$Id: testrandpower.rb,v 1.10 2007/08/19 19:25:59 bk1 Exp $-'
+  @RCS_ID='-$Id: testrandpower.rb,v 1.12 2009/04/18 05:51:14 bk1 Exp $-'
 
   # for how many seconds should this test run? change to different
   # value on demand
@@ -49,13 +51,9 @@ class TestRandomPower_class < RUNIT::TestCase
       xarr.each do |x|
         puts(" t=#{Time.new-t1} x=#{x}: ")
         next if x <= 0
-        # 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|
           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

data/version.rb

@@ -3,8 +3,10 @@
 #
 # version.rb -- extract version information from files
 #
-# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/version.rb,v 1.6 2006/04/11 19:45:40 bk1 Exp $
-# CVS-Label: $Name: ALPHA_01_03 $
+# (C) Karl Brodowsky (IT Sky Consulting GmbH) 2006-2009
+#
+# CVS-ID:    $Header: /var/cvs/long-decimal/long-decimal/version.rb,v 1.8 2009/04/18 05:51:14 bk1 Exp $
+# CVS-Label: $Name: BETA_02_01 $
 # Author:    $Author: bk1 $ (Karl Brodowsky)
 #
 
@@ -29,7 +31,6 @@ ARGV.each do |file|
     str = version
   end
 
-  # print "file=#{file}\nstr=#{str}\nversion=#{version}\nname=#{name}\n";
   print str,"\n"
 end
 

metadata

@@ -1,62 +1,69 @@
 --- !ruby/object:Gem::Specification 
-rubygems_version: 0.9.4
-specification_version: 1
 name: long-decimal
 version: !ruby/object:Gem::Version 
-  version: 0.01.03
-date: 2007-08-28 00:00:00 +02:00
-summary: LongDecimal for numbers with fixed point
-require_paths: 
-- lib
-email: http://www.velofahren.de/mail.html
-homepage: http://long-decimal.rubyforge.org/
-rubyforge_project: long-decimal
-description: ""
-autorequire: long-decimal
-default_executable: 
-bindir: bin
-has_rdoc: true
-required_ruby_version: !ruby/object:Gem::Version::Requirement 
-  requirements: 
-  - - ">"
-    - !ruby/object:Gem::Version 
-      version: 0.0.0
-  version: 
+  version: 0.02.01
 platform: ruby
-signing_key: 
-cert_chain: 
-post_install_message: 
 authors: 
 - Karl Brodowsky
+autorequire: long-decimal
+bindir: bin
+cert_chain: []
+
+date: 2009-04-21 00:00:00 +02:00
+default_executable: 
+dependencies: []
+
+description: ""
+email: http://www.velofahren.de/mail.html
+executables: []
+
+extensions: []
+
+extra_rdoc_files: []
+
 files: 
-- lib
 - test
-- Rakefile
-- README
-- install.rb
-- VERSION
-- make_doc.rb
-- version.rb
-- lib/long-decimal-extra.rb
-- lib/long-decimal.rb
-- test/testlongdecimal.rb
-- test/testrandpower.rb
-- test/testlongdecimal-extra.rb
 - test/testlongdeclib.rb
-- test/testrandom-extra.rb
 - test/testrandlib.rb
+- test/testrandom-extra.rb
+- test/testlongdecimal-extra.rb
 - test/testrandom.rb
-test_files: 
 - test/testlongdecimal.rb
+- test/testrandpower.rb
+- Rakefile
+- install.rb
+- lib
+- lib/long-decimal.rb
+- lib/long-decimal-extra.rb
+- VERSION
+- version.rb
+- make_doc.rb
+- README
+has_rdoc: true
+homepage: http://long-decimal.rubyforge.org/
+post_install_message: 
 rdoc_options: []
 
-extra_rdoc_files: []
-
-executables: []
-
-extensions: []
-
+require_paths: 
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement 
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      version: "0"
+  version: 
+required_rubygems_version: !ruby/object:Gem::Requirement 
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      version: "0"
+  version: 
 requirements: []
 
-dependencies: []
-
+rubyforge_project: long-decimal
+rubygems_version: 1.3.1
+signing_key: 
+specification_version: 2
+summary: LongDecimal for numbers with fixed point
+test_files: 
+- test/testlongdecimal.rb