The R Project SVN R

Rev

Rev 574 | Rev 5458 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 574 Rev 580
Line 33... Line 33...
33 
 
 33 
 
34 
#include "Mathlib.h"
 34 
#include "Mathlib.h"
35 
 
 35 
 
36 
double fmod(double x, double y)
 36 
double fmod(double x, double y)
37 
{
 37 
{
- 
 
 38 
    double quot;
- 
 
 39 
#ifdef IEEE_754
- 
 
 40 
    if (ISNAN(x) || ISNAN(y))
- 
 
 41 
	return x + y;
- 
 
 42 
#endif
38 
    double quot = x / y;
 43 
    quot = x / y;
39 
    return x - (quot < 0.0 ? ceil(quot) : floor(quot)) * y;
 44 
    return x - (quot < 0.0 ? ceil(quot) : floor(quot)) * y;
40 
}
 45 
}