Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]


Groups > comp.lang.python > #32964

Re: accuracy problem in calculation

References <CA+b=61CdnBbFdxsO4yz5sK_w2GxOEXTxmBj387EEKwoQzu=W6w@mail.gmail.com>
Date 2012-11-09 04:13 +1100
Subject Re: accuracy problem in calculation
From Chris Angelico <rosuav@gmail.com>
Newsgroups comp.lang.python
Message-ID <mailman.3452.1352394817.27098.python-list@python.org> (permalink)

Show all headers | View raw


On Fri, Nov 9, 2012 at 4:05 AM, Debashish Saha <silideba@gmail.com> wrote:
> (1500000000+1.00067968)-(1500000000+1.00067961)
> Out[102]: 2.384185791015625e-07
>
> 1.00067968-(1.00067961)
> Out[103]: 7.000000001866624e-08
>
> above i am showing the two different results,though the two outputs
> should be same if we do it in copy(the lass one is acceptable value).
> so my question is how to increase the accuracy(windows7(32bit)
> ,python2.7.2)

Welcome to floating point. You're working with very large and very
small numbers, and you _will_ lose accuracy.

There are a few options. It's possible that a 64-bit build of Python
will give you more accuracy, but better would be to separate your huge
numbers from your tiny ones and work with them separately.
Alternatively, switch to the Decimal or Fraction types, but be aware
that your script will probably run a lot slower.

>>> from decimal import Decimal
>>> (Decimal("1500000000")+Decimal("1.00067968"))-(Decimal("1500000000")+Decimal("1.00067961"))
Decimal('7E-8')
>>> Decimal("1.00067968")-Decimal("1.00067961")
Decimal('7E-8')

Unless something's tying you to Python 2, consider moving to Python 3.
You may find that, on Python 3.3, you can switch to Decimal without
losing too much performance.

ChrisA

Back to comp.lang.python | Previous | NextNext in thread | Find similar | Unroll thread


Thread

Re: accuracy problem in calculation Chris Angelico <rosuav@gmail.com> - 2012-11-09 04:13 +1100
  Re: accuracy problem in calculation Grant Edwards <invalid@invalid.invalid> - 2012-11-08 18:17 +0000

csiph-web