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 419520119 : equalii(GEN x, GEN y)
27 : {
28 419520119 : if ((x[1] & (LGBITS|SIGNBITS)) != (y[1] & (LGBITS|SIGNBITS))) return 0;
29 383526521 : return equaliispec(x+2, y+2, lgefint(x)-2, lgefint(y)-2);
30 : }
31 :
32 : int
33 434256085 : cmpii(GEN x, GEN y)
34 : {
35 434256085 : const long sx = signe(x), sy = signe(y);
36 434256085 : if (sx<sy) return -1;
37 427289658 : if (sx>sy) return 1;
38 416637284 : if (!sx) return 0;
39 409236029 : if (sx>0)
40 392566268 : return cmpiispec(x+2, y+2, lgefint(x)-2, lgefint(y)-2);
41 : else
42 16669761 : return -cmpiispec(x+2, y+2, lgefint(x)-2, lgefint(y)-2);
43 : }
44 :
45 : int
46 229526 : equalrr(GEN x, GEN y)
47 : {
48 : long lx, ly, i;
49 :
50 229526 : if (!signe(x)) {
51 140 : if (!signe(y)) return 1; /* all zeroes are equal */
52 28 : return expo(x) >= expo(y);
53 : }
54 229386 : if (!signe(y))
55 42 : return expo(y) >= expo(x);
56 :
57 229344 : if (x[1] != y[1]) return 0;
58 :
59 226071 : lx = lg(x);
60 226071 : ly = lg(y);
61 226071 : 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 345844 : i=2; while (i<ly && x[i]==y[i]) i++;
70 226064 : 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 664872892 : cmprr(GEN x, GEN y)
78 : {
79 664872892 : const long sx = signe(x), sy = signe(y);
80 : long ex,ey,lx,ly,lz,i;
81 :
82 664872892 : if (!sx) {
83 1768797 : if (!sy || expo(x) >= expo(y)) return 0;
84 568742 : return sy > 0? -1: 1;
85 : }
86 663104095 : if (!sy) {
87 843693 : if (expo(y) >= expo(x)) return 0;
88 843510 : return sx > 0? 1: -1;
89 : }
90 662260402 : if (sx<sy) return -1;
91 662174807 : if (sx>sy) return 1;
92 :
93 662152890 : ex=expo(x); ey=expo(y);
94 662152890 : if (ex>ey) return sx;
95 414883337 : if (ex<ey) return -sx;
96 :
97 249298955 : lx=lg(x); ly=lg(y); lz = (lx<ly)?lx:ly;
98 722288357 : i=2; while (i<lz && x[i]==y[i]) i++;
99 249298955 : if (i<lz) return ((ulong)x[i] > (ulong)y[i]) ? sx : -sx;
100 141836500 : if (lx>=ly)
101 : {
102 141842940 : while (i<lx && !x[i]) i++;
103 141842940 : 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 43620519 : absequalii(GEN x, GEN y)
112 : {
113 43620519 : if (!signe(x)) return !signe(y);
114 43525229 : if (!signe(y)) return 0;
115 43523080 : 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 1783668102 : abscmpii(GEN x, GEN y)
121 : {
122 1783668102 : if (!signe(x)) return signe(y)? -1: 0;
123 1673094429 : if (!signe(y)) return 1;
124 1554826830 : 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 1869227218 : abscmprr(GEN x, GEN y)
130 : {
131 : long ex,ey,lx,ly,lz,i;
132 :
133 1869227218 : if (!signe(x)) return signe(y)? -1: 0;
134 1866476165 : if (!signe(y)) return 1;
135 :
136 1865034637 : ex=expo(x); ey=expo(y);
137 1865034637 : if (ex>ey) return 1;
138 1238525420 : if (ex<ey) return -1;
139 :
140 835130787 : lx=lg(x); ly=lg(y); lz = (lx<ly)?lx:ly;
141 1516843913 : i=2; while (i<lz && x[i]==y[i]) i++;
142 835130787 : if (i<lz) return ((ulong)x[i] > (ulong)y[i])? 1: -1;
143 199144269 : if (lx>=ly)
144 : {
145 199144341 : while (i<lx && !x[i]) i++;
146 199144341 : return (i==lx)? 0: 1;
147 : }
148 0 : while (i<ly && !y[i]) i++;
149 0 : return (i==ly)? 0: -1;
150 : }
151 :
|
