Work towards phrases and multiple registers
[idzebra-moved-to-github.git] / index / zrpn.c
1 /*
2  * Copyright (C) 1994-1995, Index Data I/S 
3  * All rights reserved.
4  * Sebastian Hammer, Adam Dickmeiss
5  *
6  * $Log: zrpn.c,v $
7  * Revision 1.42  1996-05-09 07:28:56  quinn
8  * Work towards phrases and multiple registers
9  *
10  * Revision 1.41  1996/03/20  09:36:43  adam
11  * Function dict_lookup_grep got extra parameter, init_pos, which marks
12  * from which position in pattern approximate pattern matching should occur.
13  * Approximate pattern matching is used in relevance=re-2.
14  *
15  * Revision 1.40  1996/02/02  13:44:44  adam
16  * The public dictionary functions simply use char instead of Dict_char
17  * to represent search strings. Dict_char is used internally only.
18  *
19  * Revision 1.39  1996/01/03  16:22:13  quinn
20  * operator->roperator
21  *
22  * Revision 1.38  1995/12/11  09:12:55  adam
23  * The rec_get function returns NULL if record doesn't exist - will
24  * happen in the server if the result set records have been deleted since
25  * the creation of the set (i.e. the search).
26  * The server saves a result temporarily if it is 'volatile', i.e. the
27  * set is register dependent.
28  *
29  * Revision 1.37  1995/12/06  15:05:28  adam
30  * More verbose in count_set.
31  *
32  * Revision 1.36  1995/12/06  12:41:27  adam
33  * New command 'stat' for the index program.
34  * Filenames can be read from stdin by specifying '-'.
35  * Bug fix/enhancement of the transformation from terms to regular
36  * expressons in the search engine.
37  *
38  * Revision 1.35  1995/11/27  09:29:00  adam
39  * Bug fixes regarding conversion to regular expressions.
40  *
41  * Revision 1.34  1995/11/16  17:00:56  adam
42  * Better logging of rpn query.
43  *
44  * Revision 1.33  1995/11/01  13:58:28  quinn
45  * Moving data1 to yaz/retrieval
46  *
47  * Revision 1.32  1995/10/27  14:00:11  adam
48  * Implemented detection of database availability.
49  *
50  * Revision 1.31  1995/10/17  18:02:10  adam
51  * New feature: databases. Implemented as prefix to words in dictionary.
52  *
53  * Revision 1.30  1995/10/16  09:32:38  adam
54  * More work on relational op.
55  *
56  * Revision 1.29  1995/10/13  16:01:49  adam
57  * Work on relations.
58  *
59  * Revision 1.28  1995/10/13  12:26:43  adam
60  * Optimization of truncation.
61  *
62  * Revision 1.27  1995/10/12  17:07:22  adam
63  * Truncation works.
64  *
65  * Revision 1.26  1995/10/12  12:40:54  adam
66  * Bug fixes in rpn_prox.
67  *
68  * Revision 1.25  1995/10/10  13:59:24  adam
69  * Function rset_open changed its wflag parameter to general flags.
70  *
71  * Revision 1.24  1995/10/09  16:18:37  adam
72  * Function dict_lookup_grep got extra client data parameter.
73  *
74  * Revision 1.23  1995/10/06  16:33:37  adam
75  * Use attribute mappings.
76  *
77  * Revision 1.22  1995/10/06  15:07:39  adam
78  * Structure 'local-number' handled.
79  *
80  * Revision 1.21  1995/10/06  13:52:06  adam
81  * Bug fixes. Handler may abort further scanning.
82  *
83  * Revision 1.20  1995/10/06  11:06:33  adam
84  * Scan entries include 'occurrences' now.
85  *
86  * Revision 1.19  1995/10/06  10:43:56  adam
87  * Scan added. 'occurrences' in scan entries not set yet.
88  *
89  * Revision 1.18  1995/10/04  16:57:20  adam
90  * Key input and merge sort in one pass.
91  *
92  * Revision 1.17  1995/10/04  12:55:17  adam
93  * Bug fix in ranked search. Use=Any keys inserted.
94  *
95  * Revision 1.16  1995/10/02  16:24:40  adam
96  * Use attribute actually used in search requests.
97  *
98  * Revision 1.15  1995/10/02  15:18:52  adam
99  * New member in recRetrieveCtrl: diagnostic.
100  *
101  * Revision 1.14  1995/09/28  12:10:32  adam
102  * Bug fixes. Field prefix used in queries.
103  *
104  * Revision 1.13  1995/09/18  14:17:50  adam
105  * Minor changes.
106  *
107  * Revision 1.12  1995/09/15  14:45:21  adam
108  * Retrieve control.
109  * Work on truncation.
110  *
111  * Revision 1.11  1995/09/14  11:53:27  adam
112  * First work on regular expressions/truncations.
113  *
114  * Revision 1.10  1995/09/11  15:23:26  adam
115  * More work on relevance search.
116  *
117  * Revision 1.9  1995/09/11  13:09:35  adam
118  * More work on relevance feedback.
119  *
120  * Revision 1.8  1995/09/08  14:52:27  adam
121  * Minor changes. Dictionary is lower case now.
122  *
123  * Revision 1.7  1995/09/07  13:58:36  adam
124  * New parameter: result-set file descriptor (RSFD) to support multiple
125  * positions within the same result-set.
126  * Boolean operators: and, or, not implemented.
127  * Result-set references.
128  *
129  * Revision 1.6  1995/09/06  16:11:18  adam
130  * Option: only one word key per file.
131  *
132  * Revision 1.5  1995/09/06  10:33:04  adam
133  * More work on present. Some log messages removed.
134  *
135  * Revision 1.4  1995/09/05  15:28:40  adam
136  * More work on search engine.
137  *
138  * Revision 1.3  1995/09/04  15:20:22  adam
139  * Minor changes.
140  *
141  * Revision 1.2  1995/09/04  12:33:43  adam
142  * Various cleanup. YAZ util used instead.
143  *
144  * Revision 1.1  1995/09/04  09:10:40  adam
145  * More work on index add/del/update.
146  * Merge sort implemented.
147  * Initial work on z39 server.
148  *
149  */
150 #include <stdio.h>
151 #include <assert.h>
152 #include <unistd.h>
153 #include <ctype.h>
154
155 #include "zserver.h"
156 #include "attribute.h"
157
158 #include <rsisam.h>
159 #include <rstemp.h>
160 #include <rsnull.h>
161 #include <rsbool.h>
162 #include <rsrel.h>
163
164 int index_word_prefix_map (char *string, oid_value attrSet, int attrUse,
165                            char *basename)
166 {
167     attent *attp;
168
169     logf (LOG_DEBUG, "oid_value attrSet = %d, attrUse = %d", attrSet, attrUse);
170     attp = att_getentbyatt (attrSet, attrUse);
171     if (!attp)
172         return -1;
173     logf (LOG_DEBUG, "ord=%d", attp->attset_ordinal);
174     return index_word_prefix (string, attp->attset_ordinal,
175                               attp->local_attributes->local, basename);
176 }
177
178 typedef struct {
179     int type;
180     int major;
181     int minor;
182     Z_AttributesPlusTerm *zapt;
183 } AttrType;
184
185 static int attr_find (AttrType *src, oid_value *attributeSetP)
186 {
187     while (src->major < src->zapt->num_attributes)
188     {
189         Z_AttributeElement *element;
190
191         element = src->zapt->attributeList[src->major];
192         if (src->type == *element->attributeType)
193         {
194             switch (element->which) 
195             {
196             case Z_AttributeValue_numeric:
197                 ++(src->major);
198                 if (element->attributeSet && attributeSetP)
199                 {
200                     oident *attrset;
201
202                     attrset = oid_getentbyoid (element->attributeSet);
203                     *attributeSetP = attrset->value;
204                 }
205                 return *element->value.numeric;
206                 break;
207             case Z_AttributeValue_complex:
208                 if (src->minor >= element->value.complex->num_list ||
209                     element->value.complex->list[src->minor]->which !=  
210                     Z_StringOrNumeric_numeric)
211                     break;
212                 ++(src->minor);
213                 if (element->attributeSet && attributeSetP)
214                 {
215                     oident *attrset;
216
217                     attrset = oid_getentbyoid (element->attributeSet);
218                     *attributeSetP = attrset->value;
219                 }
220                 return *element->value.complex->list[src->minor-1]->u.numeric;
221             default:
222                 assert (0);
223             }
224         }
225         ++(src->major);
226     }
227     return -1;
228 }
229
230 static void attr_init (AttrType *src, Z_AttributesPlusTerm *zapt,
231                        int type)
232 {
233     src->zapt = zapt;
234     src->type = type;
235     src->major = 0;
236     src->minor = 0;
237 }
238
239 struct trunc_info {
240     int  *ptr;
241     int  *indx;
242     char **heap;
243     int  heapnum;
244     int  (*cmp)(const void *p1, const void *p2);
245     int  keysize;
246     char *swapbuf;
247     char *tmpbuf;
248     char *buf;
249 };
250
251 static void heap_swap (struct trunc_info *ti, int i1, int i2)
252 {
253     int swap;
254
255     swap = ti->ptr[i1];
256     ti->ptr[i1] = ti->ptr[i2];
257     ti->ptr[i2] = swap;
258 }
259
260 static void heap_delete (struct trunc_info *ti)
261 {
262     int cur = 1, child = 2;
263
264     heap_swap (ti, 1, ti->heapnum--);
265     while (child <= ti->heapnum) {
266         if (child < ti->heapnum &&
267             (*ti->cmp)(ti->heap[ti->ptr[child]],
268                        ti->heap[ti->ptr[1+child]]) > 0)
269             child++;
270         if ((*ti->cmp)(ti->heap[ti->ptr[cur]],
271                        ti->heap[ti->ptr[child]]) > 0)
272         {
273             heap_swap (ti, cur, child);
274             cur = child;
275             child = 2*cur;
276         }
277         else
278             break;
279     }
280 }
281
282 static void heap_insert (struct trunc_info *ti, const char *buf, int indx)
283 {
284     int cur, parent;
285
286     cur = ++(ti->heapnum);
287     memcpy (ti->heap[ti->ptr[cur]], buf, ti->keysize);
288     ti->indx[ti->ptr[cur]] = indx;
289     parent = cur/2;
290     while (parent && (*ti->cmp)(ti->heap[ti->ptr[parent]],
291                                 ti->heap[ti->ptr[cur]]) > 0)
292     {
293         heap_swap (ti, cur, parent);
294         cur = parent;
295         parent = cur/2;
296     }
297 }
298
299 static
300 struct trunc_info *heap_init (int size, int key_size,
301                               int (*cmp)(const void *p1, const void *p2))
302 {
303     struct trunc_info *ti = xmalloc (sizeof(*ti));
304     int i;
305
306     ++size;
307     ti->heapnum = 0;
308     ti->keysize = key_size;
309     ti->cmp = cmp;
310     ti->indx = xmalloc (size * sizeof(*ti->indx));
311     ti->heap = xmalloc (size * sizeof(*ti->heap));
312     ti->ptr = xmalloc (size * sizeof(*ti->ptr));
313     ti->swapbuf = xmalloc (ti->keysize);
314     ti->tmpbuf = xmalloc (ti->keysize);
315     ti->buf = xmalloc (size * ti->keysize);
316     for (i = size; --i >= 0; )
317     {
318         ti->ptr[i] = i;
319         ti->heap[i] = ti->buf + ti->keysize * i;
320     }
321     return ti;
322 }
323
324 static void heap_close (struct trunc_info *ti)
325 {
326     xfree (ti->ptr);
327     xfree (ti->indx);
328     xfree (ti->heap);
329     xfree (ti->swapbuf);
330     xfree (ti->tmpbuf);
331     xfree (ti);
332 }
333
334 static RSET rset_trunc_r (ISAM isam, ISAM_P *isam_p, int from, int to,
335                          int merge_chunk)
336 {
337     RSET result; 
338     RSFD result_rsfd;
339     rset_temp_parms parms;
340
341     parms.key_size = sizeof(struct it_key);
342     result = rset_create (rset_kind_temp, &parms);
343     result_rsfd = rset_open (result, RSETF_WRITE|RSETF_SORT_SYSNO);
344
345     if (to - from > merge_chunk)
346     {
347         RSFD *rsfd;
348         RSET *rset;
349         int i, i_add = (to-from)/merge_chunk + 1;
350         struct trunc_info *ti;
351         int rscur = 0;
352         int rsmax = (to-from)/i_add + 1;
353         
354         rset = xmalloc (sizeof(*rset) * rsmax);
355         rsfd = xmalloc (sizeof(*rsfd) * rsmax);
356         
357         for (i = from; i < to; i += i_add)
358         {
359             if (i_add <= to - i)
360                 rset[rscur] = rset_trunc_r (isam, isam_p, i, i+i_add,
361                                             merge_chunk);
362             else
363                 rset[rscur] = rset_trunc_r (isam, isam_p, i, to,
364                                             merge_chunk);
365             rscur++;
366         }
367         ti = heap_init (rscur, sizeof(struct it_key), key_compare);
368         for (i = rscur; --i >= 0; )
369         {
370             rsfd[i] = rset_open (rset[i], RSETF_READ|RSETF_SORT_SYSNO);
371             if (rset_read (rset[i], rsfd[i], ti->tmpbuf))
372                 heap_insert (ti, ti->tmpbuf, i);
373             else
374             {
375                 rset_close (rset[i], rsfd[i]);
376                 rset_delete (rset[i]);
377             }
378         }
379         while (ti->heapnum)
380         {
381             int n = ti->indx[ti->ptr[1]];
382
383             rset_write (result, result_rsfd, ti->heap[ti->ptr[1]]);
384
385             while (1)
386             {
387                 if (!rset_read (rset[n], rsfd[n], ti->tmpbuf))
388                 {
389                     heap_delete (ti);
390                     rset_close (rset[n], rsfd[n]);
391                     rset_delete (rset[n]);
392                     break;
393                 }
394                 if ((*ti->cmp)(ti->tmpbuf, ti->heap[ti->ptr[1]]) > 1)
395                 {
396                     heap_delete (ti);
397                     heap_insert (ti, ti->tmpbuf, n);
398                     break;
399                 }
400             }
401         }
402         xfree (rset);
403         xfree (rsfd);
404         heap_close (ti);
405     }
406     else
407     {
408         ISPT *ispt;
409         int i;
410         struct trunc_info *ti;
411
412         ispt = xmalloc (sizeof(*ispt) * (to-from));
413
414         ti = heap_init (to-from, sizeof(struct it_key),
415                         key_compare);
416         for (i = to-from; --i >= 0; )
417         {
418             ispt[i] = is_position (isam, isam_p[from+i]);
419             if (is_readkey (ispt[i], ti->tmpbuf))
420                 heap_insert (ti, ti->tmpbuf, i);
421             else
422                 is_pt_free (ispt[i]);
423         }
424         while (ti->heapnum)
425         {
426             int n = ti->indx[ti->ptr[1]];
427
428             rset_write (result, result_rsfd, ti->heap[ti->ptr[1]]);
429 #if 0
430 /* section that preserve all keys */
431             heap_delete (ti);
432             if (is_readkey (ispt[n], ti->tmpbuf))
433                 heap_insert (ti, ti->tmpbuf, n);
434             else
435                 is_pt_free (ispt[n]);
436 #else
437 /* section that preserve all keys with unique sysnos */
438             while (1)
439             {
440                 if (!is_readkey (ispt[n], ti->tmpbuf))
441                 {
442                     heap_delete (ti);
443                     is_pt_free (ispt[n]);
444                     break;
445                 }
446                 if ((*ti->cmp)(ti->tmpbuf, ti->heap[ti->ptr[1]]) > 1)
447                 {
448                     heap_delete (ti);
449                     heap_insert (ti, ti->tmpbuf, n);
450                     break;
451                 }
452             }
453 #endif
454         }
455         heap_close (ti);
456         xfree (ispt);
457     }
458     rset_close (result, result_rsfd);
459     return result;
460 }
461
462 static int isam_trunc_cmp (const void *p1, const void *p2)
463 {
464     ISAM_P i1 = *(ISAM_P*) p1;
465     ISAM_P i2 = *(ISAM_P*) p2;
466     int d;
467
468     d = is_type (i1) - is_type (i2);
469     if (d)
470         return d;
471     return is_block (i1) - is_block (i2);
472 }
473
474 static RSET rset_trunc (ISAM isam, ISAM_P *isam_p, int no)
475 {
476
477     qsort (isam_p, no, sizeof(*isam_p), isam_trunc_cmp);
478     return rset_trunc_r (isam, isam_p, 0, no, 100);
479 }
480
481 struct grep_info {
482     ISAM_P *isam_p_buf;
483     int isam_p_size;
484     int isam_p_indx;
485 };
486
487 static void add_isam_p (const char *info, struct grep_info *p)
488 {
489     if (p->isam_p_indx == p->isam_p_size)
490     {
491         ISAM_P *new_isam_p_buf;
492         
493         p->isam_p_size = 2*p->isam_p_size + 100;
494         new_isam_p_buf = xmalloc (sizeof(*new_isam_p_buf) *
495                                   p->isam_p_size);
496         if (p->isam_p_buf)
497         {
498             memcpy (new_isam_p_buf, p->isam_p_buf,
499                     p->isam_p_indx * sizeof(*p->isam_p_buf));
500             xfree (p->isam_p_buf);
501         }
502         p->isam_p_buf = new_isam_p_buf;
503     }
504     assert (*info == sizeof(*p->isam_p_buf));
505     memcpy (p->isam_p_buf + p->isam_p_indx, info+1, sizeof(*p->isam_p_buf));
506     (p->isam_p_indx)++;
507 }
508
509 static int grep_handle (char *name, const char *info, void *p)
510 {
511     logf (LOG_DEBUG, "dict name: %s", name);
512     add_isam_p (info, p);
513     return 0;
514 }
515
516 static void gen_regular_rel (char *dst, int val, int islt)
517 {
518     int dst_p = 1;
519     int w, d, i;
520     int pos = 0;
521     char numstr[20];
522
523     *dst = '(';
524     sprintf (numstr, "%d", val);
525     for (w = strlen(numstr); --w >= 0; pos++)
526     {
527         d = numstr[w];
528         if (pos > 0)
529         {
530             if (islt)
531             {
532                 if (d == '0')
533                     continue;
534                 d--;
535             } 
536             else
537             {
538                 if (d == '9')
539                     continue;
540                 d++;
541             }
542         }
543         
544         strcpy (dst + dst_p, numstr);
545         dst_p = strlen(dst) - pos - 1;
546
547         if (islt)
548         {
549             if (d != '0')
550             {
551                 dst[dst_p++] = '[';
552                 dst[dst_p++] = '0';
553                 dst[dst_p++] = '-';
554                 dst[dst_p++] = d;
555                 dst[dst_p++] = ']';
556             }
557             else
558                 dst[dst_p++] = d;
559         }
560         else
561         {
562             if (d != '9')
563             { 
564                 dst[dst_p++] = '[';
565                 dst[dst_p++] = d;
566                 dst[dst_p++] = '-';
567                 dst[dst_p++] = '9';
568                 dst[dst_p++] = ']';
569             }
570             else
571                 dst[dst_p++] = d;
572         }
573         for (i = 0; i<pos; i++)
574         {
575             dst[dst_p++] = '[';
576             dst[dst_p++] = '0';
577             dst[dst_p++] = '-';
578             dst[dst_p++] = '9';
579             dst[dst_p++] = ']';
580         }
581         dst[dst_p++] = '|';
582     }
583     dst[dst_p] = '\0';
584     if (islt)
585     {
586         for (i=1; i<pos; i++)
587             strcat (dst, "[0-9]?");
588     }
589     else
590     {
591         for (i = 0; i <= pos; i++)
592             strcat (dst, "[0-9]");
593         strcat (dst, "[0-9]*");
594     }
595     strcat (dst, ")");
596 }
597
598 static int relational_term (ZServerInfo *zi, Z_AttributesPlusTerm *zapt,
599                             const char *term_sub,
600                             char *term_dict,
601                             oid_value attributeSet,
602                             struct grep_info *grep_info,
603                             int *max_pos)
604 {
605     AttrType relation;
606     int relation_value;
607     int term_value;
608     int r;
609
610     attr_init (&relation, zapt, 2);
611     relation_value = attr_find (&relation, NULL);
612     term_value = atoi (term_sub);
613
614     switch (relation_value)
615     {
616     case 1:
617         if (term_value <= 0)
618             return 1;
619         logf (LOG_DEBUG, "Relation <");
620         gen_regular_rel (term_dict + strlen(term_dict), term_value-1, 1);
621         break;
622     case 2:
623         if (term_value < 0)
624             return 1;
625         logf (LOG_DEBUG, "Relation <=");
626         gen_regular_rel (term_dict + strlen(term_dict), term_value, 1);
627         break;
628     case 4:
629         if (term_value < 0)
630             term_value = 0;
631         logf (LOG_DEBUG, "Relation >=");
632         gen_regular_rel (term_dict + strlen(term_dict), term_value, 0);
633         break;
634     case 5:
635         if (term_value < 0)
636             term_value = 0;
637         logf (LOG_DEBUG, "Relation >");
638         gen_regular_rel (term_dict + strlen(term_dict), term_value+1, 0);
639         break;
640     default:
641         return 0;
642     }
643     logf (LOG_DEBUG, "dict_lookup_grep: %s", term_dict);
644     r = dict_lookup_grep (zi->wordDict, term_dict, 0, grep_info, max_pos,
645                           0, grep_handle);
646     if (r)
647         logf (LOG_WARN, "dict_lookup_grep fail, rel=gt: %d", r);
648     logf (LOG_DEBUG, "%d positions", grep_info->isam_p_indx);
649     return 1;
650 }
651
652 static void verbatim_char (int ch, int *indx, char *dst)
653 {
654     if (!isalnum (ch))
655         dst[(*indx)++] = '\\';
656     dst[(*indx)++] = ch;
657 }
658
659 static int trunc_term (ZServerInfo *zi, Z_AttributesPlusTerm *zapt,
660                        const char *term_sub,
661                        oid_value attributeSet, struct grep_info *grep_info,
662                        int num_bases, char **basenames)
663 {
664     char term_dict[2*IT_MAX_WORD+2];
665     int i, j, r, base_no;
666     AttrType truncation;
667     int truncation_value;
668     AttrType use;
669     int use_value;
670     oid_value curAttributeSet = attributeSet;
671
672     attr_init (&use, zapt, 1);
673     use_value = attr_find (&use, &curAttributeSet);
674     logf (LOG_DEBUG, "use value %d", use_value);
675     attr_init (&truncation, zapt, 5);
676     truncation_value = attr_find (&truncation, NULL);
677     logf (LOG_DEBUG, "truncation value %d", truncation_value);
678
679     if (use_value == -1)
680         use_value = 1016;
681
682     for (base_no = 0; base_no < num_bases; base_no++)
683     {
684         int max_pos;
685         int prefix_len = index_word_prefix_map (term_dict, curAttributeSet,
686                                                 use_value,
687                                                 basenames[base_no]);
688         if (prefix_len < 0)
689         {
690             zi->errCode = 114;
691             return -1;
692         }
693         if (!relational_term (zi, zapt, term_sub, term_dict,
694                               attributeSet, grep_info, &max_pos))
695         {
696             const char *cp;
697
698             j = prefix_len;
699             switch (truncation_value)
700             {
701             case -1:         /* not specified */
702             case 100:        /* do not truncate */
703                 term_dict[j++] = '(';
704                 for (i = 0; term_sub[i]; i++)
705                     verbatim_char (term_sub[i], &j, term_dict);
706                 strcpy (term_dict+j, ")");
707                 r = dict_lookup_grep (zi->wordDict, term_dict, 0, grep_info,
708                                       &max_pos, 0, grep_handle);
709                 if (r)
710                     logf (LOG_WARN, "dict_lookup_grep err, trunc=none:%d", r);
711                 break;
712             case 1:          /* right truncation */
713                 term_dict[j++] = '(';
714                 for (i = 0; term_sub[i]; i++)
715                     verbatim_char (term_sub[i], &j, term_dict);
716                 strcpy (term_dict+j, ".*)");
717                 dict_lookup_grep (zi->wordDict, term_dict, 0, grep_info,
718                                   &max_pos, 0, grep_handle);
719                 break;
720             case 2:          /* left truncation */
721             case 3:          /* left&right truncation */
722                 zi->errCode = 120;
723                 return -1;
724             case 101:        /* process # in term */
725                 term_dict[j++] = '(';
726                 for (i=0; term_sub[i]; i++)
727                     if (term_sub[i] == '#' && i > 2)
728                     {
729                         term_dict[j++] = '.';
730                         term_dict[j++] = '*';
731                     }
732                     else
733                         verbatim_char (term_sub[i], &j, term_dict);
734                 strcpy (term_dict+j, ")");
735                 r = dict_lookup_grep (zi->wordDict, term_dict, 0, grep_info,
736                                       &max_pos, 0, grep_handle);
737                 if (r)
738                     logf (LOG_WARN, "dict_lookup_grep err, trunc=#: %d",
739                           r);
740                 break;
741             case 102:        /* regular expression */
742                 sprintf (term_dict + j, "(%s)", term_sub);
743                 r = dict_lookup_grep (zi->wordDict, term_dict, 0, grep_info,
744                                       &max_pos, 0, grep_handle);
745                 if (r)
746                     logf (LOG_WARN, "dict_lookup_grep err, trunc=regular: %d",
747                           r);
748                 break;
749             case 103:        /* regular expression with error correction */
750                 cp = term_sub;
751                 r = 0;
752                 if (*cp == '*' && cp[1] && cp[2])
753                 {
754                     r = atoi (cp+1);
755                     cp += 2;
756                 }
757                 sprintf (term_dict + j, "(%s)", cp);
758                 r = dict_lookup_grep (zi->wordDict, term_dict, r, grep_info,
759                                       &max_pos, j, grep_handle);
760                 if (r)
761                     logf (LOG_WARN, "dict_lookup_grep err, trunc=eregular: %d",
762                           r);
763                 break;
764             }
765         }
766         if (max_pos <= strlen(basenames[base_no]))
767         {
768             zi->errCode = 109; /* Database unavailable */
769             zi->errString = basenames[base_no];
770             return -1;
771         }
772     }
773     logf (LOG_DEBUG, "%d positions", grep_info->isam_p_indx);
774     return 0;
775 }
776
777 static void trans_term (ZServerInfo *zi, Z_AttributesPlusTerm *zapt,
778                         char *termz)
779 {
780     size_t i, sizez;
781     Z_Term *term = zapt->term;
782
783     sizez = term->u.general->len;
784     if (sizez > IT_MAX_WORD)
785         sizez = IT_MAX_WORD;
786     for (i = 0; i < sizez; i++)
787         termz[i] = index_char_cvt (term->u.general->buf[i]);
788     termz[i] = '\0';
789 }
790
791 static RSET rpn_search_APT_relevance (ZServerInfo *zi, 
792                                       Z_AttributesPlusTerm *zapt,
793                                       oid_value attributeSet,
794                                       int num_bases, char **basenames)
795 {
796     rset_relevance_parms parms;
797     char termz[IT_MAX_WORD+1];
798     char term_sub[IT_MAX_WORD+1];
799     struct grep_info grep_info;
800     char *p0 = termz, *p1 = NULL;
801     RSET result;
802
803     parms.key_size = sizeof(struct it_key);
804     parms.max_rec = 100;
805     parms.cmp = key_compare;
806     parms.is = zi->wordIsam;
807
808     if (zapt->term->which != Z_Term_general)
809     {
810         zi->errCode = 124;
811         return NULL;
812     }
813     trans_term (zi, zapt, termz);
814     grep_info.isam_p_indx = 0;
815     grep_info.isam_p_size = 0;
816     grep_info.isam_p_buf = NULL;
817     while (1)
818     {
819         if ((p1 = strchr (p0, ' ')))
820         {
821             memcpy (term_sub, p0, p1-p0);
822             term_sub[p1-p0] = '\0';
823         }
824         else
825             strcpy (term_sub, p0);
826         if (trunc_term (zi, zapt, term_sub, attributeSet, &grep_info,
827                         num_bases, basenames))
828             return NULL;
829         if (!p1)
830             break;
831         p0 = p1;
832         while (*++p0 == ' ')
833             ;
834     }
835     parms.isam_positions = grep_info.isam_p_buf;
836     parms.no_isam_positions = grep_info.isam_p_indx;
837     if (grep_info.isam_p_indx > 0)
838         result = rset_create (rset_kind_relevance, &parms);
839     else
840         result = rset_create (rset_kind_null, NULL);
841     xfree (grep_info.isam_p_buf);
842     return result;
843 }
844
845 static RSET rpn_search_APT_word (ZServerInfo *zi,
846                                  Z_AttributesPlusTerm *zapt,
847                                  oid_value attributeSet,
848                                  int num_bases, char **basenames)
849 {
850     rset_isam_parms parms;
851     char termz[IT_MAX_WORD+1];
852     struct grep_info grep_info;
853     RSET result;
854
855     if (zapt->term->which != Z_Term_general)
856     {
857         zi->errCode = 124;
858         return NULL;
859     }
860     trans_term (zi, zapt, termz);
861
862     grep_info.isam_p_indx = 0;
863     grep_info.isam_p_size = 0;
864     grep_info.isam_p_buf = NULL;
865
866     if (trunc_term (zi, zapt, termz, attributeSet, &grep_info,
867                     num_bases, basenames))
868         return NULL;
869     if (grep_info.isam_p_indx < 1)
870         result = rset_create (rset_kind_null, NULL);
871     else if (grep_info.isam_p_indx == 1)
872     {
873         parms.is = zi->wordIsam;
874         parms.pos = *grep_info.isam_p_buf;
875         result = rset_create (rset_kind_isam, &parms);
876     }
877     else
878         result = rset_trunc (zi->wordIsam, grep_info.isam_p_buf,
879                              grep_info.isam_p_indx);
880     xfree (grep_info.isam_p_buf);
881     return result;
882 }
883
884 static RSET rpn_prox (RSET *rset, int rset_no)
885 {
886     int i;
887     RSFD *rsfd;
888     int  *more;
889     struct it_key **buf;
890     RSFD rsfd_result;
891     RSET result;
892     rset_temp_parms parms;
893     
894     rsfd = xmalloc (sizeof(*rsfd)*rset_no);
895     more = xmalloc (sizeof(*more)*rset_no);
896     buf = xmalloc (sizeof(*buf)*rset_no);
897
898     for (i = 0; i<rset_no; i++)
899     {
900         buf[i] = xmalloc (sizeof(**buf));
901         rsfd[i] = rset_open (rset[i], RSETF_READ|RSETF_SORT_SYSNO);
902         if (!(more[i] = rset_read (rset[i], rsfd[i], buf[i])))
903         {
904             while (i >= 0)
905             {
906                 rset_close (rset[i], rsfd[i]);
907                 xfree (buf[i]);
908                 --i;
909             }
910             xfree (rsfd);
911             xfree (more);
912             xfree (buf);
913             return rset_create (rset_kind_null, NULL);
914         }
915     }
916     parms.key_size = sizeof (struct it_key);
917     result = rset_create (rset_kind_temp, &parms);
918     rsfd_result = rset_open (result, RSETF_WRITE|RSETF_SORT_SYSNO);
919     
920     while (*more)
921     {
922         for (i = 1; i<rset_no; i++)
923         {
924             int cmp;
925             
926             if (!more[i])
927             {
928                 *more = 0;
929                 break;
930             }
931             cmp = key_compare (buf[i], buf[i-1]);
932             if (cmp > 1)
933             {
934                 more[i-1] = rset_read (rset[i-1], rsfd[i-1], buf[i-1]);
935                 break;
936             }
937             else if (cmp == 1)
938             {
939                 if (buf[i-1]->seqno+1 != buf[i]->seqno)
940                 {
941                     more[i-1] = rset_read (rset[i-1], rsfd[i-1], buf[i-1]);
942                     break;
943                 }
944             }
945             else
946             {
947                 more[i] = rset_read (rset[i], rsfd[i], buf[i]);
948                 break;
949             }
950         }
951         if (i == rset_no)
952         {
953             rset_write (result, rsfd_result, buf[0]);
954             more[0] = rset_read (*rset, *rsfd, *buf);
955         }
956     }
957     
958     for (i = 0; i<rset_no; i++)
959     {
960         rset_close (rset[i], rsfd[i]);
961         xfree (buf[i]);
962     }
963     rset_close (result, rsfd_result);
964     xfree (buf);
965     xfree (more);
966     xfree (rsfd);
967     return result;
968 }
969
970 static RSET rpn_search_APT_phrase (ZServerInfo *zi,
971                                    Z_AttributesPlusTerm *zapt,
972                                    oid_value attributeSet,
973                                    int num_bases, char **basenames)
974 {
975     char termz[IT_MAX_WORD+1];
976     char term_sub[IT_MAX_WORD+1];
977     char *p0 = termz, *p1 = NULL;
978     RSET rset[60], result;
979     int i, rset_no = 0;
980     struct grep_info grep_info;
981
982     if (zapt->term->which != Z_Term_general)
983     {
984         zi->errCode = 124;
985         return NULL;
986     }
987     trans_term (zi, zapt, termz);
988
989     grep_info.isam_p_size = 0;
990     grep_info.isam_p_buf = NULL;
991
992     while (1)
993     {
994         if ((p1 = strchr (p0, ' ')))
995         {
996             memcpy (term_sub, p0, p1-p0);
997             term_sub[p1-p0] = '\0';
998         }
999         else
1000             strcpy (term_sub, p0);
1001
1002         grep_info.isam_p_indx = 0;
1003         if (trunc_term (zi, zapt, term_sub, attributeSet, &grep_info,
1004                         num_bases, basenames))
1005             return NULL;
1006         if (grep_info.isam_p_indx == 0)
1007             rset[rset_no] = rset_create (rset_kind_null, NULL);
1008         else if (grep_info.isam_p_indx > 1)
1009             rset[rset_no] = rset_trunc (zi->wordIsam,
1010                                         grep_info.isam_p_buf,
1011                                         grep_info.isam_p_indx);
1012         else
1013         {
1014             rset_isam_parms parms;
1015             
1016             parms.is = zi->wordIsam;
1017             parms.pos = *grep_info.isam_p_buf;
1018             rset[rset_no] = rset_create (rset_kind_isam, &parms);
1019         }
1020         assert (rset[rset_no]);
1021         if (++rset_no >= sizeof(rset)/sizeof(*rset))
1022             break;
1023         if (!p1)
1024             break;
1025         p0 = p1;
1026         while (*++p0 == ' ')
1027             ;
1028     }
1029     xfree (grep_info.isam_p_buf);
1030     if (rset_no == 0)
1031         return rset_create (rset_kind_null, NULL);
1032     else if (rset_no == 1)
1033         return (rset[0]);
1034
1035     result = rpn_prox (rset, rset_no);
1036     for (i = 0; i<rset_no; i++)
1037         rset_delete (rset[i]);
1038     return result;
1039 }
1040
1041 static RSET rpn_search_APT_local (ZServerInfo *zi, Z_AttributesPlusTerm *zapt,
1042                                   oid_value attributeSet)
1043 {
1044     RSET result;
1045     RSFD rsfd;
1046     struct it_key key;
1047     rset_temp_parms parms;
1048     char termz[IT_MAX_WORD+1];
1049
1050     if (zapt->term->which != Z_Term_general)
1051     {
1052         zi->errCode = 124;
1053         return NULL;
1054     }
1055     parms.key_size = sizeof (struct it_key);
1056     result = rset_create (rset_kind_temp, &parms);
1057     rsfd = rset_open (result, RSETF_WRITE|RSETF_SORT_SYSNO);
1058
1059     trans_term (zi, zapt, termz);
1060     key.sysno = atoi (termz);
1061     if (key.sysno <= 0)
1062         key.sysno = 1;
1063     rset_write (result, rsfd, &key);
1064     rset_close (result, rsfd);
1065     return result;
1066 }
1067
1068 static RSET rpn_search_APT (ZServerInfo *zi, Z_AttributesPlusTerm *zapt,
1069                             oid_value attributeSet,
1070                             int num_bases, char **basenames)
1071 {
1072     AttrType relation;
1073     AttrType structure;
1074     int relation_value, structure_value;
1075
1076     attr_init (&relation, zapt, 2);
1077     attr_init (&structure, zapt, 4);
1078     
1079     relation_value = attr_find (&relation, NULL);
1080     structure_value = attr_find (&structure, NULL);
1081     switch (structure_value)
1082     {
1083     case -1:
1084         if (relation_value == 102) /* relevance relation */
1085             return rpn_search_APT_relevance (zi, zapt, attributeSet,
1086                                              num_bases, basenames);
1087         return rpn_search_APT_phrase (zi, zapt, attributeSet,
1088                                       num_bases, basenames);
1089     case 1: /* phrase */
1090         if (relation_value == 102) /* relevance relation */
1091             return rpn_search_APT_relevance (zi, zapt, attributeSet,
1092                                              num_bases, basenames);
1093         return rpn_search_APT_phrase (zi, zapt, attributeSet,
1094                                       num_bases, basenames);
1095         break;
1096     case 2: /* word */
1097         if (relation_value == 102) /* relevance relation */
1098             return rpn_search_APT_relevance (zi, zapt, attributeSet,
1099                                              num_bases, basenames);
1100         return rpn_search_APT_word (zi, zapt, attributeSet,
1101                                     num_bases, basenames);
1102     case 3: /* key */
1103         break;
1104     case 4: /* year */
1105         break;
1106     case 5: /* date - normalized */
1107         break;
1108     case 6: /* word list */
1109         return rpn_search_APT_relevance (zi, zapt, attributeSet,
1110                                          num_bases, basenames);
1111     case 100: /* date - un-normalized */
1112         break;
1113     case 101: /* name - normalized */
1114         break;
1115     case 102: /* date - un-normalized */
1116         break;
1117     case 103: /* structure */
1118         break;
1119     case 104: /* urx */
1120         break;
1121     case 105: /* free-form-text */
1122         return rpn_search_APT_relevance (zi, zapt, attributeSet,
1123                                          num_bases, basenames);
1124     case 106: /* document-text */
1125         return rpn_search_APT_relevance (zi, zapt, attributeSet,
1126                                          num_bases, basenames);
1127     case 107: /* local-number */
1128         return rpn_search_APT_local (zi, zapt, attributeSet);
1129     case 108: /* string */ 
1130         return rpn_search_APT_word (zi, zapt, attributeSet,
1131                                     num_bases, basenames);
1132     case 109: /* numeric string */
1133         break;
1134     }
1135     zi->errCode = 118;
1136     return NULL;
1137 }
1138
1139 static RSET rpn_search_ref (ZServerInfo *zi, Z_ResultSetId *resultSetId)
1140 {
1141     ZServerSet *s;
1142
1143     if (!(s = resultSetGet (zi, resultSetId)))
1144         return rset_create (rset_kind_null, NULL);
1145     return s->rset;
1146 }
1147
1148 static RSET rpn_search_structure (ZServerInfo *zi, Z_RPNStructure *zs,
1149                                   oid_value attributeSet,
1150                                   int num_bases, char **basenames)
1151 {
1152     RSET r = NULL;
1153     if (zs->which == Z_RPNStructure_complex)
1154     {
1155         rset_bool_parms bool_parms;
1156         int soft = 0;
1157
1158         bool_parms.rset_l = rpn_search_structure (zi, zs->u.complex->s1,
1159                                                   attributeSet,
1160                                                   num_bases, basenames);
1161         if (bool_parms.rset_l == NULL)
1162             return NULL;
1163         if (rset_is_ranked(bool_parms.rset_l))
1164             soft = 1;
1165         bool_parms.rset_r = rpn_search_structure (zi, zs->u.complex->s2,
1166                                                   attributeSet,
1167                                                   num_bases, basenames);
1168         if (bool_parms.rset_r == NULL)
1169         {
1170             rset_delete (bool_parms.rset_l);
1171             return NULL;
1172         }
1173         if (rset_is_ranked(bool_parms.rset_r))
1174             soft = 1;
1175         bool_parms.key_size = sizeof(struct it_key);
1176         bool_parms.cmp = key_compare;
1177
1178         switch (zs->u.complex->roperator->which)
1179         {
1180         case Z_Operator_and:
1181             r = rset_create (soft ? rset_kind_sand:rset_kind_and, &bool_parms);
1182             break;
1183         case Z_Operator_or:
1184             r = rset_create (soft ? rset_kind_sor:rset_kind_or, &bool_parms);
1185             break;
1186         case Z_Operator_and_not:
1187             r = rset_create (soft ? rset_kind_snot:rset_kind_not, &bool_parms);
1188             break;
1189         default:
1190             assert (0);
1191         }
1192     }
1193     else if (zs->which == Z_RPNStructure_simple)
1194     {
1195         if (zs->u.simple->which == Z_Operand_APT)
1196         {
1197             logf (LOG_DEBUG, "rpn_search_APT");
1198             r = rpn_search_APT (zi, zs->u.simple->u.attributesPlusTerm,
1199                                 attributeSet, num_bases, basenames);
1200         }
1201         else if (zs->u.simple->which == Z_Operand_resultSetId)
1202         {
1203             logf (LOG_DEBUG, "rpn_search_ref");
1204             r = rpn_search_ref (zi, zs->u.simple->u.resultSetId);
1205         }
1206         else
1207         {
1208             assert (0);
1209         }
1210     }
1211     else
1212     {
1213         assert (0);
1214     }
1215     return r;
1216 }
1217
1218 void count_set_save (RSET *r, int *count)
1219 {
1220     int psysno = 0;
1221     int kno = 0;
1222     struct it_key key;
1223     RSFD rfd, wfd;
1224     RSET w;
1225     rset_temp_parms parms;
1226
1227     logf (LOG_DEBUG, "count_set_save");
1228     *count = 0;
1229     parms.key_size = sizeof(struct it_key);
1230     w = rset_create (rset_kind_temp, &parms);
1231     wfd = rset_open (w, RSETF_WRITE|RSETF_SORT_SYSNO);
1232     rfd = rset_open (*r, RSETF_READ|RSETF_SORT_SYSNO);
1233     while (rset_read (*r, rfd, &key))
1234     {
1235         if (key.sysno != psysno)
1236         {
1237             rset_write (w, wfd, &key);
1238             psysno = key.sysno;
1239             (*count)++;
1240         }
1241         kno++;
1242     }
1243     rset_close (*r, rfd);
1244     rset_delete (*r);
1245     rset_close (w, wfd);
1246     *r = w;
1247     logf (LOG_DEBUG, "%d keys, %d distinct sysnos", kno, *count);
1248 }
1249
1250 static void count_set (RSET r, int *count)
1251 {
1252     int psysno = 0;
1253     int kno = 0;
1254     struct it_key key;
1255     RSFD rfd;
1256
1257     logf (LOG_DEBUG, "count_set");
1258     *count = 0;
1259     rfd = rset_open (r, RSETF_READ|RSETF_SORT_SYSNO);
1260     while (rset_read (r, rfd, &key))
1261     {
1262         if (key.sysno != psysno)
1263         {
1264             psysno = key.sysno;
1265             (*count)++;
1266         }
1267         kno++;
1268     }
1269     rset_close (r, rfd);
1270     logf (LOG_DEBUG, "%d keys, %d distinct sysnos", kno, *count);
1271 }
1272
1273 int rpn_search (ZServerInfo *zi,
1274                 Z_RPNQuery *rpn, int num_bases, char **basenames, 
1275                 const char *setname, int *hits)
1276 {
1277     RSET rset;
1278     oident *attrset;
1279     oid_value attributeSet;
1280
1281     zlog_rpn (rpn);
1282
1283     zi->errCode = 0;
1284     zi->errString = NULL;
1285
1286     attrset = oid_getentbyoid (rpn->attributeSetId);
1287     attributeSet = attrset->value;
1288     rset = rpn_search_structure (zi, rpn->RPNStructure, attributeSet,
1289                                  num_bases, basenames);
1290     if (!rset)
1291         return zi->errCode;
1292     if (rset_is_volatile(rset))
1293         count_set_save(&rset,hits);
1294     else
1295         count_set (rset, hits);
1296     resultSetAdd (zi, setname, 1, rset);
1297     if (zi->errCode)
1298         logf (LOG_DEBUG, "search error: %d", zi->errCode);
1299     return zi->errCode;
1300 }
1301
1302 struct scan_info {
1303     struct scan_entry *list;
1304     ODR odr;
1305     int before, after;
1306     ISAM isam;
1307     char prefix[20];
1308 };
1309
1310 static int scan_handle (char *name, const char *info, int pos, void *client)
1311 {
1312     int len_prefix, idx;
1313     ISAM_P isam_p;
1314     RSET rset;
1315     struct scan_info *scan_info = client;
1316
1317     rset_isam_parms parms;
1318
1319     len_prefix = strlen(scan_info->prefix);
1320     if (memcmp (name, scan_info->prefix, len_prefix))
1321         return 1;
1322     if (pos > 0)
1323         idx = scan_info->after - pos + scan_info->before;
1324     else
1325         idx = - pos - 1;
1326     scan_info->list[idx].term = odr_malloc (scan_info->odr,
1327                                             strlen(name + len_prefix)+1);
1328     strcpy (scan_info->list[idx].term, name + len_prefix);
1329     assert (*info == sizeof(isam_p));
1330     memcpy (&isam_p, info+1, sizeof(isam_p));
1331     parms.is = scan_info->isam;
1332     parms.pos = isam_p;
1333 #if 1
1334     rset = rset_create (rset_kind_isam, &parms);
1335     count_set (rset, &scan_info->list[idx].occurrences);
1336     rset_delete (rset);
1337 #else
1338     scan_info->list[idx].occurrences = 1;
1339 #endif
1340     logf (LOG_DEBUG, "pos=%3d idx=%3d name=%s", pos, idx, name);
1341     return 0;
1342 }
1343
1344
1345 static int dummy_handle (char *name, const char *info, void *p)
1346 {
1347     return 0;
1348 }
1349
1350 int rpn_scan (ZServerInfo *zi, Z_AttributesPlusTerm *zapt,
1351               int num_bases, char **basenames,
1352               int *position, int *num_entries, struct scan_entry **list,
1353               int *status)
1354 {
1355     int i, j, sizez, max_pos;
1356     int pos = *position;
1357     int num = *num_entries;
1358     int before;
1359     int after;
1360     char termz[IT_MAX_WORD+20];
1361     AttrType use;
1362     int use_value;
1363     Z_Term *term = zapt->term;
1364     struct scan_info scan_info;
1365
1366     logf (LOG_DEBUG, "scan, position = %d, num = %d", pos, num);
1367
1368     if (num_bases != 1)
1369         return 111;
1370     scan_info.before = before = pos-1;
1371     scan_info.after = after = 1+num-pos;
1372     scan_info.odr = zi->odr;
1373
1374     logf (LOG_DEBUG, "scan, before = %d, after = %d", before, after);
1375     
1376     scan_info.isam = zi->wordIsam;
1377     scan_info.list = odr_malloc (zi->odr, (before+after)*
1378                                  sizeof(*scan_info.list));
1379     for (j = 0; j<before+after; j++)
1380         scan_info.list[j].term = NULL;
1381     attr_init (&use, zapt, 1);
1382     use_value = attr_find (&use, NULL);
1383     logf (LOG_DEBUG, "use value %d", use_value);
1384
1385     if (use_value == -1)
1386         use_value = 1016;
1387     i = index_word_prefix (termz, 1, use_value, *basenames);
1388
1389     dict_lookup_grep (zi->wordDict, termz, 0, NULL, &max_pos, 0,
1390                       dummy_handle);
1391     if (max_pos <= strlen(*basenames))
1392     {
1393         zi->errString = *basenames;
1394         return zi->errCode = 109; /* Database unavailable */
1395     }
1396     strcpy (scan_info.prefix, termz);
1397     sizez = term->u.general->len;
1398     if (sizez > IT_MAX_WORD)
1399         sizez = IT_MAX_WORD;
1400     for (j = 0; j<sizez; j++)
1401         termz[j+i] = index_char_cvt (term->u.general->buf[j]);
1402     termz[j+i] = '\0';
1403     
1404     dict_scan (zi->wordDict, termz, &before, &after, &scan_info, scan_handle);
1405
1406     *status = BEND_SCAN_SUCCESS;
1407
1408     for (i = 0; i<scan_info.after; i++)
1409         if (scan_info.list[scan_info.before+scan_info.after-i-1].term)
1410             break;
1411     *num_entries -= i;
1412     if (i)
1413         *status = BEND_SCAN_PARTIAL;
1414
1415     for (i = 0; i<scan_info.before; i++)
1416         if (scan_info.list[i].term)
1417             break;
1418     if (i)
1419         *status = BEND_SCAN_PARTIAL;
1420     *position -= i;
1421     *num_entries -= i;
1422
1423     *list = scan_info.list+i;       /* list is set to first 'real' entry */
1424
1425     if (*num_entries == 0)          /* signal 'unsupported use-attribute' */
1426         zi->errCode = 114;          /* if no entries was found */
1427     logf (LOG_DEBUG, "position = %d, num_entries = %d",
1428           *position, *num_entries);
1429     if (zi->errCode)
1430         logf (LOG_DEBUG, "scan error: %d", zi->errCode);
1431     return 0;
1432 }
1433               
1434
1435