HoneyTable::read_item(): Read tag in one go
[xapian.git] / xapian-core / languages / steminternal.cc
blob61510e301df5520f03d336c4b1ef5da6834268c9
1 /** @file steminternal.cc
2 * @brief Base class for implementations of stemming algorithms
3 */
4 /* Derived from snowball's runtime/api.c:
6 * Copyright (c) 2001, Dr Martin Porter
7 * Copyright (c) 2004,2005, Richard Boulton
8 * Copyright (c) 2006,2007,2008,2009,2016 Olly Betts
9 * All rights reserved.
11 * Redistribution and use in source and binary forms, with or without
12 * modification, are permitted provided that the following conditions are met:
14 * * Redistributions of source code must retain the above copyright notice,
15 * this list of conditions and the following disclaimer.
17 * * Redistributions in binary form must reproduce the above copyright
18 * notice, this list of conditions and the following disclaimer in the
19 * documentation and/or other materials provided with the distribution.
21 * * Neither the name of the <ORGANIZATION> nor the names of its contributors
22 * may be used to endorse or promote products derived from this software
23 * without specific prior written permission.
25 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
26 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
28 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
29 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
30 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
31 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
32 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
33 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
34 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
35 * POSSIBILITY OF SUCH DAMAGE.
37 /* Copyright (C) 2007,2010 Olly Betts
38 * Copyright (C) 2010 Evgeny Sizikov
40 * This program is free software; you can redistribute it and/or
41 * modify it under the terms of the GNU General Public License as
42 * published by the Free Software Foundation; either version 2 of the
43 * License, or (at your option) any later version.
45 * This program is distributed in the hope that it will be useful,
46 * but WITHOUT ANY WARRANTY; without even the implied warranty of
47 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
48 * GNU General Public License for more details.
50 * You should have received a copy of the GNU General Public License
51 * along with this program; if not, write to the Free Software
52 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
55 #include <config.h>
57 #include "steminternal.h"
59 #include <xapian/error.h>
61 #include "omassert.h"
63 #include <cstdlib>
64 #include <cstring>
66 #include <string>
68 using namespace std;
70 namespace Xapian {
72 #define CREATE_SIZE 16
74 symbol *
75 SnowballStemImplementation::create_s()
77 void * mem = malloc(HEAD + (CREATE_SIZE + 1) * sizeof(symbol));
78 if (mem == NULL) throw std::bad_alloc();
79 symbol * p = reinterpret_cast<symbol*>(HEAD + static_cast<char *>(mem));
80 SET_CAPACITY(p, CREATE_SIZE);
81 SET_SIZE(p, 0);
82 return p;
86 new_c = skip_utf8(p, c, lb, l, n); skips n characters forwards from p + c
87 if n +ve, or n characters backwards from p + c - 1 if n -ve. new_c is the new
88 value of the cursor c, or -1 on failure.
90 -- used to implement hop and next in the utf8 case.
93 int
94 SnowballStemImplementation::skip_utf8(const symbol * p, int c, int lb, int l, int n)
96 if (n >= 0) {
97 for (; n > 0; --n) {
98 if (c >= l) return -1;
99 if (p[c++] >= 0xC0) { /* 1100 0000 */
100 while (c < l) {
101 /* break unless p[c] is 10------ */
102 if (p[c] >> 6 != 2) break;
103 c++;
107 } else {
108 for (; n < 0; ++n) {
109 if (c <= lb) return -1;
110 if (p[--c] >= 0x80) { /* 1000 0000 */
111 while (c > lb) {
112 if (p[c] >= 0xC0) break; /* 1100 0000 */
113 c--;
118 return c;
122 /* Increase the size of the buffer pointed to by p to at least n symbols.
123 * If insufficient memory, throw std::bad_alloc().
125 symbol *
126 SnowballStemImplementation::increase_size(symbol * p, int n)
128 int new_size = n + 20;
129 void * mem = realloc(reinterpret_cast<char *>(p) - HEAD,
130 HEAD + (new_size + 1) * sizeof(symbol));
131 if (mem == NULL) {
132 throw std::bad_alloc();
134 symbol * q = reinterpret_cast<symbol*>(HEAD + static_cast<char *>(mem));
135 SET_CAPACITY(q, new_size);
136 return q;
140 StemImplementation::~StemImplementation() { }
142 SnowballStemImplementation::~SnowballStemImplementation()
144 lose_s(p);
147 string
148 SnowballStemImplementation::operator()(const string & word)
150 const symbol * s = reinterpret_cast<const symbol *>(word.data());
151 replace_s(0, l, word.size(), s);
152 c = 0;
153 if (stem() < 0) {
154 // FIXME: Is there a better choice of exception class?
155 throw Xapian::InternalError("stemming exception!");
157 return string(reinterpret_cast<const char *>(p), l);
160 /* Code for character groupings: utf8 cases */
162 int SnowballStemImplementation::get_utf8(int * slot) {
163 int b0, b1;
164 int tmp = c;
165 if (tmp >= l) return 0;
166 b0 = p[tmp++];
167 if (b0 < 0xC0 || tmp == l) { /* 1100 0000 */
168 * slot = b0; return 1;
170 b1 = p[tmp++];
171 if (b0 < 0xE0 || tmp == l) { /* 1110 0000 */
172 * slot = (b0 & 0x1F) << 6 | (b1 & 0x3F); return 2;
174 * slot = (b0 & 0xF) << 12 | (b1 & 0x3F) << 6 | (p[tmp] & 0x3F); return 3;
177 int SnowballStemImplementation::get_b_utf8(int * slot) {
178 int b0, b1;
179 int tmp = c;
180 if (tmp <= lb) return 0;
181 b0 = p[--tmp];
182 if (b0 < 0x80 || tmp == lb) { /* 1000 0000 */
183 * slot = b0; return 1;
185 b1 = p[--tmp];
186 if (b1 >= 0xC0 || tmp == lb) { /* 1100 0000 */
187 * slot = (b1 & 0x1F) << 6 | (b0 & 0x3F); return 2;
189 * slot = (p[tmp] & 0xF) << 12 | (b1 & 0x3F) << 6 | (b0 & 0x3F); return 3;
193 SnowballStemImplementation::in_grouping_U(const unsigned char * s, int min,
194 int max, int repeat)
196 do {
197 int ch;
198 int w = get_utf8(&ch);
199 if (!w) return -1;
200 if (ch > max || (ch -= min) < 0 || (s[ch >> 3] & (0X1 << (ch & 0X7))) == 0)
201 return w;
202 c += w;
203 } while (repeat);
204 return 0;
208 SnowballStemImplementation::in_grouping_b_U(const unsigned char * s, int min,
209 int max, int repeat)
211 do {
212 int ch;
213 int w = get_b_utf8(&ch);
214 if (!w) return -1;
215 if (ch > max || (ch -= min) < 0 || (s[ch >> 3] & (0X1 << (ch & 0X7))) == 0)
216 return w;
217 c -= w;
218 } while (repeat);
219 return 0;
223 SnowballStemImplementation::out_grouping_U(const unsigned char * s, int min,
224 int max, int repeat)
226 do {
227 int ch;
228 int w = get_utf8(&ch);
229 if (!w) return -1;
230 if (!(ch > max || (ch -= min) < 0 || (s[ch >> 3] & (0X1 << (ch & 0X7))) == 0))
231 /* FIXME: try adding this so gopast in generated code is simpler: if (repeat == 2) c += w; */ return w;
232 c += w;
233 } while (repeat);
234 return 0;
238 SnowballStemImplementation::out_grouping_b_U(const unsigned char * s, int min,
239 int max, int repeat)
241 do {
242 int ch;
243 int w = get_b_utf8(&ch);
244 if (!w) return -1;
245 if (!(ch > max || (ch -= min) < 0 || (s[ch >> 3] & (0X1 << (ch & 0X7))) == 0))
246 return w;
247 c -= w;
248 } while (repeat);
249 return 0;
252 int SnowballStemImplementation::eq_s(int s_size, const symbol * s) {
253 if (l - c < s_size || memcmp(p + c, s, s_size * sizeof(symbol)) != 0)
254 return 0;
255 c += s_size;
256 return 1;
259 int SnowballStemImplementation::eq_s_b(int s_size, const symbol * s) {
260 if (c - lb < s_size || memcmp(p + c - s_size, s, s_size * sizeof(symbol)) != 0)
261 return 0;
262 c -= s_size;
263 return 1;
267 SnowballStemImplementation::find_among(const symbol * pool,
268 const struct among * v, int v_size,
269 const unsigned char * fnum,
270 const among_function * f)
272 int i = 0;
273 int j = v_size;
275 const symbol * q = p + c;
276 int c_orig = c;
278 int common_i = 0;
279 int common_j = 0;
281 int first_key_inspected = 0;
283 while (1) {
284 int k = i + ((j - i) >> 1);
285 int diff = 0;
286 int common = common_i < common_j ? common_i : common_j; /* smaller */
287 const struct among * w = v + k;
288 for (int x = common; x < w->s_size; ++x) {
289 if (c_orig + common == l) { diff = -1; break; }
290 diff = q[common] - (pool + w->s)[x];
291 if (diff != 0) break;
292 ++common;
294 if (diff < 0) { j = k; common_j = common; }
295 else { i = k; common_i = common; }
296 if (j - i <= 1) {
297 if (i > 0) break; /* v->s has been inspected */
298 if (j == i) break; /* only one item in v */
300 /* - but now we need to go round once more to get
301 v->s inspected. This looks messy, but is actually
302 the optimal approach. */
304 if (first_key_inspected) break;
305 first_key_inspected = 1;
308 while (1) {
309 const struct among * w = v + i;
310 if (common_i >= w->s_size) {
311 c = c_orig + w->s_size;
312 if (!fnum || !fnum[i]) return w->result;
314 int res = f[fnum[i] - 1](this);
315 c = c_orig + w->s_size;
316 if (res) return w->result;
319 i = w->substring_i;
320 if (i < 0) return 0;
324 /* find_among_b is for backwards processing. Same comments apply */
326 SnowballStemImplementation::find_among_b(const symbol * pool,
327 const struct among * v, int v_size,
328 const unsigned char * fnum,
329 const among_function * f)
331 int i = 0;
332 int j = v_size;
334 const symbol * q = p + c - 1;
335 int c_orig = c;
337 int common_i = 0;
338 int common_j = 0;
340 int first_key_inspected = 0;
342 while (1) {
343 int k = i + ((j - i) >> 1);
344 int diff = 0;
345 int common = common_i < common_j ? common_i : common_j;
346 const struct among * w = v + k;
347 for (int x = w->s_size - 1 - common; x >= 0; --x) {
348 if (c_orig - common == lb) { diff = -1; break; }
349 diff = q[- common] - (pool + w->s)[x];
350 if (diff != 0) break;
351 ++common;
353 if (diff < 0) { j = k; common_j = common; }
354 else { i = k; common_i = common; }
355 if (j - i <= 1) {
356 if (i > 0) break;
357 if (j == i) break;
358 if (first_key_inspected) break;
359 first_key_inspected = 1;
362 while (1) {
363 const struct among * w = v + i;
364 if (common_i >= w->s_size) {
365 c = c_orig - w->s_size;
366 if (!fnum || !fnum[i]) return w->result;
368 int res = f[fnum[i] - 1](this);
369 c = c_orig - w->s_size;
370 if (res) return w->result;
373 i = w->substring_i;
374 if (i < 0) return 0;
379 SnowballStemImplementation::replace_s(int c_bra, int c_ket, int s_size,
380 const symbol * s)
382 int adjustment;
383 int len;
384 Assert(p);
385 adjustment = s_size - (c_ket - c_bra);
386 len = SIZE(p);
387 if (adjustment != 0) {
388 if (adjustment + len > CAPACITY(p)) {
389 p = increase_size(p, adjustment + len);
391 memmove(p + c_ket + adjustment,
392 p + c_ket,
393 (len - c_ket) * sizeof(symbol));
394 SET_SIZE(p, adjustment + len);
395 l += adjustment;
396 if (c >= c_ket)
397 c += adjustment;
398 else
399 if (c > c_bra)
400 c = c_bra;
402 if (s_size != 0) memmove(p + c_bra, s, s_size * sizeof(symbol));
403 return adjustment;
406 int SnowballStemImplementation::slice_check() {
407 Assert(p);
408 if (bra < 0 || bra > ket || ket > l) {
409 #if 0
410 fprintf(stderr, "faulty slice operation:\n");
411 debug(z, -1, 0);
412 #endif
413 return -1;
415 return 0;
418 int SnowballStemImplementation::slice_from_s(int s_size, const symbol * s) {
419 if (slice_check()) return -1;
420 replace_s(bra, ket, s_size, s);
421 return 0;
424 void
425 SnowballStemImplementation::insert_s(int c_bra, int c_ket, int s_size,
426 const symbol * s)
428 int adjustment = replace_s(c_bra, c_ket, s_size, s);
429 if (c_bra <= bra) bra += adjustment;
430 if (c_bra <= ket) ket += adjustment;
433 symbol * SnowballStemImplementation::slice_to(symbol * v) {
434 if (slice_check()) return NULL;
436 int len = ket - bra;
437 if (CAPACITY(v) < len) {
438 v = increase_size(v, len);
440 memmove(v, p + bra, len * sizeof(symbol));
441 SET_SIZE(v, len);
443 return v;
446 symbol * SnowballStemImplementation::assign_to(symbol * v) {
447 int len = l;
448 if (CAPACITY(v) < len) {
449 v = increase_size(v, len);
451 memmove(v, p, len * sizeof(symbol));
452 SET_SIZE(v, len);
453 return v;
456 int SnowballStemImplementation::len_utf8(const symbol * v) {
457 int size = SIZE(v);
458 int len = 0;
459 while (size--) {
460 symbol b = *v++;
461 if (b >= 0xC0 || b < 0x80) ++len;
463 return len;
466 #if 0
467 void SnowballStemImplementation::debug(int number, int line_count) {
468 int i;
469 int limit = SIZE(p);
470 if (number >= 0) printf("%3d (line %4d): [%d]'", number, line_count, limit);
471 for (i = 0; i <= limit; ++i) {
472 if (lb == i) printf("{");
473 if (bra == i) printf("[");
474 if (c == i) printf("|");
475 if (ket == i) printf("]");
476 if (l == i) printf("}");
477 if (i < limit) {
478 int ch = p[i];
479 if (ch == 0) ch = '#';
480 printf("%c", ch);
483 printf("'\n");
485 #endif