flt 1.5.0 → 1.5.1
Sign up to get free protection for your applications and to get access to all the features.
- checksums.yaml +4 -4
- data/History.txt +8 -0
- data/README.md +765 -0
- data/expand.rb +2 -0
- data/flt.gemspec +1 -0
- data/lib/flt.rb +1 -1
- data/lib/flt/bigdecimal.rb +2 -5
- data/lib/flt/bin_num.rb +1 -4
- data/lib/flt/complex.rb +10 -9
- data/lib/flt/dec_num.rb +9 -10
- data/lib/flt/float.rb +4 -6
- data/lib/flt/math.rb +3 -4
- data/lib/flt/num.rb +54 -21
- data/lib/flt/sugar.rb +1 -2
- data/lib/flt/support.rb +5 -2
- data/lib/flt/support/flag_values.rb +1 -1
- data/lib/flt/support/rationalizer.rb +3 -3
- data/lib/flt/support/rationalizer_extra.rb +6 -6
- data/lib/flt/support/reader.rb +1 -1
- data/lib/flt/tolerance.rb +2 -7
- data/lib/flt/trigonometry.rb +0 -2
- data/lib/flt/version.rb +1 -1
- data/setup.rb +1 -2
- data/test/generate_trig_data.rb +6 -6
- data/test/helper.rb +11 -5
- data/test/reader.rb +1 -1
- data/test/test_base_digits.rb +4 -2
- data/test/test_basic.rb +7 -7
- data/test/test_big_decimal.rb +21 -21
- data/test/test_bin.rb +1 -1
- data/test/test_bin_arithmetic.rb +1 -1
- data/test/test_binfloat_conversion.rb +1 -1
- data/test/test_coercion.rb +1 -1
- data/test/test_comparisons.rb +4 -4
- data/test/test_dectest.rb +2 -2
- data/test/test_define_conversions.rb +10 -10
- data/test/test_epsilon.rb +1 -1
- data/test/test_exact.rb +12 -11
- data/test/test_flags.rb +1 -1
- data/test/test_float.rb +23 -25
- data/test/test_format.rb +1 -1
- data/test/test_formatter.rb +6 -8
- data/test/test_hex_format.rb +2 -2
- data/test/test_multithreading.rb +1 -1
- data/test/test_normalized.rb +1 -1
- data/test/test_num_constructor.rb +6 -1
- data/test/test_odd_even.rb +1 -1
- data/test/test_rationalizer.rb +1 -1
- data/test/test_round.rb +36 -1
- data/test/test_sugar.rb +6 -6
- data/test/test_to_int.rb +4 -4
- data/test/test_to_rf.rb +126 -1
- data/test/test_tol.rb +1 -1
- data/test/test_trig.rb +1 -1
- data/test/test_ulp.rb +3 -3
- metadata +18 -4
- data/README.rdoc +0 -614
data/test/test_tol.rb
CHANGED
data/test/test_trig.rb
CHANGED
data/test/test_ulp.rb
CHANGED
@@ -1,6 +1,6 @@
|
|
1
1
|
require File.expand_path(File.join(File.dirname(__FILE__),'helper.rb'))
|
2
2
|
|
3
|
-
class TestUlp < Test
|
3
|
+
class TestUlp < Minitest::Test
|
4
4
|
|
5
5
|
|
6
6
|
def setup
|
@@ -67,7 +67,7 @@ class TestUlp < Test::Unit::TestCase
|
|
67
67
|
DecNum.context.exact = true
|
68
68
|
assert DecNum(1).ulp.nan?, "No ulps can be computed in exact contexts"
|
69
69
|
DecNum.context.traps[DecNum::InvalidOperation] = true
|
70
|
-
|
70
|
+
assert_raises(DecNum::InvalidOperation) { DecNum(1).ulp }
|
71
71
|
|
72
72
|
end
|
73
73
|
|
@@ -119,7 +119,7 @@ class TestUlp < Test::Unit::TestCase
|
|
119
119
|
BinNum.context.exact = true
|
120
120
|
assert BinNum(1).ulp.nan?, "No ulps can be computed in exact contexts"
|
121
121
|
BinNum.context.traps[BinNum::InvalidOperation] = true
|
122
|
-
|
122
|
+
assert_raises(BinNum::InvalidOperation) { BinNum(1).ulp }
|
123
123
|
|
124
124
|
end
|
125
125
|
|
metadata
CHANGED
@@ -1,14 +1,14 @@
|
|
1
1
|
--- !ruby/object:Gem::Specification
|
2
2
|
name: flt
|
3
3
|
version: !ruby/object:Gem::Version
|
4
|
-
version: 1.5.
|
4
|
+
version: 1.5.1
|
5
5
|
platform: ruby
|
6
6
|
authors:
|
7
7
|
- Javier Goizueta
|
8
8
|
autorequire:
|
9
9
|
bindir: bin
|
10
10
|
cert_chain: []
|
11
|
-
date:
|
11
|
+
date: 2017-03-04 00:00:00.000000000 Z
|
12
12
|
dependencies:
|
13
13
|
- !ruby/object:Gem::Dependency
|
14
14
|
name: bundler
|
@@ -38,6 +38,20 @@ dependencies:
|
|
38
38
|
- - ">="
|
39
39
|
- !ruby/object:Gem::Version
|
40
40
|
version: '0'
|
41
|
+
- !ruby/object:Gem::Dependency
|
42
|
+
name: minitest
|
43
|
+
requirement: !ruby/object:Gem::Requirement
|
44
|
+
requirements:
|
45
|
+
- - "~>"
|
46
|
+
- !ruby/object:Gem::Version
|
47
|
+
version: '5.10'
|
48
|
+
type: :development
|
49
|
+
prerelease: false
|
50
|
+
version_requirements: !ruby/object:Gem::Requirement
|
51
|
+
requirements:
|
52
|
+
- - "~>"
|
53
|
+
- !ruby/object:Gem::Version
|
54
|
+
version: '5.10'
|
41
55
|
description: Decimal and binary arbitrary precision floating point numbers in pure
|
42
56
|
Ruby.
|
43
57
|
email:
|
@@ -50,7 +64,7 @@ files:
|
|
50
64
|
- Gemfile
|
51
65
|
- History.txt
|
52
66
|
- License.txt
|
53
|
-
- README.
|
67
|
+
- README.md
|
54
68
|
- Rakefile
|
55
69
|
- expand.rb
|
56
70
|
- flt.gemspec
|
@@ -130,7 +144,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
|
|
130
144
|
version: '0'
|
131
145
|
requirements: []
|
132
146
|
rubyforge_project:
|
133
|
-
rubygems_version: 2.
|
147
|
+
rubygems_version: 2.6.8
|
134
148
|
signing_key:
|
135
149
|
specification_version: 4
|
136
150
|
summary: Floating Point Numbers
|
data/README.rdoc
DELETED
@@ -1,614 +0,0 @@
|
|
1
|
-
= Introduction
|
2
|
-
{<img src="https://badge.fury.io/rb/flt.svg" alt="Gem Version" />}[http://badge.fury.io/rb/flt]
|
3
|
-
{<img src="https://travis-ci.org/jgoizueta/flt.svg?branch=master" alt="Build tests" />}[https://travis-ci.org/jgoizueta/flt]
|
4
|
-
|
5
|
-
This library provides arbitrary precision floating-point types for Ruby. All types and
|
6
|
-
functions are within a namespace called Flt. Decimal and Binary floating point numbers
|
7
|
-
are implemented in classes Flt::DecNum and Flt::BinNum. These types are completely
|
8
|
-
written in Ruby using the multiple precision native integers. The performance
|
9
|
-
could be improved in the future by using a C extension based on the decNumber libray.
|
10
|
-
|
11
|
-
The Flt::Tolerance classes and the Flt.Tolerance() constructor handle floating point
|
12
|
-
tolerances defined in flexible ways.
|
13
|
-
|
14
|
-
Context classes are defined in the files
|
15
|
-
flt/float.rb[link:lib/flt/float_rb.html] and
|
16
|
-
flt/bigdecimal.rb[link:lib/flt/bigdecimal_rb.html]
|
17
|
-
for Float and BigDecimal numbers that aid to the interchangeability of floating point types. This
|
18
|
-
represent the only definition of identifiers outside the Flt namespace: the methods
|
19
|
-
Float.context() and BigDecimal.context() and some contants in Float.
|
20
|
-
|
21
|
-
This library is the successor of the ruby-decimal gem, that defined the Decimal class
|
22
|
-
for decimal floating point; that class has been renamed to Flt::DecNum and support
|
23
|
-
has been added for binary floating point and tolerances.
|
24
|
-
|
25
|
-
The documentation for this package is available at http://rdoc.info/github/jgoizueta/flt/frames
|
26
|
-
|
27
|
-
The code is at http://github.com/jgoizueta/flt/
|
28
|
-
|
29
|
-
= DecNum
|
30
|
-
|
31
|
-
Flt::DecNum is a standards-compliant arbitrary precision decimal floating-point type for Ruby.
|
32
|
-
It is based on the Python Decimal class.
|
33
|
-
|
34
|
-
== Standars compliance.
|
35
|
-
|
36
|
-
DecNum is designed to be conformant to the General Decimal Arithmetic Specification
|
37
|
-
and the revised IEEE 754 standard (IEEE 754-2008).
|
38
|
-
|
39
|
-
= Examples of use
|
40
|
-
|
41
|
-
To install the library use gem from the command line:
|
42
|
-
|
43
|
-
gem install flt
|
44
|
-
|
45
|
-
Then require the library in your code (if it fails you may need to <tt>require 'rubygems'</tt> first)
|
46
|
-
require 'flt'
|
47
|
-
include Flt
|
48
|
-
|
49
|
-
Now we can use the DecNum class simply like this:
|
50
|
-
|
51
|
-
puts DecNum(1)/DecNum(3) # -> 0.3333333333333333333333333333
|
52
|
-
|
53
|
-
DecNum() is a constructor that can be used instead of DecNum.new()
|
54
|
-
|
55
|
-
== Contexts
|
56
|
-
|
57
|
-
Contexts are environments for arithmetic operations. They govern precision, set rules
|
58
|
-
for rounding, determine which signals are treated as exceptions, and limit the range
|
59
|
-
for exponents.
|
60
|
-
|
61
|
-
Each thread has an active context that can be accessed like this:
|
62
|
-
|
63
|
-
puts DecNum.context.precision # -> 28
|
64
|
-
|
65
|
-
The active context can be modified globally for the current thread:
|
66
|
-
|
67
|
-
DecNum.context.precision = 2
|
68
|
-
puts DecNum.context.precision # -> 2
|
69
|
-
puts DecNum(1)/DecNum(3) # -> 0.33
|
70
|
-
DecNum.context.precision += 7
|
71
|
-
puts DecNum.context.precision # -> 9
|
72
|
-
puts DecNum(1)/DecNum(3) # -> 0.333333333
|
73
|
-
|
74
|
-
Or it can be altered locally inside a block:
|
75
|
-
|
76
|
-
DecNum.context do
|
77
|
-
DecNum.context.precision = 5
|
78
|
-
puts DecNum.context.precision # -> 5
|
79
|
-
end
|
80
|
-
puts DecNum.context.precision # -> 9
|
81
|
-
|
82
|
-
The block for a local context can be passed the current context as an argument:
|
83
|
-
|
84
|
-
DecNum.context do |local_context|
|
85
|
-
local_context.precision = 5
|
86
|
-
puts DecNum.context.precision # -> 5
|
87
|
-
end
|
88
|
-
puts DecNum.context.precision # -> 9
|
89
|
-
|
90
|
-
A context object can also be used to define the local context:
|
91
|
-
|
92
|
-
my_context = DecNum::Context(:precision=>20)
|
93
|
-
DecNum.context(my_context) do |context|
|
94
|
-
puts context.precision
|
95
|
-
end # -> 20
|
96
|
-
|
97
|
-
And individual parameters can be assigned like this:
|
98
|
-
|
99
|
-
puts DecNum.context.precision # -> 9
|
100
|
-
puts DecNum.context.rounding # -> half_even
|
101
|
-
DecNum.context(:rounding=>:down) do |context|
|
102
|
-
puts context.precision # -> 9
|
103
|
-
puts context.rounding # -> down
|
104
|
-
end
|
105
|
-
|
106
|
-
Contexts created with the DecNum::Context() constructor
|
107
|
-
inherit from DecNum::DefaultContext.
|
108
|
-
Default context attributes can be established by modifying
|
109
|
-
that object:
|
110
|
-
|
111
|
-
DecNum::DefaultContext.precision = 10
|
112
|
-
DecNum.context = DecNum::Context(:rounding=>:half_up)
|
113
|
-
puts DecNum.context.precision # -> 10
|
114
|
-
|
115
|
-
Note that a context object assigned to DecNum.context is copied,
|
116
|
-
so it is not altered through DecNum.context:
|
117
|
-
|
118
|
-
puts my_context.precision # -> 20
|
119
|
-
DecNum.context = my_context
|
120
|
-
DecNum.context.precision = 2
|
121
|
-
puts my_context.precision # -> 20
|
122
|
-
|
123
|
-
So, DefaultContext is not altered when modifying DecNum.context.
|
124
|
-
|
125
|
-
Methods that use a context have an optional parameter to override
|
126
|
-
the active context (DecNum.context) :
|
127
|
-
|
128
|
-
DecNum.context.precision = 3
|
129
|
-
puts DecNum(1).divide(3) # -> 0.333
|
130
|
-
puts DecNum(1).divide(3, my_context) # -> 0.33333333333333333333
|
131
|
-
|
132
|
-
Individual context parameters can also be overriden:
|
133
|
-
|
134
|
-
puts DecNum(1).divide(3, :precision=>6) # -> 0.333333
|
135
|
-
|
136
|
-
There are two additional predefined contexts DecNum::ExtendedContext
|
137
|
-
and DecNum::BasicContext that are not meant to be modified; they
|
138
|
-
can be used to achieve reproducible results. We will use
|
139
|
-
DecNum::ExtendedContext in the following examples:
|
140
|
-
|
141
|
-
DecNum.context = DecNum::ExtendedContext
|
142
|
-
|
143
|
-
Most decimal operations can be executed by using either Context or DecNum methods:
|
144
|
-
|
145
|
-
puts DecNum.context.exp(1) # -> 2.71828183
|
146
|
-
puts DecNum(1).exp # -> 2.71828183
|
147
|
-
|
148
|
-
If using Context methods, values are automatically converted as if the DecNum() constructor
|
149
|
-
was used.
|
150
|
-
|
151
|
-
==Rounding
|
152
|
-
|
153
|
-
Results are normally rounded using the precision (number of significant digits)
|
154
|
-
and rounding mode defined in the context.
|
155
|
-
|
156
|
-
DecNum.context.precision = 4
|
157
|
-
puts DecNum(1)/DecNum(3) # -> 0.3333
|
158
|
-
puts DecNum('1E20')-DecNum('1E-20') # -> 1.000E+20
|
159
|
-
DecNum.context.rounding = :half_up
|
160
|
-
puts +DecNum('100.05') # -> 100.1
|
161
|
-
DecNum.context.rounding = :half_even
|
162
|
-
puts +DecNum('100.05') # -> 100.0
|
163
|
-
|
164
|
-
Note that input values are not rounded, only results; we use
|
165
|
-
the plus operator to force rounding here:
|
166
|
-
|
167
|
-
DecNum.context.precision = 4
|
168
|
-
x = DecNum('123.45678')
|
169
|
-
puts x # -> 123.45678
|
170
|
-
puts +x # -> 123.5
|
171
|
-
|
172
|
-
Precision can be also set to 'exact' to avoid rounding, by using
|
173
|
-
the exact property or using a 0 precision. In exact mode results
|
174
|
-
are never rounded and results that have an infinite number of
|
175
|
-
digits trigger the DecNum::Inexact exception.
|
176
|
-
|
177
|
-
DecNum.context.exact = true
|
178
|
-
puts DecNum('1E20')-DecNum('1E-20') # -> 99999999999999999999.99999999999999999999
|
179
|
-
puts DecNum(16).sqrt # -> 4
|
180
|
-
puts DecNum(16)/DecNum(4) # -> 4
|
181
|
-
puts DecNum(1)/DecNum(3) # -> Exception : Flt::Num::Inexact
|
182
|
-
|
183
|
-
DecNum.context.precision = 5
|
184
|
-
puts DecNum('1E20')-DecNum('1E-20') # -> 1.0000E+20
|
185
|
-
puts DecNum(16).sqrt # -> 4
|
186
|
-
puts DecNum(16)/DecNum(4) # -> 4
|
187
|
-
puts DecNum(1)/DecNum(3) # -> 0.33333
|
188
|
-
|
189
|
-
There are also some methods for explicit rounding that provide
|
190
|
-
an interface compatible with that of the Ruby Float class:
|
191
|
-
|
192
|
-
puts DecNum('101.5').round # -> 102
|
193
|
-
puts DecNum('101.5').round(0) # -> 102
|
194
|
-
puts DecNum('101.12345').round(2) # -> 101.12
|
195
|
-
puts DecNum('101.12345').round(-1) # -> 1.0E+2
|
196
|
-
puts DecNum('101.12345').round(:places=>2) # -> 101.12
|
197
|
-
puts DecNum('101.12345').round(:precision=>2) # -> 1.0E+2
|
198
|
-
puts DecNum('101.5').round(:rounding=>:half_up) # -> 102
|
199
|
-
puts DecNum('101.5').ceil # -> 102
|
200
|
-
puts DecNum('101.5').floor # -> 101
|
201
|
-
puts DecNum('101.5').truncate # -> 101
|
202
|
-
|
203
|
-
==Special values
|
204
|
-
|
205
|
-
In addition to finite numbers, a DecNum object can represent some special values:
|
206
|
-
* Infinity (+Infinity, -Infinity). The method DecNum#infinite? returns true for these to values.
|
207
|
-
DecNum.infinity DecNum.infinity(-1) can be used to get these values.
|
208
|
-
* NaN (not a number) represents undefined results. The method DecNum#nan? returns true for it and
|
209
|
-
DecNum.nan can be used to obtain it. There is a variant, sNaN (signaling NaN) that causes
|
210
|
-
an invalid operation condition if used; it can be detected with DecNum.snan?.
|
211
|
-
A NaN can also include diagnostic information in its sign and coefficient.
|
212
|
-
|
213
|
-
Any of the special values can be detected with DecNum#special?
|
214
|
-
Finite numbers can be clasified with
|
215
|
-
these methods:
|
216
|
-
* DecNum#zero? detects a zero value (note that there are two zero values: +0 and -0)
|
217
|
-
* DecNum#normal? detects normal values: those whose adjusted exponents are not less than @emin@.
|
218
|
-
* DecNum#subnormal? detects subnormal values: those whose adjusted exponents are less than @emin@.
|
219
|
-
|
220
|
-
==Exceptions
|
221
|
-
|
222
|
-
Exceptional conditions that may arise during operations have corresponding classes that represent them:
|
223
|
-
* DecNum::InvalidOperation
|
224
|
-
* DecNum::DivisionByZero
|
225
|
-
* DecNum::DivisionImpossible
|
226
|
-
* DecNum::DivisionUndefined
|
227
|
-
* DecNum::Inexact
|
228
|
-
* DecNum::Overflow
|
229
|
-
* DecNum::Underflow
|
230
|
-
* DecNum::Clamped
|
231
|
-
* DecNum::InvalidContext
|
232
|
-
* DecNum::Rounded
|
233
|
-
* DecNum::Subnormal
|
234
|
-
* DecNum::ConversionSyntax
|
235
|
-
|
236
|
-
For each condition, a flag and a trap (boolean values) exist in the context.
|
237
|
-
When a condition occurs, the corresponding flag in the context takes the value true (and remains
|
238
|
-
set until cleared) and a exception is raised if the corresponding trap has the value true.
|
239
|
-
|
240
|
-
DecNum.context.traps[DecNum::DivisionByZero] = false
|
241
|
-
DecNum.context.flags[DecNum::DivisionByZero] = false
|
242
|
-
puts DecNum(1)/DecNum(0) # -> Infinity
|
243
|
-
puts DecNum.context.flags[DecNum::DivisionByZero] # -> true
|
244
|
-
DecNum.context.traps[DecNum::DivisionByZero] = true
|
245
|
-
puts DecNum(1)/DecNum(0) # -> Exception : Flt::Num::DivisionByZero
|
246
|
-
|
247
|
-
==Numerical conversion
|
248
|
-
|
249
|
-
By default, DecNum is interoperable with Integer and Rational.
|
250
|
-
Conversion happens automatically to operands:
|
251
|
-
|
252
|
-
puts DecNum('0.1') + 1 # -> 1.1
|
253
|
-
puts 7 + DecNum('0.2') # -> 7.2
|
254
|
-
puts Rational(5,2) + DecNum('3') # -> 5.5
|
255
|
-
|
256
|
-
Conversion can also be done explicitely with
|
257
|
-
the DecNum constructor:
|
258
|
-
|
259
|
-
puts DecNum(7) # -> 7
|
260
|
-
puts DecNum(Rational(1,10)) # -> 0.1
|
261
|
-
|
262
|
-
Converting a DecNum to other numerical types can be done with specific Ruby-style methods.
|
263
|
-
|
264
|
-
puts DecNum('1.1').to_i # -> 1
|
265
|
-
puts DecNum('1.1').to_r # -> 11/10
|
266
|
-
|
267
|
-
(note the truncated result of to_i)
|
268
|
-
Or with a generic method:
|
269
|
-
puts DecNum('1.1').convert_to(Integer) # -> 1
|
270
|
-
puts DecNum('1.1').convert_to(Rational) # -> 11/10
|
271
|
-
|
272
|
-
Thera are also GDAS style conversion operations:
|
273
|
-
|
274
|
-
puts DecNum('1.1').to_integral_value # -> 1
|
275
|
-
|
276
|
-
And conversion is also possible to Float:
|
277
|
-
puts DecNum('1.1').to_f # -> 1.1
|
278
|
-
puts DecNum('1.1').convert_to(Float) # -> 1.1
|
279
|
-
puts Float(DecNum('1.1')) # -> 1.1
|
280
|
-
|
281
|
-
Types with predefined bidirectional conversion (Integer and Rational)
|
282
|
-
can be operated with DecNum on either side of an operator, and the result will be a DecNum.
|
283
|
-
For Float there is no predefined bidirectional conversion (see below how to define it)
|
284
|
-
and the result of an operation between DecNum and Float will be of type Float.
|
285
|
-
|
286
|
-
puts (DecNum('1.1') + 2.0).class # -> Float
|
287
|
-
puts (2.0 + DecNum('1.1')).class # -> Float
|
288
|
-
|
289
|
-
The conversion system is extensible. For example, we can include BigDecimal into it
|
290
|
-
by defining suitable conversion procedures:
|
291
|
-
|
292
|
-
DecNum.context.define_conversion_from(BigDecimal) do |x, context|
|
293
|
-
DecNum(x.to_s)
|
294
|
-
end
|
295
|
-
DecNum.context.define_conversion_to(BigDecimal) do |x|
|
296
|
-
BigDecimal.new(x.to_s)
|
297
|
-
end
|
298
|
-
|
299
|
-
Now we can mix BigDecimals and Decimals in expressions and convert from DecNum
|
300
|
-
to BigDecimal:
|
301
|
-
|
302
|
-
puts BigDecimal.new('1.1') + DecNum('2.2') # -> 3.3
|
303
|
-
puts DecNum('1.1').convert_to(BigDecimal) # -> 0.11E1
|
304
|
-
|
305
|
-
Note that the conversions are defined in a Context object and will be available only
|
306
|
-
when that context applies. That way we can define conversions for specific purposes
|
307
|
-
without affecting a program globally.
|
308
|
-
|
309
|
-
As another example consider conversion from Float to DecNum, which is not defined by
|
310
|
-
default because it can be defined in different ways depending on the purpose.
|
311
|
-
|
312
|
-
A Float constant such as 0.1 defines a Float object which has a numerical value close to,
|
313
|
-
but not exactly 1/10. When converting that Float to DecNum we could decide to preserve
|
314
|
-
the exact numerical value of the number or try to find a simple decimal expression within
|
315
|
-
a given tolerance. If we take the first approach we can define this conversion:
|
316
|
-
|
317
|
-
DecNum.context.define_conversion_from(Float) do |x, context|
|
318
|
-
s,e = Math.frexp(x)
|
319
|
-
s = Math.ldexp(s, Float::MANT_DIG).to_i
|
320
|
-
e -= Float::MANT_DIG
|
321
|
-
DecNum(s*(Float::RADIX**e))
|
322
|
-
end
|
323
|
-
|
324
|
-
Note that the conversion we've defined depends on the context precision:
|
325
|
-
|
326
|
-
DecNum.local_context(:precision=>20) { puts DecNum(0.1) } # -> 0.10000000000000000555
|
327
|
-
|
328
|
-
DecNum.local_context(:precision=>12) { puts DecNum(0.1) } # -> 0.100000000000
|
329
|
-
|
330
|
-
DecNum.local_context(:exact=>true) { puts DecNum(0.1) } # -> 0.1000000000000000055511151231257827021181583404541015625
|
331
|
-
|
332
|
-
A different approach for Float to DecNum conversion is to find the shortest (fewer digits) DecNum
|
333
|
-
that rounds to the Float with the binary precision that the Float has.
|
334
|
-
We will assume that the DecNum to Float conversion done with the rounding mode of the DecNum context.
|
335
|
-
The BinNum class has a method to perform this kind of conversion, so we will use it.
|
336
|
-
|
337
|
-
DecNum.context.define_conversion_from(Float) do |x, dec_context|
|
338
|
-
BinNum.context(:rounding=>dec_context.rounding) do |bin_context|
|
339
|
-
BinNum(x).to_decimal
|
340
|
-
end
|
341
|
-
end
|
342
|
-
|
343
|
-
The result is independent of the context precision.
|
344
|
-
|
345
|
-
puts DecNum(0.1) # -> 0.1
|
346
|
-
puts DecNum(1.0/3) # -> 0.3333333333333333
|
347
|
-
|
348
|
-
This conversion gives the results expected most of the time, but it must be noticed that
|
349
|
-
there must be some compromise, because different decimal literals convert to the same Float value:
|
350
|
-
|
351
|
-
puts DecNum(0.10000000000000001) # -> 0.1
|
352
|
-
|
353
|
-
There's also some uncertainty because the way the Ruby interpreter parses Float literals
|
354
|
-
may not be well specified; in the usual case (IEEE Double Floats and round-to-even)
|
355
|
-
the results will be as expected (correctly rounded Floats), but some platforms may
|
356
|
-
behave differently.
|
357
|
-
|
358
|
-
The BinNum also a instance method +to_decimal_exact+ to perform the previous 'exact' conversion, that
|
359
|
-
could have be written:
|
360
|
-
|
361
|
-
DecNum.context.define_conversion_from(Float) do |x, context|
|
362
|
-
DecNum.context(context) do
|
363
|
-
BinNum(x).to_decimal_exact
|
364
|
-
end
|
365
|
-
end
|
366
|
-
|
367
|
-
== Abbreviation
|
368
|
-
|
369
|
-
The use of DecNum can be made less verbose by requiring:
|
370
|
-
|
371
|
-
require 'flt/d'
|
372
|
-
|
373
|
-
This file defines +D+ as a synonym for +DecNum+:
|
374
|
-
|
375
|
-
D.context.precision = 3
|
376
|
-
puts +D('1.234') # -> 1.23
|
377
|
-
|
378
|
-
Some convenient methods are added to numeric classes by requiring the optional
|
379
|
-
flt/sugar.rb[link:lib/flt/sugar_rb.html]. This must be explicitely required
|
380
|
-
because it could cause conflicts with other extensions of these classes.
|
381
|
-
|
382
|
-
require 'flt/sugar'
|
383
|
-
|
384
|
-
puts 34.odd? # -> false
|
385
|
-
puts 34.even? # -> true
|
386
|
-
puts 0.1.split.inspect # -> [1, 7205759403792794, -56]
|
387
|
-
puts (-0.1).sign # -> -1
|
388
|
-
|
389
|
-
A shortcut notation for DecNum is defined in this file (based on an idea by coderrr[http://coderrr.wordpress.com]
|
390
|
-
which allows exact definitions with almost literal decimal syntax (note the underscore after the dot.)
|
391
|
-
|
392
|
-
puts 10._123456789123456789 # -> 10.123456789123456789
|
393
|
-
|
394
|
-
Additional underscores can be used to separate digits at any place except before the decimal point:
|
395
|
-
|
396
|
-
puts 100_000._000_001 # -> 100000.000001
|
397
|
-
puts 100_000._000_001.class # -> Flt::DecNum
|
398
|
-
|
399
|
-
But note that 100_000.000_001 is a valid Float (it's not a DecNum because there isn't an underscore just after the
|
400
|
-
decimal point):
|
401
|
-
|
402
|
-
puts 100_000.000_001 # -> 100000.000001
|
403
|
-
puts 100_000.000_001.class # -> Float
|
404
|
-
|
405
|
-
There's also one important caveat with this notation: negative numbers with a zero integral part must be
|
406
|
-
parenthesed (otherwise the minus has no effect because it affects only the null integer part):
|
407
|
-
|
408
|
-
puts -0._5 # -> 0.5
|
409
|
-
puts -(0._5) # -> -0.5
|
410
|
-
|
411
|
-
== Error analysis
|
412
|
-
|
413
|
-
The DecNum#ulp() method returns the value of a "unit in the last place" for a given number under
|
414
|
-
the current context.
|
415
|
-
|
416
|
-
D.context.precision = 4
|
417
|
-
puts D('1.5').ulp # -> 0.001
|
418
|
-
puts D('1.5E10').ulp # -> 1E+7
|
419
|
-
|
420
|
-
Whe can compute the error in ulps of an approximation +aprx+ to correclty rounded value +exct+ with:
|
421
|
-
|
422
|
-
def ulps(exct, aprx)
|
423
|
-
(aprx-exct).abs/exct.ulp
|
424
|
-
end
|
425
|
-
|
426
|
-
puts ulps(DecNum('0.5000'), DecNum('0.5003')) # -> 3
|
427
|
-
puts ulps(DecNum('0.5000'), DecNum('0.4997')) # -> 3
|
428
|
-
|
429
|
-
puts ulps(DecNum('0.1000'), DecNum('0.1003')) # -> 3E+1
|
430
|
-
puts ulps(DecNum('0.1000'), DecNum('0.0997')) # -> 3E+1
|
431
|
-
|
432
|
-
puts ulps(DecNum(1), DecNum(10).next_minus) # -> 8.999E+4
|
433
|
-
puts ulps(DecNum(1), DecNum(10).next_plus) # -> 9.01E+4
|
434
|
-
|
435
|
-
Note that in the definition of ulps we use exct.ulp. If we had use aprx.ulp DecNum(10).next_plus
|
436
|
-
would seem to be a better approximation to DecNum(1) than DecNum(10).next_minus. (Admittedly,
|
437
|
-
such bad approximations should not be common.)
|
438
|
-
|
439
|
-
== BinNum Input/Output
|
440
|
-
|
441
|
-
BinNum can be defined with a decimal string literal and converted to one with to_s, as DecNum,
|
442
|
-
but since this involves a change of base these are inexact operations subject to some specific precision limits.
|
443
|
-
|
444
|
-
If we define the number with a binary literal, no base conversion is involved and the result is exactly defined;
|
445
|
-
here we define a number with just one bit of precision:
|
446
|
-
|
447
|
-
x = BinNum('0.001', :base=>2)
|
448
|
-
puts x.number_of_digits # -> 1
|
449
|
-
puts x.to_s(:base=>2) # -> 0.001
|
450
|
-
|
451
|
-
Note that we could have defined it with more precision, e.g.
|
452
|
-
|
453
|
-
y = BinNum('0.001000', :base=>2)
|
454
|
-
puts y.number_of_digits # -> 4
|
455
|
-
puts y.to_s(:base=>2) # -> 0.001000
|
456
|
-
|
457
|
-
But let's get back to our one bit quantity, x, and convert it to a decimal string. Since the internal precision
|
458
|
-
is only one bit it contains very little information:
|
459
|
-
|
460
|
-
puts x # -> 0.1
|
461
|
-
|
462
|
-
We can obtain more digits with the :all_digits option which show all the decimal digits that are significative
|
463
|
-
for the given precision of 1 bit:
|
464
|
-
|
465
|
-
puts x.to_s(:all_digits=>true) # -> 0.12
|
466
|
-
|
467
|
-
We can also obtain the exact value of x by using the Num.convert_exact method to convert it to decimal (base 10):
|
468
|
-
|
469
|
-
puts Num.convert_exact(x,10) # -> 0.125
|
470
|
-
|
471
|
-
Let's convert the default decimal output back to another BinNum which will preserve its precision:
|
472
|
-
|
473
|
-
y = BinNum(x.to_s)
|
474
|
-
|
475
|
-
The result may seem ok:
|
476
|
-
|
477
|
-
puts y # -> 0.1
|
478
|
-
|
479
|
-
But is not exactly what we originally had:
|
480
|
-
|
481
|
-
puts y==x # -> false
|
482
|
-
puts y # -> 0.1
|
483
|
-
puts y.number_of_digits # -> 5
|
484
|
-
puts y.to_s(:base=>2) # -> 0.00011010
|
485
|
-
|
486
|
-
The new value y has gained some digits because of the intermediate conversion to decimal: one decimal digit
|
487
|
-
contains more information than one bit, and the result shows that.
|
488
|
-
|
489
|
-
If we wanted to preserve exactly the number we should have done this:
|
490
|
-
|
491
|
-
y = BinNum(x.to_s, :fixed, :precision=>x.number_of_digits)
|
492
|
-
puts y==x # -> true
|
493
|
-
|
494
|
-
To preserve the value we had to explicitly determine how many bits should y have.
|
495
|
-
|
496
|
-
With the :fixed options the number produced by BinNum is rounded to the context precision (which can be
|
497
|
-
overriden as in the example by other options):
|
498
|
-
|
499
|
-
puts BinNum(x.to_s, :fixed, :precision=>32).to_s(:base=>2) # -> 0.00011001100110011001100110011001101
|
500
|
-
puts BinNum(x.to_s, :fixed, :precision=>1).to_s(:base=>2) # -> 0.001
|
501
|
-
|
502
|
-
Note also that if we normalize a value we will change it's precision to that of the context:
|
503
|
-
|
504
|
-
puts x.number_of_digits # -> 1
|
505
|
-
puts x.normalize.number_of_digits # -> 53
|
506
|
-
puts x.normalize.to_s # -> 0.125
|
507
|
-
|
508
|
-
== Mathematical functions
|
509
|
-
|
510
|
-
There are two mathematical functions modules analogous to Ruby's Math for Float,
|
511
|
-
Flt::DecNum::Math and Flt::BinNum::Math.
|
512
|
-
Currently they consist of basic trigonometric functions, including hypot, logarithms and the exponential
|
513
|
-
function, and the constants e and pi.
|
514
|
-
|
515
|
-
Its functions can be accessed in a number of ways:
|
516
|
-
|
517
|
-
require 'flt/math'
|
518
|
-
DecNum.context(:precision=>10) do |context|
|
519
|
-
# As module functions, using the current context for the enclosing Num class:
|
520
|
-
puts DecNum::Math.sin(1)*DecNum::Math.pi # -> 2.643559064
|
521
|
-
# As functions of a context object:
|
522
|
-
puts context.sin(1)*context.pi # -> 2.643559064
|
523
|
-
# Through a math block:
|
524
|
-
puts DecNum.context.math{sin(1)*pi} # -> 2.643559064
|
525
|
-
puts DecNum.math{sin(1)*pi} # -> 2.643559064
|
526
|
-
# And can be *included* to be used as private instance methods:
|
527
|
-
include DecNum::Math
|
528
|
-
puts sin(1)*pi # -> 2.643559064
|
529
|
-
end
|
530
|
-
|
531
|
-
== More Information
|
532
|
-
|
533
|
-
* Decimal Floating point type: see the base Flt::Num class and the Flt::DecNum class
|
534
|
-
* Binary Floating point type: see the base Flt::Num class and the Flt::BinNum class
|
535
|
-
* Floating Point Contexts: see documentation for classes Flt::Num::ContextBase,
|
536
|
-
Flt::DecNum::Context and Flt::BinNum::Context
|
537
|
-
* Floating Point Tolerance: see the flt/tolerance.rb[link:lib/flt/tolerance_rb.html] file
|
538
|
-
and the Flt::Tolerance class
|
539
|
-
* Constructors: see Flt.DecNum(), Flt.BinNum() and Flt.Tolerance().
|
540
|
-
* Trigonometry functions: see Flt::Trigonometry.
|
541
|
-
* Complex number support: see the flt/complex.rb[link:lib/flt/complex_rb.html] file
|
542
|
-
|
543
|
-
= DecNum vs BigDecimal
|
544
|
-
|
545
|
-
--
|
546
|
-
EXPAND-
|
547
|
-
++
|
548
|
-
|
549
|
-
DecNum solves some of the difficulties of using BigDecimal.
|
550
|
-
|
551
|
-
One of the major problems with BigDecimal is that it's not easy to control the number of
|
552
|
-
significant digits of the results. While addition, subtraction and multiplication are exact (unless a limit is used),
|
553
|
-
divisions will need to be passed precision explicitly or else an indeterminate number of significant digits will be lost.
|
554
|
-
Part of the problem is that numbers don't keep track of its precision (0.1000 is not distinguishable from 0.1.)
|
555
|
-
|
556
|
-
With DecNum, Context objects are used to specify the exact number of digits to be used for all operations making
|
557
|
-
the code cleaner and the results more easily predictable.
|
558
|
-
|
559
|
-
DecNum.context.precision = 10
|
560
|
-
puts DecNum(1)/DecNum(3)
|
561
|
-
|
562
|
-
Contexts are thread-safe and can be used for individual operations:
|
563
|
-
|
564
|
-
puts DecNum(1).divide(DecNum(e), DecNum::Context(:precision=>4))
|
565
|
-
|
566
|
-
Which can be abbreviated:
|
567
|
-
|
568
|
-
puts DecNum(1).divide(DecNum(e), :precision=>4)
|
569
|
-
|
570
|
-
Or use locally in a block without affecting other code:
|
571
|
-
|
572
|
-
DecNum.context {
|
573
|
-
DecNum.context.precision = 3
|
574
|
-
puts DecNum(1)/DecNum(3)
|
575
|
-
}
|
576
|
-
puts DecNum.context.precision
|
577
|
-
|
578
|
-
Which can also be abbreviated:
|
579
|
-
|
580
|
-
DecNum.context(:precision=>3) { puts DecNum(1)/DecNum(3) }
|
581
|
-
|
582
|
-
This allows in general to write simpler code; e.g. this is an exponential function, adapted from the
|
583
|
-
'recipes' in Python's Decimal:
|
584
|
-
|
585
|
-
def exp(x, c=nil)
|
586
|
-
i, lasts, s, fact, num = 0, 0, 1, 1, 1
|
587
|
-
DecNum.context(c) do |context|
|
588
|
-
context.precision += 2
|
589
|
-
while s != lasts
|
590
|
-
lasts = s
|
591
|
-
i += 1
|
592
|
-
fact *= i
|
593
|
-
num *= x
|
594
|
-
s += num / fact
|
595
|
-
end
|
596
|
-
end
|
597
|
-
return +s
|
598
|
-
end
|
599
|
-
|
600
|
-
The final unary + applied to the result forces it to be rounded to the current precision
|
601
|
-
(because we have computed it with two extra digits)
|
602
|
-
The result of this method does not have trailing non-significant digits, as is common with BigDecimal
|
603
|
-
(e.g. in the exp implementation available in the standard Ruby library, in bigdecimal/math)
|
604
|
-
|
605
|
-
--
|
606
|
-
EXPAND+
|
607
|
-
++
|
608
|
-
|
609
|
-
= Roadmap
|
610
|
-
|
611
|
-
* Trigonometry optimizations
|
612
|
-
* Implement the missing GDA functions:
|
613
|
-
rotate, shift, trim, and, or, xor, invert,
|
614
|
-
max, min, maxmag, minmag, comparetotal, comparetotmag
|