浮动分区和隐式回合

时间:2012-10-09 17:36:59

标签: c++ casting floating-point

我在我编写的一段代码中发现了一个奇怪的行为,当我尝试将其从程序中重现时(即使用MATLAB),这种行为不会发生。

我正在对数组进行如下操作:

a[i] = a[i] / (sqrt(b[i]) * sqrt(c[i]));

其中abc是float类型的数组。 abc之前除法的值范围从20006000000

对于数组的所有元素,操作的结果为零,但是分别查看分子和分母的值,对于给定元素,操作为89509.0 / 90240,375,等于{ {1}}。

所以,我想知道为什么使用float的程序中的完整操作不起作用,并且它按预期单独工作。可能是因为几个小时的混淆调试?

1 个答案:

答案 0 :(得分:7)

当我遇到类似的情况时,我发现通常有助于将其分解为步骤,看看哪些步骤正常,哪些步骤不正常。例如,为此,它可以分解为以下内容:

float sqrtb = sqrt(b[i]);
float sqrtc = sqrt(c[i]);
float divisor = srqtb * sqrtc;
float dividend = a[i];
float quotient = dividend / divisor;
a[i] = quotient;

现在我可以浏览每个项目并查看它失败的位置(如果失败)。事实上,当我这样做时,我经常在创建步骤时发现问题,然后修复它。