Print this page
11210 libm should be cstyle(1ONBLD) clean
| 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.
|
↓ open down ↓ |
14 lines elided |
↑ open up ↑ |
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
20 20 */
21 21
22 22 /*
23 23 * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
24 24 */
25 +
25 26 /*
26 27 * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
27 28 * Use is subject to license terms.
28 29 */
29 30
30 31 #pragma weak __coshl = coshl
31 32
32 33 #include "libm.h"
33 34 #include "longdouble.h"
34 35
35 36 /*
36 37 * COSH(X)
37 38 * RETURN THE HYPERBOLIC COSINE OF X
38 39 *
39 40 * Method :
40 41 * 1. Replace x by |x| (COSH(x) = COSH(-x)).
41 42 * 2.
42 43 * [ EXP(x) - 1 ]^2
43 44 * 0 <= x <= 0.3465 : COSH(x) := 1 + -------------------
44 45 * 2*EXP(x)
45 46 *
46 47 * EXP(x) + 1/EXP(x)
47 48 * 0.3465 <= x <= thresh : COSH(x) := -------------------
48 49 * 2
49 50 * thresh <= x <= lnovft : COSH(x) := EXP(x)/2
50 51 * lnovft <= x < INF : COSH(x) := SCALBN(EXP(x-MEP1*ln2),ME)
51 52 *
52 53 *
53 54 * here
54 55 * 0.3465 a number that is near one half of ln2.
55 56 * thresh a number such that
56 57 * EXP(thresh)+EXP(-thresh)=EXP(thresh)
57 58 * lnovft logarithm of the overflow threshold
58 59 * = MEP1*ln2 chopped to machine precision.
59 60 * ME maximum exponent
60 61 * MEP1 maximum exponent plus 1
61 62 *
62 63 * Special cases:
63 64 * COSH(x) is |x| if x is +INF, -INF, or NaN.
|
↓ open down ↓ |
29 lines elided |
↑ open up ↑ |
64 65 * only COSH(0)=1 is exact for finite x.
65 66 */
66 67
67 68 static const long double C[] = {
68 69 0.5L,
69 70 1.0L,
70 71 0.3465L,
71 72 45.0L,
72 73 1.135652340629414394879149e+04L,
73 74 7.004447686242549087858985e-16L,
74 - 2.710505431213761085018632e-20L, /* 2^-65 */
75 + 2.710505431213761085018632e-20L, /* 2^-65 */
75 76 };
76 77
77 -#define half C[0]
78 -#define one C[1]
79 -#define thr1 C[2]
80 -#define thr2 C[3]
81 -#define lnovft C[4]
82 -#define lnovlo C[5]
83 -#define tinyl C[6]
78 +#define half C[0]
79 +#define one C[1]
80 +#define thr1 C[2]
81 +#define thr2 C[3]
82 +#define lnovft C[4]
83 +#define lnovlo C[5]
84 +#define tinyl C[6]
84 85
85 86 long double
86 -coshl(long double x) {
87 +coshl(long double x)
88 +{
87 89 long double w, t;
88 90
89 91 w = fabsl(x);
92 +
90 93 if (!finitel(w))
91 - return (w + w); /* x is INF or NaN */
94 + return (w + w); /* x is INF or NaN */
95 +
92 96 if (w < thr1) {
93 97 if (w < tinyl)
94 98 return (one + w); /* inexact+directed rounding */
99 +
95 100 t = expm1l(w);
96 101 w = one + t;
97 102 w = one + (t * t) / (w + w);
98 103 return (w);
99 104 }
105 +
100 106 if (w < thr2) {
101 107 t = expl(w);
102 108 return (half * (t + one / t));
103 109 }
110 +
104 111 if (w <= lnovft)
105 112 return (half * expl(w));
113 +
106 114 return (scalbnl(expl((w - lnovft) - lnovlo), 16383));
107 115 }
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX