std:: fmod, std:: fmodf, std:: fmodl

Вещественный остаток от операции деления x / y , вычисляемый данной функцией, точно соответствует значению x - iquot * y , где iquot является результатом 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 равен нулю.

std::fmod , но не std::remainder полезен для бесшумного оборачивания типов с плавающей точкой в беззнаковые целочисленные типы: ( 0.0 <= ( y = std :: fmod ( std:: rint ( x ) , 65536.0 ) ) ? y : 65536.0 + y ) находится в диапазоне [ - 0.0 , 65535.0 ] , что соответствует unsigned short , но std:: remainder ( std:: rint ( x ) , 65536.0 находится в диапазоне [ - 32767.0 , + 32768.0 ] , что выходит за пределы диапазона signed short .

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

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

Выражение x - std:: trunc ( x / y ) * y может не равняться std :: fmod ( x, y ) , когда округление x / y для инициализации аргумента std::trunc теряет слишком много точности (пример: x = 30.508474576271183309 , y = 6.1016949152542370172 ).

Дополнительные перегрузки не обязаны быть предоставлены в точности как (A) . Они должны быть лишь достаточными для обеспечения того, чтобы для их первого аргумента num1 и второго аргумента num2 :

Пример

#include <cfenv>
#include <cmath>
#include <iostream>
// #pragma STDC FENV_ACCESS ON
int main()
{
    std::cout << "fmod(+5.1, +3.0) = " << std::fmod(5.1, 3) << '\n'
              << "fmod(-5.1, +3.0) = " << std::fmod(-5.1, 3) << '\n'
              << "fmod(+5.1, -3.0) = " << std::fmod(5.1, -3) << '\n'
              << "fmod(-5.1, -3.0) = " << std::fmod(-5.1, -3) << '\n';
    // специальные значения
    std::cout << "fmod(+0.0, 1.0) = " << std::fmod(0, 1) << '\n'
              << "fmod(-0.0, 1.0) = " << std::fmod(-0.0, 1) << '\n'
              << "fmod(5.1, Inf) = " << std::fmod(5.1, INFINITY) << '\n';
    // обработка ошибок
    std::feclearexcept(FE_ALL_EXCEPT);
    std::cout << "fmod(+5.1, 0) = " << std::fmod(5.1, 0) << '\n';
    if (std::fetestexcept(FE_INVALID))
        std::cout << "    FE_INVALID raised\n";
}

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
fmod(-0.0, 1.0) = -0
fmod(5.1, Inf) = 5.1
fmod(+5.1, 0) = -nan
    FE_INVALID raised

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