Return to w_remainder.c CVS log

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

Annotation of src/lib/libm/src/w_remainder.c, Revision 1.5

1.1       jtc         1: /* @(#)w_remainder.c 5.1 93/09/24 */
                      2: /*
                      3:  * ====================================================
                      4:  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
                      5:  *
                      6:  * Developed at SunPro, a Sun Microsystems, Inc. business.
                      7:  * Permission to use, copy, modify, and distribute this
                      8:  * software is freely granted, provided that this notice
                      9:  * is preserved.
                     10:  * ====================================================
                     11:  */
1.3       jtc        12:
1.5     ! jtc        13: #if defined(LIBM_SCCS) && !defined(lint)
        !            14: static char rcsid[] = "$Id: w_remainder.c,v 1.4 1994/08/10 20:35:32 jtc Exp $";
1.3       jtc        15: #endif
1.1       jtc        16:
                     17: /*
                     18:  * wrapper remainder(x,p)
                     19:  */
                     20:
1.4       jtc        21: #include "math.h"
                     22: #include "math_private.h"
1.1       jtc        23:
                     24: #ifdef __STDC__
                     25:        double remainder(double x, double y)    /* wrapper remainder */
                     26: #else
                     27:        double remainder(x,y)                   /* wrapper remainder */
                     28:        double x,y;
                     29: #endif
                     30: {
                     31: #ifdef _IEEE_LIBM
                     32:        return __ieee754_remainder(x,y);
                     33: #else
                     34:        double z;
                     35:        z = __ieee754_remainder(x,y);
                     36:        if(_LIB_VERSION == _IEEE_ || isnan(y)) return z;
                     37:        if(y==0.0)
                     38:            return __kernel_standard(x,y,28); /* remainder(x,0) */
                     39:        else
                     40:            return z;
                     41: #endif
                     42: }