.
[glibc/history.git] / math / w_fmodf.c
blobe2e869a3ccac49c77c71bbd70bd3eed0fa9e62ff
1 /* w_fmodf.c -- float version of w_fmod.c.
2 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
3 */
5 /*
6 * ====================================================
7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
9 * Developed at SunPro, a Sun Microsystems, Inc. business.
10 * Permission to use, copy, modify, and distribute this
11 * software is freely granted, provided that this notice
12 * is preserved.
13 * ====================================================
16 #if defined(LIBM_SCCS) && !defined(lint)
17 static char rcsid[] = "$NetBSD: w_fmodf.c,v 1.3 1995/05/10 20:48:57 jtc Exp $";
18 #endif
21 * wrapper fmodf(x,y)
24 #include <math.h>
25 #include "math_private.h"
28 #ifdef __STDC__
29 float __fmodf(float x, float y) /* wrapper fmodf */
30 #else
31 float __fmodf(x,y) /* wrapper fmodf */
32 float x,y;
33 #endif
35 #ifdef _IEEE_LIBM
36 return __ieee754_fmodf(x,y);
37 #else
38 float z;
39 z = __ieee754_fmodf(x,y);
40 if(_LIB_VERSION == _IEEE_ ||__isnanf(y)||__isnanf(x)) return z;
41 if(__isinff(x)||y==(float)0.0) {
42 /* fmodf(+-Inf,y) or fmodf(x,0) */
43 return (float)__kernel_standard((double)x,(double)y,127);
44 } else
45 return z;
46 #endif
48 weak_alias (__fmodf, fmodf)