Fix sin, sincos missing underflows (bug 16526, bug 16538).
[glibc.git] / misc / hsearch_r.c
blob9f55e845cf7800e370b3d476655a1e5a86f0055d
1 /* Copyright (C) 1993-2015 Free Software Foundation, Inc.
2 This file is part of the GNU C Library.
3 Contributed by Ulrich Drepper <drepper@gnu.ai.mit.edu>, 1993.
5 The GNU C Library is free software; you can redistribute it and/or
6 modify it under the terms of the GNU Lesser General Public
7 License as published by the Free Software Foundation; either
8 version 2.1 of the License, or (at your option) any later version.
10 The GNU C Library is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 Lesser General Public License for more details.
15 You should have received a copy of the GNU Lesser General Public
16 License along with the GNU C Library; if not, see
17 <http://www.gnu.org/licenses/>. */
19 #include <errno.h>
20 #include <malloc.h>
21 #include <string.h>
23 #include <search.h>
25 /* [Aho,Sethi,Ullman] Compilers: Principles, Techniques and Tools, 1986
26 [Knuth] The Art of Computer Programming, part 3 (6.4) */
29 /* The reentrant version has no static variables to maintain the state.
30 Instead the interface of all functions is extended to take an argument
31 which describes the current status. */
32 typedef struct _ENTRY
34 unsigned int used;
35 ENTRY entry;
37 _ENTRY;
40 /* For the used double hash method the table size has to be a prime. To
41 correct the user given table size we need a prime test. This trivial
42 algorithm is adequate because
43 a) the code is (most probably) called a few times per program run and
44 b) the number is small because the table must fit in the core */
45 static int
46 isprime (unsigned int number)
48 /* no even number will be passed */
49 unsigned int div = 3;
51 while (div * div < number && number % div != 0)
52 div += 2;
54 return number % div != 0;
58 /* Before using the hash table we must allocate memory for it.
59 Test for an existing table are done. We allocate one element
60 more as the found prime number says. This is done for more effective
61 indexing as explained in the comment for the hsearch function.
62 The contents of the table is zeroed, especially the field used
63 becomes zero. */
64 int
65 __hcreate_r (nel, htab)
66 size_t nel;
67 struct hsearch_data *htab;
69 /* Test for correct arguments. */
70 if (htab == NULL)
72 __set_errno (EINVAL);
73 return 0;
76 /* There is still another table active. Return with error. */
77 if (htab->table != NULL)
78 return 0;
80 /* We need a size of at least 3. Otherwise the hash functions we
81 use will not work. */
82 if (nel < 3)
83 nel = 3;
84 /* Change nel to the first prime number not smaller as nel. */
85 nel |= 1; /* make odd */
86 while (!isprime (nel))
87 nel += 2;
89 htab->size = nel;
90 htab->filled = 0;
92 /* allocate memory and zero out */
93 htab->table = (_ENTRY *) calloc (htab->size + 1, sizeof (_ENTRY));
94 if (htab->table == NULL)
95 return 0;
97 /* everything went alright */
98 return 1;
100 libc_hidden_def (__hcreate_r)
101 weak_alias (__hcreate_r, hcreate_r)
104 /* After using the hash table it has to be destroyed. The used memory can
105 be freed and the local static variable can be marked as not used. */
106 void
107 __hdestroy_r (htab)
108 struct hsearch_data *htab;
110 /* Test for correct arguments. */
111 if (htab == NULL)
113 __set_errno (EINVAL);
114 return;
117 /* Free used memory. */
118 free (htab->table);
120 /* the sign for an existing table is an value != NULL in htable */
121 htab->table = NULL;
123 libc_hidden_def (__hdestroy_r)
124 weak_alias (__hdestroy_r, hdestroy_r)
127 /* This is the search function. It uses double hashing with open addressing.
128 The argument item.key has to be a pointer to an zero terminated, most
129 probably strings of chars. The function for generating a number of the
130 strings is simple but fast. It can be replaced by a more complex function
131 like ajw (see [Aho,Sethi,Ullman]) if the needs are shown.
133 We use an trick to speed up the lookup. The table is created by hcreate
134 with one more element available. This enables us to use the index zero
135 special. This index will never be used because we store the first hash
136 index in the field used where zero means not used. Every other value
137 means used. The used field can be used as a first fast comparison for
138 equality of the stored and the parameter value. This helps to prevent
139 unnecessary expensive calls of strcmp. */
141 __hsearch_r (item, action, retval, htab)
142 ENTRY item;
143 ACTION action;
144 ENTRY **retval;
145 struct hsearch_data *htab;
147 unsigned int hval;
148 unsigned int count;
149 unsigned int len = strlen (item.key);
150 unsigned int idx;
152 /* Compute an value for the given string. Perhaps use a better method. */
153 hval = len;
154 count = len;
155 while (count-- > 0)
157 hval <<= 4;
158 hval += item.key[count];
160 if (hval == 0)
161 ++hval;
163 /* First hash function: simply take the modul but prevent zero. */
164 idx = hval % htab->size + 1;
166 if (htab->table[idx].used)
168 /* Further action might be required according to the action value. */
169 if (htab->table[idx].used == hval
170 && strcmp (item.key, htab->table[idx].entry.key) == 0)
172 *retval = &htab->table[idx].entry;
173 return 1;
176 /* Second hash function, as suggested in [Knuth] */
177 unsigned int hval2 = 1 + hval % (htab->size - 2);
178 unsigned int first_idx = idx;
182 /* Because SIZE is prime this guarantees to step through all
183 available indeces. */
184 if (idx <= hval2)
185 idx = htab->size + idx - hval2;
186 else
187 idx -= hval2;
189 /* If we visited all entries leave the loop unsuccessfully. */
190 if (idx == first_idx)
191 break;
193 /* If entry is found use it. */
194 if (htab->table[idx].used == hval
195 && strcmp (item.key, htab->table[idx].entry.key) == 0)
197 *retval = &htab->table[idx].entry;
198 return 1;
201 while (htab->table[idx].used);
204 /* An empty bucket has been found. */
205 if (action == ENTER)
207 /* If table is full and another entry should be entered return
208 with error. */
209 if (htab->filled == htab->size)
211 __set_errno (ENOMEM);
212 *retval = NULL;
213 return 0;
216 htab->table[idx].used = hval;
217 htab->table[idx].entry = item;
219 ++htab->filled;
221 *retval = &htab->table[idx].entry;
222 return 1;
225 __set_errno (ESRCH);
226 *retval = NULL;
227 return 0;
229 libc_hidden_def (__hsearch_r)
230 weak_alias (__hsearch_r, hsearch_r)