Print this page
11210 libm should be cstyle(1ONBLD) clean
Split |
Close |
Expand all |
Collapse all |
--- old/usr/src/lib/libm/common/m9x/fminf.c
+++ new/usr/src/lib/libm/common/m9x/fminf.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 fminf = __fminf
31 32
32 -#include "libm.h" /* for islessequal macro */
33 +#include "libm.h" /* for islessequal macro */
33 34
34 35 float
35 -__fminf(float x, float y) {
36 +__fminf(float x, float y)
37 +{
36 38 /*
37 39 * On SPARC v8plus/v9, this could be implemented as follows
38 40 * (assuming %f0 = x, %f1 = y, return value left in %f0):
39 41 *
40 42 * fcmps %fcc0,%f1,%f1
41 43 * fmovsu %fcc0,%f0,%f1
42 44 * fcmps %fcc0,%f0,%f1
43 45 * fmovsug %fcc0,%f1,%f0
44 46 * st %f0,[x]
45 47 * st %f1,[y]
46 48 * ld [x],%l0
47 49 * ld [y],%l1
48 50 * or %l0,%l1,%l2
49 51 * sethi %hi(0x80000000),%l3
50 52 * and %l3,%l2,%l2
51 53 * or %l0,%l2,%l0
52 54 * st %l0,[x]
53 55 * ld [x],%f0
54 56 *
55 57 * If VIS instructions are available, use this code instead:
56 58 *
57 59 * fcmps %fcc0,%f1,%f1
58 60 * fmovsu %fcc0,%f0,%f1
59 61 * fcmps %fcc0,%f0,%f1
60 62 * fmovsug %fcc0,%f1,%f0
61 63 * fors %f0,%f1,%f2
62 64 * fzeros %f3
63 65 * fnegs %f3,%f3
64 66 * fands %f3,%f2,%f2
65 67 * fors %f0,%f2,%f0
66 68 *
↓ open down ↓ |
21 lines elided |
↑ open up ↑ |
67 69 * If VIS 3.0 instructions are available, use this:
68 70 *
69 71 * flcmps %fcc0,%f0,%f1
70 72 * fmovsge %fcc0,%f1,%f0 ! move if %fcc0 is 0 or 2
71 73 */
72 74
73 75 union {
74 76 unsigned i;
75 77 float f;
76 78 } xx, yy;
79 +
77 80 unsigned s;
78 81
79 82 /* if y is nan, replace it by x */
80 83 if (y != y)
81 84 y = x;
82 85
83 86 /* if x is nan, replace it by y */
84 87 if (x != x)
85 88 x = y;
86 89
87 90 /* At this point, x and y are either both numeric, or both NaN */
88 91 if (!isnan(x) && !islessequal(x, y))
89 92 x = y;
90 93
91 94 /*
92 95 * set the sign of the result if either x or y has its sign set
93 96 */
94 97 xx.f = x;
95 98 yy.f = y;
96 99 s = (xx.i | yy.i) & 0x80000000;
97 100 xx.i |= s;
98 101
99 102 return (xx.f);
100 103 }
↓ open down ↓ |
14 lines elided |
↑ open up ↑ |
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX