Print this page
5262 libm needs to be carefully unifdef'd
5268 libm doesn't need to hide symbols which are already local
| Split |
Close |
| Expand all |
| Collapse all |
--- old/usr/src/lib/libm/common/LD/coshl.c
+++ new/usr/src/lib/libm/common/LD/coshl.c
1 1 /*
2 2 * CDDL HEADER START
3 3 *
4 4 * The contents of this file are subject to the terms of the
5 5 * Common Development and Distribution License (the "License").
6 6 * You may not use this file except in compliance with the License.
7 7 *
8 8 * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
9 9 * or http://www.opensolaris.org/os/licensing.
10 10 * See the License for the specific language governing permissions
11 11 * and limitations under the License.
12 12 *
13 13 * When distributing Covered Code, include this CDDL HEADER in each
14 14 * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
15 15 * If applicable, add the following below this CDDL HEADER, with the
16 16 * fields enclosed by brackets "[]" replaced with your own identifying
17 17 * information: Portions Copyright [yyyy] [name of copyright owner]
18 18 *
19 19 * CDDL HEADER END
|
↓ open down ↓ |
19 lines elided |
↑ open up ↑ |
20 20 */
21 21
22 22 /*
23 23 * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
24 24 */
25 25 /*
26 26 * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
27 27 * Use is subject to license terms.
28 28 */
29 29
30 -#if defined(ELFOBJ)
31 30 #pragma weak coshl = __coshl
32 -#endif
33 31
34 32 #include "libm.h"
35 33 #include "longdouble.h"
36 34
37 35 /*
38 36 * COSH(X)
39 37 * RETURN THE HYPERBOLIC COSINE OF X
40 38 *
41 39 * Method :
42 40 * 1. Replace x by |x| (COSH(x) = COSH(-x)).
43 41 * 2.
44 42 * [ EXP(x) - 1 ]^2
45 43 * 0 <= x <= 0.3465 : COSH(x) := 1 + -------------------
46 44 * 2*EXP(x)
47 45 *
48 46 * EXP(x) + 1/EXP(x)
49 47 * 0.3465 <= x <= thresh : COSH(x) := -------------------
50 48 * 2
51 49 * thresh <= x <= lnovft : COSH(x) := EXP(x)/2
52 50 * lnovft <= x < INF : COSH(x) := SCALBN(EXP(x-MEP1*ln2),ME)
53 51 *
54 52 *
55 53 * here
56 54 * 0.3465 a number that is near one half of ln2.
57 55 * thresh a number such that
58 56 * EXP(thresh)+EXP(-thresh)=EXP(thresh)
59 57 * lnovft logarithm of the overflow threshold
60 58 * = MEP1*ln2 chopped to machine precision.
61 59 * ME maximum exponent
62 60 * MEP1 maximum exponent plus 1
63 61 *
64 62 * Special cases:
65 63 * COSH(x) is |x| if x is +INF, -INF, or NaN.
66 64 * only COSH(0)=1 is exact for finite x.
67 65 */
68 66
69 67 static const long double C[] = {
70 68 0.5L,
71 69 1.0L,
72 70 0.3465L,
73 71 45.0L,
74 72 1.135652340629414394879149e+04L,
75 73 7.004447686242549087858985e-16L,
76 74 2.710505431213761085018632e-20L, /* 2^-65 */
77 75 };
78 76
79 77 #define half C[0]
80 78 #define one C[1]
81 79 #define thr1 C[2]
82 80 #define thr2 C[3]
83 81 #define lnovft C[4]
84 82 #define lnovlo C[5]
85 83 #define tinyl C[6]
86 84
87 85 long double
88 86 coshl(long double x) {
89 87 long double w, t;
90 88
91 89 w = fabsl(x);
92 90 if (!finitel(w))
93 91 return (w + w); /* x is INF or NaN */
94 92 if (w < thr1) {
95 93 if (w < tinyl)
96 94 return (one + w); /* inexact+directed rounding */
97 95 t = expm1l(w);
98 96 w = one + t;
99 97 w = one + (t * t) / (w + w);
100 98 return (w);
101 99 }
102 100 if (w < thr2) {
103 101 t = expl(w);
104 102 return (half * (t + one / t));
105 103 }
106 104 if (w <= lnovft)
107 105 return (half * expl(w));
108 106 return (scalbnl(expl((w - lnovft) - lnovlo), 16383));
109 107 }
|
↓ open down ↓ |
67 lines elided |
↑ open up ↑ |
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX