w_logf.c revision 1.1 1 1.1 jtc /* w_logf.c -- float version of w_log.c.
2 1.1 jtc * Conversion to float by Ian Lance Taylor, Cygnus Support, ian (at) cygnus.com.
3 1.1 jtc */
4 1.1 jtc
5 1.1 jtc /*
6 1.1 jtc * ====================================================
7 1.1 jtc * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8 1.1 jtc *
9 1.1 jtc * Developed at SunPro, a Sun Microsystems, Inc. business.
10 1.1 jtc * Permission to use, copy, modify, and distribute this
11 1.1 jtc * software is freely granted, provided that this notice
12 1.1 jtc * is preserved.
13 1.1 jtc * ====================================================
14 1.1 jtc */
15 1.1 jtc
16 1.1 jtc #ifndef lint
17 1.1 jtc static char rcsid[] = "$Id: w_logf.c,v 1.1 1994/08/10 20:35:23 jtc Exp $";
18 1.1 jtc #endif
19 1.1 jtc
20 1.1 jtc /*
21 1.1 jtc * wrapper logf(x)
22 1.1 jtc */
23 1.1 jtc
24 1.1 jtc #include "math.h"
25 1.1 jtc #include "math_private.h"
26 1.1 jtc
27 1.1 jtc
28 1.1 jtc #ifdef __STDC__
29 1.1 jtc float logf(float x) /* wrapper logf */
30 1.1 jtc #else
31 1.1 jtc float logf(x) /* wrapper logf */
32 1.1 jtc float x;
33 1.1 jtc #endif
34 1.1 jtc {
35 1.1 jtc #ifdef _IEEE_LIBM
36 1.1 jtc return __ieee754_logf(x);
37 1.1 jtc #else
38 1.1 jtc float z;
39 1.1 jtc z = __ieee754_logf(x);
40 1.1 jtc if(_LIB_VERSION == _IEEE_ || isnanf(x) || x > (float)0.0) return z;
41 1.1 jtc if(x==(float)0.0)
42 1.1 jtc /* logf(0) */
43 1.1 jtc return (float)__kernel_standard((double)x,(double)x,116);
44 1.1 jtc else
45 1.1 jtc /* logf(x<0) */
46 1.1 jtc return (float)__kernel_standard((double)x,(double)x,117);
47 1.1 jtc #endif
48 1.1 jtc }
49