fmod, fmodf, fmodl

Вещественный остаток от операции деления x / y , вычисляемый данной функцией, точно равен значению x - n * y , где n является результатом x / y с усечённой дробной частью.

Возвращаемое значение имеет тот же знак, что и x , и по модулю меньше или равно y .

Параметры

Возвращаемое значение

В случае успеха возвращает остаток от деления с плавающей запятой x / y в соответствии с определением выше.

Если происходит ошибка области определения, возвращается значение, определяемое реализацией (NaN, если поддерживается).

Если происходит ошибка диапазона из-за потери значимости (underflow), возвращается корректный результат (после округления).

Обработка ошибок

Ошибки сообщаются, как указано в 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 .

double версия функции fmod ведёт себя так, как если бы была реализована следующим образом:

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 <fenv.h>
#include <math.h>
#include <stdio.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

Ссылки

Смотрите также