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 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 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