如何使用fmod并避免精度问题
我将把这个问题归结为最简单的forms:
让我们以0.05的步长从[0 .. 5.0]迭代,并且每0.25个乘数打印出’X’。
for(double d=0.0; d<=5.0; d+=0.05) { if(fmod(d,0.25) is equal 0) print 'X'; }
这当然不起作用,因为d
将是[0,0.05000000001,0.100000000002,…]导致fmod()失败。 极端的例子是当d=1.999999999998
和fmod(d,0.25) = 1
。
如何解决这个问题? 这是一个可编辑的在线示例。
我只是不以这种方式使用浮点变量来解决这个问题:
for (int i = 0; i <= 500; i += 5) { double d = i / 100.0; // in case you need to use it. if ((i % 25) == 0) print 'X'; }
它们通常存在问题,因此在避免它们方面需要额外的努力。
如果你总是有相同的固定小数,只需乘以你的循环计数器(在这种情况下乘以20)。 如果要以doubleforms访问它,除以20.实际上,您使用隐藏指数将双精度值保持为整数。
在这个例子中,我使用一个整数作为循环计数器,假设它具有所需的精度。
for(int i=0; i<=100; i+=1) { if(i % 5 == 0) print 'X'; double d = (double)i / 20.0; // use d }