xref: /onnv-gate/usr/src/lib/libbc/libc/gen/common/_Qfcompare.c (revision 0:68f95e015346)
1*0Sstevel@tonic-gate /*
2*0Sstevel@tonic-gate  * CDDL HEADER START
3*0Sstevel@tonic-gate  *
4*0Sstevel@tonic-gate  * The contents of this file are subject to the terms of the
5*0Sstevel@tonic-gate  * Common Development and Distribution License, Version 1.0 only
6*0Sstevel@tonic-gate  * (the "License").  You may not use this file except in compliance
7*0Sstevel@tonic-gate  * with the License.
8*0Sstevel@tonic-gate  *
9*0Sstevel@tonic-gate  * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
10*0Sstevel@tonic-gate  * or http://www.opensolaris.org/os/licensing.
11*0Sstevel@tonic-gate  * See the License for the specific language governing permissions
12*0Sstevel@tonic-gate  * and limitations under the License.
13*0Sstevel@tonic-gate  *
14*0Sstevel@tonic-gate  * When distributing Covered Code, include this CDDL HEADER in each
15*0Sstevel@tonic-gate  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
16*0Sstevel@tonic-gate  * If applicable, add the following below this CDDL HEADER, with the
17*0Sstevel@tonic-gate  * fields enclosed by brackets "[]" replaced with your own identifying
18*0Sstevel@tonic-gate  * information: Portions Copyright [yyyy] [name of copyright owner]
19*0Sstevel@tonic-gate  *
20*0Sstevel@tonic-gate  * CDDL HEADER END
21*0Sstevel@tonic-gate  */
22*0Sstevel@tonic-gate #pragma ident	"%Z%%M%	%I%	%E% SMI"
23*0Sstevel@tonic-gate 
24*0Sstevel@tonic-gate /*
25*0Sstevel@tonic-gate  * Copyright (c) 1988 by Sun Microsystems, Inc.
26*0Sstevel@tonic-gate  */
27*0Sstevel@tonic-gate 
28*0Sstevel@tonic-gate #include "_Qquad.h"
29*0Sstevel@tonic-gate #include "_Qglobals.h"
30*0Sstevel@tonic-gate 
31*0Sstevel@tonic-gate enum fcc_type
_fp_compare(px,py,strict)32*0Sstevel@tonic-gate _fp_compare(px, py, strict)
33*0Sstevel@tonic-gate 	unpacked       *px, *py;
34*0Sstevel@tonic-gate 	int             strict;	/* 0 if quiet NaN unexceptional, 1 if
35*0Sstevel@tonic-gate 				 * exceptional */
36*0Sstevel@tonic-gate 
37*0Sstevel@tonic-gate {
38*0Sstevel@tonic-gate 	enum fcc_type   cc;
39*0Sstevel@tonic-gate 	int  k,n;
40*0Sstevel@tonic-gate 
41*0Sstevel@tonic-gate 	if ((px->fpclass == fp_quiet) || (py->fpclass == fp_quiet) ||
42*0Sstevel@tonic-gate 	    (px->fpclass == fp_signaling) || (py->fpclass == fp_signaling)) {
43*0Sstevel@tonic-gate 		if (strict)				/* NaN */
44*0Sstevel@tonic-gate 			fpu_set_exception(fp_invalid);
45*0Sstevel@tonic-gate 		cc = fcc_unordered;
46*0Sstevel@tonic-gate 	} else if ((px->fpclass == fp_zero) && (py->fpclass == fp_zero))
47*0Sstevel@tonic-gate 		cc = fcc_equal;
48*0Sstevel@tonic-gate 	/* both zeros */
49*0Sstevel@tonic-gate 	else if (px->sign < py->sign)
50*0Sstevel@tonic-gate 		cc = fcc_greater;
51*0Sstevel@tonic-gate 	else if (px->sign > py->sign)
52*0Sstevel@tonic-gate 		cc = fcc_less;
53*0Sstevel@tonic-gate 	else {			/* signs the same, compute magnitude cc */
54*0Sstevel@tonic-gate 		if ((int) px->fpclass > (int) py->fpclass)
55*0Sstevel@tonic-gate 			cc = fcc_greater;
56*0Sstevel@tonic-gate 		else if ((int) px->fpclass < (int) py->fpclass)
57*0Sstevel@tonic-gate 			cc = fcc_less;
58*0Sstevel@tonic-gate 		else
59*0Sstevel@tonic-gate 		 /* same classes */ if (px->fpclass == fp_infinity)
60*0Sstevel@tonic-gate 			cc = fcc_equal;	/* same infinity */
61*0Sstevel@tonic-gate 		else if (px->exponent > py->exponent)
62*0Sstevel@tonic-gate 			cc = fcc_greater;
63*0Sstevel@tonic-gate 		else if (px->exponent < py->exponent)
64*0Sstevel@tonic-gate 			cc = fcc_less;
65*0Sstevel@tonic-gate 		else {	/* equal exponents */
66*0Sstevel@tonic-gate 			n = fpu_cmpli(px->significand,py->significand,4);
67*0Sstevel@tonic-gate 			if(n>0) cc = fcc_greater;
68*0Sstevel@tonic-gate 			else if(n<0) cc = fcc_less;
69*0Sstevel@tonic-gate 			else cc = fcc_equal;
70*0Sstevel@tonic-gate 		}
71*0Sstevel@tonic-gate 		if (px->sign)
72*0Sstevel@tonic-gate 			switch (cc) {	/* negative numbers */
73*0Sstevel@tonic-gate 			case fcc_less:
74*0Sstevel@tonic-gate 				cc = fcc_greater;
75*0Sstevel@tonic-gate 				break;
76*0Sstevel@tonic-gate 			case fcc_greater:
77*0Sstevel@tonic-gate 				cc = fcc_less;
78*0Sstevel@tonic-gate 				break;
79*0Sstevel@tonic-gate 			}
80*0Sstevel@tonic-gate 	}
81*0Sstevel@tonic-gate 	return (cc);
82*0Sstevel@tonic-gate }
83