/elec/propeller-clock

To get this branch, use:
bzr branch http://bzr.ed.am/elec/propeller-clock

« back to all changes in this revision

Viewing changes to src/utility/algorithm

  • Committer: edam
  • Date: 2011-11-17 11:38:11 UTC
  • Revision ID: edam@waxworlds.org-20111117113811-os1jlwxwukdxt6um
added info about RC circuits

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
 
/*      Copyright (C) 2004 Garrett A. Kajmowicz
2
 
        This file is part of the uClibc++ Library.
3
 
        This library is free software; you can redistribute it and/or
4
 
        modify it under the terms of the GNU Lesser General Public
5
 
        License as published by the Free Software Foundation; either
6
 
        version 2.1 of the License, or (at your option) any later version.
7
 
 
8
 
        This library is distributed in the hope that it will be useful,
9
 
        but WITHOUT ANY WARRANTY; without even the implied warranty of
10
 
        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
11
 
        Lesser General Public License for more details.
12
 
 
13
 
        You should have received a copy of the GNU Lesser General Public
14
 
        License along with this library; if not, write to the Free Software
15
 
        Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
16
 
*/
17
 
 
18
 
#include <cstdlib>
19
 
#include <iterator>
20
 
#include <utility.h>
21
 
#include <functional>
22
 
 
23
 
#ifndef __STD_HEADER_ALGORITHM
24
 
#define __STD_HEADER_ALGORITHM 1
25
 
 
26
 
//Elliminate any previously defined macro
27
 
#undef min
28
 
#undef max
29
 
 
30
 
#pragma GCC visibility push(default)
31
 
 
32
 
namespace std{
33
 
 
34
 
        // subclause _lib.alg.nonmodifying_, non-modifying sequence operations:
35
 
        template<class InputIterator, class Function> _UCXXEXPORT
36
 
                Function for_each(InputIterator first, InputIterator last, Function f)
37
 
        {
38
 
                while(first !=last){
39
 
                        f(*first);
40
 
                        ++first;
41
 
                }
42
 
                return f;
43
 
        }
44
 
 
45
 
 
46
 
        template<class InputIterator, class T> _UCXXEXPORT
47
 
                InputIterator find(InputIterator first, InputIterator last, const T& value)
48
 
        {
49
 
                while(first !=last && ! ( *first == value ) ){
50
 
                        ++first;
51
 
                }
52
 
                return first;
53
 
        }
54
 
 
55
 
        template<class InputIterator, class Predicate> _UCXXEXPORT
56
 
                InputIterator find_if(InputIterator first, InputIterator last, Predicate pred)
57
 
        {
58
 
                while(first !=last && !pred(*first)){
59
 
                        ++first;
60
 
                }
61
 
                return first;
62
 
        }
63
 
 
64
 
        template<class ForwardIterator1, class ForwardIterator2> _UCXXEXPORT ForwardIterator1
65
 
                find_end(ForwardIterator1 first1, ForwardIterator1 last1,
66
 
                        ForwardIterator2 first2, ForwardIterator2 last2)
67
 
        {
68
 
                ForwardIterator1 retval = last1;
69
 
                while(first1 !=last1){
70
 
                        ForwardIterator1 temp1(first1);
71
 
                        ForwardIterator2 temp2(first2);
72
 
                        while(temp2!=last2 && temp1!= last1){
73
 
                                if(*temp1 != *temp2){
74
 
                                        break;          //Exit while loop
75
 
                                }
76
 
                                ++temp1;
77
 
                                ++temp2;
78
 
                                if(temp2 == last2){     //Have successfully checked the whole sequence
79
 
                                        retval = first1;
80
 
                                }
81
 
                        }
82
 
                }
83
 
                return retval;
84
 
        }
85
 
 
86
 
        template<class ForwardIterator1, class ForwardIterator2, class BinaryPredicate> _UCXXEXPORT
87
 
                ForwardIterator1
88
 
                find_end(ForwardIterator1 first1, ForwardIterator1 last1,
89
 
                        ForwardIterator2 first2, ForwardIterator2 last2,
90
 
                        BinaryPredicate pred)
91
 
        {
92
 
                ForwardIterator1 retval = last1;
93
 
                while(first1 !=last1){
94
 
                        ForwardIterator1 temp1(first1);
95
 
                        ForwardIterator2 temp2(first2);
96
 
                        while(temp2!=last2 && temp1!= last1){
97
 
                                if( ! pred(*temp1, *temp2)){
98
 
                                        break;          //Exit while loop
99
 
                                }
100
 
                                ++temp1;
101
 
                                ++temp2;
102
 
                                if(temp2 == last2){     //Have successfully checked the whole sequence
103
 
                                        retval = first1;
104
 
                                }
105
 
                        }
106
 
                }
107
 
                return retval;
108
 
        }
109
 
 
110
 
        template<class ForwardIterator1, class ForwardIterator2> _UCXXEXPORT
111
 
                ForwardIterator1
112
 
                find_first_of(ForwardIterator1 first1, ForwardIterator1 last1,
113
 
                        ForwardIterator2 first2, ForwardIterator2 last2)
114
 
        {
115
 
                while(first1 != last1){
116
 
                        ForwardIterator2 temp2(first2);
117
 
                        while(temp2 != last2 ){
118
 
                                if(*temp2 == first1){
119
 
                                        return first1;
120
 
                                }
121
 
                                ++temp2;
122
 
                        }
123
 
                        
124
 
                }
125
 
                return last1;
126
 
        }
127
 
 
128
 
        template<class ForwardIterator1, class ForwardIterator2, class BinaryPredicate> _UCXXEXPORT
129
 
                ForwardIterator1
130
 
                find_first_of(ForwardIterator1 first1, ForwardIterator1 last1,  
131
 
                        ForwardIterator2 first2, ForwardIterator2 last2,
132
 
                        BinaryPredicate pred)
133
 
        {
134
 
                while(first1 != last1){
135
 
                        ForwardIterator2 temp2(first2);
136
 
                        while(temp2 != last2 ){
137
 
                                if(*temp2 == first1){
138
 
                                        return first1;
139
 
                                }
140
 
                                ++temp2;
141
 
                        }
142
 
                        
143
 
                }
144
 
                return last1;
145
 
        }
146
 
 
147
 
        template<class ForwardIterator> _UCXXEXPORT ForwardIterator
148
 
                adjacent_find(ForwardIterator first, ForwardIterator last)
149
 
        {
150
 
                ForwardIterator temp;
151
 
                while(first != last){
152
 
                        temp = first;
153
 
                        ++temp;
154
 
                        if(*temp == *first){
155
 
                                return first;
156
 
                        }
157
 
                        ++first;
158
 
                }
159
 
                return first;
160
 
        }
161
 
 
162
 
        template<class ForwardIterator, class BinaryPredicate> _UCXXEXPORT
163
 
                ForwardIterator
164
 
                adjacent_find(ForwardIterator first, ForwardIterator last, BinaryPredicate pred)
165
 
        {
166
 
                ForwardIterator temp;
167
 
                while(first != last){
168
 
                        temp = first;
169
 
                        ++temp;
170
 
                        if( pred(*temp, *first)){
171
 
                                return first;
172
 
                        }
173
 
                        ++first;
174
 
                }
175
 
                return first;
176
 
        }
177
 
 
178
 
        template<class InputIterator, class T> _UCXXEXPORT
179
 
                typename iterator_traits<InputIterator>::difference_type
180
 
                count(InputIterator first, InputIterator last, const T& value)
181
 
        {
182
 
                typename iterator_traits<InputIterator>::difference_type i = 0;
183
 
                while(first!=last){
184
 
                        if(*first == value){
185
 
                                ++i;
186
 
                        }
187
 
                        ++first;
188
 
                }
189
 
                return i;
190
 
        }
191
 
 
192
 
        template<class InputIterator, class Predicate> _UCXXEXPORT
193
 
                typename iterator_traits<InputIterator>::difference_type
194
 
                count_if(InputIterator first, InputIterator last, Predicate pred)
195
 
        {
196
 
                typename iterator_traits<InputIterator>::difference_type i = 0;
197
 
                while(first!=last){
198
 
                        if( pred(*first) ){
199
 
                                ++i;
200
 
                        }
201
 
                        ++first;
202
 
                }
203
 
                return i;
204
 
        }
205
 
 
206
 
        template<class InputIterator1, class InputIterator2> _UCXXEXPORT
207
 
                pair<InputIterator1, InputIterator2>    
208
 
                mismatch(InputIterator1 first1, InputIterator1 last1, InputIterator2 first2)
209
 
        {
210
 
                while(first1 != last1){
211
 
                        if(*first1 != *first2){
212
 
                                break;
213
 
                        }
214
 
                        ++first1;
215
 
                        ++first2;
216
 
                }
217
 
 
218
 
                pair<InputIterator1, InputIterator2> retval;
219
 
                retval.first = first1;
220
 
                retval.second = first2;
221
 
                return retval;
222
 
        }
223
 
 
224
 
        template<class InputIterator1, class InputIterator2, class BinaryPredicate> _UCXXEXPORT
225
 
                pair<InputIterator1, InputIterator2>
226
 
                mismatch(InputIterator1 first1, InputIterator1 last1,
227
 
                        InputIterator2 first2, BinaryPredicate pred)
228
 
        {
229
 
                while(first1 != last1){
230
 
                        if( !pred(*first1, *first2) ){
231
 
                                break;
232
 
                        }
233
 
                        ++first1;
234
 
                        ++first2;
235
 
                }
236
 
 
237
 
                pair<InputIterator1, InputIterator2> retval;
238
 
                retval.first = first1;
239
 
                retval.second = first2;
240
 
                return retval;
241
 
        }
242
 
 
243
 
        template<class InputIterator1, class InputIterator2> _UCXXEXPORT
244
 
                bool
245
 
                equal(InputIterator1 first1, InputIterator1 last1, InputIterator2 first2)
246
 
        {
247
 
                while(first1 !=last1){
248
 
                        if(*first1 != *first2){
249
 
                                return false;
250
 
                        }
251
 
                        ++first1;
252
 
                        ++first2;
253
 
                }
254
 
                return true;
255
 
        }
256
 
 
257
 
        template<class InputIterator1, class InputIterator2, class BinaryPredicate> _UCXXEXPORT
258
 
                bool
259
 
                equal(InputIterator1 first1, InputIterator1 last1,
260
 
                        InputIterator2 first2, BinaryPredicate pred)
261
 
        {
262
 
                while(first1 !=last1){
263
 
                        if( !pred(*first1, *first2) ){
264
 
                                return false;
265
 
                        }
266
 
                        ++first1;
267
 
                        ++first2;
268
 
                }
269
 
                return true;
270
 
        }
271
 
 
272
 
        template<class ForwardIterator1, class ForwardIterator2> _UCXXEXPORT
273
 
                ForwardIterator1 search(ForwardIterator1 first1, ForwardIterator1 last1,
274
 
                        ForwardIterator2 first2, ForwardIterator2 last2)
275
 
        {
276
 
                equal_to<typename iterator_traits<ForwardIterator1>::value_type> c;
277
 
                return search(first1, last1, first2, last2, c);
278
 
        }
279
 
 
280
 
 
281
 
        template<class ForwardIterator1, class ForwardIterator2, class BinaryPredicate> _UCXXEXPORT
282
 
                ForwardIterator1
283
 
                search(ForwardIterator1 first1, ForwardIterator1 last1,
284
 
                        ForwardIterator2 first2, ForwardIterator2 last2,
285
 
                        BinaryPredicate pred)
286
 
        {
287
 
                while(first1 != last1){
288
 
                        ForwardIterator1 temp1(first1);
289
 
                        ForwardIterator2 temp2(first2);
290
 
                        while(temp2 != last2 && temp1 != last1){
291
 
                                if( !pred(*temp2, *temp1) ){
292
 
                                        break;
293
 
                                }
294
 
                                ++temp1;
295
 
                                ++temp2;
296
 
                                if(temp2 == last2){
297
 
                                        return first1;
298
 
                                }
299
 
                        }
300
 
                        ++first1;
301
 
                }
302
 
                return last1;
303
 
        }
304
 
 
305
 
 
306
 
        template<class ForwardIterator, class Size, class T> _UCXXEXPORT
307
 
                ForwardIterator
308
 
                search_n(ForwardIterator first, ForwardIterator last, Size count, const T& value)
309
 
        {
310
 
                while( first != last ){
311
 
                        if(*first == value){
312
 
                                ForwardIterator temp(first);
313
 
                                Size i = 1;
314
 
                                ++first;        //Avoid doing the same comparison over again
315
 
                                while(i < count && *temp == value){
316
 
                                        ++first;
317
 
                                        ++i;
318
 
                                }
319
 
                                if(i == count){
320
 
                                        return first;
321
 
                                }
322
 
                        }
323
 
                        ++first;
324
 
                }
325
 
                return first;           
326
 
        }
327
 
 
328
 
 
329
 
        template<class ForwardIterator, class Size, class T, class BinaryPredicate> _UCXXEXPORT
330
 
                ForwardIterator
331
 
                search_n(ForwardIterator first, ForwardIterator last,
332
 
                        Size count, const T& value, BinaryPredicate pred)
333
 
        {
334
 
                while( first != last ){
335
 
                        if( pred(*first, value) ){
336
 
                                ForwardIterator temp(first);
337
 
                                Size i = 1;
338
 
                                ++first;        //Avoid doing the same comparison over again
339
 
                                while(i < count && pred(*temp, value) ){
340
 
                                        ++first;
341
 
                                        ++i;
342
 
                                }
343
 
                                if(i == count){
344
 
                                        return first;
345
 
                                }
346
 
                        }
347
 
                        ++first;
348
 
                }
349
 
                return first;           
350
 
 
351
 
        }
352
 
 
353
 
        // subclause _lib.alg.modifying.operations_, modifying sequence operations:
354
 
 
355
 
        template<class InputIterator, class OutputIterator> _UCXXEXPORT
356
 
                OutputIterator
357
 
                copy(InputIterator first, InputIterator last, OutputIterator result)
358
 
        {
359
 
                while(first != last){
360
 
                        *result = *first;
361
 
                        ++first;
362
 
                        ++result;
363
 
                }
364
 
                return result;
365
 
        }
366
 
 
367
 
        template<class BidirectionalIterator1, class BidirectionalIterator2> _UCXXEXPORT
368
 
                BidirectionalIterator2
369
 
                copy_backward(BidirectionalIterator1 first, BidirectionalIterator1 last,
370
 
                        BidirectionalIterator2 result)
371
 
        {
372
 
                while(first !=last ){
373
 
                        --result;
374
 
                        --last;
375
 
                        *result = *last;
376
 
                }
377
 
                return result;
378
 
        }
379
 
 
380
 
        template<class T> _UCXXEXPORT void swap(T& a, T& b){
381
 
                T temp(a);
382
 
                a = b;
383
 
                b = temp;
384
 
        }
385
 
 
386
 
        template<class ForwardIterator1, class ForwardIterator2> _UCXXEXPORT
387
 
                ForwardIterator2
388
 
                swap_ranges(ForwardIterator1 first1, ForwardIterator1 last1, ForwardIterator2 first2)
389
 
        {
390
 
                while(first1 != last1){
391
 
                        iter_swap(first1, first2);
392
 
                        ++first1;
393
 
                        ++first2;
394
 
                }
395
 
                return first2;
396
 
        }
397
 
 
398
 
 
399
 
        template<class ForwardIterator1, class ForwardIterator2> _UCXXEXPORT
400
 
                void
401
 
                iter_swap(ForwardIterator1 a, ForwardIterator2 b)
402
 
        {
403
 
                typename iterator_traits<ForwardIterator1>::value_type temp(*a);
404
 
                *a = *b;
405
 
                *b = temp;
406
 
        }
407
 
 
408
 
 
409
 
        template<class InputIterator, class OutputIterator, class UnaryOperation> _UCXXEXPORT
410
 
                OutputIterator
411
 
                transform(InputIterator first, InputIterator last,
412
 
                        OutputIterator result, UnaryOperation op)
413
 
        {
414
 
                while(first != last){
415
 
                        *result = op(*first);
416
 
                        ++first;
417
 
                        ++result;
418
 
                }
419
 
                return result;
420
 
        }
421
 
 
422
 
 
423
 
        template<class InputIterator1, class InputIterator2, class OutputIterator, class BinaryOperation> _UCXXEXPORT
424
 
                OutputIterator transform(InputIterator1 first1, InputIterator1 last1,
425
 
                        InputIterator2 first2, OutputIterator result,
426
 
                        BinaryOperation binary_op)
427
 
        {
428
 
                while(first1 != last1){
429
 
                        *result = binary_op(*first1, *first2);
430
 
                        ++first1;
431
 
                        ++first2;
432
 
                        ++result;
433
 
                }
434
 
                return result;
435
 
        }
436
 
 
437
 
 
438
 
        template<class ForwardIterator, class T> _UCXXEXPORT
439
 
                void
440
 
                replace(ForwardIterator first, ForwardIterator last,
441
 
                        const T& old_value, const T& new_value)
442
 
        {
443
 
                while(first != last){
444
 
                        if(*first == old_value){
445
 
                                *first = new_value;
446
 
                        }
447
 
                        ++first;
448
 
                }
449
 
        }
450
 
 
451
 
        template<class ForwardIterator, class Predicate, class T> _UCXXEXPORT
452
 
                void
453
 
                replace_if(ForwardIterator first, ForwardIterator last,
454
 
                        Predicate pred, const T& new_value)
455
 
        {
456
 
                while(first != last){
457
 
                        if( pred(*first) ){
458
 
                                *first = new_value;
459
 
                        }
460
 
                        ++first;
461
 
                }
462
 
        }
463
 
 
464
 
 
465
 
        template<class InputIterator, class OutputIterator, class T> _UCXXEXPORT
466
 
                OutputIterator
467
 
                replace_copy(InputIterator first, InputIterator last,
468
 
                        OutputIterator result, const T& old_value, const T& new_value)
469
 
        {
470
 
                while(first != last){
471
 
                        if(*first == old_value){
472
 
                                *result = new_value;
473
 
                        }else{
474
 
                                *result = *first;
475
 
                        }
476
 
                        ++first;
477
 
                        ++result;
478
 
                }
479
 
                return result;
480
 
        }
481
 
 
482
 
        template<class Iterator, class OutputIterator, class Predicate, class T> _UCXXEXPORT
483
 
                OutputIterator
484
 
                replace_copy_if(Iterator first, Iterator last,
485
 
                        OutputIterator result, Predicate pred, const T& new_value)
486
 
        {
487
 
                while(first != last){
488
 
                        if( pred(*first) ){
489
 
                                *result = new_value;
490
 
                        }else{
491
 
                                *result = *first;
492
 
                        }
493
 
                        ++first;
494
 
                        ++result;
495
 
                }
496
 
                return result;
497
 
        }
498
 
 
499
 
        template<class ForwardIterator, class T> _UCXXEXPORT
500
 
                void
501
 
                fill(ForwardIterator first, ForwardIterator last, const T& value)
502
 
        {
503
 
                while(first != last){
504
 
                        *first = value;
505
 
                        ++first;
506
 
                }
507
 
        }
508
 
 
509
 
        template<class OutputIterator, class Size, class T> _UCXXEXPORT
510
 
                void
511
 
                fill_n(OutputIterator first, Size n, const T& value)
512
 
        {
513
 
                while(n != 0){
514
 
                        *first = value;
515
 
                        --n;
516
 
                        ++first;
517
 
                }
518
 
        }
519
 
 
520
 
        template<class ForwardIterator, class Generator> _UCXXEXPORT
521
 
                void
522
 
                generate(ForwardIterator first, ForwardIterator last, Generator gen)
523
 
        {
524
 
                while(first != last){
525
 
                        *first = gen();
526
 
                        ++first;
527
 
                }
528
 
        }
529
 
 
530
 
        template<class OutputIterator, class Size, class Generator> _UCXXEXPORT
531
 
                void
532
 
                generate_n(OutputIterator first, Size n, Generator gen)
533
 
        {
534
 
                while(n !=0){
535
 
                        *first = gen();
536
 
                        --n;
537
 
                        ++first;
538
 
                }
539
 
        }
540
 
 
541
 
        template<class ForwardIterator, class T> _UCXXEXPORT
542
 
                ForwardIterator
543
 
                remove(ForwardIterator first, ForwardIterator last, const T& value)
544
 
        {
545
 
                ForwardIterator temp(first);
546
 
                while(temp !=last){
547
 
                        if(*temp == value){
548
 
 
549
 
                        }else{
550
 
                                *first = *temp;
551
 
                                ++first;
552
 
                        }
553
 
                        ++temp;
554
 
                }
555
 
                return first;
556
 
        }
557
 
 
558
 
        template<class ForwardIterator, class Predicate> _UCXXEXPORT
559
 
                ForwardIterator
560
 
                remove_if(ForwardIterator first, ForwardIterator last, Predicate pred)
561
 
        {
562
 
                ForwardIterator temp(first);
563
 
                while(temp !=last){
564
 
                        if( pred(*temp) ){
565
 
 
566
 
                        }else{
567
 
                                *first = *temp;
568
 
                                ++first;
569
 
                        }
570
 
                        ++temp;
571
 
                }
572
 
                return first;
573
 
        }
574
 
 
575
 
 
576
 
        template<class InputIterator, class OutputIterator, class T> _UCXXEXPORT
577
 
                OutputIterator
578
 
                remove_copy(InputIterator first, InputIterator last,
579
 
                        OutputIterator result, const T& value)
580
 
        {
581
 
                while(first !=last){
582
 
                        if(*first != value){
583
 
                                *result = *first;
584
 
                                ++result;
585
 
                        }
586
 
                        ++first;
587
 
                }
588
 
                return result;
589
 
        }
590
 
 
591
 
        template<class InputIterator, class OutputIterator, class Predicate> _UCXXEXPORT
592
 
                OutputIterator
593
 
                remove_copy_if(InputIterator first, InputIterator last,
594
 
                        OutputIterator result, Predicate pred)
595
 
        {
596
 
                while(first !=last){
597
 
                        if( !pred(*first) ){
598
 
                                *result = *first;
599
 
                                ++result;
600
 
                        }
601
 
                        ++first;
602
 
                }
603
 
                return result;
604
 
        }
605
 
 
606
 
        template<class ForwardIterator> _UCXXEXPORT
607
 
                ForwardIterator
608
 
                unique(ForwardIterator first, ForwardIterator last)
609
 
        {
610
 
                ForwardIterator new_val(first);
611
 
                ForwardIterator old_val(first);
612
 
 
613
 
                while(new_val !=last){
614
 
                        if(*new_val == *old_val && new_val != old_val){
615
 
 
616
 
                        }else{
617
 
                                *first = *new_val;
618
 
                                old_val = new_val;
619
 
                                ++first;
620
 
                        }
621
 
                        ++new_val;
622
 
                }
623
 
                return first;
624
 
        }
625
 
 
626
 
        template<class ForwardIterator, class BinaryPredicate> _UCXXEXPORT
627
 
                ForwardIterator
628
 
                unique(ForwardIterator first, ForwardIterator last, BinaryPredicate pred)
629
 
        {
630
 
                ForwardIterator new_val(first);
631
 
                ForwardIterator old_val(first);
632
 
 
633
 
                while(new_val !=last){
634
 
                        if( pred(*new_val, *old_val) && new_val != old_val){
635
 
 
636
 
                        }else{
637
 
                                *first = *new_val;
638
 
                                old_val = new_val;
639
 
                                ++first;
640
 
                        }
641
 
                        ++new_val;
642
 
                }
643
 
                return first;
644
 
        }
645
 
 
646
 
        template<class InputIterator, class OutputIterator> _UCXXEXPORT
647
 
                OutputIterator
648
 
                unique_copy(InputIterator first, InputIterator last, OutputIterator result)
649
 
        {
650
 
                InputIterator temp(first);
651
 
                while(first !=last ){
652
 
                        if(*first == *temp && first != temp){
653
 
 
654
 
                        }else{
655
 
                                *result = *first;
656
 
                                temp = first;
657
 
                                ++result;
658
 
                        }
659
 
                        ++first;
660
 
                }
661
 
                return result;
662
 
        }
663
 
 
664
 
        template<class InputIterator, class OutputIterator, class BinaryPredicate> _UCXXEXPORT
665
 
                OutputIterator
666
 
                unique_copy(InputIterator first, InputIterator last,
667
 
                        OutputIterator result, BinaryPredicate pred)
668
 
        {
669
 
                InputIterator temp(first);
670
 
                while(first !=last ){
671
 
                        if( pred(*first, *temp) && first != temp){
672
 
 
673
 
                        }else{
674
 
                                *result = *first;
675
 
                                temp = first;
676
 
                                ++result;
677
 
                        }
678
 
                        ++first;
679
 
                }
680
 
                return result;
681
 
        }
682
 
 
683
 
        template<class BidirectionalIterator> _UCXXEXPORT
684
 
                void
685
 
                reverse(BidirectionalIterator first, BidirectionalIterator last)
686
 
        {
687
 
                --last;         //Don't work with one past the end element
688
 
                while(first < last){
689
 
                        iter_swap(first, last);
690
 
                        ++first;
691
 
                        --last;
692
 
                }
693
 
        }
694
 
 
695
 
        template<class BidirectionalIterator, class OutputIterator> _UCXXEXPORT
696
 
                OutputIterator
697
 
                reverse_copy(BidirectionalIterator first, BidirectionalIterator last,
698
 
                OutputIterator result)
699
 
        {
700
 
                while(last > first){
701
 
                        --last;
702
 
                        *result = *last;
703
 
                        ++result;
704
 
                }
705
 
                return result;
706
 
        }
707
 
 
708
 
        template<class ForwardIterator> _UCXXEXPORT
709
 
                void
710
 
                rotate(ForwardIterator first, ForwardIterator middle, ForwardIterator last)
711
 
        {
712
 
                if ((first == middle) || (last  == middle)){
713
 
                        return;
714
 
                }
715
 
 
716
 
                ForwardIterator first2 = middle;
717
 
 
718
 
                do {
719
 
                        swap(*first, *first2);
720
 
                        first++;
721
 
                        first2++;
722
 
                        if(first == middle){
723
 
                                middle = first2;
724
 
                        }
725
 
                } while(first2 != last);
726
 
 
727
 
                first2 = middle;
728
 
 
729
 
                while (first2 != last) {
730
 
                        swap(*first, *first2);
731
 
                        first++;
732
 
                        first2++;
733
 
                        if (first == middle){
734
 
                                middle = first2;
735
 
                        }else if (first2 == last){
736
 
                                first2 = middle;
737
 
                        }
738
 
                }
739
 
        }
740
 
 
741
 
        template<class ForwardIterator, class OutputIterator> _UCXXEXPORT
742
 
                OutputIterator
743
 
                rotate_copy(ForwardIterator first, ForwardIterator middle,
744
 
                        ForwardIterator last, OutputIterator result)
745
 
        {
746
 
                ForwardIterator temp(middle);
747
 
                while(temp !=last){
748
 
                        *result = *temp;
749
 
                        ++temp;
750
 
                        ++result;
751
 
                }
752
 
                while(first != middle){
753
 
                        *result = *first;
754
 
                        ++first;
755
 
                        ++result;
756
 
                }
757
 
                return result;
758
 
        }
759
 
 
760
 
 
761
 
        template<class RandomAccessIterator> _UCXXEXPORT
762
 
                void
763
 
                random_shuffle(RandomAccessIterator first, RandomAccessIterator last)
764
 
        {
765
 
                --last;
766
 
                while(first != last){
767
 
                        iter_swap(first, (first + (rand() % (last - first) ) ) );
768
 
                        ++first;
769
 
                }
770
 
        }
771
 
 
772
 
        template<class RandomAccessIterator, class RandomNumberGenerator> _UCXXEXPORT
773
 
                void
774
 
                random_shuffle(RandomAccessIterator first, RandomAccessIterator last,
775
 
                RandomNumberGenerator& rand)
776
 
        {
777
 
                --last;
778
 
                while(first != last){
779
 
                        iter_swap(first, (first + (rand(last - first) ) ) );
780
 
                        ++first;
781
 
                }
782
 
        }
783
 
 
784
 
        // _lib.alg.partitions_, partitions:
785
 
        template<class BidirectionalIterator, class Predicate> _UCXXEXPORT BidirectionalIterator
786
 
                partition(BidirectionalIterator first, BidirectionalIterator last, Predicate pred)
787
 
        {
788
 
                return stable_partition(first, last, pred);
789
 
        }
790
 
 
791
 
        template<class BidirectionalIterator, class Predicate> _UCXXEXPORT BidirectionalIterator
792
 
                stable_partition(BidirectionalIterator first, BidirectionalIterator last, Predicate pred)
793
 
        {
794
 
                //first now points to the first non-desired element
795
 
                while( first != last && pred(*first) ){
796
 
                        ++first;
797
 
                }
798
 
 
799
 
                BidirectionalIterator tempb;
800
 
                BidirectionalIterator tempe = first;
801
 
 
802
 
                while( tempe != last ){
803
 
                        //Find the next desired element
804
 
                        while( !pred(*tempe) ){
805
 
                                ++tempe;
806
 
 
807
 
                                //See if we should exit
808
 
                                if(tempe == last){
809
 
                                        return first;
810
 
                                }
811
 
                        }
812
 
 
813
 
                        //Rotate the element back to the begining
814
 
                        tempb = tempe;
815
 
                        while(tempb != first){
816
 
                                iter_swap(tempb, tempb-1 );
817
 
                                --tempb;
818
 
                        }
819
 
                        //First now has a desired element
820
 
                        ++first;
821
 
                }
822
 
 
823
 
                return first;
824
 
        }
825
 
 
826
 
        template<class RandomAccessIterator> _UCXXEXPORT
827
 
                void sort(RandomAccessIterator first, RandomAccessIterator last)
828
 
        {
829
 
                less<typename iterator_traits<RandomAccessIterator>::value_type> c;
830
 
                sort(first, last, c );
831
 
        }
832
 
 
833
 
        template<class RandomAccessIterator, class Compare> _UCXXEXPORT
834
 
                void sort(RandomAccessIterator first, RandomAccessIterator last, Compare comp)
835
 
        {
836
 
                stable_sort(first, last, comp);
837
 
        }
838
 
 
839
 
        template<class RandomAccessIterator> _UCXXEXPORT
840
 
                void stable_sort(RandomAccessIterator first, RandomAccessIterator last)
841
 
        {
842
 
                less<typename iterator_traits<RandomAccessIterator>::value_type> c;
843
 
                stable_sort(first, last, c);
844
 
        }
845
 
 
846
 
        template<class RandomAccessIterator, class Compare> _UCXXEXPORT
847
 
                void stable_sort(RandomAccessIterator first, RandomAccessIterator last, Compare comp)
848
 
        {
849
 
                //FIXME - bubble sort
850
 
                RandomAccessIterator temp;
851
 
                --last;
852
 
                while(last - first > 0){
853
 
                        temp = last;
854
 
                        while(temp != first){
855
 
                                if( comp( *temp, *(temp-1) ) ){
856
 
                                        iter_swap( temp-1, temp);
857
 
                                }
858
 
                                --temp;
859
 
                        }
860
 
                        ++first;
861
 
                }
862
 
        }
863
 
 
864
 
        template<class RandomAccessIterator> _UCXXEXPORT
865
 
                void partial_sort(RandomAccessIterator first, RandomAccessIterator middle, RandomAccessIterator last)
866
 
        {
867
 
                less<typename iterator_traits<RandomAccessIterator>::value_type> c;
868
 
                partial_sort(first, middle, last, c);
869
 
        }
870
 
        template<class RandomAccessIterator, class Compare> _UCXXEXPORT
871
 
                void partial_sort(RandomAccessIterator first, RandomAccessIterator middle,
872
 
                        RandomAccessIterator last, Compare comp)
873
 
        {
874
 
                //Fixme - partial bubble sort
875
 
                RandomAccessIterator temp;
876
 
                --last;
877
 
                while(first < middle){
878
 
                        temp = last;
879
 
                        while(temp != first){
880
 
                                if( comp(*temp, *(temp -1 ) ) ) {
881
 
                                        iter_swap( temp-1, temp);
882
 
                                }
883
 
                                --temp;
884
 
                        }
885
 
                        ++first;
886
 
                }
887
 
        }
888
 
        template<class InputIterator, class RandomAccessIterator> _UCXXEXPORT
889
 
                RandomAccessIterator
890
 
                partial_sort_copy(InputIterator first, InputIterator last,
891
 
                        RandomAccessIterator result_first, RandomAccessIterator result_last)
892
 
        {
893
 
                less<typename iterator_traits<RandomAccessIterator>::value_type> c;
894
 
                return partial_sort_copy(first, last, result_first, result_last, c);
895
 
        }
896
 
        template<class InputIterator, class RandomAccessIterator, class Compare> _UCXXEXPORT
897
 
                RandomAccessIterator
898
 
                partial_sort_copy(InputIterator first, InputIterator last,
899
 
                        RandomAccessIterator result_first, RandomAccessIterator result_last, Compare comp)
900
 
        {
901
 
                size_t output_size = last-first;
902
 
                size_t temp_size = result_last - result_first;
903
 
                if(temp_size < output_size){
904
 
                        output_size = temp_size;
905
 
                }
906
 
 
907
 
                RandomAccessIterator middle = result_first + output_size;
908
 
                RandomAccessIterator temp = result_first;
909
 
 
910
 
                while(temp != middle){
911
 
                        *temp = *first;
912
 
                        ++temp;
913
 
                        ++first;
914
 
                }
915
 
                sort(result_first, middle, comp);
916
 
 
917
 
                while( first != last){
918
 
                        if( comp( *first, *(middle-1) ) ){
919
 
                                *(middle-1) = *first;
920
 
                                sort(result_first, middle, comp);
921
 
                        }
922
 
                        ++first;
923
 
                }
924
 
 
925
 
                return middle;
926
 
        }
927
 
        template<class RandomAccessIterator> _UCXXEXPORT
928
 
                void nth_element(RandomAccessIterator first, RandomAccessIterator nth, RandomAccessIterator last)
929
 
        {
930
 
                less<typename iterator_traits<RandomAccessIterator>::value_type> c;
931
 
                nth_element(first, nth, last, c);
932
 
        }
933
 
        template<class RandomAccessIterator, class Compare> _UCXXEXPORT
934
 
                void nth_element(RandomAccessIterator first, RandomAccessIterator nth,
935
 
                        RandomAccessIterator last, Compare comp)
936
 
        {
937
 
                partial_sort(first, nth, last, comp);
938
 
        }
939
 
 
940
 
        template<class ForwardIterator, class T> _UCXXEXPORT
941
 
                ForwardIterator lower_bound(ForwardIterator first, ForwardIterator last,
942
 
                        const T& value)
943
 
        {
944
 
                less<typename iterator_traits<ForwardIterator>::value_type> c;
945
 
                return lower_bound(first, last, value, c);
946
 
        }
947
 
 
948
 
        template<class ForwardIterator, class T, class Compare> _UCXXEXPORT
949
 
                ForwardIterator lower_bound(ForwardIterator first, ForwardIterator last,
950
 
                        const T& value, Compare comp)
951
 
        {
952
 
                if(first == last){
953
 
                        return last;
954
 
                }
955
 
                //If below or equal to the first element
956
 
                if( comp(*first, value) == false){
957
 
                        return first;
958
 
                }
959
 
                ForwardIterator middle;
960
 
 
961
 
                //If greater than the top element
962
 
                middle = first;
963
 
                advance(middle, distance(first, last) - 1);
964
 
                if( comp(*middle, value) ){
965
 
                        return last;
966
 
                }
967
 
 
968
 
                //Now begin the actual search for the begining
969
 
                while( distance(first, last) > 1 ){
970
 
                        //Find middle
971
 
                        middle = first;
972
 
                        advance(middle, (distance(first, last)/2) );
973
 
                        if( !comp(*middle, value) ){
974
 
                                last = middle;
975
 
                        }else{
976
 
                                first = middle;
977
 
                        }
978
 
                }
979
 
        
980
 
                if( !comp(*first, value) ){
981
 
                        return first;
982
 
                }
983
 
                return last;
984
 
        }
985
 
 
986
 
        template<class ForwardIterator, class T> _UCXXEXPORT
987
 
                ForwardIterator upper_bound(ForwardIterator first, ForwardIterator last,
988
 
                        const T& value)
989
 
        {
990
 
                less<typename iterator_traits<ForwardIterator>::value_type> c;
991
 
                return upper_bound(first, last, value, c);
992
 
        }
993
 
 
994
 
 
995
 
        template<class ForwardIterator, class T, class Compare> _UCXXEXPORT
996
 
                ForwardIterator upper_bound(ForwardIterator first, ForwardIterator last,
997
 
                        const T& value, Compare comp)
998
 
        {
999
 
                if(first == last){
1000
 
                        return last;
1001
 
                }
1002
 
                //If below the first element
1003
 
                if( comp(value, *first) == true){
1004
 
                        return first;
1005
 
                }
1006
 
 
1007
 
                ForwardIterator middle;
1008
 
 
1009
 
                //If greater than the top element
1010
 
                middle = first;
1011
 
                advance(middle, distance(first, last) - 1);
1012
 
                if( comp(*middle, value) ){
1013
 
                        return last;
1014
 
                }
1015
 
 
1016
 
                //Now begin the actual search for the end
1017
 
                while( distance(first, last) > 1 ){
1018
 
                        //Find middle
1019
 
                        middle = first;
1020
 
                        advance(middle, (distance(first, last)/2) );
1021
 
                        if( comp(value, *middle) ){
1022
 
                                last = middle;
1023
 
                        }else{
1024
 
                                first = middle;
1025
 
                        }
1026
 
                }
1027
 
        
1028
 
                if( comp(value, *first) ){
1029
 
                        return first;
1030
 
                }
1031
 
                return last;
1032
 
        }
1033
 
 
1034
 
 
1035
 
        template<class ForwardIterator, class T> _UCXXEXPORT
1036
 
                pair<ForwardIterator, ForwardIterator>
1037
 
                equal_range(ForwardIterator first, ForwardIterator last, const T& value)
1038
 
        {
1039
 
                less<typename iterator_traits<ForwardIterator>::value_type> c;
1040
 
                return equal_range(first, last, value, c);
1041
 
        }
1042
 
 
1043
 
        template<class ForwardIterator, class T, class Compare> _UCXXEXPORT
1044
 
                pair<ForwardIterator, ForwardIterator>
1045
 
                equal_range(ForwardIterator first, ForwardIterator last, const T& value, Compare comp)
1046
 
        {
1047
 
                pair<ForwardIterator, ForwardIterator> retval;
1048
 
                retval.first = lower_bound(first, last, value, comp);
1049
 
                //Reuse retval.first in order to (possibly) save a few comparisons
1050
 
                retval.second = upper_bound(retval.first, last, value, comp);
1051
 
                return retval;
1052
 
 
1053
 
        }
1054
 
 
1055
 
        template<class ForwardIterator, class T> _UCXXEXPORT
1056
 
                bool binary_search(ForwardIterator first, ForwardIterator last,
1057
 
                        const T& value)
1058
 
        {
1059
 
                less<typename iterator_traits<ForwardIterator>::value_type> c;
1060
 
                return binary_search(first, last, value, c);
1061
 
        }
1062
 
 
1063
 
        template<class ForwardIterator, class T, class Compare> _UCXXEXPORT
1064
 
                bool binary_search(ForwardIterator first, ForwardIterator last,
1065
 
                        const T& value, Compare comp)
1066
 
        {
1067
 
                if( distance(first, last) == 0){
1068
 
                        return false;
1069
 
                }
1070
 
 
1071
 
                ForwardIterator middle;
1072
 
 
1073
 
                while( distance(first, last) > 1 ){
1074
 
                        //Set middle between first and last
1075
 
                        middle = first;
1076
 
                        advance(middle, distance(first, last)/2 );
1077
 
 
1078
 
                        if( comp(*middle, value ) == true){
1079
 
                                first = middle;
1080
 
                        }else{
1081
 
                                last = middle;
1082
 
                        }
1083
 
                }
1084
 
 
1085
 
                if( !comp(*first, value) && !comp(value, *first) ){
1086
 
                        return true;
1087
 
                }       
1088
 
                if( !comp(*last, value) && !comp(value, *last) ){
1089
 
                        return true;
1090
 
                }       
1091
 
 
1092
 
                return false;
1093
 
        }
1094
 
 
1095
 
        // _lib.alg.merge_, merge:
1096
 
        template<class InputIterator1, class InputIterator2, class OutputIterator> _UCXXEXPORT
1097
 
                OutputIterator
1098
 
                merge(InputIterator1 first1, InputIterator1 last1,
1099
 
                        InputIterator2 first2, InputIterator2 last2, OutputIterator result)
1100
 
        {
1101
 
                less<typename iterator_traits<InputIterator1>::value_type> c;
1102
 
                return merge(first1, last1, first2, last2, result, c);
1103
 
        }
1104
 
        template<class InputIterator1, class InputIterator2, class OutputIterator, class Compare> _UCXXEXPORT
1105
 
                OutputIterator
1106
 
                merge(InputIterator1 first1, InputIterator1 last1,
1107
 
                        InputIterator2 first2, InputIterator2 last2, OutputIterator result, Compare comp)
1108
 
        {
1109
 
                while( first1 != last1 && first2 != last2){
1110
 
                        //If in this order so first1 elements which are equal come first
1111
 
                        if( comp(*first2, *first1) ){
1112
 
                                *result = *first2;
1113
 
                                ++first2;
1114
 
                        }else{
1115
 
                                *result = *first1;
1116
 
                                ++first1;
1117
 
                        }
1118
 
                        ++result;
1119
 
                }
1120
 
                //Clean up remaining elements
1121
 
                while(first1 != last1){
1122
 
                        *result = *first1;
1123
 
                        ++first1;
1124
 
                        ++result;
1125
 
                }
1126
 
                while(first2 != last2){
1127
 
                        *result = *first2;
1128
 
                        ++first2;
1129
 
                        ++result;
1130
 
                }
1131
 
                return result;
1132
 
        }
1133
 
        template<class BidirectionalIterator> _UCXXEXPORT
1134
 
                void inplace_merge(BidirectionalIterator first,
1135
 
                        BidirectionalIterator middle, BidirectionalIterator last)
1136
 
        {
1137
 
                less<typename iterator_traits<BidirectionalIterator>::value_type> c;
1138
 
                inplace_merge(first, middle, last, c);
1139
 
        }
1140
 
        template<class BidirectionalIterator, class Compare> _UCXXEXPORT
1141
 
                void inplace_merge(BidirectionalIterator first,
1142
 
                        BidirectionalIterator middle, BidirectionalIterator last, Compare comp)
1143
 
        {
1144
 
                //FIXME - using a bubble exchange method
1145
 
                while(first != middle && middle !=last){
1146
 
                        if( comp(*middle, *first) ){
1147
 
                                BidirectionalIterator temp(middle);
1148
 
                                while( temp != first){
1149
 
                                        iter_swap(temp, temp-1);
1150
 
                                        --temp;
1151
 
                                }
1152
 
                                ++middle;
1153
 
                        }
1154
 
                        ++first;
1155
 
                }
1156
 
        }
1157
 
 
1158
 
        // _lib.alg.set.operations_, set operations:
1159
 
        template<class InputIterator1, class InputIterator2> _UCXXEXPORT
1160
 
                bool includes(InputIterator1 first1, InputIterator1 last1,
1161
 
                        InputIterator2 first2, InputIterator2 last2)
1162
 
        {
1163
 
                less<typename iterator_traits<InputIterator1>::value_type> c;
1164
 
                return includes(first1, last1, first2, last2, c );
1165
 
        }
1166
 
 
1167
 
        template<class InputIterator1, class InputIterator2, class Compare> _UCXXEXPORT
1168
 
                bool includes(InputIterator1 first1, InputIterator1 last1,
1169
 
                        InputIterator2 first2, InputIterator2 last2, Compare comp)
1170
 
        {
1171
 
                while(first2 != last2){
1172
 
                        //Advance the large set until no longer smaller than the element we are looking for
1173
 
                        while( comp(*first1, *first2) ){
1174
 
                                ++first1;
1175
 
                                //If we are at the end of the list, we don't have the desired element
1176
 
                                if(first1 == last1){
1177
 
                                        return false;
1178
 
                                }
1179
 
                        }
1180
 
                        //If we are past the element we want, it isn't here
1181
 
                        if( comp(*first2, *first1) ){
1182
 
                                return false;
1183
 
                        }
1184
 
                        ++first2;       //Move to next element
1185
 
                }
1186
 
                return true;
1187
 
        }
1188
 
 
1189
 
        template<class InputIterator1, class InputIterator2, class OutputIterator> _UCXXEXPORT
1190
 
                OutputIterator set_union(InputIterator1 first1, InputIterator1 last1,
1191
 
                        InputIterator2 first2, InputIterator2 last2, OutputIterator result)
1192
 
        {
1193
 
                less<typename iterator_traits<InputIterator1>::value_type> c;
1194
 
                return set_union(first1, last1, first2, last2, result, c);
1195
 
        }
1196
 
 
1197
 
        template<class InputIterator1, class InputIterator2, class OutputIterator, class Compare> _UCXXEXPORT
1198
 
                OutputIterator set_union(InputIterator1 first1, InputIterator1 last1,
1199
 
                        InputIterator2 first2, InputIterator2 last2, OutputIterator result, Compare comp)
1200
 
        {
1201
 
                while( first1 != last1 && first2 != last2){
1202
 
                        if( comp(*first2, *first1) ){
1203
 
                                *result = *first2;
1204
 
 
1205
 
                                //Elliminate duplicates
1206
 
                                InputIterator2 temp = first2;
1207
 
                                while( first1 != last1 && !comp( *temp, *first1) ){
1208
 
                                        ++first1;
1209
 
                                }
1210
 
                                while( first2 != last2 && !comp( *temp, *first2) ){
1211
 
                                        ++first2;
1212
 
                                }
1213
 
                        }else{
1214
 
                                *result = *first1;
1215
 
                                        //Elliminate duplicates
1216
 
                                InputIterator1 temp = first1;
1217
 
                                while( first1 != last1 && !comp( *temp, *first1) ){
1218
 
                                        ++first1;
1219
 
                                }
1220
 
                                while( first2 != last2 && !comp( *temp, *first2) ){
1221
 
                                        ++first2;
1222
 
                                }
1223
 
                        }
1224
 
                        ++result;
1225
 
                }
1226
 
 
1227
 
                //Clean up remaining elements
1228
 
                while(first1 != last1){
1229
 
                        *result = *first1;
1230
 
                        ++result;
1231
 
                        InputIterator1 temp = first1;
1232
 
                        while( first1 != last1 && !comp( *temp, *first1) ){
1233
 
                                ++first1;
1234
 
                        }
1235
 
                }
1236
 
 
1237
 
                while(first2 != last2){
1238
 
                        *result = *first2;
1239
 
                        ++result;
1240
 
                        InputIterator2 temp = first2;
1241
 
                        while( first2 != last2 && !comp( *temp, *first2) ){
1242
 
                                ++first2;
1243
 
                        }
1244
 
                }
1245
 
                return result;
1246
 
        }
1247
 
 
1248
 
        template<class InputIterator1, class InputIterator2, class OutputIterator> _UCXXEXPORT
1249
 
                OutputIterator set_intersection(InputIterator1 first1, InputIterator1 last1,
1250
 
                        InputIterator2 first2, InputIterator2 last2, OutputIterator result)
1251
 
        {
1252
 
                less<typename iterator_traits<InputIterator1>::value_type> c;
1253
 
                return set_intersection(first1, last1, first2, last2, result, c);
1254
 
        }
1255
 
        template<class InputIterator1, class InputIterator2, class OutputIterator, class Compare> _UCXXEXPORT
1256
 
                OutputIterator set_intersection(InputIterator1 first1, InputIterator1 last1,
1257
 
                        InputIterator2 first2, InputIterator2 last2, OutputIterator result, Compare comp)
1258
 
        {
1259
 
                while( first1 != last1 && first2 != last2){
1260
 
                        if( comp(*first2, *first1) ){
1261
 
                                ++first2;
1262
 
                        }else if( comp(*first1, *first2) ) {
1263
 
                                ++first1;
1264
 
                        }else{
1265
 
                                *result = *first1;
1266
 
                                ++result;
1267
 
                                InputIterator1 temp = first1;
1268
 
                                while( first1 != last1 && !comp( *temp, *first1) ){
1269
 
                                        ++first1;
1270
 
                                }
1271
 
                                ++first2;
1272
 
                        }
1273
 
                }
1274
 
                return result;
1275
 
        }
1276
 
        template<class InputIterator1, class InputIterator2, class OutputIterator> _UCXXEXPORT
1277
 
                OutputIterator set_difference(InputIterator1 first1, InputIterator1 last1,
1278
 
                        InputIterator2 first2, InputIterator2 last2, OutputIterator result)
1279
 
        {
1280
 
                less<typename iterator_traits<InputIterator1>::value_type> c;
1281
 
                return set_difference(first1, last1, first2, last2, result, c);
1282
 
        }
1283
 
        template<class InputIterator1, class InputIterator2, class OutputIterator, class Compare> _UCXXEXPORT
1284
 
                OutputIterator set_difference(InputIterator1 first1, InputIterator1 last1,
1285
 
                        InputIterator2 first2, InputIterator2 last2, OutputIterator result, Compare comp)
1286
 
        {
1287
 
                while( first1 != last1 && first2 != last2){
1288
 
                        if( comp(*first1, *first2) ){
1289
 
                                *result = *first1;
1290
 
                                ++result;
1291
 
 
1292
 
                                //Elliminate duplicates
1293
 
                                InputIterator1 temp = first1;
1294
 
                                while( first1 != last1 && !comp( *temp, *first1) ){
1295
 
                                        ++first1;
1296
 
                                }
1297
 
                        }else if( comp(*first2, *first1) ){
1298
 
 
1299
 
                                //Elliminate duplicates
1300
 
                                InputIterator2 temp = first2;
1301
 
                                while( first2 != last2 && !comp( *temp, *first2) ){
1302
 
                                        ++first2;
1303
 
                                }
1304
 
                
1305
 
                        }else{  //They are identical - skip
1306
 
                                //Elliminate duplicates
1307
 
                                InputIterator1 temp = first1;
1308
 
                                while( first1 != last1 && !comp( *temp, *first1) ){
1309
 
                                        ++first1;
1310
 
                                }
1311
 
                                while( first2 != last2 && !comp( *temp, *first2) ){
1312
 
                                        ++first2;
1313
 
                                }
1314
 
                        }
1315
 
                }
1316
 
 
1317
 
                //Clean up remaining elements
1318
 
                while(first1 != last1){
1319
 
                        *result = *first1;
1320
 
                        ++result;
1321
 
                        InputIterator1 temp = first1;
1322
 
                        while( first1 != last1 && !comp( *temp, *first1) ){
1323
 
                                ++first1;
1324
 
                        }
1325
 
                }
1326
 
 
1327
 
                return result;
1328
 
        }
1329
 
        template<class InputIterator1, class InputIterator2, class OutputIterator> _UCXXEXPORT
1330
 
                OutputIterator set_symmetric_difference(InputIterator1 first1, InputIterator1 last1,
1331
 
                        InputIterator2 first2, InputIterator2 last2, OutputIterator result)
1332
 
        {
1333
 
                less<typename iterator_traits<InputIterator1>::value_type> c;
1334
 
                return set_symmetric_difference(first1, last1, first2, last2, result, c);
1335
 
        }
1336
 
        template<class InputIterator1, class InputIterator2, class OutputIterator, class Compare> _UCXXEXPORT
1337
 
                OutputIterator set_symmetric_difference(InputIterator1 first1, InputIterator1 last1,
1338
 
                        InputIterator2 first2, InputIterator2 last2, OutputIterator result, Compare comp)
1339
 
        {
1340
 
                while( first1 != last1 && first2 != last2){
1341
 
                        if( comp(*first1, *first2) ){
1342
 
                                *result = *first1;
1343
 
                                ++result;
1344
 
 
1345
 
                                //Elliminate duplicates
1346
 
                                InputIterator1 temp = first1;
1347
 
                                while( first1 != last1 && !comp( *temp, *first1) ){
1348
 
                                        ++first1;
1349
 
                                }
1350
 
                        }else if( comp(*first2, *first1) ){
1351
 
                                *result = *first2;
1352
 
                                ++result;
1353
 
 
1354
 
                                //Elliminate duplicates
1355
 
                                InputIterator2 temp = first2;
1356
 
                                while( first2 != last2 && !comp( *temp, *first2) ){
1357
 
                                        ++first2;
1358
 
                                }
1359
 
                
1360
 
                        }else{  //They are identical - skip
1361
 
                                //Elliminate duplicates
1362
 
                                InputIterator1 temp = first1;
1363
 
                                while( first1 != last1 && !comp( *temp, *first1) ){
1364
 
                                        ++first1;
1365
 
                                }
1366
 
                                while( first2 != last2 && !comp( *temp, *first2) ){
1367
 
                                        ++first2;
1368
 
                                }
1369
 
                        }
1370
 
                }
1371
 
 
1372
 
                //Clean up remaining elements
1373
 
                while(first1 != last1){
1374
 
                        *result = *first1;
1375
 
                        ++result;
1376
 
                        InputIterator1 temp = first1;
1377
 
                        while( first1 != last1 && !comp( *temp, *first1) ){
1378
 
                                ++first1;
1379
 
                        }
1380
 
                }
1381
 
 
1382
 
                while(first2 != last2){
1383
 
                        *result = *first2;
1384
 
                        ++result;
1385
 
                        InputIterator2 temp = first2;
1386
 
                        while( first2 != last2 && !comp( *temp, *first2) ){
1387
 
                                ++first2;
1388
 
                        }
1389
 
                }
1390
 
 
1391
 
                return result;
1392
 
        }
1393
 
 
1394
 
        // _lib.alg.heap.operations_, heap operations:
1395
 
 
1396
 
        template<class RandomAccessIterator> _UCXXEXPORT
1397
 
                void push_heap(RandomAccessIterator first, RandomAccessIterator last)
1398
 
        {
1399
 
                less<typename iterator_traits<RandomAccessIterator>::value_type> c;
1400
 
                push_heap(first, last, c);
1401
 
        }
1402
 
 
1403
 
        template<class RandomAccessIterator, class Compare> _UCXXEXPORT
1404
 
                void push_heap(RandomAccessIterator first, RandomAccessIterator last, Compare comp)
1405
 
        {
1406
 
                // *(last - 1) is the last element
1407
 
                RandomAccessIterator begin, middle, end;
1408
 
                begin = first;
1409
 
                end = last - 2;
1410
 
                if(last - first < 2){           //Empty heap
1411
 
                        return;
1412
 
                }
1413
 
                if ( comp(*(last-1), *end) ){   //Belongs past the end - already in place
1414
 
                        return;
1415
 
                }
1416
 
                while(end - begin > 1){
1417
 
                        middle = begin + ((end - begin)/2);
1418
 
                        if( comp(*middle, *(last-1) ) ){
1419
 
                                end = middle;
1420
 
                        }else{
1421
 
                                begin = middle;
1422
 
                        }
1423
 
                }
1424
 
 
1425
 
                if( comp(*begin, *(last-1)) ){
1426
 
                        middle = begin;
1427
 
                }else{
1428
 
                        middle = end;
1429
 
                }
1430
 
 
1431
 
                //Now all we do is swap the elements up to the begining
1432
 
                --last;
1433
 
 
1434
 
                while(last > middle){
1435
 
                        iter_swap(last, last-1);
1436
 
                        --last;
1437
 
                }
1438
 
        }
1439
 
 
1440
 
        template<class RandomAccessIterator> _UCXXEXPORT
1441
 
                void pop_heap(RandomAccessIterator first, RandomAccessIterator last)
1442
 
        {
1443
 
                less<typename iterator_traits<RandomAccessIterator>::value_type> c;
1444
 
                pop_heap(first, last, c);
1445
 
        }
1446
 
        template<class RandomAccessIterator, class Compare> _UCXXEXPORT
1447
 
                void pop_heap(RandomAccessIterator first, RandomAccessIterator last, Compare)
1448
 
        {
1449
 
                --last;
1450
 
                while(first < last){
1451
 
                        iter_swap( first, first+1);
1452
 
                        ++first;
1453
 
                }
1454
 
        }
1455
 
 
1456
 
        template<class RandomAccessIterator> _UCXXEXPORT
1457
 
                void make_heap(RandomAccessIterator first, RandomAccessIterator last)
1458
 
        {
1459
 
                less<typename iterator_traits<RandomAccessIterator>::value_type> c;
1460
 
                make_heap(first, last, c);
1461
 
        }
1462
 
        template<class RandomAccessIterator, class Compare> _UCXXEXPORT
1463
 
                void make_heap(RandomAccessIterator first, RandomAccessIterator last, Compare comp)
1464
 
        {
1465
 
                sort( reverse_iterator<RandomAccessIterator>(last), reverse_iterator<RandomAccessIterator>(first), comp);
1466
 
        }
1467
 
        template<class RandomAccessIterator> _UCXXEXPORT
1468
 
                void sort_heap(RandomAccessIterator first, RandomAccessIterator last)
1469
 
        {
1470
 
                less<typename iterator_traits<RandomAccessIterator>::value_type> c;
1471
 
                sort_heap(first, last, c);
1472
 
        }
1473
 
        template<class RandomAccessIterator, class Compare> _UCXXEXPORT
1474
 
                void sort_heap(RandomAccessIterator first, RandomAccessIterator last, Compare comp)
1475
 
        {
1476
 
                sort( first, last, comp);
1477
 
        }
1478
 
 
1479
 
 
1480
 
        // _lib.alg.min.max_, minimum and maximum:
1481
 
        template<class T> _UCXXEXPORT
1482
 
                const T& min(const T& a, const T& b)
1483
 
        {
1484
 
                if(b < a){
1485
 
                        return b;
1486
 
                }
1487
 
                return a;
1488
 
        }
1489
 
 
1490
 
        template<class T, class Compare> _UCXXEXPORT
1491
 
                const T& min(const T& a, const T& b, Compare comp)
1492
 
        {
1493
 
                if( comp(b, a) == true){
1494
 
                        return b;
1495
 
                }
1496
 
                return a;
1497
 
        }
1498
 
 
1499
 
        template<class T> _UCXXEXPORT
1500
 
                const T& max(const T& a, const T& b)
1501
 
        {
1502
 
                if( a < b){
1503
 
                        return b;
1504
 
                }
1505
 
                return a;
1506
 
        }
1507
 
 
1508
 
        template<class T, class Compare> _UCXXEXPORT
1509
 
                const T& max(const T& a, const T& b, Compare comp)
1510
 
        {
1511
 
                if( comp(a, b) ){
1512
 
                        return b;
1513
 
                }
1514
 
                return a;
1515
 
        }
1516
 
 
1517
 
        template<class ForwardIterator> _UCXXEXPORT
1518
 
                ForwardIterator min_element(ForwardIterator first, ForwardIterator last)
1519
 
        {
1520
 
                less<typename iterator_traits<ForwardIterator>::value_type> c;
1521
 
                return min_element(first, last, c);
1522
 
        }
1523
 
 
1524
 
        template<class ForwardIterator, class Compare> _UCXXEXPORT
1525
 
                ForwardIterator min_element(ForwardIterator first, ForwardIterator last, Compare comp)
1526
 
        {
1527
 
                ForwardIterator retval = first;
1528
 
                ++first;
1529
 
                while(first != last){
1530
 
                        if( comp( *first, *retval) ){
1531
 
                                retval = first;
1532
 
                        }
1533
 
                        ++first;
1534
 
                }
1535
 
                return retval;
1536
 
        }
1537
 
 
1538
 
        template<class ForwardIterator> _UCXXEXPORT
1539
 
                ForwardIterator max_element(ForwardIterator first, ForwardIterator last)
1540
 
        {
1541
 
                less<typename iterator_traits<ForwardIterator>::value_type> c;
1542
 
                return max_element(first, last, c);
1543
 
        }
1544
 
 
1545
 
        template<class ForwardIterator, class Compare> _UCXXEXPORT
1546
 
                ForwardIterator max_element(ForwardIterator first, ForwardIterator last, Compare comp)
1547
 
        {
1548
 
                ForwardIterator retval = first;
1549
 
                ++first;
1550
 
                while(first != last){
1551
 
                        if( comp( *retval, *first ) ){
1552
 
                                retval = first;
1553
 
                        }
1554
 
                        ++first;
1555
 
                }
1556
 
                return retval;
1557
 
        }
1558
 
 
1559
 
        template<class InputIterator1, class InputIterator2> _UCXXEXPORT
1560
 
                bool lexicographical_compare(InputIterator1 first1, InputIterator1 last1,
1561
 
                        InputIterator2 first2, InputIterator2 last2)
1562
 
        {
1563
 
                less<typename iterator_traits<InputIterator1>::value_type> c;
1564
 
                return lexicographical_compare(first1, last1, first2, last2, c);
1565
 
        }
1566
 
 
1567
 
        template<class InputIterator1, class InputIterator2, class Compare> _UCXXEXPORT
1568
 
                bool lexicographical_compare(InputIterator1 first1, InputIterator1 last1,
1569
 
                        InputIterator2 first2, InputIterator2 last2, Compare comp)
1570
 
        {
1571
 
                while(first1 != last1 && first2 != last2){
1572
 
                        if( comp(*first1, *first2) ){
1573
 
                                return true;
1574
 
                        }
1575
 
                        if( comp(*first2, *first1) ){
1576
 
                                return false;
1577
 
                        }
1578
 
                        ++first1;
1579
 
                        ++first2;
1580
 
                }
1581
 
                return first1==last1 && first2 != last2;
1582
 
        }
1583
 
 
1584
 
        // _lib.alg.permutation.generators_, permutations
1585
 
        template<class BidirectionalIterator> _UCXXEXPORT
1586
 
                bool next_permutation(BidirectionalIterator first, BidirectionalIterator last)
1587
 
        {
1588
 
                less<typename iterator_traits<BidirectionalIterator>::value_type> c;
1589
 
                return next_permutation(first, last, c);
1590
 
        }
1591
 
 
1592
 
        template<class BidirectionalIterator, class Compare> _UCXXEXPORT
1593
 
                bool next_permutation(BidirectionalIterator first, BidirectionalIterator last, Compare comp)
1594
 
        {
1595
 
                if(first == last){
1596
 
                        return false;   // No data
1597
 
                }
1598
 
                BidirectionalIterator i = last;
1599
 
                --i;
1600
 
                if(i == first){
1601
 
                        return false;   //Only one element
1602
 
                }
1603
 
                BidirectionalIterator ii = i;
1604
 
                --ii;
1605
 
                //If the last two items are in order, swap them and call it done
1606
 
                if( comp(*ii, *i) ){
1607
 
                        iter_swap(ii, i);
1608
 
                        return true;
1609
 
                }
1610
 
 
1611
 
 
1612
 
                //This part of the algorithm copied from G++ header files ver. 3.4.2
1613
 
                i = last;
1614
 
                --i;
1615
 
                for(;;){
1616
 
                        ii = i;
1617
 
                        --i;
1618
 
                        if ( comp(*i, *ii) ){
1619
 
                                BidirectionalIterator j = last;
1620
 
                                --j;
1621
 
                                while ( !comp(*i, *j)){
1622
 
                                        --j;
1623
 
                                }
1624
 
                                iter_swap(i, j);
1625
 
                                reverse(ii, last);
1626
 
                                return true;
1627
 
                        }
1628
 
                        if (i == first){
1629
 
                                reverse(first, last);
1630
 
                                return false;
1631
 
                        }
1632
 
                }
1633
 
                
1634
 
 
1635
 
        }
1636
 
 
1637
 
        template<class BidirectionalIterator> _UCXXEXPORT
1638
 
                bool prev_permutation(BidirectionalIterator first, BidirectionalIterator last)
1639
 
        {
1640
 
                less<typename iterator_traits<BidirectionalIterator>::value_type> c;
1641
 
                return prev_permutation(first, last, c);
1642
 
        }
1643
 
        
1644
 
        template<class BidirectionalIterator, class Compare> _UCXXEXPORT
1645
 
                bool prev_permutation(BidirectionalIterator first, BidirectionalIterator last, Compare comp)
1646
 
        {
1647
 
                if(first == last){
1648
 
                        return false;   // No data
1649
 
                }
1650
 
                BidirectionalIterator i = last;
1651
 
                --i;
1652
 
                if(i == first){
1653
 
                        return false;   //Only one element
1654
 
                }
1655
 
                BidirectionalIterator ii = i;
1656
 
                --ii;
1657
 
                //If the last two items are in reverse order, swap them and call it done
1658
 
                if( comp(*i, *ii) ){
1659
 
                        iter_swap(ii, i);
1660
 
                        return true;
1661
 
                }
1662
 
 
1663
 
                //Copied from GNU GCC header files version 3.4.2
1664
 
                i = last;
1665
 
                --i;
1666
 
 
1667
 
                for(;;){
1668
 
                        ii = i;
1669
 
                        --i;
1670
 
                        if ( comp(*ii, *i)){
1671
 
                                BidirectionalIterator j = last;
1672
 
                                --j;
1673
 
                                while ( !comp(*j, *i)){
1674
 
                                        --j;
1675
 
                                }
1676
 
                                iter_swap(i, j);
1677
 
                                reverse(ii, last);
1678
 
                                return true;
1679
 
                        }
1680
 
                        if (i == first){
1681
 
                                reverse(first, last);
1682
 
                                return false;
1683
 
                        }
1684
 
                }
1685
 
 
1686
 
        }
1687
 
 
1688
 
}
1689
 
 
1690
 
#pragma GCC visibility pop
1691
 
 
1692
 
#endif
1693
 
 
1694
 
 
1695