/* * 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.2 07/08/85 */ float r_mod(x,y) float *x, *y; { double floor(), quotient = *x / *y; if (quotient >= 0.0) quotient = floor(quotient); else quotient = -floor(-quotient); return(*x - (*y) * quotient ); }