math - 浮点算术四舍五入到最接近

Question

我对浮点运算中的四舍五入感到有些困惑。令 a、b 和 c 为归一化的双精度浮点数。a+b=b+a 是否正确，其中 + 正确舍入到最接近的浮点加法？我最初的猜测是肯定的，这总是正确的，但我不完全理解四舍五入到最近。有人可以举一个例子，当 a+b != b+a 使用浮点加法并四舍五入到最近？

Answer 1

无论舍入模式如何，正确实现的 IEEE-754 浮点加法都是可交换的（a+b 等于 b+a）。

舍入模式会影响精确的数学结果如何舍入以适合目标格式。由于 a+b 和 b+a 的精确数学结果相同，因此它们的舍入相同。

Answer 2

如上所述，加法是可交换的，但不是结合的。通过运行以下 (MS Visual Studio) C++ 代码可以看出舍入模式的差异：

#include <iostream>
#include <float.h>

#pragma fenv_access(on)

using namespace std;

int main(int argc, char* argv[])
{
    float a = 1.75f, b = 1e-6f;

    cout.setf(ios::fixed,ios::floatfield);
    cout.precision(7);

    cout << "a = " << a << ", b = " << b << endl;

    _controlfp_s(NULL, _RC_DOWN,_MCW_RC);

    cout << "Result of a + b rounded down: " << a+b << endl;
    cout << "Result of a - b rounded down: " << a-b << endl;

    _controlfp_s(NULL, _RC_UP,_MCW_RC);
    cout << "Result of a + b rounded up: " << a+b << endl;
    cout << "Result of a - b rounded up: " << a-b << endl;

    _controlfp_s(NULL, _RC_NEAR,_MCW_RC);
    cout << "Result of a + b rounded to nearest: " << a+b << endl;
    cout << "Result of a - b rounded to nearest: " << a-b << endl;

    return 0;
}

输出：

a = 1.7500000, b = 0.0000010
Result of a + b rounded down: 1.7500010
Result of a - b rounded down: 1.7499989
Result of a + b rounded up: 1.7500011
Result of a - b rounded up: 1.7499990
Result of a + b rounded to nearest: 1.7500010
Result of a - b rounded to nearest: 1.7499990