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