forked from nothings/stb
-
Notifications
You must be signed in to change notification settings - Fork 0
/
stb_ds.h
1860 lines (1594 loc) · 65.9 KB
/
stb_ds.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
/* stb_ds.h - v0.64 - public domain data structures - Sean Barrett 2019
This is a single-header-file library that provides easy-to-use
dynamic arrays and hash tables for C (also works in C++).
For a gentle introduction:
http://nothings.org/stb_ds
To use this library, do this in *one* C or C++ file:
#define STB_DS_IMPLEMENTATION
#include "stb_ds.h"
TABLE OF CONTENTS
Table of Contents
Compile-time options
License
Documentation
Notes
Notes - Dynamic arrays
Notes - Hash maps
Credits
COMPILE-TIME OPTIONS
#define STBDS_NO_SHORT_NAMES
This flag needs to be set globally.
By default stb_ds exposes shorter function names that are not qualified
with the "stbds_" prefix. If these names conflict with the names in your
code, define this flag.
#define STBDS_SIPHASH_2_4
This flag only needs to be set in the file containing #define STB_DS_IMPLEMENTATION.
By default stb_ds.h hashes using a weaker variant of SipHash and a custom hash for
4- and 8-byte keys. On 64-bit platforms, you can define the above flag to force
stb_ds.h to use specification-compliant SipHash-2-4 for all keys. Doing so makes
hash table insertion about 20% slower on 4- and 8-byte keys, 5% slower on
64-byte keys, and 10% slower on 256-byte keys on my test computer.
#define STBDS_REALLOC(context,ptr,size) better_realloc
#define STBDS_FREE(context,ptr) better_free
These defines only need to be set in the file containing #define STB_DS_IMPLEMENTATION.
By default stb_ds uses stdlib realloc() and free() for memory management. You can
substitute your own functions instead by defining these symbols. You must either
define both, or neither. Note that at the moment, 'context' will always be NULL.
@TODO add an array/hash initialization function that takes a memory context pointer.
#define STBDS_UNIT_TESTS
Defines a function stbds_unit_tests() that checks the functioning of the data structures.
Note that on older versions of gcc (e.g. 5.x.x) you may need to build with '-std=c++0x'
(or equivalentally '-std=c++11') when using anonymous structures as seen on the web
page or in STBDS_UNIT_TESTS.
LICENSE
Placed in the public domain and also MIT licensed.
See end of file for detailed license information.
DOCUMENTATION
Dynamic Arrays
Non-function interface:
Declare an empty dynamic array of type T
T* foo = NULL;
Access the i'th item of a dynamic array 'foo' of type T, T* foo:
foo[i]
Functions (actually macros)
arrfree:
void arrfree(T*);
Frees the array.
arrlen:
ptrdiff_t arrlen(T*);
Returns the number of elements in the array.
arrlenu:
size_t arrlenu(T*);
Returns the number of elements in the array as an unsigned type.
arrpop:
T arrpop(T* a)
Removes the final element of the array and returns it.
arrput:
T arrput(T* a, T b);
Appends the item b to the end of array a. Returns b.
arrins:
T arrins(T* a, int p, T b);
Inserts the item b into the middle of array a, into a[p],
moving the rest of the array over. Returns b.
arrinsn:
void arrins(T* a, int p, int n);
Inserts n uninitialized items into array a starting at a[p],
moving the rest of the array over.
arraddn:
T* arraddn(T* a, int n)
Appends n uninitialized items onto array at the end.
Returns a pointer to the first uninitialized item added.
arrdel:
void arrdel(T* a, int p);
Deletes the element at a[p], moving the rest of the array over.
arrdeln:
void arrdel(T* a, int p, int n);
Deletes n elements starting at a[p], moving the rest of the array over.
arrdelswap:
void arrdelswap(T* a, int p);
Deletes the element at a[p], replacing it with the element from
the end of the array. O(1) performance.
arrsetlen:
void arrsetlen(T* a, int n);
Changes the length of the array to n. Allocates uninitialized
slots at the end if necessary.
arrsetcap:
size_t arrsetcap(T* a, int n);
Sets the length of allocated storage to at least n. It will not
change the length of the array.
arrcap:
size_t arrcap(T* a);
Returns the number of total elements the array can contain without
needing to be reallocated.
Hash maps & String hash maps
Given T is a structure type: struct { TK key; TV value; }. Note that some
functions do not require TV value and can have other fields. For string
hash maps, TK must be 'char *'.
Special interface:
stbds_rand_seed:
void stbds_rand_seed(size_t seed);
For security against adversarially chosen data, you should seed the
library with a strong random number. Or at least seed it with time().
stbds_hash_string:
size_t stbds_hash_string(char *str, size_t seed);
Returns a hash value for a string.
stbds_hash_bytes:
size_t stbds_hash_bytes(void *p, size_t len, size_t seed);
These functions hash an arbitrary number of bytes. The function
uses a custom hash for 4- and 8-byte data, and a weakened version
of SipHash for everything else. On 64-bit platforms you can get
specification-compliant SipHash-2-4 on all data by defining
STBDS_SIPHASH_2_4, at a significant cost in speed.
Non-function interface:
Declare an empty hash map of type T
T* foo = NULL;
Access the i'th entry in a hash table T* foo:
foo[i]
Function interface (actually macros):
hmfree
shfree
void hmfree(T*);
void shfree(T*);
Frees the hashmap and sets the pointer to NULL.
hmlen
shlen
ptrdiff_t hmlen(T*)
ptrdiff_t shlen(T*)
Returns the number of elements in the hashmap.
hmlenu
shlenu
size_t hmlenu(T*)
size_t shlenu(T*)
Returns the number of elements in the hashmap.
hmgeti
shgeti
hmgeti_ts
ptrdiff_t hmgeti(T*, TK key)
ptrdiff_t shgeti(T*, char* key)
ptrdiff_t hmgeti_ts(T*, TK key, ptrdiff_t tempvar)
Returns the index in the hashmap which has the key 'key', or -1
if the key is not present.
hmget
hmget_ts
shget
TV hmget(T*, TK key)
TV shget(T*, char* key)
TV hmget_ts(T*, TK key, ptrdiff_t tempvar)
Returns the value corresponding to 'key' in the hashmap.
The structure must have a 'value' field
hmgets
shgets
T hmgets(T*, TK key)
T shgets(T*, char* key)
Returns the structure corresponding to 'key' in the hashmap.
hmgetp
shgetp
hmgetp_ts
hmgetp_null
shgetp_null
T* hmgetp(T*, TK key)
T* shgetp(T*, char* key)
T* hmgetp_ts(T*, TK key, ptrdiff_t tempvar)
T* hmgetp_null(T*, TK key)
T* shgetp_null(T*, char *key)
Returns a pointer to the structure corresponding to 'key' in
the hashmap. Functions ending in "_null" return NULL if the key
is not present in the hashmap; the others return a pointer to a
structure holding the default value (but not the searched-for key).
hmdefault
shdefault
TV hmdefault(T*, TV value)
TV shdefault(T*, TV value)
Sets the default value for the hashmap, the value which will be
returned by hmget/shget if the key is not present.
hmdefaults
shdefaults
TV hmdefaults(T*, T item)
TV shdefaults(T*, T item)
Sets the default struct for the hashmap, the contents which will be
returned by hmgets/shgets if the key is not present.
hmput
shput
TV hmput(T*, TK key, TV value)
TV shput(T*, char* key, TV value)
Inserts a <key,value> pair into the hashmap. If the key is already
present in the hashmap, updates its value.
hmputs
shputs
T hmputs(T*, T item)
T shputs(T*, T item)
Inserts a struct with T.key into the hashmap. If the struct is already
present in the hashmap, updates it.
hmdel
shdel
int hmdel(T*, TK key)
int shdel(T*, char* key)
If 'key' is in the hashmap, deletes its entry and returns 1.
Otherwise returns 0.
Function interface (actually macros) for strings only:
sh_new_strdup
void sh_new_strdup(T*);
Overwrites the existing pointer with a newly allocated
string hashmap which will automatically allocate and free
each string key using realloc/free
sh_new_arena
void sh_new_arena(T*);
Overwrites the existing pointer with a newly allocated
string hashmap which will automatically allocate each string
key to a string arena. Every string key ever used by this
hash table remains in the arena until the arena is freed.
Additionally, any key which is deleted and reinserted will
be allocated multiple times in the string arena.
NOTES
* These data structures are realloc'd when they grow, and the macro
"functions" write to the provided pointer. This means: (a) the pointer
must be an lvalue, and (b) the pointer to the data structure is not
stable, and you must maintain it the same as you would a realloc'd
pointer. For example, if you pass a pointer to a dynamic array to a
function which updates it, the function must return back the new
pointer to the caller. This is the price of trying to do this in C.
* The following are the only functions that are thread-safe on a single data
structure, i.e. can be run in multiple threads simultaneously on the same
data structure
hmlen shlen
hmlenu shlenu
hmget_ts shget_ts
hmgeti_ts shgeti_ts
hmgets_ts shgets_ts
* You iterate over the contents of a dynamic array and a hashmap in exactly
the same way, using arrlen/hmlen/shlen:
for (i=0; i < arrlen(foo); ++i)
... foo[i] ...
* All operations except arrins/arrdel are O(1) amortized, but individual
operations can be slow, so these data structures may not be suitable
for real time use. Dynamic arrays double in capacity as needed, so
elements are copied an average of once. Hash tables double/halve
their size as needed, with appropriate hysteresis to maintain O(1)
performance.
NOTES - DYNAMIC ARRAY
* If you know how long a dynamic array is going to be in advance, you can avoid
extra memory allocations by using arrsetlen to allocate it to that length in
advance and use foo[n] while filling it out, or arrsetcap to allocate the memory
for that length and use arrput/arrpush as normal.
* Unlike some other versions of the dynamic array, this version should
be safe to use with strict-aliasing optimizations.
NOTES - HASH MAP
* For compilers other than GCC and clang (e.g. Visual Studio), for hmput/hmget/hmdel
and variants, the key must be an lvalue (so the macro can take the address of it).
Extensions are used that eliminate this requirement if you're using C99 and later
in GCC or clang, or if you're using C++ in GCC. But note that this can make your
code less portable.
* To test for presence of a key in a hashmap, just do 'hmgeti(foo,key) >= 0'.
* The iteration order of your data in the hashmap is determined solely by the
order of insertions and deletions. In particular, if you never delete, new
keys are always added at the end of the array. This will be consistent
across all platforms and versions of the library. However, you should not
attempt to serialize the internal hash table, as the hash is not consistent
between different platforms, and may change with future versions of the library.
* Use sh_new_arena() for string hashmaps that you never delete from. Initialize
with NULL if you're managing the memory for your strings, or your strings are
never freed (at least until the hashmap is freed). Otherwise, use sh_new_strdup().
@TODO: make an arena variant that garbage collects the strings with a trivial
copy collector into a new arena whenever the table shrinks / rebuilds. Since
current arena recommendation is to only use arena if it never deletes, then
this can just replace current arena implementation.
* If adversarial input is a serious concern and you're on a 64-bit platform,
enable STBDS_SIPHASH_2_4 (see the 'Compile-time options' section), and pass
a strong random number to stbds_rand_seed.
* The default value for the hash table is stored in foo[-1], so if you
use code like 'hmget(T,k)->value = 5' you can accidentally overwrite
the value stored by hmdefault if 'k' is not present.
CREDITS
Sean Barrett -- library, idea for dynamic array API/implementation
Per Vognsen -- idea for hash table API/implementation
Rafael Sachetto -- arrpop()
Bugfixes:
Andy Durdin
Shane Liesegang
Vinh Truong
*/
#ifdef STBDS_UNIT_TESTS
#define _CRT_SECURE_NO_WARNINGS
#endif
#ifndef INCLUDE_STB_DS_H
#define INCLUDE_STB_DS_H
#include <stddef.h>
#include <string.h>
#ifndef STBDS_NO_SHORT_NAMES
#define arrlen stbds_arrlen
#define arrlenu stbds_arrlenu
#define arrput stbds_arrput
#define arrpush stbds_arrput
#define arrpop stbds_arrpop
#define arrfree stbds_arrfree
#define arraddn stbds_arraddn
#define arrsetlen stbds_arrsetlen
#define arrlast stbds_arrlast
#define arrins stbds_arrins
#define arrinsn stbds_arrinsn
#define arrdel stbds_arrdel
#define arrdeln stbds_arrdeln
#define arrdelswap stbds_arrdelswap
#define arrcap stbds_arrcap
#define arrsetcap stbds_arrsetcap
#define hmput stbds_hmput
#define hmputs stbds_hmputs
#define hmget stbds_hmget
#define hmget_ts stbds_hmget_ts
#define hmgets stbds_hmgets
#define hmgetp stbds_hmgetp
#define hmgetp_ts stbds_hmgetp_ts
#define hmgetp_null stbds_hmgetp_null
#define hmgeti stbds_hmgeti
#define hmgeti_ts stbds_hmgeti_ts
#define hmdel stbds_hmdel
#define hmlen stbds_hmlen
#define hmlenu stbds_hmlenu
#define hmfree stbds_hmfree
#define hmdefault stbds_hmdefault
#define hmdefaults stbds_hmdefaults
#define shput stbds_shput
#define shputi stbds_shputi
#define shputs stbds_shputs
#define shget stbds_shget
#define shgeti stbds_shgeti
#define shgets stbds_shgets
#define shgetp stbds_shgetp
#define shgetp_null stbds_shgetp_null
#define shdel stbds_shdel
#define shlen stbds_shlen
#define shlenu stbds_shlenu
#define shfree stbds_shfree
#define shdefault stbds_shdefault
#define shdefaults stbds_shdefaults
#define sh_new_arena stbds_sh_new_arena
#define sh_new_strdup stbds_sh_new_strdup
#define stralloc stbds_stralloc
#define strreset stbds_strreset
#endif
#if defined(STBDS_REALLOC) && !defined(STBDS_FREE) || !defined(STBDS_REALLOC) && defined(STBDS_FREE)
#error "You must define both STBDS_REALLOC and STBDS_FREE, or neither."
#endif
#if !defined(STBDS_REALLOC) && !defined(STBDS_FREE)
#include <stdlib.h>
#define STBDS_REALLOC(c,p,s) realloc(p,s)
#define STBDS_FREE(c,p) free(p)
#endif
#ifdef __cplusplus
extern "C" {
#endif
// for security against attackers, seed the library with a random number, at least time() but stronger is better
extern void stbds_rand_seed(size_t seed);
// these are the hash functions used internally if you want to test them or use them for other purposes
extern size_t stbds_hash_bytes(void *p, size_t len, size_t seed);
extern size_t stbds_hash_string(char *str, size_t seed);
// this is a simple string arena allocator, initialize with e.g. 'stbds_string_arena my_arena={0}'.
typedef struct stbds_string_arena stbds_string_arena;
extern char * stbds_stralloc(stbds_string_arena *a, char *str);
extern void stbds_strreset(stbds_string_arena *a);
// have to #define STBDS_UNIT_TESTS to call this
extern void stbds_unit_tests(void);
///////////////
//
// Everything below here is implementation details
//
extern void * stbds_arrgrowf(void *a, size_t elemsize, size_t addlen, size_t min_cap);
extern void stbds_hmfree_func(void *p, size_t elemsize);
extern void * stbds_hmget_key(void *a, size_t elemsize, void *key, size_t keysize, int mode);
extern void * stbds_hmget_key_ts(void *a, size_t elemsize, void *key, size_t keysize, ptrdiff_t *temp, int mode);
extern void * stbds_hmput_default(void *a, size_t elemsize);
extern void * stbds_hmput_key(void *a, size_t elemsize, void *key, size_t keysize, int mode);
extern void * stbds_hmdel_key(void *a, size_t elemsize, void *key, size_t keysize, size_t keyoffset, int mode);
extern void * stbds_shmode_func(size_t elemsize, int mode);
#ifdef __cplusplus
}
#endif
#if defined(__GNUC__) || defined(__clang__)
#define STBDS_HAS_TYPEOF
#ifdef __cplusplus
//#define STBDS_HAS_LITERAL_ARRAY // this is currently broken for clang
#endif
#endif
#if !defined(__cplusplus)
#if defined(__STDC_VERSION__) && __STDC_VERSION__ >= 199901L
#define STBDS_HAS_LITERAL_ARRAY
#endif
#endif
// this macro takes the address of the argument, but on gcc/clang can accept rvalues
#if defined(STBDS_HAS_LITERAL_ARRAY) && defined(STBDS_HAS_TYPEOF)
#if __clang__
#define STBDS_ADDRESSOF(typevar, value) ((__typeof__(typevar)[1]){value}) // literal array decays to pointer to value
#else
#define STBDS_ADDRESSOF(typevar, value) ((typeof(typevar)[1]){value}) // literal array decays to pointer to value
#endif
#else
#define STBDS_ADDRESSOF(typevar, value) &(value)
#endif
#define STBDS_OFFSETOF(var,field) ((char *) &(var)->field - (char *) (var))
#define stbds_header(t) ((stbds_array_header *) (t) - 1)
#define stbds_temp(t) stbds_header(t)->temp
#define stbds_temp_key(t) stbds_header(t)->temp_key
#define stbds_arrsetcap(a,n) (stbds_arrgrow(a,0,n))
#define stbds_arrsetlen(a,n) ((stbds_arrcap(a) < (size_t) (n) ? stbds_arrsetcap((a),(size_t)(n)),0 : 0), (a) ? stbds_header(a)->length = (size_t) (n) : 0)
#define stbds_arrcap(a) ((a) ? stbds_header(a)->capacity : 0)
#define stbds_arrlen(a) ((a) ? (ptrdiff_t) stbds_header(a)->length : 0)
#define stbds_arrlenu(a) ((a) ? stbds_header(a)->length : 0)
#define stbds_arrput(a,v) (stbds_arrmaybegrow(a,1), (a)[stbds_header(a)->length++] = (v))
#define stbds_arrpush stbds_arrput // synonym
#define stbds_arrpop(a) (stbds_header(a)->length--, (a)[stbds_header(a)->length])
#define stbds_arraddn(a,n) (stbds_arrmaybegrow(a,n), stbds_header(a)->length += (n), stbds_header(a)->length-(n))
#define stbds_arrlast(a) ((a)[stbds_header(a)->length-1])
#define stbds_arrfree(a) ((void) ((a) ? STBDS_FREE(NULL,stbds_header(a)) : (void)0), (a)=NULL)
#define stbds_arrdel(a,i) stbds_arrdeln(a,i,1)
#define stbds_arrdeln(a,i,n) (memmove(&(a)[i], &(a)[(i)+(n)], sizeof *(a) * (stbds_header(a)->length-(n)-(i))), stbds_header(a)->length -= (n))
#define stbds_arrdelswap(a,i) ((a)[i] = stbds_arrlast(a), stbds_header(a)->length -= 1)
#define stbds_arrinsn(a,i,n) (stbds_arraddn((a),(n)), memmove(&(a)[(i)+(n)], &(a)[i], sizeof *(a) * (stbds_header(a)->length-(n)-(i))))
#define stbds_arrins(a,i,v) (stbds_arrinsn((a),(i),1), (a)[i]=(v))
#define stbds_arrmaybegrow(a,n) ((!(a) || stbds_header(a)->length + (n) > stbds_header(a)->capacity) \
? (stbds_arrgrow(a,n,0),0) : 0)
#define stbds_arrgrow(a,b,c) ((a) = stbds_arrgrowf_wrapper((a), sizeof *(a), (b), (c)))
#define stbds_hmput(t, k, v) \
((t) = stbds_hmput_key_wrapper((t), sizeof *(t), (void*) STBDS_ADDRESSOF((t)->key, (k)), sizeof (t)->key, 0), \
(t)[stbds_temp((t)-1)].key = (k), \
(t)[stbds_temp((t)-1)].value = (v))
#define stbds_hmputs(t, s) \
((t) = stbds_hmput_key_wrapper((t), sizeof *(t), &(s).key, sizeof (s).key, STBDS_HM_BINARY), \
(t)[stbds_temp((t)-1)] = (s))
#define stbds_hmgeti(t,k) \
((t) = stbds_hmget_key_wrapper((t), sizeof *(t), (void*) STBDS_ADDRESSOF((t)->key, (k)), sizeof (t)->key, STBDS_HM_BINARY), \
stbds_temp((t)-1))
#define stbds_hmgeti_ts(t,k,temp) \
((t) = stbds_hmget_key_ts_wrapper((t), sizeof *(t), (void*) STBDS_ADDRESSOF((t)->key, (k)), sizeof (t)->key, &(temp), STBDS_HM_BINARY), \
(temp))
#define stbds_hmgetp(t, k) \
((void) stbds_hmgeti(t,k), &(t)[stbds_temp((t)-1)])
#define stbds_hmgetp_ts(t, k, temp) \
((void) stbds_hmgeti_ts(t,k,temp), &(t)[temp])
#define stbds_hmdel(t,k) \
(((t) = stbds_hmdel_key_wrapper((t),sizeof *(t), (void*) STBDS_ADDRESSOF((t)->key, (k)), sizeof (t)->key, STBDS_OFFSETOF((t),key), STBDS_HM_BINARY)),(t)?stbds_temp((t)-1):0)
#define stbds_hmdefault(t, v) \
((t) = stbds_hmput_default_wrapper((t), sizeof *(t)), (t)[-1].value = (v))
#define stbds_hmdefaults(t, s) \
((t) = stbds_hmput_default_wrapper((t), sizeof *(t)), (t)[-1] = (s))
#define stbds_hmfree(p) \
((void) ((p) != NULL ? stbds_hmfree_func((p)-1,sizeof*(p)),0 : 0),(p)=NULL)
#define stbds_hmgets(t, k) (*stbds_hmgetp(t,k))
#define stbds_hmget(t, k) (stbds_hmgetp(t,k)->value)
#define stbds_hmget_ts(t, k, temp) (stbds_hmgetp_ts(t,k,temp)->value)
#define stbds_hmlen(t) ((t) ? (ptrdiff_t) stbds_header((t)-1)->length-1 : 0)
#define stbds_hmlenu(t) ((t) ? stbds_header((t)-1)->length-1 : 0)
#define stbds_hmgetp_null(t,k) (stbds_hmgeti(t,k) == -1 ? NULL : &(t)[stbds_temp(t)-1])
#define stbds_shput(t, k, v) \
((t) = stbds_hmput_key_wrapper((t), sizeof *(t), (void*) (k), sizeof (t)->key, STBDS_HM_STRING), \
(t)[stbds_temp((t)-1)].value = (v))
#define stbds_shputi(t, k, v) \
((t) = stbds_hmput_key_wrapper((t), sizeof *(t), (void*) (k), sizeof (t)->key, STBDS_HM_STRING), \
(t)[stbds_temp((t)-1)].value = (v), stbds_temp((t)-1))
#define stbds_shputs(t, s) \
((t) = stbds_hmput_key_wrapper((t), sizeof *(t), (void*) (s).key, sizeof (s).key, STBDS_HM_STRING), \
(t)[stbds_temp((t)-1)] = (s), \
(t)[stbds_temp((t)-1)].key = stbds_temp_key((t)-1))
#define stbds_pshput(t, p) \
((t) = stbds_hmput_key_wrapper((t), sizeof *(t), (void*) (p)->key, sizeof (p)->key, STBDS_HM_PTR_TO_STRING), \
(t)[stbds_temp((t)-1)] = (p))
#define stbds_shgeti(t,k) \
((t) = stbds_hmget_key_wrapper((t), sizeof *(t), (void*) (k), sizeof (t)->key, STBDS_HM_STRING), \
stbds_temp((t)-1))
#define stbds_pshgeti(t,k) \
((t) = stbds_hmget_key_wrapper((t), sizeof *(t), (void*) (k), sizeof (*(t))->key, STBDS_HM_PTR_TO_STRING), \
stbds_temp((t)-1))
#define stbds_shgetp(t, k) \
((void) stbds_shgeti(t,k), &(t)[stbds_temp((t)-1)])
#define stbds_pshget(t, k) \
((void) stbds_pshgeti(t,k), (t)[stbds_temp((t)-1)])
#define stbds_shdel(t,k) \
(((t) = stbds_hmdel_key_wrapper((t),sizeof *(t), (void*) (k), sizeof (t)->key, STBDS_OFFSETOF((t),key), STBDS_HM_STRING)),(t)?stbds_temp((t)-1):0)
#define stbds_pshdel(t,k) \
(((t) = stbds_hmdel_key_wrapper((t),sizeof *(t), (void*) (k), sizeof (*(t))->key, STBDS_OFFSETOF(*(t),key), STBDS_HM_PTR_TO_STRING)),(t)?stbds_temp((t)-1):0)
#define stbds_sh_new_arena(t) \
((t) = stbds_shmode_func_wrapper(t, sizeof *(t), STBDS_SH_ARENA))
#define stbds_sh_new_strdup(t) \
((t) = stbds_shmode_func_wrapper(t, sizeof *(t), STBDS_SH_STRDUP))
#define stbds_shdefault(t, v) stbds_hmdefault(t,v)
#define stbds_shdefaults(t, s) stbds_hmdefaults(t,s)
#define stbds_shfree stbds_hmfree
#define stbds_shlenu stbds_hmlenu
#define stbds_shgets(t, k) (*stbds_shgetp(t,k))
#define stbds_shget(t, k) (stbds_shgetp(t,k)->value)
#define stbds_shgetp_null(t,k) (stbds_shgeti(t,k) == -1 ? NULL : &(t)[stbds_temp(t)-1])
#define stbds_shlen stbds_hmlen
typedef struct
{
size_t length;
size_t capacity;
void * hash_table;
ptrdiff_t temp;
char * temp_key;
} stbds_array_header;
typedef struct stbds_string_block
{
struct stbds_string_block *next;
char storage[8];
} stbds_string_block;
struct stbds_string_arena
{
stbds_string_block *storage;
size_t remaining;
unsigned char block;
unsigned char mode; // this isn't used by the string arena itself
};
#define STBDS_HM_BINARY 0
#define STBDS_HM_STRING 1
enum
{
STBDS_SH_NONE,
STBDS_SH_DEFAULT,
STBDS_SH_STRDUP,
STBDS_SH_ARENA
};
#ifdef __cplusplus
// in C we use implicit assignment from these void*-returning functions to T*.
// in C++ these templates make the same code work
template<class T> static T * stbds_arrgrowf_wrapper(T *a, size_t elemsize, size_t addlen, size_t min_cap) {
return (T*)stbds_arrgrowf((void *)a, elemsize, addlen, min_cap);
}
template<class T> static T * stbds_hmget_key_wrapper(T *a, size_t elemsize, void *key, size_t keysize, int mode) {
return (T*)stbds_hmget_key((void*)a, elemsize, key, keysize, mode);
}
template<class T> static T * stbds_hmget_key_ts_wrapper(T *a, size_t elemsize, void *key, size_t keysize, ptrdiff_t *temp, int mode) {
return (T*)stbds_hmget_key_ts((void*)a, elemsize, key, keysize, temp, mode);
}
template<class T> static T * stbds_hmput_default_wrapper(T *a, size_t elemsize) {
return (T*)stbds_hmput_default((void *)a, elemsize);
}
template<class T> static T * stbds_hmput_key_wrapper(T *a, size_t elemsize, void *key, size_t keysize, int mode) {
return (T*)stbds_hmput_key((void*)a, elemsize, key, keysize, mode);
}
template<class T> static T * stbds_hmdel_key_wrapper(T *a, size_t elemsize, void *key, size_t keysize, size_t keyoffset, int mode){
return (T*)stbds_hmdel_key((void*)a, elemsize, key, keysize, keyoffset, mode);
}
template<class T> static T * stbds_shmode_func_wrapper(T *, size_t elemsize, int mode) {
return (T*)stbds_shmode_func(elemsize, mode);
}
#else
#define stbds_arrgrowf_wrapper stbds_arrgrowf
#define stbds_hmget_key_wrapper stbds_hmget_key
#define stbds_hmget_key_ts_wrapper stbds_hmget_key_ts
#define stbds_hmput_default_wrapper stbds_hmput_default
#define stbds_hmput_key_wrapper stbds_hmput_key
#define stbds_hmdel_key_wrapper stbds_hmdel_key
#define stbds_shmode_func_wrapper(t,e,m) stbds_shmode_func(e,m)
#endif
#endif // INCLUDE_STB_DS_H
//////////////////////////////////////////////////////////////////////////////
//
// IMPLEMENTATION
//
#ifdef STB_DS_IMPLEMENTATION
#include <assert.h>
#include <string.h>
#ifndef STBDS_ASSERT
#define STBDS_ASSERT_WAS_UNDEFINED
#define STBDS_ASSERT(x) ((void) 0)
#endif
#ifdef STBDS_STATISTICS
#define STBDS_STATS(x) x
size_t stbds_array_grow;
size_t stbds_hash_grow;
size_t stbds_hash_shrink;
size_t stbds_hash_rebuild;
size_t stbds_hash_probes;
size_t stbds_hash_alloc;
size_t stbds_rehash_probes;
size_t stbds_rehash_items;
#else
#define STBDS_STATS(x)
#endif
//
// stbds_arr implementation
//
//int *prev_allocs[65536];
//int num_prev;
void *stbds_arrgrowf(void *a, size_t elemsize, size_t addlen, size_t min_cap)
{
void *b;
size_t min_len = stbds_arrlen(a) + addlen;
// compute the minimum capacity needed
if (min_len > min_cap)
min_cap = min_len;
if (min_cap <= stbds_arrcap(a))
return a;
// increase needed capacity to guarantee O(1) amortized
if (min_cap < 2 * stbds_arrcap(a))
min_cap = 2 * stbds_arrcap(a);
else if (min_cap < 4)
min_cap = 4;
//if (num_prev < 65536) if (a) prev_allocs[num_prev++] = (int *) ((char *) a+1);
//if (num_prev == 2201)
// num_prev = num_prev;
b = STBDS_REALLOC(NULL, (a) ? stbds_header(a) : 0, elemsize * min_cap + sizeof(stbds_array_header));
//if (num_prev < 65536) prev_allocs[num_prev++] = (int *) (char *) b;
b = (char *) b + sizeof(stbds_array_header);
if (a == NULL) {
stbds_header(b)->length = 0;
stbds_header(b)->hash_table = 0;
} else {
STBDS_STATS(++stbds_array_grow);
}
stbds_header(b)->capacity = min_cap;
return b;
}
//
// stbds_hm hash table implementation
//
#ifdef STBDS_INTERNAL_SMALL_BUCKET
#define STBDS_BUCKET_LENGTH 4
#else
#define STBDS_BUCKET_LENGTH 8
#endif
#define STBDS_BUCKET_SHIFT (STBDS_BUCKET_LENGTH == 8 ? 3 : 2)
#define STBDS_BUCKET_MASK (STBDS_BUCKET_LENGTH-1)
#define STBDS_CACHE_LINE_SIZE 64
#define STBDS_ALIGN_FWD(n,a) (((n) + (a) - 1) & ~((a)-1))
typedef struct
{
size_t hash [STBDS_BUCKET_LENGTH];
ptrdiff_t index[STBDS_BUCKET_LENGTH];
} stbds_hash_bucket; // in 32-bit, this is one 64-byte cache line; in 64-bit, each array is one 64-byte cache line
typedef struct
{
size_t slot_count;
size_t used_count;
size_t used_count_threshold;
size_t used_count_shrink_threshold;
size_t tombstone_count;
size_t tombstone_count_threshold;
size_t seed;
size_t slot_count_log2;
stbds_string_arena string;
stbds_hash_bucket *storage; // not a separate allocation, just 64-byte aligned storage after this struct
} stbds_hash_index;
#define STBDS_INDEX_EMPTY -1
#define STBDS_INDEX_DELETED -2
#define STBDS_INDEX_IN_USE(x) ((x) >= 0)
#define STBDS_HASH_EMPTY 0
#define STBDS_HASH_DELETED 1
static size_t stbds_hash_seed=0x31415926;
void stbds_rand_seed(size_t seed)
{
stbds_hash_seed = seed;
}
#define stbds_load_32_or_64(var, temp, v32, v64_hi, v64_lo) \
temp = v64_lo ^ v32, temp <<= 16, temp <<= 16, temp >>= 16, temp >>= 16, /* discard if 32-bit */ \
var = v64_hi, var <<= 16, var <<= 16, /* discard if 32-bit */ \
var ^= temp ^ v32
#define STBDS_SIZE_T_BITS ((sizeof (size_t)) * 8)
static size_t stbds_probe_position(size_t hash, size_t slot_count, size_t slot_log2)
{
size_t pos;
pos = hash & (slot_count-1);
#ifdef STBDS_INTERNAL_BUCKET_START
pos &= ~STBDS_BUCKET_MASK;
#endif
return pos;
}
static size_t stbds_log2(size_t slot_count)
{
size_t n=0;
while (slot_count > 1) {
slot_count >>= 1;
++n;
}
return n;
}
static stbds_hash_index *stbds_make_hash_index(size_t slot_count, stbds_hash_index *ot)
{
stbds_hash_index *t;
t = (stbds_hash_index *) STBDS_REALLOC(NULL,0,(slot_count >> STBDS_BUCKET_SHIFT) * sizeof(stbds_hash_bucket) + sizeof(stbds_hash_index) + STBDS_CACHE_LINE_SIZE-1);
t->storage = (stbds_hash_bucket *) STBDS_ALIGN_FWD((size_t) (t+1), STBDS_CACHE_LINE_SIZE);
t->slot_count = slot_count;
t->slot_count_log2 = stbds_log2(slot_count);
t->tombstone_count = 0;
t->used_count = 0;
#if 0 // A1
t->used_count_threshold = slot_count*12/16; // if 12/16th of table is occupied, grow
t->tombstone_count_threshold = slot_count* 2/16; // if tombstones are 2/16th of table, rebuild
t->used_count_shrink_threshold = slot_count* 4/16; // if table is only 4/16th full, shrink
#elif 1 // A2
//t->used_count_threshold = slot_count*12/16; // if 12/16th of table is occupied, grow
//t->tombstone_count_threshold = slot_count* 3/16; // if tombstones are 3/16th of table, rebuild
//t->used_count_shrink_threshold = slot_count* 4/16; // if table is only 4/16th full, shrink
// compute without overflowing
t->used_count_threshold = slot_count - (slot_count>>2);
t->tombstone_count_threshold = (slot_count>>3) + (slot_count>>4);
t->used_count_shrink_threshold = slot_count >> 2;
#elif 0 // B1
t->used_count_threshold = slot_count*13/16; // if 13/16th of table is occupied, grow
t->tombstone_count_threshold = slot_count* 2/16; // if tombstones are 2/16th of table, rebuild
t->used_count_shrink_threshold = slot_count* 5/16; // if table is only 5/16th full, shrink
#else // C1
t->used_count_threshold = slot_count*14/16; // if 14/16th of table is occupied, grow
t->tombstone_count_threshold = slot_count* 2/16; // if tombstones are 2/16th of table, rebuild
t->used_count_shrink_threshold = slot_count* 6/16; // if table is only 6/16th full, shrink
#endif
// Following statistics were measured on a Core i7-6700 @ 4.00Ghz, compiled with clang 7.0.1 -O2
// Note that the larger tables have high variance as they were run fewer times
// A1 A2 B1 C1
// 0.10ms : 0.10ms : 0.10ms : 0.11ms : 2,000 inserts creating 2K table
// 0.96ms : 0.95ms : 0.97ms : 1.04ms : 20,000 inserts creating 20K table
// 14.48ms : 14.46ms : 10.63ms : 11.00ms : 200,000 inserts creating 200K table
// 195.74ms : 196.35ms : 203.69ms : 214.92ms : 2,000,000 inserts creating 2M table
// 2193.88ms : 2209.22ms : 2285.54ms : 2437.17ms : 20,000,000 inserts creating 20M table
// 65.27ms : 53.77ms : 65.33ms : 65.47ms : 500,000 inserts & deletes in 2K table
// 72.78ms : 62.45ms : 71.95ms : 72.85ms : 500,000 inserts & deletes in 20K table
// 89.47ms : 77.72ms : 96.49ms : 96.75ms : 500,000 inserts & deletes in 200K table
// 97.58ms : 98.14ms : 97.18ms : 97.53ms : 500,000 inserts & deletes in 2M table
// 118.61ms : 119.62ms : 120.16ms : 118.86ms : 500,000 inserts & deletes in 20M table
// 192.11ms : 194.39ms : 196.38ms : 195.73ms : 500,000 inserts & deletes in 200M table
if (slot_count <= STBDS_BUCKET_LENGTH)
t->used_count_shrink_threshold = 0;
// to avoid infinite loop, we need to guarantee that at least one slot is empty and will terminate probes
STBDS_ASSERT(t->used_count_threshold + t->tombstone_count_threshold < t->slot_count);
STBDS_STATS(++stbds_hash_alloc);
if (ot) {
t->string = ot->string;
// reuse old seed so we can reuse old hashes so below "copy out old data" doesn't do any hashing
t->seed = ot->seed;
} else {
size_t a,b,temp;
memset(&t->string, 0, sizeof(t->string));
t->seed = stbds_hash_seed;
// LCG
// in 32-bit, a = 2147001325 b = 715136305
// in 64-bit, a = 2862933555777941757 b = 3037000493
stbds_load_32_or_64(a,temp, 2147001325, 0x27bb2ee6, 0x87b0b0fd);
stbds_load_32_or_64(b,temp, 715136305, 0, 0xb504f32d);
stbds_hash_seed = stbds_hash_seed * a + b;
}
{
size_t i,j;
for (i=0; i < slot_count >> STBDS_BUCKET_SHIFT; ++i) {
stbds_hash_bucket *b = &t->storage[i];
for (j=0; j < STBDS_BUCKET_LENGTH; ++j)
b->hash[j] = STBDS_HASH_EMPTY;
for (j=0; j < STBDS_BUCKET_LENGTH; ++j)
b->index[j] = STBDS_INDEX_EMPTY;
}
}
// copy out the old data, if any
if (ot) {
size_t i,j;
t->used_count = ot->used_count;
for (i=0; i < ot->slot_count >> STBDS_BUCKET_SHIFT; ++i) {
stbds_hash_bucket *ob = &ot->storage[i];
for (j=0; j < STBDS_BUCKET_LENGTH; ++j) {
if (STBDS_INDEX_IN_USE(ob->index[j])) {
size_t hash = ob->hash[j];
size_t pos = stbds_probe_position(hash, t->slot_count, t->slot_count_log2);
size_t step = STBDS_BUCKET_LENGTH;
STBDS_STATS(++stbds_rehash_items);
for (;;) {
size_t limit,z;
stbds_hash_bucket *bucket;
bucket = &t->storage[pos >> STBDS_BUCKET_SHIFT];
STBDS_STATS(++stbds_rehash_probes);
for (z=pos & STBDS_BUCKET_MASK; z < STBDS_BUCKET_LENGTH; ++z) {
if (bucket->hash[z] == 0) {
bucket->hash[z] = hash;
bucket->index[z] = ob->index[j];
goto done;
}
}
limit = pos & STBDS_BUCKET_MASK;
for (z = 0; z < limit; ++z) {
if (bucket->hash[z] == 0) {
bucket->hash[z] = hash;
bucket->index[z] = ob->index[j];
goto done;
}
}
pos += step; // quadratic probing
step += STBDS_BUCKET_LENGTH;
pos &= (t->slot_count-1);
}
}
done:
;
}
}
}
return t;
}
#define STBDS_ROTATE_LEFT(val, n) (((val) << (n)) | ((val) >> (STBDS_SIZE_T_BITS - (n))))
#define STBDS_ROTATE_RIGHT(val, n) (((val) >> (n)) | ((val) << (STBDS_SIZE_T_BITS - (n))))
size_t stbds_hash_string(char *str, size_t seed)
{
size_t hash = seed;
while (*str)
hash = STBDS_ROTATE_LEFT(hash, 9) + (unsigned char) *str++;
// Thomas Wang 64-to-32 bit mix function, hopefully also works in 32 bits
hash ^= seed;
hash = (~hash) + (hash << 18);
hash ^= hash ^ STBDS_ROTATE_RIGHT(hash,31);
hash = hash * 21;
hash ^= hash ^ STBDS_ROTATE_RIGHT(hash,11);
hash += (hash << 6);
hash ^= STBDS_ROTATE_RIGHT(hash,22);
return hash+seed;
}
#ifdef STBDS_SIPHASH_2_4