home *** CD-ROM | disk | FTP | other *** search
/ Mac-Source 1994 July / Mac-Source_July_1994.iso / C and C++ / Libraries / stringsearch / bmsource / uf.om.inc.c < prev    next >
Encoding:
C/C++ Source or Header  |  1991-05-06  |  2.8 KB  |  147 lines  |  [TEXT/MPS ]
  1. /*
  2.     search routine generated by gen.
  3.     skip=uf, match=om (using omr), shift=inc
  4. */
  5. #include    "freq.h"
  6. /*
  7.  * The authors of this software are Andrew Hume and Daniel Sunday.
  8.  * 
  9.  * Copyright (c) 1991 by AT&T and Daniel Sunday.
  10.  * 
  11.  * Permission to use, copy, modify, and distribute this software for any
  12.  * purpose without fee is hereby granted, provided that this entire notice
  13.  * is included in all copies of any software which is or includes a copy
  14.  * or modification of this software and in all copies of the supporting
  15.  * documentation for such software.
  16.  * 
  17.  * THIS SOFTWARE IS BEING PROVIDED "AS IS", WITHOUT ANY EXPRESS OR IMPLIED
  18.  * WARRANTY.  IN PARTICULAR, NEITHER THE AUTHORS NOR AT&T MAKE ANY
  19.  * REPRESENTATION OR WARRANTY OF ANY KIND CONCERNING THE MERCHANTABILITY
  20.  * OF THIS SOFTWARE OR ITS FITNESS FOR ANY PARTICULAR PURPOSE.
  21.  */
  22.  
  23. #ifndef    CHARTYPE
  24. #define    CHARTYPE    unsigned char
  25. #endif
  26. #define    MAXPAT    256
  27.  
  28. static struct om
  29. {
  30.     int loc;
  31.     CHARTYPE c;
  32. };
  33.  
  34. om_pcmp(p1, p2)
  35.     register struct om *p1, *p2;
  36. {
  37.     register float f = freq[p1->c] - freq[p2->c];
  38.  
  39.     if(f < 0)
  40.         return(-1);
  41.     else if(f > 0)
  42.         return(1);
  43.     return(p2->loc - p1->loc);
  44. }
  45. #include    "stats.h"
  46.  
  47. #ifndef    TABTYPE
  48. #define    TABTYPE    long
  49. #endif
  50. typedef TABTYPE Tab;
  51.  
  52. static struct
  53. {
  54.     int patlen;
  55.     CHARTYPE pat[MAXPAT];
  56.     Tab delta[256];
  57.     struct om om[1024];
  58.     int nom;
  59. } pat;
  60.  
  61. prep(base, m)
  62.     CHARTYPE *base;
  63.     register m;
  64. {
  65.     CHARTYPE *skipc;
  66.     register CHARTYPE *pe, *p;
  67.     register int j;
  68.     register Tab *d;
  69.     register struct om *op;
  70.     register CHARTYPE *opb, *ope;
  71.  
  72.     pat.patlen = m;
  73.     if(m > MAXPAT)
  74.         abort();
  75.     memcpy(pat.pat, base, m);
  76.     skipc = 0;
  77.     stats.len = m;
  78.     d = pat.delta;
  79.     for(j = 0; j < 256; j++)
  80.         d[j] = pat.patlen;
  81.     for(p = pat.pat, pe = p+m-1; p < pe; p++)
  82.         d[*p] = pe-p;
  83.     d[*p] = 0;
  84.     skipc = (CHARTYPE *)p;
  85.     for(op = pat.om, opb = pat.pat, ope = opb+m-1; opb < ope; opb++){
  86.         op->loc = opb - ope;
  87.         op->c = *opb;
  88.         op++;
  89.     }
  90.     qsort(pat.om, pat.nom = op-pat.om, sizeof(pat.om[0]), om_pcmp);
  91. }
  92.  
  93. exec(base, n)
  94.     CHARTYPE *base;
  95. {
  96.     int nmatch = 0;
  97.     register CHARTYPE *e, *s;
  98.     register Tab *d0 = pat.delta;
  99.     register k;
  100.     register struct om *op, *oep;
  101.  
  102.     s = base+pat.patlen-1;
  103.     e = base+n;
  104.     memset(e, pat.pat[pat.patlen-1], pat.patlen);
  105.     oep = pat.om+pat.nom;
  106.     while(s < e){
  107. #ifdef    STATS
  108.         k = d0[*s];
  109.         stats.jump++;
  110.         while(k){
  111.             stats.jump++; stats.step[k]++;
  112.             k = d0[*(s += k)];
  113.             stats.jump++; stats.step[k]++;
  114.             k = d0[*(s += k)];
  115.             stats.jump++; stats.step[k]++;
  116.             k = d0[*(s += k)];
  117.         }
  118. #else
  119.         k = d0[*s];
  120.         while(k){
  121.             k = d0[*(s += k)];
  122.             k = d0[*(s += k)];
  123.             k = d0[*(s += k)];
  124.         }
  125. #endif
  126.         if(s >= e)
  127.             break;
  128. #ifdef    STATS
  129.         stats.slow++;
  130. #endif
  131.         for(op = pat.om; op < oep; ++op){
  132. #ifdef    STATS
  133.             stats.cmp++;
  134. #endif
  135.             if(op->c != s[op->loc])
  136.                 goto mismatch;
  137.         }
  138.         nmatch++;
  139.     mismatch:
  140.         s++;
  141. #ifdef    STATS
  142.         stats.step[1]++;
  143. #endif
  144.     }
  145.     return(nmatch);
  146. }
  147.