flt 1.2.1 → 1.3.0

data/lib/flt/version.rb

@@ -1,8 +1,8 @@
 module Flt
   module VERSION #:nodoc:
     MAJOR = 1
-    MINOR = 2
-    TINY  = 1
+    MINOR = 3
+    TINY  = 0
 
     STRING = [MAJOR, MINOR, TINY].join('.')
   end

data/test/generate_trig_data.rb

@@ -0,0 +1,169 @@
+# Generate test data for trigonometry tests (test_trig.rbG)
+
+require File.expand_path(File.join(File.dirname(__FILE__),'helper.rb'))
+require File.dirname(__FILE__) + '/../lib/flt/math'
+include Flt
+
+def around(x)
+  [x,x.next_plus,x.next_minus,x.next_plus.next_plus,x.next_minus.next_minus]
+end
+
+def near(x)
+  ulp = x.ulp
+  prec = x.num_class.context.precision
+  xs = around(x)
+  [2,3,4].each do |k|
+    d = (prec/k)*ulp
+    xs += around(x-d)
+    xs += around(x+d)
+  end
+  xs.uniq
+end
+
+# random angles in radians for sin,cos,tan
+def angle_data(num_class)
+  xs = []
+  pi = num_class::Math.send(:half_cycle)
+  (-8..8).each do |k|
+    x = k*pi/4
+    xs += near(x)
+  end
+  pi2 = num_class::Math.send(:quarter_cycle)
+  50.times{ xs << random_num_one(num_class)*pi2}
+  base = xs.dup
+  (2...10).each do |k|
+    xs += base.map{|x| x+num_class.int_radix_power(k)}
+  end
+  xs.uniq
+end
+
+# random data in [-1,-1] for asin, acos
+def one_data(num_class)
+  half = num_class.Num(1)/2
+  zero = num_class.Num(0)
+  one = num_class.Num(1)
+  xs = [-one, -half, zero, half, one].map{|x| near(x)}.flatten
+  50.times{ xs << random_num_one(num_class)}
+  xs.uniq.reject{|x| x<-1 || x>1}
+end
+
+def random_num_one(num_class=DecNum)
+  context = num_class.context
+  if false # rand(20)==0
+    # generate 5% of subnormals
+    f = rand(context.radix**(context.precision-1))
+    e = context.etiny
+  elsif rand(20)==0
+    # and some singular values too
+    if rand(1) == 0
+      f = context.radix**context.precision - 1
+      f -= rand(3)
+    else
+      f = context.radix**(context.precision - 1)
+      f += rand(3)
+    end
+    # e = random_integer(context.etiny, context.etop)
+    if rand(1)==0
+      e = random_integer(-context.precision-1,-context.precision)
+     else
+      e = random_integer(-context.precision-10,-context.precision)
+    end
+  else
+    f = rand(context.radix**context.precision)
+    # e = random_integer(context.etiny, context.etop)
+    if rand(1)==0
+      e = random_integer(-context.precision-1,-context.precision)
+     else
+      e = random_integer(-context.precision-10,-context.precision)
+    end
+  end
+  # f = -f if rand(1)==0
+  context.Num(f, e)
+end
+
+
+def gen_test_data(num_class, prec, angle_units=:rad)
+  dir = File.dirname(__FILE__)+'/trigtest'
+  num_class.context(:precision=>prec) do
+    num_class.context.angle = angle_units
+    DecNum.context.traps[DecNum::DivisionByZero] = false
+    radix = "#{num_class.context.radix}_"
+    units = "_#{num_class.context.angle}"
+    angles = angle_data(num_class)
+    extra_prec = prec + 100
+    extra_elimit = [10000, num_class.context.emax.abs, num_class.context.emin.abs].max
+    File.open("#{dir}/sin#{radix}#{prec}#{units}.txt","w") do |out|
+      angles.each do |angle|
+        result = +num_class.context(:precision=>extra_prec, :elimit=>extra_elimit) {
+          num_class::Math.sin(angle)
+        }
+        out.puts "#{angle.to_s(:base=>num_class.radix)}\t#{result.to_s(:base=>num_class.radix)}"
+      end
+    end
+    File.open("#{dir}/cos#{radix}#{prec}#{units}.txt","w") do |out|
+      angles.each do |angle|
+        result = +num_class.context(:precision=>extra_prec, :elimit=>extra_elimit) {
+          num_class::Math.cos(angle)
+        }
+        out.puts "#{angle.to_s(:base=>num_class.radix)}\t#{result.to_s(:base=>num_class.radix)}"
+      end
+    end
+    File.open("#{dir}/tan#{radix}#{prec}#{units}.txt","w") do |out|
+      angles.each do |angle|
+        result = +num_class.context(:precision=>extra_prec, :elimit=>extra_elimit) {
+          num_class::Math.tan(angle)
+        }
+        out.puts "#{angle.to_s(:base=>num_class.radix)}\t#{result.to_s(:base=>num_class.radix)}"
+      end
+    end
+    xs = one_data(num_class)
+    File.open("#{dir}/asin#{radix}#{prec}#{units}.txt","w") do |out|
+      xs.each do |x|
+        result = +num_class.context(:precision=>extra_prec, :elimit=>extra_elimit) {
+          num_class::Math.asin(x)
+        }
+        out.puts "#{x.to_s(:base=>num_class.radix)}\t#{result.to_s(:base=>num_class.radix)}"
+      end
+    end
+    File.open("#{dir}/acos#{radix}#{prec}#{units}.txt","w") do |out|
+      xs.each do |x|
+        result = +num_class.context(:precision=>extra_prec, :elimit=>extra_elimit) {
+          num_class::Math.acos(x)
+        }
+        out.puts "#{x.to_s(:base=>num_class.radix)}\t#{result.to_s(:base=>num_class.radix)}"
+      end
+    end
+    xs += Array.new(100){random_num(num_class)}
+    File.open("#{dir}/atan#{radix}#{prec}#{units}.txt","w") do |out|
+      xs.each do |x|
+        result = +num_class.context(:precision=>extra_prec, :elimit=>extra_elimit) {
+          num_class::Math.atan(x)
+        }
+        out.puts "#{x.to_s(:base=>num_class.radix)}\t#{result.to_s(:base=>num_class.radix)}"
+      end
+    end
+  end
+end
+
+srand 12322
+gen_test_data DecNum, 12, :rad
+
+srand 12322
+gen_test_data DecNum, 12, :deg
+
+srand 12322
+BinNum.context = BinNum::IEEEDoubleContext
+gen_test_data BinNum, 53, :rad
+
+
+# TODO:
+# 1. prepare test data for 12, 15, 25, 50 digits with gen_test_data
+# 2. check with RPL, Mathematica:
+#    programs that process test data and prepare RPL/Mathematica programs that yield results
+#    then check results
+#    RPL: (to be executed in emulator for speed/ease of I/O)
+#     \<< [angle1, angle2, angle3, ...] N \-> D N \<< 1 N FOR I D I GET SIN NEXT N \->LIST \>> |>>
+#    Mathematica:
+#      F[x_]:=N[Sin[x],prec]
+#      Map[F, {angle1.to_r.to_s, angle2.to_r.to_s, ...}]
+# 3. add special number tests (nans, infinities)

data/test/test_basic.rb

@@ -101,6 +101,11 @@ class TestBasic < Test::Unit::TestCase
     }
     assert_equal 10, DecNum.context.precision
 
+    DecNum.context(:extra_precision=>4) {
+      assert_equal 14, DecNum.context.precision
+    }
+    assert_equal 10, DecNum.context.precision
+
     DecNum.local_context(DecNum::BasicContext) {
       assert_equal :half_up, DecNum.context.rounding
       assert_equal 9, DecNum.context.precision

metadata

@@ -1,12 +1,13 @@
 --- !ruby/object:Gem::Specification 
 name: flt
 version: !ruby/object:Gem::Version 
+  hash: 27
   prerelease: false
   segments: 
   - 1
-  - 2
-  - 1
-  version: 1.2.1
+  - 3
+  - 0
+  version: 1.3.0
 platform: ruby
 authors: 
 - Javier Goizueta
@@ -14,7 +15,7 @@ autorequire:
 bindir: bin
 cert_chain: []
 
-date: 2010-06-16 00:00:00 +02:00
+date: 2010-06-22 00:00:00 +02:00
 default_executable: 
 dependencies: 
 - !ruby/object:Gem::Dependency 
@@ -25,6 +26,7 @@ dependencies:
     requirements: 
     - - ">="
       - !ruby/object:Gem::Version 
+        hash: 9
         segments: 
         - 2
         - 1
@@ -52,6 +54,7 @@ files:
 - lib/flt/b.rb
 - lib/flt/bigdecimal.rb
 - lib/flt/bin_num.rb
+- lib/flt/complex.rb
 - lib/flt/d.rb
 - lib/flt/dec_num.rb
 - lib/flt/float.rb
@@ -61,6 +64,7 @@ files:
 - lib/flt/support.rb
 - lib/flt/tolerance.rb
 - lib/flt/tolerance/sugar.rb
+- lib/flt/trigonometry.rb
 - lib/flt/version.rb
 - setup.rb
 - tasks/ann.rake
@@ -77,6 +81,7 @@ files:
 - tasks/svn.rake
 - tasks/test.rake
 - test/all_tests.rb
+- test/generate_trig_data.rb
 - test/helper.rb
 - test/reader.rb
 - test/test_basic.rb
@@ -123,6 +128,7 @@ required_ruby_version: !ruby/object:Gem::Requirement
   requirements: 
   - - ">="
     - !ruby/object:Gem::Version 
+      hash: 3
       segments: 
       - 0
       version: "0"
@@ -131,6 +137,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
   requirements: 
   - - ">="
     - !ruby/object:Gem::Version 
+      hash: 3
       segments: 
       - 0
       version: "0"