home *** CD-ROM | disk | FTP | other *** search
/ Programmer 7500 / MAX_PROGRAMMERS.iso / INFO / FORTRAN / F2C_SRC.ZIP / F2C / LIBF77 / R_MOD.C < prev    next >
Encoding:
C/C++ Source or Header  |  1991-06-10  |  443 b   |  27 lines
  1. #include "f2c.h"
  2.  
  3. double r_mod(x,y)
  4. real *x, *y;
  5. {
  6. #ifdef IEEE_drem
  7.     double drem(), xa, ya, z;
  8.     if ((ya = *y) < 0.)
  9.         ya = -ya;
  10.     z = drem(xa = *x, ya);
  11.     if (xa > 0) {
  12.         if (z < 0)
  13.             z += ya;
  14.         }
  15.     else if (z > 0)
  16.         z -= ya;
  17.     return z;
  18. #else
  19.     double floor(), quotient;
  20.     if( (quotient = (double)*x / *y) >= 0)
  21.         quotient = floor(quotient);
  22.     else
  23.         quotient = -floor(-quotient);
  24.     return(*x - (*y) * quotient );
  25. #endif
  26. }
  27.