Home | History | Annotate | Line # | Download | only in bigram
locate.bigram.c revision 1.7
      1  1.7  jdolecek /*	$NetBSD: locate.bigram.c,v 1.7 2000/03/20 19:17:35 jdolecek Exp $	*/
      2  1.3       jtc 
      3  1.1       cgd /*
      4  1.3       jtc  * Copyright (c) 1989, 1993
      5  1.3       jtc  *	The Regents of the University of California.  All rights reserved.
      6  1.1       cgd  *
      7  1.1       cgd  * This code is derived from software contributed to Berkeley by
      8  1.1       cgd  * James A. Woods.
      9  1.1       cgd  *
     10  1.1       cgd  * Redistribution and use in source and binary forms, with or without
     11  1.1       cgd  * modification, are permitted provided that the following conditions
     12  1.1       cgd  * are met:
     13  1.1       cgd  * 1. Redistributions of source code must retain the above copyright
     14  1.1       cgd  *    notice, this list of conditions and the following disclaimer.
     15  1.1       cgd  * 2. Redistributions in binary form must reproduce the above copyright
     16  1.1       cgd  *    notice, this list of conditions and the following disclaimer in the
     17  1.1       cgd  *    documentation and/or other materials provided with the distribution.
     18  1.1       cgd  * 3. All advertising materials mentioning features or use of this software
     19  1.1       cgd  *    must display the following acknowledgement:
     20  1.1       cgd  *	This product includes software developed by the University of
     21  1.1       cgd  *	California, Berkeley and its contributors.
     22  1.1       cgd  * 4. Neither the name of the University nor the names of its contributors
     23  1.1       cgd  *    may be used to endorse or promote products derived from this software
     24  1.1       cgd  *    without specific prior written permission.
     25  1.1       cgd  *
     26  1.1       cgd  * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
     27  1.1       cgd  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
     28  1.1       cgd  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
     29  1.1       cgd  * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
     30  1.1       cgd  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
     31  1.1       cgd  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
     32  1.1       cgd  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
     33  1.1       cgd  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
     34  1.1       cgd  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
     35  1.1       cgd  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
     36  1.1       cgd  * SUCH DAMAGE.
     37  1.1       cgd  */
     38  1.1       cgd 
     39  1.6     lukem #include <sys/cdefs.h>
     40  1.1       cgd #ifndef lint
     41  1.6     lukem __COPYRIGHT("@(#) Copyright (c) 1989, 1993\n\
     42  1.6     lukem 	The Regents of the University of California.  All rights reserved.\n");
     43  1.1       cgd #endif /* not lint */
     44  1.1       cgd 
     45  1.1       cgd #ifndef lint
     46  1.3       jtc #if 0
     47  1.4       jtc static char sccsid[] = "@(#)locate.bigram.c	8.2 (Berkeley) 4/28/95";
     48  1.3       jtc #endif
     49  1.7  jdolecek __RCSID("$NetBSD: locate.bigram.c,v 1.7 2000/03/20 19:17:35 jdolecek Exp $");
     50  1.1       cgd #endif /* not lint */
     51  1.1       cgd 
     52  1.1       cgd /*
     53  1.1       cgd  *  bigram < text > bigrams
     54  1.1       cgd  *
     55  1.1       cgd  * List bigrams for 'updatedb' script.
     56  1.1       cgd  * Use 'code' to encode a file using this output.
     57  1.1       cgd  */
     58  1.1       cgd 
     59  1.1       cgd #include <stdio.h>
     60  1.7  jdolecek #include <stdlib.h>
     61  1.1       cgd #include <sys/param.h>			/* for MAXPATHLEN */
     62  1.1       cgd 
     63  1.7  jdolecek int	main __P((int, char **));
     64  1.7  jdolecek static int compare_bigrams __P((const void *, const void *));
     65  1.7  jdolecek static void add_bigram __P((u_char, u_char));
     66  1.7  jdolecek 
     67  1.7  jdolecek static char buf1[MAXPATHLEN] = " ";
     68  1.7  jdolecek static char buf2[MAXPATHLEN];
     69  1.7  jdolecek 
     70  1.7  jdolecek struct bigram {
     71  1.7  jdolecek 	int count;
     72  1.7  jdolecek 	u_char b1, b2;	/* needed for final sorting */
     73  1.7  jdolecek };
     74  1.7  jdolecek 
     75  1.7  jdolecek struct bigram bigrams[256 * 256];
     76  1.7  jdolecek 
     77  1.7  jdolecek static void
     78  1.7  jdolecek add_bigram(i1, i2)
     79  1.7  jdolecek 	u_char i1, i2;
     80  1.7  jdolecek {
     81  1.7  jdolecek 	if (i1 == '\n') {
     82  1.7  jdolecek 		bigrams[0].count++;
     83  1.7  jdolecek 		i1 = 0;
     84  1.7  jdolecek 	}
     85  1.7  jdolecek 	if (i2 == '\n') {
     86  1.7  jdolecek 		bigrams[0].count++;
     87  1.7  jdolecek 		i2 = i1;
     88  1.7  jdolecek 		i1 = 0;
     89  1.7  jdolecek 	}
     90  1.7  jdolecek 
     91  1.7  jdolecek 	bigrams[(i1<<8)+i2].count++;
     92  1.7  jdolecek }
     93  1.7  jdolecek 
     94  1.7  jdolecek static int
     95  1.7  jdolecek compare_bigrams(item1, item2)
     96  1.7  jdolecek 	const void *item1, *item2;
     97  1.7  jdolecek {
     98  1.7  jdolecek 	const struct bigram *it1=item1, *it2=item2;
     99  1.1       cgd 
    100  1.7  jdolecek 	if (it1->count != it2->count)
    101  1.7  jdolecek 		return it2->count - it1->count;
    102  1.7  jdolecek 	else if (it1->b1 != it2->b1)
    103  1.7  jdolecek 		return it2->b1 - it1->b1;
    104  1.7  jdolecek 	else
    105  1.7  jdolecek 		return it2->b2 - it2->b2;
    106  1.7  jdolecek }
    107  1.6     lukem 
    108  1.6     lukem int
    109  1.6     lukem main(argc, argv)
    110  1.6     lukem 	int argc;
    111  1.6     lukem 	char *argv[];
    112  1.1       cgd {
    113  1.6     lukem   	char *cp;
    114  1.6     lukem 	char *oldpath = buf1, *path = buf2;
    115  1.7  jdolecek 	struct bigram *bg;
    116  1.7  jdolecek 	int i;
    117  1.7  jdolecek 
    118  1.7  jdolecek 	/* initialize bigram array */
    119  1.7  jdolecek 	memset(bigrams, 0, sizeof(bigrams));
    120  1.7  jdolecek 	for(i=0; i < 65536; i++) {
    121  1.7  jdolecek 		bigrams[i].b1 = i / 256;
    122  1.7  jdolecek 		bigrams[i].b2 = i % 256;
    123  1.7  jdolecek 	}
    124  1.1       cgd 
    125  1.1       cgd      	while ( fgets ( path, sizeof(buf2), stdin ) != NULL ) {
    126  1.1       cgd 
    127  1.1       cgd 		/* skip longest common prefix */
    128  1.1       cgd 		for ( cp = path; *cp == *oldpath; cp++, oldpath++ )
    129  1.4       jtc 			if ( *oldpath == '\0' )
    130  1.1       cgd 				break;
    131  1.7  jdolecek 
    132  1.1       cgd 		/*
    133  1.1       cgd 		 * output post-residue bigrams only
    134  1.1       cgd 		 */
    135  1.7  jdolecek 		for(; cp[0] != '\0' && cp[1] != '\0'; cp += 2)
    136  1.7  jdolecek 			add_bigram((u_char)cp[0], (u_char)cp[1]);
    137  1.7  jdolecek 
    138  1.7  jdolecek 		if (path == buf1)		/* swap pointers */
    139  1.1       cgd 			path = buf2, oldpath = buf1;
    140  1.1       cgd 		else
    141  1.1       cgd 			path = buf1, oldpath = buf2;
    142  1.1       cgd    	}
    143  1.7  jdolecek 
    144  1.7  jdolecek 	/* sort the bigrams by how many times it appeared and their value */
    145  1.7  jdolecek 	heapsort((void *)bigrams, 256 * 256, sizeof(struct bigram),
    146  1.7  jdolecek 			compare_bigrams);
    147  1.7  jdolecek 
    148  1.7  jdolecek 	/* write 128 most frequent bigrams out */
    149  1.7  jdolecek 	bg = bigrams;
    150  1.7  jdolecek 	for (i = 0; i < 128 && bg->count > 0; i++, bg++) {
    151  1.7  jdolecek 		if (bg->b1 != '\0')
    152  1.7  jdolecek 			fputc(bg->b1, stdout);
    153  1.7  jdolecek 		if (bg->b2 != '\0')
    154  1.7  jdolecek 			fputc(bg->b2, stdout);
    155  1.7  jdolecek 	}
    156  1.7  jdolecek 
    157  1.4       jtc 	return (0);
    158  1.1       cgd }
    159