xref: /minix3/lib/libm/src/e_scalbf.c (revision 2fe8fb192fe7e8720e3e7a77f928da545e872a6a)
1*2fe8fb19SBen Gras /* e_scalbf.c -- float version of e_scalb.c.
2*2fe8fb19SBen Gras  * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
3*2fe8fb19SBen Gras  */
4*2fe8fb19SBen Gras 
5*2fe8fb19SBen Gras /*
6*2fe8fb19SBen Gras  * ====================================================
7*2fe8fb19SBen Gras  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8*2fe8fb19SBen Gras  *
9*2fe8fb19SBen Gras  * Developed at SunPro, a Sun Microsystems, Inc. business.
10*2fe8fb19SBen Gras  * Permission to use, copy, modify, and distribute this
11*2fe8fb19SBen Gras  * software is freely granted, provided that this notice
12*2fe8fb19SBen Gras  * is preserved.
13*2fe8fb19SBen Gras  * ====================================================
14*2fe8fb19SBen Gras  */
15*2fe8fb19SBen Gras 
16*2fe8fb19SBen Gras #include <sys/cdefs.h>
17*2fe8fb19SBen Gras #if defined(LIBM_SCCS) && !defined(lint)
18*2fe8fb19SBen Gras __RCSID("$NetBSD: e_scalbf.c,v 1.7 2010/04/23 19:17:07 drochner Exp $");
19*2fe8fb19SBen Gras #endif
20*2fe8fb19SBen Gras 
21*2fe8fb19SBen Gras #include "namespace.h"
22*2fe8fb19SBen Gras #include "math.h"
23*2fe8fb19SBen Gras #include "math_private.h"
24*2fe8fb19SBen Gras 
25*2fe8fb19SBen Gras #ifdef _SCALB_INT
26*2fe8fb19SBen Gras float
__ieee754_scalbf(float x,int fn)27*2fe8fb19SBen Gras __ieee754_scalbf(float x, int fn)
28*2fe8fb19SBen Gras #else
29*2fe8fb19SBen Gras float
30*2fe8fb19SBen Gras __ieee754_scalbf(float x, float fn)
31*2fe8fb19SBen Gras #endif
32*2fe8fb19SBen Gras {
33*2fe8fb19SBen Gras #ifdef _SCALB_INT
34*2fe8fb19SBen Gras 	return scalbnf(x,fn);
35*2fe8fb19SBen Gras #else
36*2fe8fb19SBen Gras 	if (isnanf(x)||isnanf(fn)) return x*fn;
37*2fe8fb19SBen Gras 	if (!finitef(fn)) {
38*2fe8fb19SBen Gras 	    if(fn>(float)0.0) return x*fn;
39*2fe8fb19SBen Gras 	    else       return x/(-fn);
40*2fe8fb19SBen Gras 	}
41*2fe8fb19SBen Gras 	if (rintf(fn)!=fn) return (fn-fn)/(fn-fn);
42*2fe8fb19SBen Gras 	if ( fn > (float)65000.0) return scalbnf(x, 65000);
43*2fe8fb19SBen Gras 	if (-fn > (float)65000.0) return scalbnf(x,-65000);
44*2fe8fb19SBen Gras 	return scalbnf(x,(int)fn);
45*2fe8fb19SBen Gras #endif
46*2fe8fb19SBen Gras }
47