Line data Source code
1 : /*
2 : * contrib/intarray/_int_tool.c
3 : */
4 : #include "postgres.h"
5 :
6 : #include <limits.h>
7 :
8 : #include "_int.h"
9 : #include "catalog/pg_type.h"
10 : #include "common/int.h"
11 : #include "lib/qunique.h"
12 :
13 : /* arguments are assumed sorted & unique-ified */
14 : bool
15 188186 : inner_int_contains(ArrayType *a, ArrayType *b)
16 : {
17 : int na,
18 : nb;
19 : int i,
20 : j,
21 : n;
22 : int *da,
23 : *db;
24 :
25 188186 : na = ARRNELEMS(a);
26 188186 : nb = ARRNELEMS(b);
27 188186 : da = ARRPTR(a);
28 188186 : db = ARRPTR(b);
29 :
30 188186 : i = j = n = 0;
31 460562 : while (i < na && j < nb)
32 : {
33 447086 : if (da[i] < db[j])
34 259468 : i++;
35 187618 : else if (da[i] == db[j])
36 : {
37 12908 : n++;
38 12908 : i++;
39 12908 : j++;
40 : }
41 : else
42 174710 : break; /* db[j] is not in da */
43 : }
44 :
45 188186 : return (n == nb);
46 : }
47 :
48 : /* arguments are assumed sorted */
49 : bool
50 48834 : inner_int_overlap(ArrayType *a, ArrayType *b)
51 : {
52 : int na,
53 : nb;
54 : int i,
55 : j;
56 : int *da,
57 : *db;
58 :
59 48834 : na = ARRNELEMS(a);
60 48834 : nb = ARRNELEMS(b);
61 48834 : da = ARRPTR(a);
62 48834 : db = ARRPTR(b);
63 :
64 48834 : i = j = 0;
65 228300 : while (i < na && j < nb)
66 : {
67 185276 : if (da[i] < db[j])
68 93378 : i++;
69 91898 : else if (da[i] == db[j])
70 5810 : return true;
71 : else
72 86088 : j++;
73 : }
74 :
75 43024 : return false;
76 : }
77 :
78 : ArrayType *
79 4303696 : inner_int_union(ArrayType *a, ArrayType *b)
80 : {
81 4303696 : ArrayType *r = NULL;
82 :
83 4303696 : CHECKARRVALID(a);
84 4303696 : CHECKARRVALID(b);
85 :
86 4303696 : if (ARRISEMPTY(a) && ARRISEMPTY(b))
87 376 : return new_intArrayType(0);
88 4303320 : if (ARRISEMPTY(a))
89 27280 : r = copy_intArrayType(b);
90 4303320 : if (ARRISEMPTY(b))
91 8694 : r = copy_intArrayType(a);
92 :
93 4303320 : if (!r)
94 : {
95 4267346 : int na = ARRNELEMS(a),
96 4267346 : nb = ARRNELEMS(b);
97 4267346 : int *da = ARRPTR(a),
98 4267346 : *db = ARRPTR(b);
99 : int i,
100 : j,
101 : *dr;
102 :
103 4267346 : r = new_intArrayType(na + nb);
104 4267346 : dr = ARRPTR(r);
105 :
106 : /* union */
107 4267346 : i = j = 0;
108 162302468 : while (i < na && j < nb)
109 : {
110 158035122 : if (da[i] == db[j])
111 : {
112 10116628 : *dr++ = da[i++];
113 10116628 : j++;
114 : }
115 147918494 : else if (da[i] < db[j])
116 125214588 : *dr++ = da[i++];
117 : else
118 22703906 : *dr++ = db[j++];
119 : }
120 :
121 15152916 : while (i < na)
122 10885570 : *dr++ = da[i++];
123 8069238 : while (j < nb)
124 3801892 : *dr++ = db[j++];
125 :
126 4267346 : r = resize_intArrayType(r, dr - ARRPTR(r));
127 : }
128 :
129 4303320 : if (ARRNELEMS(r) > 1)
130 4303312 : r = _int_unique(r);
131 :
132 4303320 : return r;
133 : }
134 :
135 : ArrayType *
136 3431794 : inner_int_inter(ArrayType *a, ArrayType *b)
137 : {
138 : ArrayType *r;
139 : int na,
140 : nb;
141 : int *da,
142 : *db,
143 : *dr;
144 : int i,
145 : j,
146 : k;
147 :
148 3431794 : if (ARRISEMPTY(a) || ARRISEMPTY(b))
149 35136 : return new_intArrayType(0);
150 :
151 3396658 : na = ARRNELEMS(a);
152 3396658 : nb = ARRNELEMS(b);
153 3396658 : da = ARRPTR(a);
154 3396658 : db = ARRPTR(b);
155 3396658 : r = new_intArrayType(Min(na, nb));
156 3396658 : dr = ARRPTR(r);
157 :
158 3396658 : i = j = k = 0;
159 42477906 : while (i < na && j < nb)
160 : {
161 39081248 : if (da[i] < db[j])
162 17979048 : i++;
163 21102200 : else if (da[i] == db[j])
164 : {
165 3488024 : if (k == 0 || dr[k - 1] != db[j])
166 3488024 : dr[k++] = db[j];
167 3488024 : i++;
168 3488024 : j++;
169 : }
170 : else
171 17614176 : j++;
172 : }
173 :
174 3396658 : if (k == 0)
175 : {
176 2690726 : pfree(r);
177 2690726 : return new_intArrayType(0);
178 : }
179 : else
180 705932 : return resize_intArrayType(r, k);
181 : }
182 :
183 : void
184 8356136 : rt__int_size(ArrayType *a, float *size)
185 : {
186 8356136 : *size = (float) ARRNELEMS(a);
187 8356136 : }
188 :
189 : /* comparison function for isort() and _int_unique() */
190 : static inline int
191 677732574 : isort_cmp(const void *a, const void *b, void *arg)
192 : {
193 677732574 : int32 aval = *((const int32 *) a);
194 677732574 : int32 bval = *((const int32 *) b);
195 :
196 677732574 : if (*((bool *) arg))
197 : {
198 : /* compare for ascending order */
199 677732562 : if (aval < bval)
200 294269848 : return -1;
201 383462714 : if (aval > bval)
202 381250380 : return 1;
203 : }
204 : else
205 : {
206 12 : if (aval > bval)
207 8 : return -1;
208 4 : if (aval < bval)
209 4 : return 1;
210 : }
211 2212334 : return 0;
212 : }
213 :
214 : #define ST_SORT isort
215 : #define ST_ELEMENT_TYPE int32
216 : #define ST_COMPARE(a, b, ascending) isort_cmp(a, b, ascending)
217 : #define ST_COMPARE_ARG_TYPE void
218 : #define ST_SCOPE
219 : #define ST_DEFINE
220 : #include "lib/sort_template.h"
221 :
222 : /* Create a new int array with room for "num" elements */
223 : ArrayType *
224 10595618 : new_intArrayType(int num)
225 : {
226 : ArrayType *r;
227 : int nbytes;
228 :
229 : /* if no elements, return a zero-dimensional array */
230 10595618 : if (num <= 0)
231 : {
232 : Assert(num == 0);
233 2726238 : r = construct_empty_array(INT4OID);
234 2726238 : return r;
235 : }
236 :
237 7869380 : nbytes = ARR_OVERHEAD_NONULLS(1) + sizeof(int) * num;
238 :
239 7869380 : r = (ArrayType *) palloc0(nbytes);
240 :
241 7869380 : SET_VARSIZE(r, nbytes);
242 7869380 : ARR_NDIM(r) = 1;
243 7869380 : r->dataoffset = 0; /* marker for no null bitmap */
244 7869380 : ARR_ELEMTYPE(r) = INT4OID;
245 7869380 : ARR_DIMS(r)[0] = num;
246 7869380 : ARR_LBOUND(r)[0] = 1;
247 :
248 7869380 : return r;
249 : }
250 :
251 : ArrayType *
252 9966262 : resize_intArrayType(ArrayType *a, int num)
253 : {
254 : int nbytes;
255 : int i;
256 :
257 : /* if no elements, return a zero-dimensional array */
258 9966262 : if (num <= 0)
259 : {
260 : Assert(num == 0);
261 468 : a = construct_empty_array(INT4OID);
262 468 : return a;
263 : }
264 :
265 9965794 : if (num == ARRNELEMS(a))
266 7715538 : return a;
267 :
268 2250256 : nbytes = ARR_DATA_OFFSET(a) + sizeof(int) * num;
269 :
270 2250256 : a = (ArrayType *) repalloc(a, nbytes);
271 :
272 2250256 : SET_VARSIZE(a, nbytes);
273 : /* usually the array should be 1-D already, but just in case ... */
274 4500512 : for (i = 0; i < ARR_NDIM(a); i++)
275 : {
276 2250256 : ARR_DIMS(a)[i] = num;
277 2250256 : num = 1;
278 : }
279 2250256 : return a;
280 : }
281 :
282 : ArrayType *
283 38386 : copy_intArrayType(ArrayType *a)
284 : {
285 : ArrayType *r;
286 38386 : int n = ARRNELEMS(a);
287 :
288 38386 : r = new_intArrayType(n);
289 38386 : memcpy(ARRPTR(r), ARRPTR(a), n * sizeof(int32));
290 38386 : return r;
291 : }
292 :
293 : /* num for compressed key */
294 : int
295 55160 : internal_size(int *a, int len)
296 : {
297 : int i;
298 55160 : int64 size = 0;
299 :
300 13198216 : for (i = 0; i < len; i += 2)
301 : {
302 13143056 : if (!i || a[i] != a[i - 1]) /* do not count repeated range */
303 13143056 : size += (int64) (a[i + 1]) - (int64) (a[i]) + 1;
304 : }
305 :
306 55160 : if (size > (int64) INT_MAX || size < (int64) INT_MIN)
307 0 : return -1; /* overflow */
308 55160 : return (int) size;
309 : }
310 :
311 : /* unique-ify elements of r in-place ... r must be sorted already */
312 : ArrayType *
313 4987340 : _int_unique(ArrayType *r)
314 : {
315 4987340 : int num = ARRNELEMS(r);
316 4987340 : bool ascending = true;
317 :
318 4987340 : num = qunique_arg(ARRPTR(r), num, sizeof(int), isort_cmp,
319 : &ascending);
320 :
321 4987340 : return resize_intArrayType(r, num);
322 : }
323 :
324 : void
325 0 : gensign(BITVECP sign, int *a, int len, int siglen)
326 : {
327 : int i;
328 :
329 : /* we assume that the sign vector is previously zeroed */
330 0 : for (i = 0; i < len; i++)
331 : {
332 0 : HASH(sign, *a, siglen);
333 0 : a++;
334 : }
335 0 : }
336 :
337 : int32
338 2 : intarray_match_first(ArrayType *a, int32 elem)
339 : {
340 : int32 *aa,
341 : c,
342 : i;
343 :
344 2 : CHECKARRVALID(a);
345 2 : c = ARRNELEMS(a);
346 2 : aa = ARRPTR(a);
347 4 : for (i = 0; i < c; i++)
348 4 : if (aa[i] == elem)
349 2 : return (i + 1);
350 0 : return 0;
351 : }
352 :
353 : ArrayType *
354 8 : intarray_add_elem(ArrayType *a, int32 elem)
355 : {
356 : ArrayType *result;
357 : int32 *r;
358 : int32 c;
359 :
360 8 : CHECKARRVALID(a);
361 8 : c = ARRNELEMS(a);
362 8 : result = new_intArrayType(c + 1);
363 8 : r = ARRPTR(result);
364 8 : if (c > 0)
365 8 : memcpy(r, ARRPTR(a), c * sizeof(int32));
366 8 : r[c] = elem;
367 8 : return result;
368 : }
369 :
370 : ArrayType *
371 2 : intarray_concat_arrays(ArrayType *a, ArrayType *b)
372 : {
373 : ArrayType *result;
374 2 : int32 ac = ARRNELEMS(a);
375 2 : int32 bc = ARRNELEMS(b);
376 :
377 2 : CHECKARRVALID(a);
378 2 : CHECKARRVALID(b);
379 2 : result = new_intArrayType(ac + bc);
380 2 : if (ac)
381 2 : memcpy(ARRPTR(result), ARRPTR(a), ac * sizeof(int32));
382 2 : if (bc)
383 2 : memcpy(ARRPTR(result) + ac, ARRPTR(b), bc * sizeof(int32));
384 2 : return result;
385 : }
386 :
387 : ArrayType *
388 2 : int_to_intset(int32 elem)
389 : {
390 : ArrayType *result;
391 : int32 *aa;
392 :
393 2 : result = new_intArrayType(1);
394 2 : aa = ARRPTR(result);
395 2 : aa[0] = elem;
396 2 : return result;
397 : }
|
