1 /* w_log10f.c -- float version of w_log10.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 
16 #include <sys/cdefs.h>
17 #if defined(LIBM_SCCS) && !defined(lint)
18 __RCSID("$NetBSD: w_log2f.c,v 1.2 2024/07/16 14:52:50 riastradh Exp $");
19 #endif
20 
21 /*
22  * wrapper log2f(X)
23  */
24 
25 #include "namespace.h"
26 
27 #include "math.h"
28 #include "math_private.h"
29 
__weak_alias(log2f,_log2f)30 __weak_alias(log2f, _log2f)
31 float
32 log2f(float x)                /* wrapper log2f */
33 {
34 #ifdef _IEEE_LIBM
35           return __ieee754_log2f(x);
36 #else
37           float z;
38           z = __ieee754_log2f(x);
39           if(_LIB_VERSION == _IEEE_ || isnanf(x)) return z;
40           if(x<=(float)0.0) {
41               if(x==(float)0.0)
42                   /* log2(0) */
43                   return (float)__kernel_standard((double)x,(double)x,148);
44               else
45                   /* log2(x<0) */
46                   return (float)__kernel_standard((double)x,(double)x,149);
47           } else
48               return z;
49 #endif
50 }
51