关于python浮点mod运算的解释。

luffy.deng

问题来自http://bbs.chinaunix.net/thread-1708133-1-1.html 这个帖子
但是我感觉帖子里对python浮点 % 运算并没有彻底讨论清楚，因此这里说一下我个人的理解。
关于python浮点mod运算的结果 <<core python>>里提到这样一个公式，x-(math.floor(x/y)*y)
这实际上是说python对浮点商是向下取整的，这个可以通过float_divmod的实现得到求证。
......
mod = fmod(vx, wx);
        /* fmod is typically exact, so vx-mod is *mathematically* an
           exact multiple of wx.  But this is fp arithmetic, and fp
           vx - mod is an approximation; the result is that div may
           not be an exact integral value after the division, although
           it will always be very close to one.
        */
        div = (vx - mod) / wx;  //注意这里
       
        if (mod) {
                /* ensure the remainder has the same sign as the denominator */
                if ((wx < 0) != (mod < 0)) {
                        mod += wx;
                        div -= 1.0;
                }
        }
......

2gua

其实是，关键要理解地板除的含义。

t6760915

我测试了下确实是，学习了

jjj137

嗯，感谢楼主分享，以前每次用之前先自己试一遍，汗……

可可熊

math.fmod(x, y)
Return fmod(x, y), as defined by the platform C library. Note that the Python expression x % y may not return the same result. The intent of the C standard is that fmod(x, y) be exactly (mathematically; to infinite precision) equal to x - n*y for some integer n such that the result has the same sign as x and magnitude less than abs(y). Python’s x % y returns a result with the sign of y instead, and may not be exactly computable for float arguments. For example, fmod(-1e-100, 1e100) is -1e-100, but the result of Python’s -1e-100 % 1e100 is 1e100-1e-100, which cannot be represented exactly as a float, and rounds to the surprising 1e100. For this reason, function fmod() is generally preferred when working with floats, while Python’s x % y is preferred when working with integers.

看下fmod的定义以及%的定义会更清楚些：
http://python.org/doc/current/ref/binary.html#l2h-411

luffy.deng

回复 5# 可可熊
从文档上看还真是按绝对值来的，我理解的有偏差。不过实际效果还是向下取整数。

zhuzhu316

回复 6# luffy.deng


    。。我还是想多问句，有中文手册吗？

t6760915

凡事google下先

jjj137

回复  luffy.deng


    。。我还是想多问句，有中文手册吗？
zhuzhu316 发表于 2010-05-25 11:33

    有也都是很旧的，而且手册这种东西，一般翻译了也会晦涩难懂吧……往往看起来比英文难得多。

daybreakcx

其实可以这么看，我们就是要确定k令：0 <= a + k * b < b
原来a==-25.5，b==2.25
当k==12时-25.5 + 12 * 2.25 == 1.5
这个值就是我们要求的
其他的也是一样的…………之所以会有地板除，是在找k的时候方便