boyermoore.c

/**
 * This is an implementation of the (turbo) bayer/moore string search algorithm.
 * Based upon the code presented at http://www-igm.univ-mlv.fr/~lecroq/string/.
 *
 * @file boyermoore.c
 */

#include "boyermoore.h"

#include <limits.h>
#include <stdlib.h>


#define MAX(a,b) ((a)>(b)?(a):(b))

static void suffixes(char *x, int m, unsigned short *suff)
{
   int g, i;
   int f = 0; /* This is to silent the compiler */

   suff[m - 1] = m;
   g = m - 1;
   for (i = m - 2; i >= 0; --i)
   {
      if (i > g && suff[i + m - 1 - f] < i - g)
         suff[i] = suff[i + m - 1 - f];
      else {
         if (i < g)
            g = i;
         f = i;
         while (g >= 0 && x[g] == x[g + m - 1 - f])
            --g;
         suff[i] = f - g;
      }
   }
}

/**
 * Computes the good suffix skip table.
 *
 * @param x
 * @param m
 * @param suff
 * @param bmGs
 */
static void preBmGs(char *x, int m, unsigned short suff[], unsigned short bmGs[])
{
   int i, j;

   suffixes(x, m, suff);

   for (i = 0; i < m; ++i)
      bmGs[i] = m;
   j = 0;
   for (i = m - 1; i >= 0; --i)
      if (suff[i] == i + 1)
         for (; j < m - 1 - i; ++j)
            if (bmGs[j] == m)
               bmGs[j] = m - 1 - i;
   for (i = 0; i <= m - 2; ++i)
      bmGs[m - 1 - suff[i]] = m - 1 - i;
}


/**
 * @brief The context data structure of the boyermoore search
 */
struct boyermoore_context
{
	/** @brief The bad character skip table */
	unsigned short skip_table[UCHAR_MAX];

	/** @brief The good suffix table */
	unsigned short *good_suffix_table;

	char *pat; /**! @brief pattern to be searched for */
	int plen; /**! @brief length of the pattern */
};

/*****************************************************************************/

struct boyermoore_context *boyermoore_create_context(char *p, int plen)
{
	struct boyermoore_context *context;
	unsigned short *next;
	unsigned short *suff;
	int i;

	if (!(context = (struct boyermoore_context*)malloc(sizeof(*context))))
		return NULL;

	if (!(context->good_suffix_table = next = (unsigned short*)malloc((plen + 1)*sizeof(context->good_suffix_table[0]))))
	{
		free(context);
		return NULL;
	}

	if (!(suff = (unsigned short*)(unsigned short*)malloc((plen + 1)*sizeof(suff[0]))))
	{
		free(next);
		free(context);
		return NULL;
	}

	context->pat = p;
	context->plen = plen;

	/* Prepare bad character skip table */
	for (i = 0; i < (int)(sizeof(context->skip_table)/sizeof(context->skip_table[0])); i++)
		context->skip_table[i] = plen;
	for (i = 0; i < plen - 1; i++)
		context->skip_table[(unsigned char)p[i]] = plen - i - 1;

    preBmGs(p, plen, suff, next);
    free(suff);
	return context;
}

/*****************************************************************************/

void boyermoore_delete_context(struct boyermoore_context *context)
{
	if (context)
	{
		free(context->good_suffix_table);
		free(context);
	}
}

/*****************************************************************************/

int boyermoore(struct boyermoore_context *context, char *str, int n, bm_callback callback, void *user_data)
{
	int i, j;
	int rc = -1;

	unsigned short *bmBc = context->skip_table;
	unsigned short *bmGs = context->good_suffix_table;
	int m = context->plen;
	char *substr = context->pat;

	/* Searching */
	j = 0;
	while (j <= n - m)
	{
		for (i = m - 1; i >= 0 && substr[i] == str[i + j]; --i);
		if (i < 0)
		{
			rc = j;
			if (!callback || !callback(substr,j,user_data))
				return rc;
			j += bmGs[0];
		}  else
		{
			j += MAX(bmGs[i], bmBc[(unsigned char)str[i + j]] - m + 1 + i);
		}
	}

   return rc;
}