forked from RcppCore/RcppParallel
-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathconcurrent_vector.h
More file actions
1374 lines (1177 loc) · 57.7 KB
/
concurrent_vector.h
File metadata and controls
1374 lines (1177 loc) · 57.7 KB
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
/*
Copyright (c) 2005-2017 Intel Corporation
Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at
http://www.apache.org/licenses/LICENSE-2.0
Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.
*/
#ifndef __TBB_concurrent_vector_H
#define __TBB_concurrent_vector_H
#include "tbb_stddef.h"
#include "tbb_exception.h"
#include "atomic.h"
#include "cache_aligned_allocator.h"
#include "blocked_range.h"
#include "tbb_machine.h"
#include "tbb_profiling.h"
#include <new>
#include <cstring> // for memset()
#include __TBB_STD_SWAP_HEADER
#if !TBB_USE_EXCEPTIONS && _MSC_VER
// Suppress "C++ exception handler used, but unwind semantics are not enabled" warning in STL headers
#pragma warning (push)
#pragma warning (disable: 4530)
#endif
#include <algorithm>
#include <iterator>
#if !TBB_USE_EXCEPTIONS && _MSC_VER
#pragma warning (pop)
#endif
#if _MSC_VER==1500 && !__INTEL_COMPILER
// VS2008/VC9 seems to have an issue; limits pull in math.h
#pragma warning( push )
#pragma warning( disable: 4985 )
#endif
#include <limits> /* std::numeric_limits */
#if _MSC_VER==1500 && !__INTEL_COMPILER
#pragma warning( pop )
#endif
#if __TBB_INITIALIZER_LISTS_PRESENT
#include <initializer_list>
#endif
#if defined(_MSC_VER) && !defined(__INTEL_COMPILER)
// Workaround for overzealous compiler warnings in /Wp64 mode
#pragma warning (push)
#if defined(_Wp64)
#pragma warning (disable: 4267)
#endif
#pragma warning (disable: 4127) //warning C4127: conditional expression is constant
#endif
namespace tbb {
template<typename T, class A = cache_aligned_allocator<T> >
class concurrent_vector;
//! @cond INTERNAL
namespace internal {
template<typename Container, typename Value>
class vector_iterator;
//! Bad allocation marker
static void *const vector_allocation_error_flag = reinterpret_cast<void*>(size_t(63));
//! Exception helper function
template<typename T>
void handle_unconstructed_elements(T* array, size_t n_of_elements){
std::memset( (void*) array, 0, n_of_elements * sizeof( T ) );
}
//! Base class of concurrent vector implementation.
/** @ingroup containers */
class concurrent_vector_base_v3 {
protected:
// Basic types declarations
typedef size_t segment_index_t;
typedef size_t size_type;
// Using enumerations due to Mac linking problems of static const variables
enum {
// Size constants
default_initial_segments = 1, // 2 initial items
//! Number of slots for segment pointers inside the class
pointers_per_short_table = 3, // to fit into 8 words of entire structure
pointers_per_long_table = sizeof(segment_index_t) * 8 // one segment per bit
};
public:
struct segment_not_used {};
struct segment_allocated {};
struct segment_allocation_failed {};
protected:
class segment_t;
class segment_value_t {
void* array;
private:
//TODO: More elegant way to grant access to selected functions _only_?
friend class segment_t;
explicit segment_value_t(void* an_array):array(an_array) {}
public:
friend bool operator==(segment_value_t const& lhs, segment_not_used ) { return lhs.array == 0;}
friend bool operator==(segment_value_t const& lhs, segment_allocated) { return lhs.array > internal::vector_allocation_error_flag;}
friend bool operator==(segment_value_t const& lhs, segment_allocation_failed) { return lhs.array == internal::vector_allocation_error_flag;}
template<typename argument_type>
friend bool operator!=(segment_value_t const& lhs, argument_type arg) { return ! (lhs == arg);}
template<typename T>
T* pointer() const { return static_cast<T*>(const_cast<void*>(array)); }
};
friend void enforce_segment_allocated(segment_value_t const& s, internal::exception_id exception = eid_bad_last_alloc){
if(s != segment_allocated()){
internal::throw_exception(exception);
}
}
// Segment pointer.
class segment_t {
atomic<void*> array;
public:
segment_t(){ store<relaxed>(segment_not_used());}
//Copy ctor and assignment operator are defined to ease using of stl algorithms.
//These algorithms usually not a synchronization point, so, semantic is
//intentionally relaxed here.
segment_t(segment_t const& rhs ){ array.store<relaxed>(rhs.array.load<relaxed>());}
void swap(segment_t & rhs ){
tbb::internal::swap<relaxed>(array, rhs.array);
}
segment_t& operator=(segment_t const& rhs ){
array.store<relaxed>(rhs.array.load<relaxed>());
return *this;
}
template<memory_semantics M>
segment_value_t load() const { return segment_value_t(array.load<M>());}
template<memory_semantics M>
void store(segment_not_used) {
array.store<M>(0);
}
template<memory_semantics M>
void store(segment_allocation_failed) {
__TBB_ASSERT(load<relaxed>() != segment_allocated(),"transition from \"allocated\" to \"allocation failed\" state looks non-logical");
array.store<M>(internal::vector_allocation_error_flag);
}
template<memory_semantics M>
void store(void* allocated_segment_pointer) __TBB_NOEXCEPT(true) {
__TBB_ASSERT(segment_value_t(allocated_segment_pointer) == segment_allocated(),
"other overloads of store should be used for marking segment as not_used or allocation_failed" );
array.store<M>(allocated_segment_pointer);
}
#if TBB_USE_ASSERT
~segment_t() {
__TBB_ASSERT(load<relaxed>() != segment_allocated(), "should have been freed by clear" );
}
#endif /* TBB_USE_ASSERT */
};
friend void swap(segment_t & , segment_t & ) __TBB_NOEXCEPT(true);
// Data fields
//! allocator function pointer
void* (*vector_allocator_ptr)(concurrent_vector_base_v3 &, size_t);
//! count of segments in the first block
atomic<size_type> my_first_block;
//! Requested size of vector
atomic<size_type> my_early_size;
//! Pointer to the segments table
atomic<segment_t*> my_segment;
//! embedded storage of segment pointers
segment_t my_storage[pointers_per_short_table];
// Methods
concurrent_vector_base_v3() {
//Here the semantic is intentionally relaxed.
//The reason this is next:
//Object that is in middle of construction (i.e. its constructor is not yet finished)
//cannot be used concurrently until the construction is finished.
//Thus to flag other threads that construction is finished, some synchronization with
//acquire-release semantic should be done by the (external) code that uses the vector.
//So, no need to do the synchronization inside the vector.
my_early_size.store<relaxed>(0);
my_first_block.store<relaxed>(0); // here is not default_initial_segments
my_segment.store<relaxed>(my_storage);
}
__TBB_EXPORTED_METHOD ~concurrent_vector_base_v3();
//these helpers methods use the fact that segments are allocated so
//that every segment size is a (increasing) power of 2.
//with one exception 0 segment has size of 2 as well segment 1;
//e.g. size of segment with index of 3 is 2^3=8;
static segment_index_t segment_index_of( size_type index ) {
return segment_index_t( __TBB_Log2( index|1 ) );
}
static segment_index_t segment_base( segment_index_t k ) {
return (segment_index_t(1)<<k & ~segment_index_t(1));
}
static inline segment_index_t segment_base_index_of( segment_index_t &index ) {
segment_index_t k = segment_index_of( index );
index -= segment_base(k);
return k;
}
static size_type segment_size( segment_index_t k ) {
return segment_index_t(1)<<k; // fake value for k==0
}
static bool is_first_element_in_segment(size_type element_index){
//check if element_index is a power of 2 that is at least 2.
//The idea is to detect if the iterator crosses a segment boundary,
//and 2 is the minimal index for which it's true
__TBB_ASSERT(element_index, "there should be no need to call "
"is_first_element_in_segment for 0th element" );
return is_power_of_two_at_least( element_index, 2 );
}
//! An operation on an n-element array starting at begin.
typedef void (__TBB_EXPORTED_FUNC *internal_array_op1)(void* begin, size_type n );
//! An operation on n-element destination array and n-element source array.
typedef void (__TBB_EXPORTED_FUNC *internal_array_op2)(void* dst, const void* src, size_type n );
//! Internal structure for compact()
struct internal_segments_table {
segment_index_t first_block;
segment_t table[pointers_per_long_table];
};
void __TBB_EXPORTED_METHOD internal_reserve( size_type n, size_type element_size, size_type max_size );
size_type __TBB_EXPORTED_METHOD internal_capacity() const;
void internal_grow( size_type start, size_type finish, size_type element_size, internal_array_op2 init, const void *src );
size_type __TBB_EXPORTED_METHOD internal_grow_by( size_type delta, size_type element_size, internal_array_op2 init, const void *src );
void* __TBB_EXPORTED_METHOD internal_push_back( size_type element_size, size_type& index );
segment_index_t __TBB_EXPORTED_METHOD internal_clear( internal_array_op1 destroy );
void* __TBB_EXPORTED_METHOD internal_compact( size_type element_size, void *table, internal_array_op1 destroy, internal_array_op2 copy );
void __TBB_EXPORTED_METHOD internal_copy( const concurrent_vector_base_v3& src, size_type element_size, internal_array_op2 copy );
void __TBB_EXPORTED_METHOD internal_assign( const concurrent_vector_base_v3& src, size_type element_size,
internal_array_op1 destroy, internal_array_op2 assign, internal_array_op2 copy );
//! Obsolete
void __TBB_EXPORTED_METHOD internal_throw_exception(size_type) const;
void __TBB_EXPORTED_METHOD internal_swap(concurrent_vector_base_v3& v);
void __TBB_EXPORTED_METHOD internal_resize( size_type n, size_type element_size, size_type max_size, const void *src,
internal_array_op1 destroy, internal_array_op2 init );
size_type __TBB_EXPORTED_METHOD internal_grow_to_at_least_with_result( size_type new_size, size_type element_size, internal_array_op2 init, const void *src );
//! Deprecated entry point for backwards compatibility to TBB 2.1.
void __TBB_EXPORTED_METHOD internal_grow_to_at_least( size_type new_size, size_type element_size, internal_array_op2 init, const void *src );
private:
//! Private functionality
class helper;
friend class helper;
template<typename Container, typename Value>
friend class vector_iterator;
};
inline void swap(concurrent_vector_base_v3::segment_t & lhs, concurrent_vector_base_v3::segment_t & rhs) __TBB_NOEXCEPT(true) {
lhs.swap(rhs);
}
typedef concurrent_vector_base_v3 concurrent_vector_base;
//! Meets requirements of a forward iterator for STL and a Value for a blocked_range.*/
/** Value is either the T or const T type of the container.
@ingroup containers */
template<typename Container, typename Value>
class vector_iterator
{
//! concurrent_vector over which we are iterating.
Container* my_vector;
//! Index into the vector
size_t my_index;
//! Caches my_vector->internal_subscript(my_index)
/** NULL if cached value is not available */
mutable Value* my_item;
template<typename C, typename T>
friend vector_iterator<C,T> operator+( ptrdiff_t offset, const vector_iterator<C,T>& v );
template<typename C, typename T, typename U>
friend bool operator==( const vector_iterator<C,T>& i, const vector_iterator<C,U>& j );
template<typename C, typename T, typename U>
friend bool operator<( const vector_iterator<C,T>& i, const vector_iterator<C,U>& j );
template<typename C, typename T, typename U>
friend ptrdiff_t operator-( const vector_iterator<C,T>& i, const vector_iterator<C,U>& j );
template<typename C, typename U>
friend class internal::vector_iterator;
#if !__TBB_TEMPLATE_FRIENDS_BROKEN
template<typename T, class A>
friend class tbb::concurrent_vector;
#else
public:
#endif
vector_iterator( const Container& vector, size_t index, void *ptr = 0 ) :
my_vector(const_cast<Container*>(&vector)),
my_index(index),
my_item(static_cast<Value*>(ptr))
{}
public:
//! Default constructor
vector_iterator() : my_vector(NULL), my_index(~size_t(0)), my_item(NULL) {}
vector_iterator( const vector_iterator<Container,typename Container::value_type>& other ) :
my_vector(other.my_vector),
my_index(other.my_index),
my_item(other.my_item)
{}
vector_iterator operator+( ptrdiff_t offset ) const {
return vector_iterator( *my_vector, my_index+offset );
}
vector_iterator &operator+=( ptrdiff_t offset ) {
my_index+=offset;
my_item = NULL;
return *this;
}
vector_iterator operator-( ptrdiff_t offset ) const {
return vector_iterator( *my_vector, my_index-offset );
}
vector_iterator &operator-=( ptrdiff_t offset ) {
my_index-=offset;
my_item = NULL;
return *this;
}
Value& operator*() const {
Value* item = my_item;
if( !item ) {
item = my_item = &my_vector->internal_subscript(my_index);
}
__TBB_ASSERT( item==&my_vector->internal_subscript(my_index), "corrupt cache" );
return *item;
}
Value& operator[]( ptrdiff_t k ) const {
return my_vector->internal_subscript(my_index+k);
}
Value* operator->() const {return &operator*();}
//! Pre increment
vector_iterator& operator++() {
size_t element_index = ++my_index;
if( my_item ) {
//TODO: consider using of knowledge about "first_block optimization" here as well?
if( concurrent_vector_base::is_first_element_in_segment(element_index)) {
//if the iterator crosses a segment boundary, the pointer become invalid
//as possibly next segment is in another memory location
my_item= NULL;
} else {
++my_item;
}
}
return *this;
}
//! Pre decrement
vector_iterator& operator--() {
__TBB_ASSERT( my_index>0, "operator--() applied to iterator already at beginning of concurrent_vector" );
size_t element_index = my_index--;
if( my_item ) {
if(concurrent_vector_base::is_first_element_in_segment(element_index)) {
//if the iterator crosses a segment boundary, the pointer become invalid
//as possibly next segment is in another memory location
my_item= NULL;
} else {
--my_item;
}
}
return *this;
}
//! Post increment
vector_iterator operator++(int) {
vector_iterator result = *this;
operator++();
return result;
}
//! Post decrement
vector_iterator operator--(int) {
vector_iterator result = *this;
operator--();
return result;
}
// STL support
typedef ptrdiff_t difference_type;
typedef Value value_type;
typedef Value* pointer;
typedef Value& reference;
typedef std::random_access_iterator_tag iterator_category;
};
template<typename Container, typename T>
vector_iterator<Container,T> operator+( ptrdiff_t offset, const vector_iterator<Container,T>& v ) {
return vector_iterator<Container,T>( *v.my_vector, v.my_index+offset );
}
template<typename Container, typename T, typename U>
bool operator==( const vector_iterator<Container,T>& i, const vector_iterator<Container,U>& j ) {
return i.my_index==j.my_index && i.my_vector == j.my_vector;
}
template<typename Container, typename T, typename U>
bool operator!=( const vector_iterator<Container,T>& i, const vector_iterator<Container,U>& j ) {
return !(i==j);
}
template<typename Container, typename T, typename U>
bool operator<( const vector_iterator<Container,T>& i, const vector_iterator<Container,U>& j ) {
return i.my_index<j.my_index;
}
template<typename Container, typename T, typename U>
bool operator>( const vector_iterator<Container,T>& i, const vector_iterator<Container,U>& j ) {
return j<i;
}
template<typename Container, typename T, typename U>
bool operator>=( const vector_iterator<Container,T>& i, const vector_iterator<Container,U>& j ) {
return !(i<j);
}
template<typename Container, typename T, typename U>
bool operator<=( const vector_iterator<Container,T>& i, const vector_iterator<Container,U>& j ) {
return !(j<i);
}
template<typename Container, typename T, typename U>
ptrdiff_t operator-( const vector_iterator<Container,T>& i, const vector_iterator<Container,U>& j ) {
return ptrdiff_t(i.my_index)-ptrdiff_t(j.my_index);
}
template<typename T, class A>
class allocator_base {
public:
typedef typename A::template
rebind<T>::other allocator_type;
allocator_type my_allocator;
allocator_base(const allocator_type &a = allocator_type() ) : my_allocator(a) {}
};
} // namespace internal
//! @endcond
//! Concurrent vector container
/** concurrent_vector is a container having the following main properties:
- It provides random indexed access to its elements. The index of the first element is 0.
- It ensures safe concurrent growing its size (different threads can safely append new elements).
- Adding new elements does not invalidate existing iterators and does not change indices of existing items.
@par Compatibility
The class meets all Container Requirements and Reversible Container Requirements from
C++ Standard (See ISO/IEC 14882:2003(E), clause 23.1). But it doesn't meet
Sequence Requirements due to absence of insert() and erase() methods.
@par Exception Safety
Methods working with memory allocation and/or new elements construction can throw an
exception if allocator fails to allocate memory or element's default constructor throws one.
Concurrent vector's element of type T must conform to the following requirements:
- Throwing an exception is forbidden for destructor of T.
- Default constructor of T must not throw an exception OR its non-virtual destructor must safely work when its object memory is zero-initialized.
.
Otherwise, the program's behavior is undefined.
@par
If an exception happens inside growth or assignment operation, an instance of the vector becomes invalid unless it is stated otherwise in the method documentation.
Invalid state means:
- There are no guarantees that all items were initialized by a constructor. The rest of items is zero-filled, including item where exception happens.
- An invalid vector instance cannot be repaired; it is unable to grow anymore.
- Size and capacity reported by the vector are incorrect, and calculated as if the failed operation were successful.
- Attempt to access not allocated elements using operator[] or iterators results in access violation or segmentation fault exception, and in case of using at() method a C++ exception is thrown.
.
If a concurrent grow operation successfully completes, all the elements it has added to the vector will remain valid and accessible even if one of subsequent grow operations fails.
@par Fragmentation
Unlike an STL vector, a concurrent_vector does not move existing elements if it needs
to allocate more memory. The container is divided into a series of contiguous arrays of
elements. The first reservation, growth, or assignment operation determines the size of
the first array. Using small number of elements as initial size incurs fragmentation that
may increase element access time. Internal layout can be optimized by method compact() that
merges several smaller arrays into one solid.
@par Changes since TBB 2.1
- Fixed guarantees of concurrent_vector::size() and grow_to_at_least() methods to assure elements are allocated.
- Methods end()/rbegin()/back() are partly thread-safe since they use size() to get the end of vector
- Added resize() methods (not thread-safe)
- Added cbegin/cend/crbegin/crend methods
- Changed return type of methods grow* and push_back to iterator
@par Changes since TBB 2.0
- Implemented exception-safety guarantees
- Added template argument for allocator
- Added allocator argument in constructors
- Faster index calculation
- First growth call specifies a number of segments to be merged in the first allocation.
- Fixed memory blow up for swarm of vector's instances of small size
- Added grow_by(size_type n, const_reference t) growth using copying constructor to init new items.
- Added STL-like constructors.
- Added operators ==, < and derivatives
- Added at() method, approved for using after an exception was thrown inside the vector
- Added get_allocator() method.
- Added assign() methods
- Added compact() method to defragment first segments
- Added swap() method
- range() defaults on grainsize = 1 supporting auto grainsize algorithms.
@ingroup containers */
template<typename T, class A>
class concurrent_vector: protected internal::allocator_base<T, A>,
private internal::concurrent_vector_base {
private:
template<typename I>
class generic_range_type: public blocked_range<I> {
public:
typedef T value_type;
typedef T& reference;
typedef const T& const_reference;
typedef I iterator;
typedef ptrdiff_t difference_type;
generic_range_type( I begin_, I end_, size_t grainsize_ = 1) : blocked_range<I>(begin_,end_,grainsize_) {}
template<typename U>
generic_range_type( const generic_range_type<U>& r) : blocked_range<I>(r.begin(),r.end(),r.grainsize()) {}
generic_range_type( generic_range_type& r, split ) : blocked_range<I>(r,split()) {}
};
template<typename C, typename U>
friend class internal::vector_iterator;
public:
//------------------------------------------------------------------------
// STL compatible types
//------------------------------------------------------------------------
typedef internal::concurrent_vector_base_v3::size_type size_type;
typedef typename internal::allocator_base<T, A>::allocator_type allocator_type;
typedef T value_type;
typedef ptrdiff_t difference_type;
typedef T& reference;
typedef const T& const_reference;
typedef T *pointer;
typedef const T *const_pointer;
typedef internal::vector_iterator<concurrent_vector,T> iterator;
typedef internal::vector_iterator<concurrent_vector,const T> const_iterator;
#if !defined(_MSC_VER) || _CPPLIB_VER>=300
// Assume ISO standard definition of std::reverse_iterator
typedef std::reverse_iterator<iterator> reverse_iterator;
typedef std::reverse_iterator<const_iterator> const_reverse_iterator;
#else
// Use non-standard std::reverse_iterator
typedef std::reverse_iterator<iterator,T,T&,T*> reverse_iterator;
typedef std::reverse_iterator<const_iterator,T,const T&,const T*> const_reverse_iterator;
#endif /* defined(_MSC_VER) && (_MSC_VER<1300) */
//------------------------------------------------------------------------
// Parallel algorithm support
//------------------------------------------------------------------------
typedef generic_range_type<iterator> range_type;
typedef generic_range_type<const_iterator> const_range_type;
//------------------------------------------------------------------------
// STL compatible constructors & destructors
//------------------------------------------------------------------------
//! Construct empty vector.
explicit concurrent_vector(const allocator_type &a = allocator_type())
: internal::allocator_base<T, A>(a), internal::concurrent_vector_base()
{
vector_allocator_ptr = &internal_allocator;
}
//Constructors are not required to have synchronization
//(for more details see comment in the concurrent_vector_base constructor).
#if __TBB_INITIALIZER_LISTS_PRESENT
//! Constructor from initializer_list
concurrent_vector(std::initializer_list<T> init_list, const allocator_type &a = allocator_type())
: internal::allocator_base<T, A>(a), internal::concurrent_vector_base()
{
vector_allocator_ptr = &internal_allocator;
__TBB_TRY {
internal_assign_iterators(init_list.begin(), init_list.end());
} __TBB_CATCH(...) {
segment_t *table = my_segment.load<relaxed>();;
internal_free_segments( table, internal_clear(&destroy_array), my_first_block.load<relaxed>());
__TBB_RETHROW();
}
}
#endif //# __TBB_INITIALIZER_LISTS_PRESENT
//! Copying constructor
concurrent_vector( const concurrent_vector& vector, const allocator_type& a = allocator_type() )
: internal::allocator_base<T, A>(a), internal::concurrent_vector_base()
{
vector_allocator_ptr = &internal_allocator;
__TBB_TRY {
internal_copy(vector, sizeof(T), ©_array);
} __TBB_CATCH(...) {
segment_t *table = my_segment.load<relaxed>();
internal_free_segments( table, internal_clear(&destroy_array), my_first_block.load<relaxed>());
__TBB_RETHROW();
}
}
#if __TBB_CPP11_RVALUE_REF_PRESENT
//! Move constructor
//TODO add __TBB_NOEXCEPT(true) and static_assert(std::has_nothrow_move_constructor<A>::value)
concurrent_vector( concurrent_vector&& source)
: internal::allocator_base<T, A>(std::move(source)), internal::concurrent_vector_base()
{
vector_allocator_ptr = &internal_allocator;
concurrent_vector_base_v3::internal_swap(source);
}
concurrent_vector( concurrent_vector&& source, const allocator_type& a)
: internal::allocator_base<T, A>(a), internal::concurrent_vector_base()
{
vector_allocator_ptr = &internal_allocator;
//C++ standard requires instances of an allocator being compared for equality,
//which means that memory allocated by one instance is possible to deallocate with the other one.
if (a == source.my_allocator) {
concurrent_vector_base_v3::internal_swap(source);
} else {
__TBB_TRY {
internal_copy(source, sizeof(T), &move_array);
} __TBB_CATCH(...) {
segment_t *table = my_segment.load<relaxed>();
internal_free_segments( table, internal_clear(&destroy_array), my_first_block.load<relaxed>());
__TBB_RETHROW();
}
}
}
#endif
//! Copying constructor for vector with different allocator type
template<class M>
concurrent_vector( const concurrent_vector<T, M>& vector, const allocator_type& a = allocator_type() )
: internal::allocator_base<T, A>(a), internal::concurrent_vector_base()
{
vector_allocator_ptr = &internal_allocator;
__TBB_TRY {
internal_copy(vector.internal_vector_base(), sizeof(T), ©_array);
} __TBB_CATCH(...) {
segment_t *table = my_segment.load<relaxed>();
internal_free_segments( table, internal_clear(&destroy_array), my_first_block.load<relaxed>() );
__TBB_RETHROW();
}
}
//! Construction with initial size specified by argument n
explicit concurrent_vector(size_type n)
{
vector_allocator_ptr = &internal_allocator;
__TBB_TRY {
internal_resize( n, sizeof(T), max_size(), NULL, &destroy_array, &initialize_array );
} __TBB_CATCH(...) {
segment_t *table = my_segment.load<relaxed>();
internal_free_segments( table, internal_clear(&destroy_array), my_first_block.load<relaxed>() );
__TBB_RETHROW();
}
}
//! Construction with initial size specified by argument n, initialization by copying of t, and given allocator instance
concurrent_vector(size_type n, const_reference t, const allocator_type& a = allocator_type())
: internal::allocator_base<T, A>(a)
{
vector_allocator_ptr = &internal_allocator;
__TBB_TRY {
internal_resize( n, sizeof(T), max_size(), static_cast<const void*>(&t), &destroy_array, &initialize_array_by );
} __TBB_CATCH(...) {
segment_t *table = my_segment.load<relaxed>();
internal_free_segments( table, internal_clear(&destroy_array), my_first_block.load<relaxed>() );
__TBB_RETHROW();
}
}
//! Construction with copying iteration range and given allocator instance
template<class I>
concurrent_vector(I first, I last, const allocator_type &a = allocator_type())
: internal::allocator_base<T, A>(a)
{
vector_allocator_ptr = &internal_allocator;
__TBB_TRY {
internal_assign_range(first, last, static_cast<is_integer_tag<std::numeric_limits<I>::is_integer> *>(0) );
} __TBB_CATCH(...) {
segment_t *table = my_segment.load<relaxed>();
internal_free_segments( table, internal_clear(&destroy_array), my_first_block.load<relaxed>() );
__TBB_RETHROW();
}
}
//! Assignment
concurrent_vector& operator=( const concurrent_vector& vector ) {
if( this != &vector )
internal_assign(vector, sizeof(T), &destroy_array, &assign_array, ©_array);
return *this;
}
#if __TBB_CPP11_RVALUE_REF_PRESENT
//TODO: add __TBB_NOEXCEPT()
//! Move assignment
concurrent_vector& operator=( concurrent_vector&& other ) {
__TBB_ASSERT(this != &other, "Move assignment to itself is prohibited ");
typedef typename tbb::internal::allocator_traits<A>::propagate_on_container_move_assignment pocma_t;
if(pocma_t::value || this->my_allocator == other.my_allocator) {
concurrent_vector trash (std::move(*this));
internal_swap(other);
if (pocma_t::value) {
this->my_allocator = std::move(other.my_allocator);
}
} else {
internal_assign(other, sizeof(T), &destroy_array, &move_assign_array, &move_array);
}
return *this;
}
#endif
//TODO: add an template assignment operator? (i.e. with different element type)
//! Assignment for vector with different allocator type
template<class M>
concurrent_vector& operator=( const concurrent_vector<T, M>& vector ) {
if( static_cast<void*>( this ) != static_cast<const void*>( &vector ) )
internal_assign(vector.internal_vector_base(),
sizeof(T), &destroy_array, &assign_array, ©_array);
return *this;
}
#if __TBB_INITIALIZER_LISTS_PRESENT
//! Assignment for initializer_list
concurrent_vector& operator=( std::initializer_list<T> init_list ) {
internal_clear(&destroy_array);
internal_assign_iterators(init_list.begin(), init_list.end());
return *this;
}
#endif //#if __TBB_INITIALIZER_LISTS_PRESENT
//------------------------------------------------------------------------
// Concurrent operations
//------------------------------------------------------------------------
//! Grow by "delta" elements.
/** Returns iterator pointing to the first new element. */
iterator grow_by( size_type delta ) {
return iterator(*this, delta ? internal_grow_by( delta, sizeof(T), &initialize_array, NULL ) : my_early_size.load());
}
//! Grow by "delta" elements using copying constructor.
/** Returns iterator pointing to the first new element. */
iterator grow_by( size_type delta, const_reference t ) {
return iterator(*this, delta ? internal_grow_by( delta, sizeof(T), &initialize_array_by, static_cast<const void*>(&t) ) : my_early_size.load());
}
/** Returns iterator pointing to the first new element. */
template<typename I>
iterator grow_by( I first, I last ) {
typename std::iterator_traits<I>::difference_type delta = std::distance(first, last);
__TBB_ASSERT( delta >= 0, NULL);
return iterator(*this, delta ? internal_grow_by(delta, sizeof(T), ©_range<I>, static_cast<const void*>(&first)) : my_early_size.load());
}
#if __TBB_INITIALIZER_LISTS_PRESENT
/** Returns iterator pointing to the first new element. */
iterator grow_by( std::initializer_list<T> init_list ) {
return grow_by( init_list.begin(), init_list.end() );
}
#endif //#if __TBB_INITIALIZER_LISTS_PRESENT
//! Append minimal sequence of elements such that size()>=n.
/** The new elements are default constructed. Blocks until all elements in range [0..n) are allocated.
May return while other elements are being constructed by other threads.
Returns iterator that points to beginning of appended sequence.
If no elements were appended, returns iterator pointing to nth element. */
iterator grow_to_at_least( size_type n ) {
size_type m=0;
if( n ) {
m = internal_grow_to_at_least_with_result( n, sizeof(T), &initialize_array, NULL );
if( m>n ) m=n;
}
return iterator(*this, m);
};
/** Analogous to grow_to_at_least( size_type n ) with exception that the new
elements are initialized by copying of t instead of default construction. */
iterator grow_to_at_least( size_type n, const_reference t ) {
size_type m=0;
if( n ) {
m = internal_grow_to_at_least_with_result( n, sizeof(T), &initialize_array_by, &t);
if( m>n ) m=n;
}
return iterator(*this, m);
};
//! Push item
/** Returns iterator pointing to the new element. */
iterator push_back( const_reference item )
{
push_back_helper prolog(*this);
new(prolog.internal_push_back_result()) T(item);
return prolog.return_iterator_and_dismiss();
}
#if __TBB_CPP11_RVALUE_REF_PRESENT
//! Push item, move-aware
/** Returns iterator pointing to the new element. */
iterator push_back( T&& item )
{
push_back_helper prolog(*this);
new(prolog.internal_push_back_result()) T(std::move(item));
return prolog.return_iterator_and_dismiss();
}
#if __TBB_CPP11_VARIADIC_TEMPLATES_PRESENT
//! Push item, create item "in place" with provided arguments
/** Returns iterator pointing to the new element. */
template<typename... Args>
iterator emplace_back( Args&&... args )
{
push_back_helper prolog(*this);
new(prolog.internal_push_back_result()) T(std::forward<Args>(args)...);
return prolog.return_iterator_and_dismiss();
}
#endif //__TBB_CPP11_VARIADIC_TEMPLATES_PRESENT
#endif //__TBB_CPP11_RVALUE_REF_PRESENT
//! Get reference to element at given index.
/** This method is thread-safe for concurrent reads, and also while growing the vector,
as long as the calling thread has checked that index < size(). */
reference operator[]( size_type index ) {
return internal_subscript(index);
}
//! Get const reference to element at given index.
const_reference operator[]( size_type index ) const {
return internal_subscript(index);
}
//! Get reference to element at given index. Throws exceptions on errors.
reference at( size_type index ) {
return internal_subscript_with_exceptions(index);
}
//! Get const reference to element at given index. Throws exceptions on errors.
const_reference at( size_type index ) const {
return internal_subscript_with_exceptions(index);
}
//! Get range for iterating with parallel algorithms
range_type range( size_t grainsize = 1 ) {
return range_type( begin(), end(), grainsize );
}
//! Get const range for iterating with parallel algorithms
const_range_type range( size_t grainsize = 1 ) const {
return const_range_type( begin(), end(), grainsize );
}
//------------------------------------------------------------------------
// Capacity
//------------------------------------------------------------------------
//! Return size of vector. It may include elements under construction
size_type size() const {
size_type sz = my_early_size, cp = internal_capacity();
return cp < sz ? cp : sz;
}
//! Return false if vector is not empty or has elements under construction at least.
bool empty() const {return !my_early_size;}
//! Maximum size to which array can grow without allocating more memory. Concurrent allocations are not included in the value.
size_type capacity() const {return internal_capacity();}
//! Allocate enough space to grow to size n without having to allocate more memory later.
/** Like most of the methods provided for STL compatibility, this method is *not* thread safe.
The capacity afterwards may be bigger than the requested reservation. */
void reserve( size_type n ) {
if( n )
internal_reserve(n, sizeof(T), max_size());
}
//! Resize the vector. Not thread-safe.
void resize( size_type n ) {
internal_resize( n, sizeof(T), max_size(), NULL, &destroy_array, &initialize_array );
}
//! Resize the vector, copy t for new elements. Not thread-safe.
void resize( size_type n, const_reference t ) {
internal_resize( n, sizeof(T), max_size(), static_cast<const void*>(&t), &destroy_array, &initialize_array_by );
}
//! Optimize memory usage and fragmentation.
void shrink_to_fit();
//! Upper bound on argument to reserve.
size_type max_size() const {return (~size_type(0))/sizeof(T);}
//------------------------------------------------------------------------
// STL support
//------------------------------------------------------------------------
//! start iterator
iterator begin() {return iterator(*this,0);}
//! end iterator
iterator end() {return iterator(*this,size());}
//! start const iterator
const_iterator begin() const {return const_iterator(*this,0);}
//! end const iterator
const_iterator end() const {return const_iterator(*this,size());}
//! start const iterator
const_iterator cbegin() const {return const_iterator(*this,0);}
//! end const iterator
const_iterator cend() const {return const_iterator(*this,size());}
//! reverse start iterator
reverse_iterator rbegin() {return reverse_iterator(end());}
//! reverse end iterator
reverse_iterator rend() {return reverse_iterator(begin());}
//! reverse start const iterator
const_reverse_iterator rbegin() const {return const_reverse_iterator(end());}
//! reverse end const iterator
const_reverse_iterator rend() const {return const_reverse_iterator(begin());}
//! reverse start const iterator
const_reverse_iterator crbegin() const {return const_reverse_iterator(end());}
//! reverse end const iterator
const_reverse_iterator crend() const {return const_reverse_iterator(begin());}
//! the first item
reference front() {
__TBB_ASSERT( size()>0, NULL);
const segment_value_t& segment_value = my_segment[0].template load<relaxed>();
return (segment_value.template pointer<T>())[0];
}
//! the first item const
const_reference front() const {
__TBB_ASSERT( size()>0, NULL);
const segment_value_t& segment_value = my_segment[0].template load<relaxed>();
return (segment_value.template pointer<const T>())[0];
}
//! the last item
reference back() {
__TBB_ASSERT( size()>0, NULL);
return internal_subscript( size()-1 );
}
//! the last item const
const_reference back() const {
__TBB_ASSERT( size()>0, NULL);
return internal_subscript( size()-1 );
}
//! return allocator object
allocator_type get_allocator() const { return this->my_allocator; }
//! assign n items by copying t item
void assign(size_type n, const_reference t) {
clear();
internal_resize( n, sizeof(T), max_size(), static_cast<const void*>(&t), &destroy_array, &initialize_array_by );
}