Home | History | Annotate | Line # | Download | only in stdlib
bsearch.c revision 1.5
      1  1.5  thorpej /*	$NetBSD: bsearch.c,v 1.5 1995/12/28 08:52:10 thorpej Exp $	*/
      2  1.5  thorpej 
      3  1.1      cgd /*
      4  1.1      cgd  * Copyright (c) 1990 Regents of the University of California.
      5  1.1      cgd  * All rights reserved.
      6  1.1      cgd  *
      7  1.1      cgd  * Redistribution and use in source and binary forms, with or without
      8  1.1      cgd  * modification, are permitted provided that the following conditions
      9  1.1      cgd  * are met:
     10  1.1      cgd  * 1. Redistributions of source code must retain the above copyright
     11  1.1      cgd  *    notice, this list of conditions and the following disclaimer.
     12  1.1      cgd  * 2. Redistributions in binary form must reproduce the above copyright
     13  1.1      cgd  *    notice, this list of conditions and the following disclaimer in the
     14  1.1      cgd  *    documentation and/or other materials provided with the distribution.
     15  1.1      cgd  * 3. All advertising materials mentioning features or use of this software
     16  1.1      cgd  *    must display the following acknowledgement:
     17  1.1      cgd  *	This product includes software developed by the University of
     18  1.1      cgd  *	California, Berkeley and its contributors.
     19  1.1      cgd  * 4. Neither the name of the University nor the names of its contributors
     20  1.1      cgd  *    may be used to endorse or promote products derived from this software
     21  1.1      cgd  *    without specific prior written permission.
     22  1.1      cgd  *
     23  1.1      cgd  * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
     24  1.1      cgd  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
     25  1.1      cgd  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
     26  1.1      cgd  * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
     27  1.1      cgd  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
     28  1.1      cgd  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
     29  1.1      cgd  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
     30  1.1      cgd  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
     31  1.1      cgd  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
     32  1.1      cgd  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
     33  1.1      cgd  * SUCH DAMAGE.
     34  1.1      cgd  */
     35  1.1      cgd 
     36  1.1      cgd #if defined(LIBC_SCCS) && !defined(lint)
     37  1.5  thorpej #if 0
     38  1.5  thorpej static char *sccsid = "from: @(#)bsearch.c	5.4 (Berkeley) 2/23/91";
     39  1.5  thorpej #else
     40  1.5  thorpej static char *rcsid = "$NetBSD: bsearch.c,v 1.5 1995/12/28 08:52:10 thorpej Exp $";
     41  1.5  thorpej #endif
     42  1.1      cgd #endif /* LIBC_SCCS and not lint */
     43  1.1      cgd 
     44  1.1      cgd #include <stdlib.h>
     45  1.1      cgd 
     46  1.1      cgd /*
     47  1.1      cgd  * Perform a binary search.
     48  1.1      cgd  *
     49  1.1      cgd  * The code below is a bit sneaky.  After a comparison fails, we
     50  1.1      cgd  * divide the work in half by moving either left or right. If lim
     51  1.1      cgd  * is odd, moving left simply involves halving lim: e.g., when lim
     52  1.1      cgd  * is 5 we look at item 2, so we change lim to 2 so that we will
     53  1.1      cgd  * look at items 0 & 1.  If lim is even, the same applies.  If lim
     54  1.1      cgd  * is odd, moving right again involes halving lim, this time moving
     55  1.1      cgd  * the base up one item past p: e.g., when lim is 5 we change base
     56  1.1      cgd  * to item 3 and make lim 2 so that we will look at items 3 and 4.
     57  1.1      cgd  * If lim is even, however, we have to shrink it by one before
     58  1.1      cgd  * halving: e.g., when lim is 4, we still looked at item 2, so we
     59  1.1      cgd  * have to make lim 3, then halve, obtaining 1, so that we will only
     60  1.1      cgd  * look at item 3.
     61  1.1      cgd  */
     62  1.1      cgd void *
     63  1.1      cgd bsearch(key, base0, nmemb, size, compar)
     64  1.1      cgd 	register const void *key;
     65  1.1      cgd 	const void *base0;
     66  1.1      cgd 	size_t nmemb;
     67  1.1      cgd 	register size_t size;
     68  1.1      cgd 	register int (*compar) __P((const void *, const void *));
     69  1.1      cgd {
     70  1.1      cgd 	register const char *base = base0;
     71  1.1      cgd 	register int lim, cmp;
     72  1.1      cgd 	register const void *p;
     73  1.1      cgd 
     74  1.1      cgd 	for (lim = nmemb; lim != 0; lim >>= 1) {
     75  1.1      cgd 		p = base + (lim >> 1) * size;
     76  1.1      cgd 		cmp = (*compar)(key, p);
     77  1.1      cgd 		if (cmp == 0)
     78  1.1      cgd 			return ((void *)p);
     79  1.1      cgd 		if (cmp > 0) {	/* key > p: move right */
     80  1.1      cgd 			base = (char *)p + size;
     81  1.1      cgd 			lim--;
     82  1.1      cgd 		} /* else move left */
     83  1.1      cgd 	}
     84  1.1      cgd 	return (NULL);
     85  1.1      cgd }
     86