intervals 0.3.63 → 0.4.75

data/test/test_struct_float.rb

@@ -0,0 +1,72 @@
+#!/usr/bin/ruby
+
+# Tests for the Struct::Float class.
+#
+# Copyright (c) 2006 Stefano Taschini <taschini@ieee.org>.
+# Licensed under the terms of LGPL: http://www.gnu.org/copyleft/lesser.html
+
+require 'test/unit'
+require File.dirname(__FILE__) + '/../lib/struct_float.rb' if __FILE__ == $0
+
+# Tests Struct::Float.
+class Struct::Float::Test < Test::Unit::TestCase
+
+  def assert_consistent(f,s)
+    assert_equal(f,s.to_f)
+    assert_equal(s,f.to_struct)
+    assert_equal(s, eval(s.inspect))
+    assert(f.ulp.infinite? || (f/f.ulp).to_i == f/f.ulp, [f,f.ulp].inspect)
+  end
+
+  def assert_nan(s)
+    # NaN does not have a unique representation, and therefore
+    # it does not make sense to have it as a constant as for
+    # Infinity. In particular a NaN has 2047 as biased exponent and
+    # non-zero fraction
+    assert(2047, s.biased_exp)
+    assert(!s.fraction.zero?, s.fraction)
+    # You cannot test for NaN with equality. You MUST use nan?
+    assert(s.to_f.nan?, s.to_f)
+    assert_equal(s, eval(s.inspect))
+  end
+
+  def test_consistency
+    assert_consistent Infinity, Struct::Float[1,2047,0]
+    assert_consistent -Infinity, Struct::Float[-1,2047,0]
+    assert_consistent 2.0 ** -1074, Struct::Float[1,0,1]
+    assert_consistent 0.0, Struct::Float[1,0,0]
+    assert_consistent -0.0, Struct::Float[-1,0,0]
+    assert_consistent 1.0, Struct::Float[1,1023, 0]
+    assert_consistent 1/3.0, Struct::Float[1,1021,0x5555555555555]
+    assert_consistent 1/5.0, Struct::Float[1,1020,0x999999999999a]
+  end
+
+  def test_nan
+    assert_nan( (0.0/0.0).to_struct)
+    assert_nan( Struct::Float[-1,2047,34])
+    assert_nan( Struct::Float[+1,2047,34])
+  end
+
+  def test_denormalized
+    assert_consistent 34 * 2.0 ** -(1023+51), Struct::Float[1,0,34]
+  end
+
+  def test_normalized
+    assert_consistent -(2**52+34) * 2.0**(1043-1023-52), Struct::Float[-1,1043,34]
+  end
+
+  def test_samples
+    n = 1000
+    srand 12345
+    n.times{
+      s = Struct::Float[2*rand(2) -1, rand(2048), rand(2**52)]
+      f = s.to_f
+      if f.nan?
+        assert_nan(s)
+      else
+        assert_consistent(f,s)
+      end
+    }
+  end
+
+end

data/var/clean-nan.rb

@@ -0,0 +1,59 @@
+#!/usr/bin/ruby
+#
+# Alternative implementation of Interval::Simple.multiply that does
+# not rely on catching exceptions to detect NaNs
+#
+# Copyright (c) 2006 Stefano Taschini <taschini@ieee.org>.
+# Licensed under the terms of LGPL: http://www.gnu.org/copyleft/lesser.html
+
+require File.dirname(__FILE__) + '/../lib/interval'
+
+module Enumerable
+
+  # Enumerable#min fails if the collection includes a NaN. Instead,
+  # this methods returns a NaN.
+  def ieee_min
+    find( proc { min } ){|x| x.respond_to?(:nan?) && x.nan?}
+  end
+
+  # Enumerable#max fails if the collection includes a NaN. Instead,
+  # this methods returns a NaN.
+  def ieee_max
+    find( proc { max } ){|x| x.respond_to?(:nan?) && x.nan?}
+  end
+
+end
+
+class Interval::Simple
+
+  # Used to implement Interval's multiplicative operator (\*).
+  def multiply (other)
+    self.class.new(
+      FPU.down {
+        [inf * other.inf, inf * other.sup,
+         sup * other.inf, sup * other.sup]}.ieee_min,
+      FPU.up {
+        [inf * other.inf, inf * other.sup,
+         sup * other.inf, sup * other.sup]}.ieee_max)
+  end
+
+end
+
+# Tests exception-free multiplication.
+class CleanNanTest < Test::Unit::TestCase
+
+  def assert_nan(s)
+    assert(s.nan?, "#{s} is not nan")
+  end
+
+  def test_ieee_min_and_max
+    assert_equal(1,[1,2,3].ieee_min)
+    assert_equal(3,[1,2,3].ieee_max)
+    assert_nan([1,0.0/0.0,3].ieee_min)
+    assert_nan([1,0.0/0.0,3].ieee_max)
+  end
+
+end
+
+Test::Unit.run = (__FILE__ != $0)
+

metadata

@@ -3,8 +3,8 @@ rubygems_version: 0.8.11
 specification_version: 1
 name: intervals
 version: !ruby/object:Gem::Version 
-  version: 0.3.63
-date: 2006-03-31 00:00:00 +02:00
+  version: 0.4.75
+date: 2006-04-09 00:00:00 +02:00
 summary: Interval arithmetic in Ruby.
 require_paths: 
 - lib
@@ -28,14 +28,18 @@ cert_chain:
 authors: 
 - Stefano Taschini
 files: 
+- test/runall.rb
 - test/data_sin.txt
 - test/data_cosh.txt
 - test/data_exp.txt
+- test/test_interval.rb
 - test/data_log.txt
 - test/data_sinh.txt
 - test/data_tan.txt
 - test/data_atan.txt
+- test/test_struct_float.rb
 - test/data_cos.txt
+- test/test_fpu.rb
 - ext/jamis-mod.rb
 - ext/extconf.rb
 - ext/fpu.c
@@ -162,10 +166,11 @@ files:
 - lib/interval.rb
 - lib/struct_float.rb
 - lib/fpu.rb
+- var/clean-nan.rb
 - README.txt
 - VERSION.txt
 test_files: 
-- lib/interval.rb
+- test/runall.rb
 rdoc_options: 
 - --inline-source
 - --template