1 2/* 3 * Copyright (c) 1998, 2001, Oracle and/or its affiliates. All rights reserved. 4 * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER. 5 * 6 * This code is free software; you can redistribute it and/or modify it 7 * under the terms of the GNU General Public License version 2 only, as 8 * published by the Free Software Foundation. Oracle designates this 9 * particular file as subject to the "Classpath" exception as provided 10 * by Oracle in the LICENSE file that accompanied this code. 11 * 12 * This code is distributed in the hope that it will be useful, but WITHOUT 13 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 14 * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 15 * version 2 for more details (a copy is included in the LICENSE file that 16 * accompanied this code). 17 * 18 * You should have received a copy of the GNU General Public License version 19 * 2 along with this work; if not, write to the Free Software Foundation, 20 * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA. 21 * 22 * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA 23 * or visit www.oracle.com if you need additional information or have any 24 * questions. 25 */ 26 27/* 28 * wrapper scalb(double x, double fn) is provide for 29 * passing various standard test suite. One 30 * should use scalbn() instead. 31 */ 32 33#include "fdlibm.h" 34 35#include <errno.h> 36 37#ifdef __STDC__ 38#ifdef _SCALB_INT 39 double scalb(double x, int fn) /* wrapper scalb */ 40#else 41 double scalb(double x, double fn) /* wrapper scalb */ 42#endif 43#else 44 double scalb(x,fn) /* wrapper scalb */ 45#ifdef _SCALB_INT 46 double x; int fn; 47#else 48 double x,fn; 49#endif 50#endif 51{ 52#ifdef _IEEE_LIBM 53 return __ieee754_scalb(x,fn); 54#else 55 double z; 56 z = __ieee754_scalb(x,fn); 57 if(_LIB_VERSION == _IEEE_) return z; 58 if(!(finite(z)||isnan(z))&&finite(x)) { 59 return __kernel_standard(x,(double)fn,32); /* scalb overflow */ 60 } 61 if(z==0.0&&z!=x) { 62 return __kernel_standard(x,(double)fn,33); /* scalb underflow */ 63 } 64#ifndef _SCALB_INT 65 if(!finite(fn)) errno = ERANGE; 66#endif 67 return z; 68#endif 69} 70