如何修复C++中负浮点数的取模计算函数？

我编写了一段用于计算两个浮点数取模的C++代码，希望结果始终落在[0, denom)区间（已知denom为大于0的整数）：

auto mod(float x, float denom)
{
    return x>= 0 ? std::fmod(x, denom) : denom + std::fmod(x + 1.0f, denom) - 1.0f;
}

但处理负浮点数输入时，部分场景结果不符合预期：比如输入x=-4.75、denom=4时输出-0.75（预期3.25），x=-4.5时输出-0.5（预期3.5）。测试代码如下：

for(size_t k = 0; k != 65; ++k)
{
    auto x = 0.25f*(static_cast<float>(k) - 32);
    printf("%.8g %.8g\n", x, mod(x, 4));
}

错误输出示例（关键错误行）：

-5 3
-4.75 -0.75 <== should be 3.25
-4.5 -0.5   <== should be 3.5
-4.25 -0.25 <== should be 3.75
-4 0

问题分析 #

原代码对负数的处理逻辑存在漏洞：当x是denom整数倍附近的负数时（比如x=-4.75），x+1.0f得到-3.75，std::fmod(-3.75,4)返回-3.75，代入公式denom + std::fmod(...) -1.0f后计算得4 + (-3.75) -1 = -0.75，与预期不符。本质是错误地试图通过调整x来修正std::fmod的符号问题，逻辑绕弯且覆盖不全。

修复方案 #

方案1：直接修正std::fmod的负结果 #

std::fmod的结果符号与输入的x一致，因此只需判断结果是否为负，若为负则加上denom即可得到[0, denom)区间的正确值：

#include <cmath>

auto mod(float x, float denom)
{
    float result = std::fmod(x, denom);
    return result >= 0 ? result : result + denom;
}

测试验证：x=-4.75时，std::fmod(-4.75,4)返回-0.75，加上4后得到3.25；x=-4时，std::fmod返回0，直接返回即可，完全符合预期。

方案2：基于数学定义的通用实现 #

取模的数学定义为：余数 = 被除数 - 商 × 除数，其中商是对x/denom向下取整的结果（std::floor实现）。这种方法无需依赖std::fmod的特性，逻辑更直观：

#include <cmath>

auto mod(float x, float denom)
{
    float quotient = std::floor(x / denom);
    return x - quotient * denom;
}

测试验证：x=-4.75时，x/denom = -1.1875，std::floor后得到-2，计算-4.75 - (-2)*4 = 3.25；x=-4时，x/denom=-1，std::floor后为-1，计算-4 - (-1)*4 = 0，结果正确。

总结 #

两种方案均能正确处理所有正负浮点数输入，方案1更简洁高效，方案2更贴合数学定义，可根据需求选择。

内容的提问来源于stack exchange，提问作者user877329