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