00001 
00002 
00003 
00004 
00005 
00006 
00007 
00008 
00009 
00010 
00011 
00012 
00013 #if defined(LIBM_SCCS) && !defined(lint)
00014 static char rcsid[] = "$NetBSD: w_fmod.c,v 1.6 1995/05/10 20:48:55 jtc Exp $";
00015 #endif
00016 
00017 
00018 
00019 
00020 
00021 #include "math.h"
00022 #include "mathP.h"
00023 
00024 
00025 #ifdef __STDC__
00026     double fmod(double x, double y) 
00027 #else
00028     double fmod(x,y)        
00029     double x,y;
00030 #endif
00031 {
00032 #ifdef _IEEE_LIBM
00033     return __ieee754_fmod(x,y);
00034 #else
00035     double z;
00036     z = __ieee754_fmod(x,y);
00037     if(_LIB_VERSION == _IEEE_ ||isnan(y)||isnan(x)) return z;
00038     if(y==0.0) {
00039             return __kernel_standard(x,y,27); 
00040     } else
00041         return z;
00042 #endif
00043 }