Return to w_asinf.c CVS log

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

Annotation of src/lib/libm/src/w_asinf.c, Revision 1.2

1.1       jtc         1: /* w_asinf.c -- float version of w_asin.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.2     ! jtc        16: #if defined(LIBM_SCCS) && !defined(lint)
        !            17: static char rcsid[] = "$Id: w_asinf.c,v 1.1 1994/08/10 20:33:34 jtc Exp $";
1.1       jtc        18: #endif
                     19:
                     20: /*
                     21:  * wrapper asinf(x)
                     22:  */
                     23:
                     24:
                     25: #include "math.h"
                     26: #include "math_private.h"
                     27:
                     28:
                     29: #ifdef __STDC__
                     30:        float asinf(float x)            /* wrapper asinf */
                     31: #else
                     32:        float asinf(x)                  /* wrapper asinf */
                     33:        float x;
                     34: #endif
                     35: {
                     36: #ifdef _IEEE_LIBM
                     37:        return __ieee754_asinf(x);
                     38: #else
                     39:        float z;
                     40:        z = __ieee754_asinf(x);
                     41:        if(_LIB_VERSION == _IEEE_ || isnanf(x)) return z;
                     42:        if(fabsf(x)>(float)1.0) {
                     43:            /* asinf(|x|>1) */
                     44:            return (float)__kernel_standard((double)x,(double)x,102);
                     45:        } else
                     46:            return z;
                     47: #endif
                     48: }