我在我编写的一段代码中发现了一个奇怪的行为,当我尝试将其从程序中重现时(即使用MATLAB),这种行为不会发生。
我正在对数组进行如下操作:
a[i] = a[i] / (sqrt(b[i]) * sqrt(c[i]));
其中a
,b
和c
是float类型的数组。 a
,b
和c
之前除法的值范围从2000
到6000000
对于数组的所有元素,操作的结果为零,但是分别查看分子和分母的值,对于给定元素,操作为89509.0 / 90240,375
,等于{ {1}}。
所以,我想知道为什么使用float的程序中的完整操作不起作用,并且它按预期单独工作。可能是因为几个小时的混淆调试?
答案 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;
现在我可以浏览每个项目并查看它失败的位置(如果失败)。事实上,当我这样做时,我经常在创建步骤时发现问题,然后修复它。