2

作为更大模拟的一部分,我用 Python 编写了以下函数:

#!/usr/bin/python
counter = 1
while (counter < 10000):
    oldpa = .5
    t = 1
    while (t < counter):
        newpa = ((oldpa * t) + 1) / (t + 1)
        t = t + 1
        oldpa = newpa
    counter = counter + 1
    print str(counter) + "\t" + str(oldpa)

然后,我开始用 C 重写模拟,以便它运行得更快(同时也给自己一个花时间学习 C 的借口)。这是我的上述函数的 C 版本。

#include <stdio.h>
main()
{
  int counter, t;
  float oldpa, newpa;
  counter = 1;
  while ( counter < 10000 )
    {
      oldpa = .5;
      t = 1;
      while ( t < counter )
        {
          newpa = ((oldpa * t) + 1) / (t + 1);
          t = t + 1;
          oldpa = newpa;
        }
      counter = counter + 1;
      printf("%d\t%f\n", counter, oldpa);
   }
}

现在,有趣的事情来了。当我运行 Python 函数时,结果收敛到 0.999950,但是当我运行 C 函数时,它收敛到 0.999883。就我的模拟而言,这种差异实际上可以忽略不计,但我仍然想知道为什么我得到不同的结果

4

1 回答 1

3

Python 中的浮点值几乎总是 IEEE-754双精度,对应于 C 或 C++ double。如果您想要更高的精度,请查看小数模块

于 2013-03-23T12:54:52.343 回答