double d_mod(x
,y
) doublereal
*x
, *y
;
double drem(double, double);
double d_mod(doublereal
*x
, doublereal
*y
)
if( (quotient
= *x
/ *y
) >= 0)
quotient
= floor(quotient
);
quotient
= -floor(-quotient
);
return(*x
- (*y
) * quotient
);