将浮点数输入到字符串流时出现“浮点无效操作”

Question

我有一段简单的代码，它从 FORTRAN 生成的 REAL 数组中提取浮点数，然后将其插入到流中以进行日志记录。虽然这适用于前 30 个案例，但在 31 日它会因“浮点无效操作”而崩溃。

代码是：

int FunctionDeclaration(float* mrSwap)
{
...
float swap_float;
stringstream message_stream;
...
swap_float = *(mrSwap+30-1);
...
message_stream.clear();
message_stream <<  30 << "\t" << swap_float << "\tblah blah blah \t";

在调试时，实例崩溃前的 swap_float 的值（在上面的最后一行）是 1711696.3 - 除了这比到目前为止的大多数值大得多之外，它并没有什么特别的地方。

我也试过用 cerr 替换 message_stream，也遇到了同样的问题。迄今为止，我一直认为 cerr 几乎是坚不可摧的 - 一个简单的浮子怎么能摧毁它？

编辑：

感谢 cmets：我添加了 mrSwap 的声明。 mrSwap 大约有 200 长，所以我离结束还有很长的路要走。它在我的控制之外填充，并且单个条目可能不会被填充 - 但据我所知，这只是意味着 swap_float 将被设置为随机浮点数？

Answer 1

我只是添加此答案以突出显示上述 TonyK 答案中的正确解决方案 - 因为我们做了几个循环，答案已被编辑，并且由于 cmets 中有几个要点，实际答案可能不会立即出现明显。所有功劳都应归功于 TonyK 的解决方案。

“另一种可能是延迟浮点陷阱。浮点协处理器（这些天内置在 CPU 中）由于某些非法操作而产生浮点中断，但中断并没有在尝试下一个浮点运算之前会被注意到。所以这个崩溃可能是上一个浮点运算的结果，它可能在任何地方。" - TonyK

这确实是问题所在：在我使用 IsSame 进行的比较中，另一个值是 NaN（在这种情况下这是一个有效值），尽管它很高兴地从 swap_float 中减去它，但它设置了一个标志来表示报告下一个操作为错误。我不得不说我完全没有意识到这是可能的——我认为如果它有效，它就会有效。

Answer 2

个人条目可能不会被填充 - 但尽我所能 理解，这只是意味着 swap_float 将被设置为 随机浮动？

绝对不是。 IEEE 浮点数中的某些位模式表示无效的数字——例如，溢出算术运算的结果，或无效的结果（例如 0.0/0.0）。这里令人费解的是，调试器显然接受该数字为有效，而cout 不接受。

尝试获取swap_float 的位布局。在 32 位系统上：

int i = *(int*)&swap_float;

然后以十六进制打印i，让我们知道你看到了什么。

更新为添加：来自 Mike 的评论，i=1238430338，即 49D0F282 十六进制。这是一个有效的浮点数，正好等于 1711696.25。所以我不知道发生了什么事，我很害怕。我唯一可以建议的是，编译器可能正在将无效的浮点数直接从mrSwap 数组加载到浮点寄存器组中，而不通过swapFloat。因此，swapFloat 的真实值对调试器根本不可用。要检查这一点，请尝试

int j = *(int*)(mrSwap+30-1);

告诉我们你看到了什么。

再次更新添加：另一种可能性是延迟浮点陷阱。浮点协处理器（这些天内置在 CPU 中）由于某些非法操作而产生浮点中断，但直到 下一个 浮点操作完成后才会注意到该中断尝试过。所以这个崩溃可能是 previous 浮点运算的结果，它可能在任何地方。祝你好运……