]>
Commit | Line | Data |
---|---|---|
f7eac6eb RM |
1 | /* @(#)e_scalb.c 5.1 93/09/24 */ |
2 | /* | |
3 | * ==================================================== | |
4 | * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. | |
5 | * | |
6 | * Developed at SunPro, a Sun Microsystems, Inc. business. | |
7 | * Permission to use, copy, modify, and distribute this | |
d210ca02 | 8 | * software is freely granted, provided that this notice |
f7eac6eb RM |
9 | * is preserved. |
10 | * ==================================================== | |
11 | */ | |
12 | ||
13 | #if defined(LIBM_SCCS) && !defined(lint) | |
14 | static char rcsid[] = "$NetBSD: e_scalb.c,v 1.6 1995/05/10 20:46:09 jtc Exp $"; | |
15 | #endif | |
16 | ||
17 | /* | |
18 | * __ieee754_scalb(x, fn) is provide for | |
d210ca02 | 19 | * passing various standard test suite. One |
f7eac6eb RM |
20 | * should use scalbn() instead. |
21 | */ | |
22 | ||
d210ca02 | 23 | #include <fenv.h> |
9d13fb24 | 24 | #include <math.h> |
9277c064 | 25 | #include <math_private.h> |
f7eac6eb RM |
26 | |
27 | #ifdef _SCALB_INT | |
28 | #ifdef __STDC__ | |
29 | double __ieee754_scalb(double x, int fn) | |
30 | #else | |
31 | double __ieee754_scalb(x,fn) | |
32 | double x; int fn; | |
33 | #endif | |
34 | #else | |
35 | #ifdef __STDC__ | |
36 | double __ieee754_scalb(double x, double fn) | |
37 | #else | |
38 | double __ieee754_scalb(x,fn) | |
39 | double x, fn; | |
40 | #endif | |
41 | #endif | |
42 | { | |
43 | #ifdef _SCALB_INT | |
b4012b75 | 44 | return __scalbn(x,fn); |
f7eac6eb | 45 | #else |
b4012b75 UD |
46 | if (__isnan(x)||__isnan(fn)) return x*fn; |
47 | if (!__finite(fn)) { | |
f7eac6eb | 48 | if(fn>0.0) return x*fn; |
478b92f0 UD |
49 | else if (x == 0) |
50 | return x; | |
51 | else if (!__finite (x)) | |
d210ca02 | 52 | { |
bb769ab6 | 53 | # ifdef FE_INVALID |
d210ca02 | 54 | feraiseexcept (FE_INVALID); |
bb769ab6 | 55 | # endif |
d210ca02 UD |
56 | return __nan (""); |
57 | } | |
f7eac6eb RM |
58 | else return x/(-fn); |
59 | } | |
d210ca02 UD |
60 | if (__rint(fn)!=fn) |
61 | { | |
bb769ab6 | 62 | # ifdef FE_INVALID |
d210ca02 | 63 | feraiseexcept (FE_INVALID); |
bb769ab6 | 64 | # endif |
d210ca02 UD |
65 | return __nan (""); |
66 | } | |
b4012b75 UD |
67 | if ( fn > 65000.0) return __scalbn(x, 65000); |
68 | if (-fn > 65000.0) return __scalbn(x,-65000); | |
69 | return __scalbn(x,(int)fn); | |
f7eac6eb RM |
70 | #endif |
71 | } |