lib/libm/fmod.c

   1 #include "mconf.h"
   2 #include "mathf.h"
   3
   4 double fmod(double x, double y) {
   5         if (y == 0) {
   6                 mtherr("fmod", DOMAIN);
   7                 return 0.0;
   8         } else {
   9                 double div;
  10                 div = x/y;
  11                 if (div < 0)
  12                         div = -floor(-div);
  13                 else
  14                         div = floor(div);
  15                 return (x - div * y);
  16         }
  17 }