fmod, fmodf, fmodl

El resto en punto flotante de la operación de división x / y calculado por esta función es exactamente el valor x - n * y , donde n es x / y con su parte fraccionaria truncada.

El valor devuelto tiene el mismo signo que x y es menor o igual que y en magnitud.

Parámetros

Valor de retorno

Si tiene éxito, devuelve el resto en punto flotante de la división x / y como se definió anteriormente.

Si ocurre un error de dominio, se devuelve un valor definido por la implementación (NaN donde esté soportado).

Si ocurre un error de rango debido a desbordamiento inferior, se devuelve el resultado correcto (después del redondeo).

Manejo de errores

Los errores se reportan como se especifica en math_errhandling .

Puede ocurrir un error de dominio si y es cero.

Si la implementación soporta aritmética de punto flotante IEEE (IEC 60559):

• Si x es ±0 y y no es cero, se retorna ±0.
• Si x es ±∞ y y no es NaN, se retorna NaN y se genera FE_INVALID .
• Si y es ±0 y x no es NaN, se retorna NaN y se genera FE_INVALID .
• Si y es ±∞ y x es finito, x es retornado.
• Si cualquier argumento es NaN, se retorna NaN.

Notas

POSIX requiere que ocurra un error de dominio si x es infinito o y es cero.

fmod , pero no remainder es útil para realizar el ajuste silencioso de tipos de punto flotante a tipos enteros sin signo: ( 0.0 <= ( y = fmod ( rint ( x ) , 65536.0 ) ) ? y : 65536.0 + y ) está en el rango [ - 0.0 , 65535.0 ] , que corresponde a unsigned short , pero remainder ( rint ( x ) , 65536.0 ) está en el rango [ - 32767.0 , + 32768.0 ] , que está fuera del rango de signed short .

La versión double de fmod se comporta como si estuviera implementada de la siguiente manera:

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);
}

Ejemplo

#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));
    // special values
    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));
    // error handling
    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

Referencias

Véase también