Line data Source code
1 : #line 2 "../src/kernel/none/cmp.c"
2 : /* Copyright (C) 2002-2003 The PARI group.
3 :
4 : This file is part of the PARI/GP package.
5 :
6 : PARI/GP is free software; you can redistribute it and/or modify it under the
7 : terms of the GNU General Public License as published by the Free Software
8 : Foundation; either version 2 of the License, or (at your option) any later
9 : version. It is distributed in the hope that it will be useful, but WITHOUT
10 : ANY WARRANTY WHATSOEVER.
11 :
12 : Check the License for details. You should have received a copy of it, along
13 : with the package; see the file 'COPYING'. If not, write to the Free Software
14 : Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. */
15 :
16 :
17 : /********************************************************************/
18 : /** **/
19 : /** Comparison routines **/
20 : /** **/
21 : /********************************************************************/
22 :
23 : /*They depend on cmpiispec and equaliispec in mp.c*/
24 :
25 : int
26 426695894 : equalii(GEN x, GEN y)
27 : {
28 426695894 : if ((x[1] & (LGBITS|SIGNBITS)) != (y[1] & (LGBITS|SIGNBITS))) return 0;
29 390550870 : return equaliispec(x+2, y+2, lgefint(x)-2, lgefint(y)-2);
30 : }
31 :
32 : int
33 465002012 : cmpii(GEN x, GEN y)
34 : {
35 465002012 : const long sx = signe(x), sy = signe(y);
36 465002012 : if (sx<sy) return -1;
37 457189387 : if (sx>sy) return 1;
38 445765319 : if (!sx) return 0;
39 436089263 : if (sx>0)
40 417343961 : return cmpiispec(x+2, y+2, lgefint(x)-2, lgefint(y)-2);
41 : else
42 18745302 : return -cmpiispec(x+2, y+2, lgefint(x)-2, lgefint(y)-2);
43 : }
44 :
45 : int
46 229610 : equalrr(GEN x, GEN y)
47 : {
48 : long lx, ly, i;
49 :
50 229610 : if (!signe(x)) {
51 140 : if (!signe(y)) return 1; /* all zeroes are equal */
52 28 : return expo(x) >= expo(y);
53 : }
54 229470 : if (!signe(y))
55 42 : return expo(y) >= expo(x);
56 :
57 229428 : if (x[1] != y[1]) return 0;
58 :
59 226155 : lx = lg(x);
60 226155 : ly = lg(y);
61 226155 : if (lx < ly)
62 : {
63 7 : i=2; while (i<lx && x[i]==y[i]) i++;
64 7 : if (i<lx) return 0;
65 0 : for (; i < ly; i++) if (y[i]) return 0;
66 : }
67 : else
68 : {
69 346020 : i=2; while (i<ly && x[i]==y[i]) i++;
70 226148 : if (i<ly) return 0;
71 1162 : for (; i < lx; i++) if (x[i]) return 0;
72 : }
73 1162 : return 1;
74 : }
75 :
76 : int
77 568407610 : cmprr(GEN x, GEN y)
78 : {
79 568407610 : const long sx = signe(x), sy = signe(y);
80 : long ex,ey,lx,ly,lz,i;
81 :
82 568407610 : if (!sx) {
83 1770252 : if (!sy || expo(x) >= expo(y)) return 0;
84 570220 : return sy > 0? -1: 1;
85 : }
86 566637358 : if (!sy) {
87 846902 : if (expo(y) >= expo(x)) return 0;
88 846752 : return sx > 0? 1: -1;
89 : }
90 565790456 : if (sx<sy) return -1;
91 565695798 : if (sx>sy) return 1;
92 :
93 565673740 : ex=expo(x); ey=expo(y);
94 565673740 : if (ex>ey) return sx;
95 370836600 : if (ex<ey) return -sx;
96 :
97 233069543 : lx=lg(x); ly=lg(y); lz = (lx<ly)?lx:ly;
98 663965561 : i=2; while (i<lz && x[i]==y[i]) i++;
99 233069543 : if (i<lz) return ((ulong)x[i] > (ulong)y[i]) ? sx : -sx;
100 126578249 : if (lx>=ly)
101 : {
102 126584150 : while (i<lx && !x[i]) i++;
103 126584150 : return (i==lx) ? 0 : sx;
104 : }
105 6 : while (i<ly && !y[i]) i++;
106 6 : return (i==ly) ? 0 : -sx;
107 : }
108 :
109 : /* x and y are integers. Return 1 if |x| == |y|, 0 otherwise */
110 : int
111 44727169 : absequalii(GEN x, GEN y)
112 : {
113 44727169 : if (!signe(x)) return !signe(y);
114 44631523 : if (!signe(y)) return 0;
115 44629374 : return equaliispec(x+2, y+2, lgefint(x)-2, lgefint(y)-2);
116 : }
117 :
118 : /* x and y are integers. Return sign(|x| - |y|) */
119 : int
120 1937102504 : abscmpii(GEN x, GEN y)
121 : {
122 1937102504 : if (!signe(x)) return signe(y)? -1: 0;
123 1822016897 : if (!signe(y)) return 1;
124 1701602996 : return cmpiispec(x+2, y+2, lgefint(x)-2, lgefint(y)-2);
125 : }
126 :
127 : /* x and y are reals. Return sign(|x| - |y|) */
128 : int
129 1536213138 : abscmprr(GEN x, GEN y)
130 : {
131 : long ex,ey,lx,ly,lz,i;
132 :
133 1536213138 : if (!signe(x)) return signe(y)? -1: 0;
134 1533457681 : if (!signe(y)) return 1;
135 :
136 1532016542 : ex=expo(x); ey=expo(y);
137 1532016542 : if (ex>ey) return 1;
138 1030032743 : if (ex<ey) return -1;
139 :
140 704541945 : lx=lg(x); ly=lg(y); lz = (lx<ly)?lx:ly;
141 1317600264 : i=2; while (i<lz && x[i]==y[i]) i++;
142 704541945 : if (i<lz) return ((ulong)x[i] > (ulong)y[i])? 1: -1;
143 175281854 : if (lx>=ly)
144 : {
145 175281943 : while (i<lx && !x[i]) i++;
146 175281943 : return (i==lx)? 0: 1;
147 : }
148 0 : while (i<ly && !y[i]) i++;
149 0 : return (i==ly)? 0: -1;
150 : }
151 :
|