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