索引运行时间间隔错误

Question

我不明白为什么会发生这种情况。ngrid将是10的倍数。对于ngrid=10，比如说，当我打印出0.7*ngrid时，我得到7，这当然是正确的。但当我打印出来时，它会打印出6以及7&8.为什么会发生这种情况？提前谢谢。

    for(int i=(0.7*ngrid); i<0.8*ngrid; i++)
    {
    cout << i;
    operator()(0.2*ngrid,i) = -3.0;
    }

ngrid是int类型，在构造函数中声明

class Field
{
  private:
    double* data;
    int     n;
    int index(int x, int y) const { return x + (n+1)*y; }
  public:

    Field(int ngrid)  : n(ngrid)
    {
      data = new double[(ngrid+1)*(ngrid+1)];
      int l=(ngrid+1)*(ngrid+1);
      for (int i=0;i<l;i++) data[i]=0.0;
    }

Answer 1

你在这里做了一些危险的事情。首先，期望浮点运算的结果是一个精确的整数。事实证明，即使ngrid是int类型，0.7也是不能用IEEE浮点格式表示的。如果使用此页面，您可以看到0.7实际上表示为0.699999988070971。乘以10，然后将其转换为int，现在得到6。

类似地，0.8表示为0.800000011920929。因此，当将8与8.00000011920929进行比较时，可以正确地确定它仍然较小，因此也打印了8。

因此，当for循环开始执行时，i被初始化为6而不是7。这就是为什么6被打印出来的原因。这个故事的寓意是，不要指望浮点运算能很好地计算出整数。。。即使你只处理小数点后一位。这台计算机不能用以10为基数的数字工作，所以并不是所有的数字都能很容易地表示出来。

与其进行浮点运算（假设ngrid是int），不如将其乘以7，然后除以7。这将执行整数运算，而不是浮点运算。既然你说ngrid总是10的倍数，那么结果应该总是一个整数，你不必担心截断。

Answer 2

不要混合整数和浮点/双精度。for仅用于数组索引或其他整数walkthru。使用while。

或者，使用固定点，将所有内容乘以，比如100。

编辑：当程序员在代码中发现for时，他们会认为"这是一种对索引或整数坐标的遍历"。如果你有这样的问题，比如当x*0.7达到一个水平时，循环中断，最好使用while，这意味着"我们正在重复一些事情，直到发生某种传导"。有些语言甚至不允许在for中使用非整数。

Answer 3

我的谦虚意见：

将0.7*ngrid更改为（7*ngrid）/10，并将0.8*ngridd更改为（8*ng里德）/10

所以你只有整数

int i=0.7*ngrid由于0.7而将ngrid强制转换为浮点/双精度，然后再次强制转换为int，因为i是延迟整数，将乘积四舍五入到较低的整数。