1*05a0b428SJohn Marino /* e_scalbf.c -- float version of e_scalb.c.
2*05a0b428SJohn Marino * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
3*05a0b428SJohn Marino */
4*05a0b428SJohn Marino
5*05a0b428SJohn Marino /*
6*05a0b428SJohn Marino * ====================================================
7*05a0b428SJohn Marino * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8*05a0b428SJohn Marino *
9*05a0b428SJohn Marino * Developed at SunPro, a Sun Microsystems, Inc. business.
10*05a0b428SJohn Marino * Permission to use, copy, modify, and distribute this
11*05a0b428SJohn Marino * software is freely granted, provided that this notice
12*05a0b428SJohn Marino * is preserved.
13*05a0b428SJohn Marino * ====================================================
14*05a0b428SJohn Marino */
15*05a0b428SJohn Marino
16*05a0b428SJohn Marino #include "math.h"
17*05a0b428SJohn Marino #include "math_private.h"
18*05a0b428SJohn Marino
19*05a0b428SJohn Marino #ifdef _SCALB_INT
20*05a0b428SJohn Marino float
scalbf(float x,int fn)21*05a0b428SJohn Marino scalbf(float x, int fn)
22*05a0b428SJohn Marino {
23*05a0b428SJohn Marino return scalbnf(x,fn);
24*05a0b428SJohn Marino }
25*05a0b428SJohn Marino
26*05a0b428SJohn Marino #else
27*05a0b428SJohn Marino
28*05a0b428SJohn Marino float
scalbf(float x,float fn)29*05a0b428SJohn Marino scalbf(float x, float fn)
30*05a0b428SJohn Marino {
31*05a0b428SJohn Marino if (isnanf(x)||isnanf(fn)) return x*fn;
32*05a0b428SJohn Marino if (!finitef(fn)) {
33*05a0b428SJohn Marino if(fn>(float)0.0) return x*fn;
34*05a0b428SJohn Marino else return x/(-fn);
35*05a0b428SJohn Marino }
36*05a0b428SJohn Marino if (rintf(fn)!=fn) return (fn-fn)/(fn-fn);
37*05a0b428SJohn Marino if ( fn > (float)65000.0) return scalbnf(x, 65000);
38*05a0b428SJohn Marino if (-fn > (float)65000.0) return scalbnf(x,-65000);
39*05a0b428SJohn Marino return scalbnf(x,(int)fn);
40*05a0b428SJohn Marino }
41*05a0b428SJohn Marino #endif
42