Project

General

Profile

Actions

Bug #8542

closed

BigMath::exp modifies its first argument

Added by GSnyder (Garth Snyder) almost 11 years ago. Updated over 4 years ago.

Status:
Closed
Target version:
-
ruby -v:
ruby 2.0.0p195 (2013-05-14 revision 40734) [i386-solaris2.11]
Backport:
[ruby-core:55550]

Description

=begin
I noticed this when creating the patch at https://github.com/ruby/ruby/pull/332. I believe it affects everything from ruby 1.8 up.

BigMath::exp is implemented for negative numbers according to the identity E ** -x == 1 / (E ** x). At bigdecimal.c:2742 (in GitHub commit 258acf3, in the implementation of BigMath_s_exp), the code flips the sign on the input argument for negative numbers. Later, it returns the reciprocal of the result.

Problem: When the first argument is a BigDecimal (or, presumably, Bignum), the original argument is modified, so its sign bit remains flipped. Hence:

$ irb
2.0.0-p195 :001 > require 'bigdecimal/math'
=> true
2.0.0-p195 :002 > x = BigDecimal(-1)
=> #BigDecimal:827d960,'-0.1E1',9(36)
2.0.0-p195 :003 > BigMath.exp(x, 20)
=> #<BigDecimal:822dca8,'0.3678794411 714423216E0',27(54)>
2.0.0-p195 :004 > x
=> #BigDecimal:827d960,'0.1E1',9(36)

Here, x is modified when used as an argument to BigMath.exp AND the answer is wrong. I have already submitted the patch mentioned above for the latter problem, but I'm not sure what the appropriate fix would be for the sign modification. Just resetting the sign bit on exit would be easy but not thread safe; BigMath.exp really shouldn't be modifying the argument at all. But copying the whole argument is potentially wasteful if the precision is high.

I suspect that the special calculation track for negative values is actually not needed at all. Without the patch I just submitted, BigMath.exp is reliably returning the reciprocal of the correct answer, which means that it's properly calculating the correct answer by using the (negative) VALUE x passed in as the original argument -- at least for immediate values. So perhaps the basic iteration loop is just as valid for negative exponents as it is already implemented.
=end


Files

bigdecimal-exp-no-mutate-8542.patch (1.63 KB) bigdecimal-exp-no-mutate-8542.patch jeremyevans0 (Jeremy Evans), 08/10/2019 05:29 AM
Actions

Also available in: Atom PDF

Like0
Like0Like0Like0Like0Like0Like0Like0Like0