xref: /illumos-gate/usr/src/lib/libc/sparc/fp/_Q_cmpe.c (revision 1da57d551424de5a9d469760be7c4b4d4f10a755)
1 /*
2  * CDDL HEADER START
3  *
4  * The contents of this file are subject to the terms of the
5  * Common Development and Distribution License, Version 1.0 only
6  * (the "License").  You may not use this file except in compliance
7  * with the License.
8  *
9  * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
10  * or http://www.opensolaris.org/os/licensing.
11  * See the License for the specific language governing permissions
12  * and limitations under the License.
13  *
14  * When distributing Covered Code, include this CDDL HEADER in each
15  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
16  * If applicable, add the following below this CDDL HEADER, with the
17  * fields enclosed by brackets "[]" replaced with your own identifying
18  * information: Portions Copyright [yyyy] [name of copyright owner]
19  *
20  * CDDL HEADER END
21  */
22 /*
23  * Copyright (c) 1994-1997, by Sun Microsystems, Inc.
24  * All rights reserved.
25  */
26 
27 #include "quad.h"
28 
29 #ifdef __sparcv9
30 #define	_Q_cmpe	_Qp_cmpe
31 #endif
32 
33 /*
34  * _Q_cmpe(x, y) returns the condition code that would result from
35  * fcmpeq *x, *y (and raises the same exceptions).
36  */
37 enum fcc_type
_Q_cmpe(const union longdouble * x,const union longdouble * y)38 _Q_cmpe(const union longdouble *x, const union longdouble *y)
39 {
40 	unsigned int	xm, ym, fsr;
41 
42 	if (QUAD_ISNAN(*x) || QUAD_ISNAN(*y)) {
43 		/* nan, signal invalid */
44 		__quad_getfsrp(&fsr);
45 		if (fsr & FSR_NVM) {
46 			__quad_fcmpeq(x, y, &fsr);
47 			return ((fsr >> 10) & 3);
48 		} else {
49 			fsr = (fsr & ~FSR_CEXC) | FSR_NVA | FSR_NVC;
50 			__quad_setfsrp(&fsr);
51 		}
52 		return (fcc_unordered);
53 	}
54 
55 	/* ignore sign of zero */
56 	xm = x->l.msw;
57 	if (QUAD_ISZERO(*x))
58 		xm &= 0x7fffffff;
59 	ym = y->l.msw;
60 	if (QUAD_ISZERO(*y))
61 		ym &= 0x7fffffff;
62 
63 	if ((xm ^ ym) & 0x80000000)	/* x and y have opposite signs */
64 		return ((ym & 0x80000000)? fcc_greater : fcc_less);
65 
66 	if (xm & 0x80000000) {
67 		if (xm > ym)
68 			return (fcc_less);
69 		if (xm < ym)
70 			return (fcc_greater);
71 		if (x->l.frac2 > y->l.frac2)
72 			return (fcc_less);
73 		if (x->l.frac2 < y->l.frac2)
74 			return (fcc_greater);
75 		if (x->l.frac3 > y->l.frac3)
76 			return (fcc_less);
77 		if (x->l.frac3 < y->l.frac3)
78 			return (fcc_greater);
79 		if (x->l.frac4 > y->l.frac4)
80 			return (fcc_less);
81 		if (x->l.frac4 < y->l.frac4)
82 			return (fcc_greater);
83 		return (fcc_equal);
84 	}
85 	if (xm < ym)
86 		return (fcc_less);
87 	if (xm > ym)
88 		return (fcc_greater);
89 	if (x->l.frac2 < y->l.frac2)
90 		return (fcc_less);
91 	if (x->l.frac2 > y->l.frac2)
92 		return (fcc_greater);
93 	if (x->l.frac3 < y->l.frac3)
94 		return (fcc_less);
95 	if (x->l.frac3 > y->l.frac3)
96 		return (fcc_greater);
97 	if (x->l.frac4 < y->l.frac4)
98 		return (fcc_less);
99 	if (x->l.frac4 > y->l.frac4)
100 		return (fcc_greater);
101 	return (fcc_equal);
102 }
103