Clarification on the Decimal type in Python

by

The Decimal class is best for financial type addition, subtraction multiplication, division type problems:

>>> (1.1+2.2-3.3)*10000000000000000000
4440.892098500626                            # relevant for government invoices...
>>> import decimal
>>> D=decimal.Decimal
>>> (D('1.1')+D('2.2')-D('3.3'))*10000000000000000000
Decimal('0.0')

The Fraction module works well with the rational number problem domain you describe:

>>> from fractions import Fraction
>>> f = Fraction(1) / Fraction(3)
>>> f
Fraction(1, 3)
>>> f * 3 < 1
False
>>> f * 3 == 1
True

For pure multi precision floating point for scientific work, consider mpmath.

If your problem can be held to the symbolic realm, consider sympy. Here is how you would handle the 1/3 issue:

>>> sympy.sympify('1/3')*3
1
>>> (sympy.sympify('1/3')*3) == 1
True

Sympy uses mpmath for arbitrary precision floating point, includes the ability to handle rational numbers and irrational numbers symbolically.

Consider the pure floating point representation of the irrational value of √2: 

>>> math.sqrt(2)
1.4142135623730951
>>> math.sqrt(2)*math.sqrt(2)
2.0000000000000004
>>> math.sqrt(2)*math.sqrt(2)==2
False

Compare to sympy:

>>> sympy.sqrt(2)
sqrt(2)                              # treated symbolically
>>> sympy.sqrt(2)*sympy.sqrt(2)==2
True

You can also reduce values:

>>> import sympy
>>> sympy.sqrt(8)
2*sqrt(2)                            # √8 == √(4 x 2) == 2*√2...

However, you can see issues with Sympy similar to straight floating point if not careful:

>>> 1.1+2.2-3.3
4.440892098500626e-16
>>> sympy.sympify('1.1+2.2-3.3')
4.44089209850063e-16                   # :-(

This is better done with Decimal:

>>> D('1.1')+D('2.2')-D('3.3')
Decimal('0.0')

Or using Fractions or Sympy and keeping values such as 1.1 as ratios:

>>> sympy.sympify('11/10+22/10-33/10')==0
True
>>> Fraction('1.1')+Fraction('2.2')-Fraction('3.3')==0
True

Or use Rational in sympy:

>>> frac=sympy.Rational
>>> frac('1.1')+frac('2.2')-frac('3.3')==0
True
>>> frac('1/3')*3
1

You can play with sympy live.

More Related Contents:

  1. Python JSON serialize a Decimal object
  2. How to avoid floating point errors? [duplicate]
  3. How to properly round-up half float numbers?
  4. How to get numbers after decimal point?
  5. Convert hex to float
  6. Add zeros to a float after the decimal point in Python
  7. rounding errors in Python floor division
  8. Rounding a number in Python but keeping ending zeros
  9. Convert float to string in positional format (without scientific notation and false precision)
  10. What is the difference between int() and floor() in Python 3?
  11. Decimal place issues with floats and decimal.Decimal
  12. How to get largest possible precision? (Python – Decimal)
  13. Rounding floats so that they sum to precisely 1
  14. Print a float number in normal form, not exponential form / scientific notation [duplicate]
  15. ‘is’ operator behaves unexpectedly with floats
  16. code that count the amount of yes and no in a list [closed]
  17. Relative imports – ModuleNotFoundError: No module named x
  18. How to get a random number between a float range?
  19. python ctypes issue on different OSes
  20. Threads and tkinter
  21. tkinter.TclError: image “pyimage3” doesn’t exist
  22. Simplest async/await example possible in Python
  23. Object of type ‘map’ has no len() in Python 3
  24. What is the purpose of Python’s itertools.repeat?
  25. How to use threading in Python?
  26. Nested list comprehension scope
  27. Identifying the data type of an input
  28. “Failed to load HostKeys” warning while connecting to SFTP server with pysftp
  29. Using 100% of all cores with the multiprocessing module
  30. What should I put in the body of an abstract method?

Leave a Comment