cdd7d8f818e839a0c4a40ff08e70a0d921e893cb
[sip-router] / hash_func.c
1 /*
2  * $Id$
3  *
4  * Copyright (C) 2001-2003 Fhg Fokus
5  *
6  * This file is part of ser, a free SIP server.
7  *
8  * ser is free software; you can redistribute it and/or modify
9  * it under the terms of the GNU General Public License as published by
10  * the Free Software Foundation; either version 2 of the License, or
11  * (at your option) any later version
12  *
13  * For a license to use the ser software under conditions
14  * other than those described here, or to purchase support for this
15  * software, please contact iptel.org by e-mail at the following addresses:
16  *    info@iptel.org
17  *
18  * ser is distributed in the hope that it will be useful,
19  * but WITHOUT ANY WARRANTY; without even the implied warranty of
20  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
21  * GNU General Public License for more details.
22  *
23  * You should have received a copy of the GNU General Public License 
24  * along with this program; if not, write to the Free Software 
25  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
26  */
27
28
29
30 #ifndef _CRC_H_
31 #define _CRC_H_
32
33 extern unsigned long int crc_32_tab[];
34 extern unsigned short int ccitt_tab[];
35 extern unsigned short int crc_16_tab[];
36
37 #endif
38
39 #include <stdio.h>
40 #include <string.h>
41 #include <stdlib.h>
42 #include "hash_func.h"
43 #include "dprint.h"
44 #include "crc.h"
45 #include "ut.h"
46
47
48 int new_hash( str call_id, str cseq_nr )
49 {
50         int hash_code = 0;
51         int i,j, k, third;
52         int ci_len, cs_len;
53         char *ci, *cs;
54
55         /* trim EoLs */
56 /*
57         ci_len = call_id.len;
58         while (ci_len && ((c=call_id.s[ci_len-1])==0 || c=='\r' || c=='\n'))
59                 ci_len--;
60         cs_len = cseq_nr.len;
61         while (cs_len && ((c=cseq_nr.s[cs_len-1])==0 || c=='\r' || c=='\n'))
62                 cs_len--;
63 */
64         trim_len( ci_len, ci, call_id );
65         trim_len( cs_len, cs, cseq_nr );
66
67         /* run the cycle from the end ... we are interested in the
68            most-right digits ... and just take the %10 value of it
69         */
70         third=(ci_len-1)/3;
71         for ( i=ci_len-1, j=2*third, k=third;
72                 k>0 ; i--, j--, k-- ) {
73                 hash_code+=crc_16_tab[(unsigned char)(*(ci+i)) /*+7*/ ]+
74                         ccitt_tab[*(ci+k)+63]+
75                         ccitt_tab[*(ci+j)+13];
76         }
77         for( i=0 ; i<cs_len ; i++ )
78                 //hash_code+=crc_32_tab[(cseq_nr.s[i]+hash_code)%243];
79                 hash_code+=ccitt_tab[*(cs+i)+123];
80
81         /* hash_code conditioning */
82 #ifdef _BUG
83         /* not flat ... % 111b has shorter period than
84        & 111b by one and results in different distribution;
85            ( 7 % 7 == 0, 7 %7 == 1 )
86            % is used as a part of the hash function too, not only
87            for rounding; & is not flat; whoever comes up with
88            a nicer flat hash function which does not take
89            costly division is welcome; feel free to verify
90            distribution using hashtest()
91     */
92         hash_code &= (TABLE_ENTRIES-1); /* TABLE_ENTRIES = 2^k */
93 #endif
94         hash_code=hash_code%(TABLE_ENTRIES-1)+1;
95         return hash_code;
96 }
97
98
99
100 int new_hash2( str call_id, str cseq_nr )
101 {
102 #define h_inc h+=v^(v>>3)
103         char* p;
104         register unsigned v;
105         register unsigned h;
106         
107         h=0;
108         
109         
110         for (p=call_id.s; p<=(call_id.s+call_id.len-4); p+=4){
111                 v=(*p<<24)+(p[1]<<16)+(p[2]<<8)+p[3];
112                 h_inc;
113         }
114         v=0;
115         for (;p<(call_id.s+call_id.len); p++){ v<<=8; v+=*p;}
116         h_inc;
117         
118         for (p=cseq_nr.s; p<=(cseq_nr.s+cseq_nr.len-4); p+=4){
119                 v=(*p<<24)+(p[1]<<16)+(p[2]<<8)+p[3];
120                 h_inc;
121         }
122         v=0;
123         for (;p<(cseq_nr.s+cseq_nr.len); p++){ v<<=8; v+=*p;}
124         h_inc;
125         
126         h=((h)+(h>>11))+((h>>13)+(h>>23));
127         return (h)&(TABLE_ENTRIES-1);
128 }
129
130
131
132 void hashtest_cycle( int hits[TABLE_ENTRIES+5], char *ip )
133 {
134         long int i,j,k, l;
135         int  hashv;
136         static char buf1[1024];
137         static char buf2[1024];
138         str call_id; 
139         str cseq;
140
141         call_id.s=buf1;
142         cseq.s=buf2;
143
144         for (i=987654328;i<987654328+10;i++)
145                 for (j=85296341;j<85296341+10;j++)
146                         for (k=987654;k<=987654+10;k++)
147                                 for (l=101;l<201;l++) {
148                                         call_id.len=sprintf( buf1, "%d-%d-%d@%s",(int)i,(int)j,
149                                                 (int)k, ip );
150                                         cseq.len=sprintf( buf2, "%d", (int)l );
151                                         /* printf("%s\t%s\n", buf1, buf2 ); */
152                                         hashv=hash( call_id, cseq );
153                                         hits[ hashv ]++;
154                                 }
155 }
156
157 void hashtest()
158 {
159         int hits[TABLE_ENTRIES+5];
160         int i;
161
162         memset( hits, 0, sizeof hits );
163         hashtest_cycle( hits, "192.168.99.100" );
164         hashtest_cycle( hits, "172.168.99.100" );
165         hashtest_cycle( hits, "142.168.99.100" );
166         for (i=0; i<TABLE_ENTRIES+5; i++)
167                 printf("[%d. %d]\n", i, hits[i] );
168         exit(0);
169 }
170