usr/src/usr.bin/f77/libF77/r_mod.c

/*
 * Copyright (c) 1980 Regents of the University of California.
 * All rights reserved.  The Berkeley software License Agreement
 * specifies the terms and conditions for redistribution.
 *
 *      @(#)r_mod.c     5.3     %G%
 */

float flt_retval;

float r_mod(x,y)
float *x, *y;
{
double floor(), quotient = *x / *y;
if (quotient >= 0.0)
        quotient = floor(quotient);
else
        quotient = -floor(-quotient);
flt_retval = *x - (*y) * quotient ;
return(flt_retval);
}