Return to w_remainderf.c CVS log

Up to [cvs.NetBSD.org] / src / lib / libm / src

Annotation of src/lib/libm/src/w_remainderf.c, Revision 1.4

1.1       jtc         1: /* w_remainderf.c -- float version of w_remainder.c.
                      2:  * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
                      3:  */
                      4:
                      5: /*
                      6:  * ====================================================
                      7:  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
                      8:  *
                      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:  * ====================================================
                     14:  */
                     15:
1.4     ! lukem      16: #include <sys/cdefs.h>
1.2       jtc        17: #if defined(LIBM_SCCS) && !defined(lint)
1.4     ! lukem      18: __RCSID("$NetBSD$");
1.1       jtc        19: #endif
                     20:
                     21: /*
                     22:  * wrapper remainderf(x,p)
                     23:  */
                     24:
                     25: #include "math.h"
                     26: #include "math_private.h"
                     27:
                     28: #ifdef __STDC__
                     29:        float remainderf(float x, float y)      /* wrapper remainder */
                     30: #else
                     31:        float remainderf(x,y)                   /* wrapper remainder */
                     32:        float x,y;
                     33: #endif
                     34: {
                     35: #ifdef _IEEE_LIBM
                     36:        return __ieee754_remainderf(x,y);
                     37: #else
                     38:        float z;
                     39:        z = __ieee754_remainderf(x,y);
                     40:        if(_LIB_VERSION == _IEEE_ || isnanf(y)) return z;
                     41:        if(y==(float)0.0)
                     42:            /* remainder(x,0) */
                     43:            return (float)__kernel_standard((double)x,(double)y,128);
                     44:        else
                     45:            return z;
                     46: #endif
                     47: }