diff options
author | Linus Torvalds <torvalds@ppc970.osdl.org> | 2005-04-16 18:20:36 -0400 |
---|---|---|
committer | Linus Torvalds <torvalds@ppc970.osdl.org> | 2005-04-16 18:20:36 -0400 |
commit | 1da177e4c3f41524e886b7f1b8a0c1fc7321cac2 (patch) | |
tree | 0bba044c4ce775e45a88a51686b5d9f90697ea9d /arch/parisc/math-emu/sfcmp.c |
Linux-2.6.12-rc2v2.6.12-rc2
Initial git repository build. I'm not bothering with the full history,
even though we have it. We can create a separate "historical" git
archive of that later if we want to, and in the meantime it's about
3.2GB when imported into git - space that would just make the early
git days unnecessarily complicated, when we don't have a lot of good
infrastructure for it.
Let it rip!
Diffstat (limited to 'arch/parisc/math-emu/sfcmp.c')
-rw-r--r-- | arch/parisc/math-emu/sfcmp.c | 155 |
1 files changed, 155 insertions, 0 deletions
diff --git a/arch/parisc/math-emu/sfcmp.c b/arch/parisc/math-emu/sfcmp.c new file mode 100644 index 000000000000..1466fb46e4c0 --- /dev/null +++ b/arch/parisc/math-emu/sfcmp.c | |||
@@ -0,0 +1,155 @@ | |||
1 | /* | ||
2 | * Linux/PA-RISC Project (http://www.parisc-linux.org/) | ||
3 | * | ||
4 | * Floating-point emulation code | ||
5 | * Copyright (C) 2001 Hewlett-Packard (Paul Bame) <bame@debian.org> | ||
6 | * | ||
7 | * This program is free software; you can redistribute it and/or modify | ||
8 | * it under the terms of the GNU General Public License as published by | ||
9 | * the Free Software Foundation; either version 2, or (at your option) | ||
10 | * any later version. | ||
11 | * | ||
12 | * This program is distributed in the hope that it will be useful, | ||
13 | * but WITHOUT ANY WARRANTY; without even the implied warranty of | ||
14 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | ||
15 | * GNU General Public License for more details. | ||
16 | * | ||
17 | * You should have received a copy of the GNU General Public License | ||
18 | * along with this program; if not, write to the Free Software | ||
19 | * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA | ||
20 | */ | ||
21 | /* | ||
22 | * BEGIN_DESC | ||
23 | * | ||
24 | * File: | ||
25 | * @(#) pa/spmath/sfcmp.c $Revision: 1.1 $ | ||
26 | * | ||
27 | * Purpose: | ||
28 | * sgl_cmp: compare two values | ||
29 | * | ||
30 | * External Interfaces: | ||
31 | * sgl_fcmp(leftptr, rightptr, cond, status) | ||
32 | * | ||
33 | * Internal Interfaces: | ||
34 | * | ||
35 | * Theory: | ||
36 | * <<please update with a overview of the operation of this file>> | ||
37 | * | ||
38 | * END_DESC | ||
39 | */ | ||
40 | |||
41 | |||
42 | #include "float.h" | ||
43 | #include "sgl_float.h" | ||
44 | |||
45 | /* | ||
46 | * sgl_cmp: compare two values | ||
47 | */ | ||
48 | int | ||
49 | sgl_fcmp (sgl_floating_point * leftptr, sgl_floating_point * rightptr, | ||
50 | unsigned int cond, unsigned int *status) | ||
51 | |||
52 | /* The predicate to be tested */ | ||
53 | |||
54 | { | ||
55 | register unsigned int left, right; | ||
56 | register int xorresult; | ||
57 | |||
58 | /* Create local copies of the numbers */ | ||
59 | left = *leftptr; | ||
60 | right = *rightptr; | ||
61 | |||
62 | /* | ||
63 | * Test for NaN | ||
64 | */ | ||
65 | if( (Sgl_exponent(left) == SGL_INFINITY_EXPONENT) | ||
66 | || (Sgl_exponent(right) == SGL_INFINITY_EXPONENT) ) | ||
67 | { | ||
68 | /* Check if a NaN is involved. Signal an invalid exception when | ||
69 | * comparing a signaling NaN or when comparing quiet NaNs and the | ||
70 | * low bit of the condition is set */ | ||
71 | if( ( (Sgl_exponent(left) == SGL_INFINITY_EXPONENT) | ||
72 | && Sgl_isnotzero_mantissa(left) | ||
73 | && (Exception(cond) || Sgl_isone_signaling(left))) | ||
74 | || | ||
75 | ( (Sgl_exponent(right) == SGL_INFINITY_EXPONENT) | ||
76 | && Sgl_isnotzero_mantissa(right) | ||
77 | && (Exception(cond) || Sgl_isone_signaling(right)) ) ) | ||
78 | { | ||
79 | if( Is_invalidtrap_enabled() ) { | ||
80 | Set_status_cbit(Unordered(cond)); | ||
81 | return(INVALIDEXCEPTION); | ||
82 | } | ||
83 | else Set_invalidflag(); | ||
84 | Set_status_cbit(Unordered(cond)); | ||
85 | return(NOEXCEPTION); | ||
86 | } | ||
87 | /* All the exceptional conditions are handled, now special case | ||
88 | NaN compares */ | ||
89 | else if( ((Sgl_exponent(left) == SGL_INFINITY_EXPONENT) | ||
90 | && Sgl_isnotzero_mantissa(left)) | ||
91 | || | ||
92 | ((Sgl_exponent(right) == SGL_INFINITY_EXPONENT) | ||
93 | && Sgl_isnotzero_mantissa(right)) ) | ||
94 | { | ||
95 | /* NaNs always compare unordered. */ | ||
96 | Set_status_cbit(Unordered(cond)); | ||
97 | return(NOEXCEPTION); | ||
98 | } | ||
99 | /* infinities will drop down to the normal compare mechanisms */ | ||
100 | } | ||
101 | /* First compare for unequal signs => less or greater or | ||
102 | * special equal case */ | ||
103 | Sgl_xortointp1(left,right,xorresult); | ||
104 | if( xorresult < 0 ) | ||
105 | { | ||
106 | /* left negative => less, left positive => greater. | ||
107 | * equal is possible if both operands are zeros. */ | ||
108 | if( Sgl_iszero_exponentmantissa(left) | ||
109 | && Sgl_iszero_exponentmantissa(right) ) | ||
110 | { | ||
111 | Set_status_cbit(Equal(cond)); | ||
112 | } | ||
113 | else if( Sgl_isone_sign(left) ) | ||
114 | { | ||
115 | Set_status_cbit(Lessthan(cond)); | ||
116 | } | ||
117 | else | ||
118 | { | ||
119 | Set_status_cbit(Greaterthan(cond)); | ||
120 | } | ||
121 | } | ||
122 | /* Signs are the same. Treat negative numbers separately | ||
123 | * from the positives because of the reversed sense. */ | ||
124 | else if( Sgl_all(left) == Sgl_all(right) ) | ||
125 | { | ||
126 | Set_status_cbit(Equal(cond)); | ||
127 | } | ||
128 | else if( Sgl_iszero_sign(left) ) | ||
129 | { | ||
130 | /* Positive compare */ | ||
131 | if( Sgl_all(left) < Sgl_all(right) ) | ||
132 | { | ||
133 | Set_status_cbit(Lessthan(cond)); | ||
134 | } | ||
135 | else | ||
136 | { | ||
137 | Set_status_cbit(Greaterthan(cond)); | ||
138 | } | ||
139 | } | ||
140 | else | ||
141 | { | ||
142 | /* Negative compare. Signed or unsigned compares | ||
143 | * both work the same. That distinction is only | ||
144 | * important when the sign bits differ. */ | ||
145 | if( Sgl_all(left) > Sgl_all(right) ) | ||
146 | { | ||
147 | Set_status_cbit(Lessthan(cond)); | ||
148 | } | ||
149 | else | ||
150 | { | ||
151 | Set_status_cbit(Greaterthan(cond)); | ||
152 | } | ||
153 | } | ||
154 | return(NOEXCEPTION); | ||
155 | } | ||