How to avoid floating point errors? [duplicate]

by

This really has nothing to do with Python – you’d see the same behavior in any language using your hardware’s binary floating-point arithmetic. First read the docs.

After you read that, you’ll better understand that you’re not adding one one-hundredth in your code. This is exactly what you’re adding:

>>> from decimal import Decimal
>>> Decimal(.01)
Decimal('0.01000000000000000020816681711721685132943093776702880859375')

That string shows the exact decimal value of the binary floating (“double precision” in C) approximation to the exact decimal value 0.01. The thing you’re really adding is a little bigger than 1/100.

Controlling floating-point numeric errors is the field called “numerical analysis”, and is a very large and complex topic. So long as you’re startled by the fact that floats are just approximations to decimal values, use the decimal module. That will take away a world of “shallow” problems for you. For example, given this small modification to your function:

from decimal import Decimal as D

def sqrt(num):
    root = D(0)
    while root * root < num:
        root += D("0.01")
    return root

then:

>>> sqrt(4)
Decimal('2.00')
>>> sqrt(9)
Decimal('3.00')

It’s not really more accurate, but may be less surprising in simple examples because now it’s adding exactly one one-hundredth.

An alternative is to stick to floats and add something that is exactly representable as a binary float: values of the form I/2**J. For example, instead of adding 0.01, add 0.125 (1/8) or 0.0625 (1/16).

Then look up “Newton’s method” for computing square roots 😉

More Related Contents:

  1. How to properly round-up half float numbers?
  2. rounding errors in Python floor division
  3. Convert float to string in positional format (without scientific notation and false precision)
  4. What is the difference between int() and floor() in Python 3?
  5. Clarification on the Decimal type in Python
  6. Rounding floats so that they sum to precisely 1
  7. Print a float number in normal form, not exponential form / scientific notation [duplicate]
  8. ‘is’ operator behaves unexpectedly with floats
  9. a pythonic way of doubling consecutive binary bits [closed]
  10. How to display position of an element in an array? [closed]
  11. How to make Firefox headless programmatically in Selenium with Python?
  12. NameError: global name ‘unicode’ is not defined – in Python 3
  13. Installing scipy in Python 3.5 on 32-bit Windows 7 Machine
  14. How to run pip from different versions of python using the python command?
  15. How to melt 2 columns at the same time?
  16. i am sending commands through serial port in python but they are sent multiple times instead of one
  17. Why is there no xrange function in Python3?
  18. Get all keys of a nested dictionary [duplicate]
  19. Type hints when unpacking a tuple?
  20. Sorting by absolute value without changing to absolute value
  21. Why isn’t .ico file defined when setting window’s icon?
  22. Printing result of the zip() function in Python 3 gives “zip object at ” [duplicate]
  23. FutureWarning: The default value of regex will change from True to False in a future version
  24. Multiple keys per value
  25. Python making combined bar and line plot with secondary y-axis
  26. Is there a way to use asyncio.Queue in multiple threads?
  27. match dates using python regular expressions
  28. How to get value of variable entered from user input?
  29. Add custom form fields that are not part of the model (Django)
  30. Writing to a new file if it doesn’t exist, and appending to a file if it does

Leave a Comment