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