fmod, fmodf, fmodl

此函数计算的除法 x/y 的浮点余数是 x - n*y 的准确值，其中 n 是截断小数部分的 x/y 。

返回值与 x 拥有相同符号，且绝对值小于 y 。

参数

返回值

若成功，则返回定义于上的除法 x/y 的浮点余数。

若出现定义域错误，则返回实现定义值（受支持的平台上为 NaN ）。

若出现下溢所指定值域错误，则返回（舍入后的）正确结果。

错误处理

报告 math_errhandling 中指定的错误。

若 y 为零则可能发生定义域错误。

若实现支持 IEEE 浮点算术（ IEC 60559 ），

• 若 x 为 ±0 且 y 非零，则返回 ±0
• 若 x 为 ±∞ 且 y 非 NaN ，则返回 NaN 并引发 FE_INVALID
• 若 y 为 ±0 且 x 非 NaN ，则返回 NaN 并引发 FE_INVALID
• 若 y 为 ±∞ 且 x 有限，则返回 x 。
• 若任一参数为 NaN ，则返回 NaN 。

注意

POSIX 要求若 x 为无穷大或 y 为零则出现定义域错误。

fmod ，但不是 remainder ，适于安静地包装浮点类型到无符号整数类型： (0.0 <= (y = fmod( rint(x), 65536.0 )) ? y : 65536.0 + y) 在范围 [-0.0 .. 65535.0] 内，它对应 unsigned short ，但 remainder(rint(x), 65536.0) 在范围 [-32767.0, +32768.0] 内，它在 signed short 的范围外。

fmod 的 double 版本表现为如同实现如下：

double fmod(double x, double y)
{
#pragma STDC FENV_ACCESS ON
    double result = remainder(fabs(x), (y = fabs(y)));
    if (signbit(result)) result += y;
    return copysign(result, x);
}

示例

#include <stdio.h>
#include <math.h>
#include <fenv.h>
 
#pragma STDC FENV_ACCESS ON
int main(void)
{
    printf("fmod(+5.1, +3.0) = %.1f\n", fmod(5.1,3));
    printf("fmod(-5.1, +3.0) = %.1f\n", fmod(-5.1,3));
    printf("fmod(+5.1, -3.0) = %.1f\n", fmod(5.1,-3));
    printf("fmod(-5.1, -3.0) = %.1f\n", fmod(-5.1,-3));
 
    // 特殊值
    printf("fmod(+0.0, 1.0) = %.1f\n", fmod(0, 1));
    printf("fmod(-0.0, 1.0) = %.1f\n", fmod(-0.0, 1));
    printf("fmod(+5.1, Inf) = %.1f\n", fmod(5.1, INFINITY));
 
    // 错误处理
    feclearexcept(FE_ALL_EXCEPT);
    printf("fmod(+5.1, 0) = %.1f\n", fmod(5.1, 0));
    if(fetestexcept(FE_INVALID)) puts("    FE_INVALID raised");
}

fmod(+5.1, +3.0) = 2.1
fmod(-5.1, +3.0) = -2.1
fmod(+5.1, -3.0) = 2.1
fmod(-5.1, -3.0) = -2.1
fmod(+0.0, 1.0) = 0.0
fmod(-0.0, 1.0) = -0.0
fmod(+5.1, Inf) = 5.1
fmod(+5.1, 0) = nan
    FE_INVALID raised

引用