1*05a0b428SJohn Marino /* @(#)e_scalb.c 5.1 93/09/24 */
2*05a0b428SJohn Marino /*
3*05a0b428SJohn Marino * ====================================================
4*05a0b428SJohn Marino * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
5*05a0b428SJohn Marino *
6*05a0b428SJohn Marino * Developed at SunPro, a Sun Microsystems, Inc. business.
7*05a0b428SJohn Marino * Permission to use, copy, modify, and distribute this
8*05a0b428SJohn Marino * software is freely granted, provided that this notice
9*05a0b428SJohn Marino * is preserved.
10*05a0b428SJohn Marino * ====================================================
11*05a0b428SJohn Marino */
12*05a0b428SJohn Marino
13*05a0b428SJohn Marino /*
14*05a0b428SJohn Marino * scalb(x, fn) is provide for
15*05a0b428SJohn Marino * passing various standard test suite. One
16*05a0b428SJohn Marino * should use scalbn() instead.
17*05a0b428SJohn Marino */
18*05a0b428SJohn Marino
19*05a0b428SJohn Marino #include "math.h"
20*05a0b428SJohn Marino #include "math_private.h"
21*05a0b428SJohn Marino
22*05a0b428SJohn Marino #ifdef _SCALB_INT
23*05a0b428SJohn Marino double
scalb(double x,int fn)24*05a0b428SJohn Marino scalb(double x, int fn)
25*05a0b428SJohn Marino {
26*05a0b428SJohn Marino return scalbn(x, fn);
27*05a0b428SJohn Marino }
28*05a0b428SJohn Marino
29*05a0b428SJohn Marino #else
30*05a0b428SJohn Marino
31*05a0b428SJohn Marino double
scalb(double x,double fn)32*05a0b428SJohn Marino scalb(double x, double fn)
33*05a0b428SJohn Marino {
34*05a0b428SJohn Marino if (isnan(x)||isnan(fn)) return x*fn;
35*05a0b428SJohn Marino if (!finite(fn)) {
36*05a0b428SJohn Marino if(fn>0.0) return x*fn;
37*05a0b428SJohn Marino else return x/(-fn);
38*05a0b428SJohn Marino }
39*05a0b428SJohn Marino if (rint(fn)!=fn) return (fn-fn)/(fn-fn);
40*05a0b428SJohn Marino if ( fn > 65000.0) return scalbn(x, 65000);
41*05a0b428SJohn Marino if (-fn > 65000.0) return scalbn(x,-65000);
42*05a0b428SJohn Marino return scalbn(x,(int)fn);
43*05a0b428SJohn Marino }
44*05a0b428SJohn Marino #endif
45