1 /** @file steminternal.cc
2 * @brief Base class for implementations of stemming algorithms
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
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
57 #include "steminternal.h"
59 #include <xapian/error.h>
72 #define CREATE_SIZE 16
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
);
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.
94 SnowballStemImplementation::skip_utf8(const symbol
* p
, int c
, int lb
, int l
, int n
)
98 if (c
>= l
) return -1;
99 if (p
[c
++] >= 0xC0) { /* 1100 0000 */
101 /* break unless p[c] is 10------ */
102 if (p
[c
] >> 6 != 2) break;
109 if (c
<= lb
) return -1;
110 if (p
[--c
] >= 0x80) { /* 1000 0000 */
112 if (p
[c
] >= 0xC0) break; /* 1100 0000 */
122 /* Increase the size of the buffer pointed to by p to at least n symbols.
123 * If insufficient memory, throw std::bad_alloc().
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
));
132 throw std::bad_alloc();
134 symbol
* q
= reinterpret_cast<symbol
*>(HEAD
+ static_cast<char *>(mem
));
135 SET_CAPACITY(q
, new_size
);
140 StemImplementation::~StemImplementation() { }
142 SnowballStemImplementation::~SnowballStemImplementation()
148 SnowballStemImplementation::operator()(const string
& word
)
150 const symbol
* s
= reinterpret_cast<const symbol
*>(word
.data());
151 replace_s(0, l
, word
.size(), s
);
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
) {
165 if (tmp
>= l
) return 0;
167 if (b0
< 0xC0 || tmp
== l
) { /* 1100 0000 */
168 * slot
= b0
; return 1;
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
) {
180 if (tmp
<= lb
) return 0;
182 if (b0
< 0x80 || tmp
== lb
) { /* 1000 0000 */
183 * slot
= b0
; return 1;
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
,
198 int w
= get_utf8(&ch
);
200 if (ch
> max
|| (ch
-= min
) < 0 || (s
[ch
>> 3] & (0X1 << (ch
& 0X7))) == 0)
208 SnowballStemImplementation::in_grouping_b_U(const unsigned char * s
, int min
,
213 int w
= get_b_utf8(&ch
);
215 if (ch
> max
|| (ch
-= min
) < 0 || (s
[ch
>> 3] & (0X1 << (ch
& 0X7))) == 0)
223 SnowballStemImplementation::out_grouping_U(const unsigned char * s
, int min
,
228 int w
= get_utf8(&ch
);
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
;
238 SnowballStemImplementation::out_grouping_b_U(const unsigned char * s
, int min
,
243 int w
= get_b_utf8(&ch
);
245 if (!(ch
> max
|| (ch
-= min
) < 0 || (s
[ch
>> 3] & (0X1 << (ch
& 0X7))) == 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)
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)
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
)
275 const symbol
* q
= p
+ c
;
281 int first_key_inspected
= 0;
284 int k
= i
+ ((j
- i
) >> 1);
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;
294 if (diff
< 0) { j
= k
; common_j
= common
; }
295 else { i
= k
; common_i
= common
; }
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;
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
;
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
)
334 const symbol
* q
= p
+ c
- 1;
340 int first_key_inspected
= 0;
343 int k
= i
+ ((j
- i
) >> 1);
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;
353 if (diff
< 0) { j
= k
; common_j
= common
; }
354 else { i
= k
; common_i
= common
; }
358 if (first_key_inspected
) break;
359 first_key_inspected
= 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
;
379 SnowballStemImplementation::replace_s(int c_bra
, int c_ket
, int s_size
,
385 adjustment
= s_size
- (c_ket
- c_bra
);
387 if (adjustment
!= 0) {
388 if (adjustment
+ len
> CAPACITY(p
)) {
389 p
= increase_size(p
, adjustment
+ len
);
391 memmove(p
+ c_ket
+ adjustment
,
393 (len
- c_ket
) * sizeof(symbol
));
394 SET_SIZE(p
, adjustment
+ len
);
402 if (s_size
!= 0) memmove(p
+ c_bra
, s
, s_size
* sizeof(symbol
));
406 int SnowballStemImplementation::slice_check() {
408 if (bra
< 0 || bra
> ket
|| ket
> l
) {
410 fprintf(stderr
, "faulty slice operation:\n");
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
);
425 SnowballStemImplementation::insert_s(int c_bra
, int c_ket
, int s_size
,
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
;
437 if (CAPACITY(v
) < len
) {
438 v
= increase_size(v
, len
);
440 memmove(v
, p
+ bra
, len
* sizeof(symbol
));
446 symbol
* SnowballStemImplementation::assign_to(symbol
* v
) {
448 if (CAPACITY(v
) < len
) {
449 v
= increase_size(v
, len
);
451 memmove(v
, p
, len
* sizeof(symbol
));
456 int SnowballStemImplementation::len_utf8(const symbol
* v
) {
461 if (b
>= 0xC0 || b
< 0x80) ++len
;
467 void SnowballStemImplementation::debug(int number
, int line_count
) {
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("}");
479 if (ch
== 0) ch
= '#';