flt 1.5.3 → 1.5.4

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: ec3f00c70c9eb39227e1b5acf490afd57b9dcedfa3d4a1fe5b935debc6cafcd9
-  data.tar.gz: 94808cc1dfe07687d93ff0c98b247373223683cad4a603ee7bda1ba33534ed79
+  metadata.gz: 884a9616a37e6a5ff27ccd739aad4ad825b04c89b00fc9ff888111b9e386f96a
+  data.tar.gz: 00f414c95fb04488487b1ffc020b35cc0fa9d0e162ffa8cffd7f9ae76cac936b
 SHA512:
-  metadata.gz: b6f9272c9d849c240ff3697730e2465f4fbcbda9ebed66d75c6d36e866a02b5963cd664c0c0583e439747685c34562ed2ab83038f305ed12aa7a708848ec37c7
-  data.tar.gz: 1f1c4783c08ea9c324495d341884aae192e79e657b51eb7ac69ce74ab55fabbd2646c892194463e60041d8273687e57552a43f652973346ade4319b403b41652
+  metadata.gz: 80fd95e11624c2e89789a9bd67cb530b9cde332d2f57b21603dfaa78ce43317032bf0aaf40bda33c2783cb2ef53dc381531d5daf1f3815b57b43511fe72f6895
+  data.tar.gz: 3d1c0e0c41a4b71e875eb0e4af7e285049b5775cec8537c7d06a60ef7ad22aed4f6810a99d5b22e3b0211c47881f65137a6e000e43c205f386c08149c7923802

data/History.txt

@@ -1,3 +1,10 @@
+== 1.5.4 2020-04-08
+
+* Bugfix
+  - Allow scaleb when result is valid, e.g. now these are not invalid results (using float-formats):
+    + IEEE_binary256.context.minimum_normal.log10
+    + IEEE_binary256.context.maximum_finite.log10
+
 == 1.5.3 2020-04-07
 
 * Bugfix

data/lib/flt/num.rb

@@ -2470,21 +2470,25 @@ class Num < Numeric
 
   # Adds a value to the exponent.
   def scaleb(other, context=nil)
-
     context = define_context(context)
     other = _convert(other)
     ans = _check_nans(context, other)
     return ans if ans
     return context.exception(InvalidOperation) if other.infinite? || other.exponent != 0
+    exp_inc = other.to_i
     unless context.exact?
       liminf = -2 * (context.emax + context.precision)
       limsup =  2 * (context.emax + context.precision)
-      i = other.to_i
-      return context.exception(InvalidOperation) if !((liminf <= i) && (i <= limsup))
+      if finite?
+        # Allow scaling a finite number which is originally out of the range
+        # of the context if the end result is within it
+        liminf = [context.etiny - @exp, liminf].min
+        limsup = [context.etop - @exp, limsup].max
+      end
+      return context.exception(InvalidOperation) if !((liminf <= exp_inc) && (exp_inc <= limsup))
     end
     return Num(self) if infinite?
-    return Num(@sign, @coeff, @exp+i)._fix(context)
-
+    return Num(@sign, @coeff, @exp + exp_inc)._fix(context)
   end
 
   # Naive implementation of exponential and logarithm functions; should be replaced
@@ -2706,7 +2710,13 @@ class Num < Numeric
       @sign*0.0
     else
       f = nil
-      f ||= to_s.to_f if num_class.radix == 10 # very precise, but slow
+      if num_class.radix == 10
+        # to_f can emit verbose warnings on overflow/underflow
+        old_verbose, $VERBOSE = $VERBOSE, nil
+        # very precise, but slow
+        f ||= to_s.to_f
+        $VERBOSE = old_verbose
+      end
       unless f
         c = @coeff.to_f
         if c.finite?

data/lib/flt/trigonometry.rb

@@ -134,8 +134,9 @@ module Flt
       rev_sign = false
       s = nil
       num_class.context(self) do |local_context|
-        local_context.precision += 3 # extra digits for intermediate steps
+        local_context.precision += 2 # extra digits for intermediate steps
         x,k,pi_2 = local_context.reduce_angle2(x,2)
+        local_context.precision += 2
         rev_sign = true if k>1
         if k % 2 == 0
           x = pi_2 - x
@@ -161,9 +162,10 @@ module Flt
       sign = x.sign
       s = nil
       num_class.context(self) do |local_context|
-        local_context.precision += 3 # extra digits for intermediate steps
+        local_context.precision += 2 # extra digits for intermediate steps
         x = x.copy_sign(+1) if sign<0
         x,k,pi_2 = local_context.reduce_angle2(x,2)
+        local_context.precision += 2
         sign = -sign if k>1
         x = pi_2 - x if k % 2 == 1
         x = local_context.to_rad(x)

data/lib/flt/version.rb

@@ -1,3 +1,3 @@
 module Flt
-  VERSION = "1.5.3"
+  VERSION = "1.5.4"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: flt
 version: !ruby/object:Gem::Version
-  version: 1.5.3
+  version: 1.5.4
 platform: ruby
 authors:
 - Javier Goizueta
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2020-04-07 00:00:00.000000000 Z
+date: 2020-04-08 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler