4BSD/usr/src/lib/libF77/d_mod.c

Compare this file to the similar file:
Show the results in this format:

double d_mod(x,y)
double *x, *y;
{
double floor(), quotient;
if( (quotient = *x / *y) >= 0)
	quotient = floor(quotient);
else
	quotient = -floor(-quotient);
return(*x - (*y) * quotient );
}