Path: csiph.com!feeder.erje.net!2.eu.feeder.erje.net!newsfeed.kamp.net!newsfeed.kamp.net!fu-berlin.de!uni-berlin.de!not-for-mail
From: Chris Angelico <rosuav@gmail.com>
Newsgroups: comp.lang.python
Subject: Re: What exactly is "exact" (was Clean Singleton Docstrings)
Date: Mon, 18 Jul 2016 20:37:23 +1000
Lines: 50
Message-ID: <mailman.72.1468838251.2307.python-list@python.org>
References: <nlmpl5$9l4$1@dont-email.me> <nlnl9m$j3b$1@ger.gmane.org> <mailman.137.1467963521.2295.python-list@python.org> <df9aaf24-394f-451a-9518-7e3a8dc67318@googlegroups.com> <nm6kbb$t6b$1@ger.gmane.org> <mailman.51.1468454075.21009.python-list@python.org> <7046c74e-e5ea-4dde-8847-8c556756a563@googlegroups.com> <5789B601.9000604@stoneleaf.us> <mailman.26.1468642801.2307.python-list@python.org> <9d2a0934-bf26-42a8-9bb2-e4e75b2d2ad3@googlegroups.com> <CAPTjJmq+VKdu7x_nzfP8fb=0EiAqWJO6_qKNsabQSnbBStyMsA@mail.gmail.com> <1468652038.1427719.667878897.460955C0@webmail.messagingengine.com> <mailman.28.1468652041.2307.python-list@python.org> <578a027c$0$1616$c3e8da3$5496439d@news.astraweb.com> <69d8d886-c18a-4c05-8bca-9a05afacd9a9@googlegroups.com> <578ca188$0$1505$c3e8da3$5496439d@news.astraweb.com> <877fcjcllf.fsf@elektro.pacujo.net> <CAPTjJmqYuhfR68g4EWWfj=GQSQMBmo3r3DRCej39xLtjPfR2ow@mail.gmail.com> <mailman.71.1468836914.2307.python-list@python.org> <402d8e8a-6acf-44e8-99e2-fa09da10d7bc@googlegroups.com> <CAPTjJmoz_neVM0kgWqnman_ecp4-Hjg_re8RuhJB1OeHAMRorQ@mail.gmail.com>
Mime-Version: 1.0
Content-Type: text/plain; charset=UTF-8
In-Reply-To: <402d8e8a-6acf-44e8-99e2-fa09da10d7bc@googlegroups.com>
Precedence: list
Xref: csiph.com comp.lang.python:111602

On Mon, Jul 18, 2016 at 8:24 PM, Rustom Mody <rustompmody@gmail.com> wrote:
> I dont know what point you are trying to make
> Here is behavior.  Should one use == ??
>
> Python 2.7.11+ (default, Apr 17 2016, 14:00:29)
> [GCC 5.3.1 20160413] on linux2
> Type "help", "copyright", "credits" or "license" for more information.
>>>> .1+.1+.1 == .3
> False
>
> EPSILON = 0.000001 # Adjust to control numeric accuracy
> def is_equal(f1, f2, epsilon=EPSILON):
>     if abs(f1) > abs(f2):
>         f1, f2 = f2, f1
>     return abs(f2-f1) < f1*epsilon
>>>> ... ... ... ...
>
>>>> is_equal(.1+.1+.1, .3)
> True
>>>>

Sure, simple equality hasn't given you an "intuitive" result - but
that's because your source code is in decimal. So you need to properly
understand it. If you work with something that can be represented
precisely in binary, you have no problems with equality:

>>> 7/4 + 9/4 + 11/4 + 13/4 == 10
True

It's only the repeating values that have that problem. And if you
disassemble your example, it's obvious why you get False:

>>> dis.dis(lambda: .1+.1+.1 == .3)
  1           0 LOAD_CONST               4 (0.30000000000000004)
              2 LOAD_CONST               2 (0.3)
              4 COMPARE_OP               2 (==)
              6 RETURN_VALUE

Of course those two values aren't the same. And it'd be just as
obvious if you looked at them in binary. It'd look like this in
decimal:

0.6667 + 0.6667 + 0.6667 == 2.0

Well, duh, all those intermediate values are rounded up, so you're
getting the sum of the rounding, and of course it's not equal. That's
why you need to get some comprehension of floating-point and how it
*really* functions.

ChrisA