rdxsort 0.3.0

Fast Radix Sort
Documentation 
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

use template::RdxSortTemplate;

use std::cmp;

impl RdxSortTemplate for () {
    #[inline]
    fn cfg_nbuckets() -> usize {
        0
    }

    #[inline]
    fn cfg_nrounds() -> usize {
        0
    }

    #[inline]
    fn get_bucket(&self, _round: usize) -> usize {
        unreachable!()
    }

    #[inline]
    fn reverse(_round: usize, _bucket: usize) -> bool {
        unreachable!()
    }
}

impl<A> RdxSortTemplate for (A,) where A: RdxSortTemplate
{
    #[inline]
    fn cfg_nbuckets() -> usize {
        A::cfg_nbuckets()
    }

    #[inline]
    fn cfg_nrounds() -> usize {
        A::cfg_nrounds()
    }

    #[inline]
    fn get_bucket(&self, round: usize) -> usize {
        self.0.get_bucket(round)
    }

    #[inline]
    fn reverse(round: usize, bucket: usize) -> bool {
        A::reverse(round, bucket)
    }
}

impl<A, B> RdxSortTemplate for (A, B)
    where A: RdxSortTemplate,
          B: RdxSortTemplate
{
    #[inline]
    fn cfg_nbuckets() -> usize {
        cmp::max(A::cfg_nbuckets(), B::cfg_nbuckets())
    }

    #[inline]
    fn cfg_nrounds() -> usize {
        A::cfg_nrounds() + B::cfg_nrounds()
    }

    #[inline]
    fn get_bucket(&self, round: usize) -> usize {
        if round < B::cfg_nrounds() {
            self.1.get_bucket(round)
        } else {
            self.0.get_bucket(round - B::cfg_nrounds())
        }
    }

    #[inline]
    fn reverse(round: usize, bucket: usize) -> bool {
        if round < B::cfg_nrounds() {
            B::reverse(round, bucket)
        } else {
            A::reverse(round - B::cfg_nrounds(), bucket)
        }
    }
}

impl<A, B, C> RdxSortTemplate for (A, B, C)
    where A: RdxSortTemplate,
          B: RdxSortTemplate,
          C: RdxSortTemplate
{
    #[inline]
    fn cfg_nbuckets() -> usize {
        cmp::max(A::cfg_nbuckets(),
                 cmp::max(B::cfg_nbuckets(), C::cfg_nbuckets()))
    }

    #[inline]
    fn cfg_nrounds() -> usize {
        A::cfg_nrounds() + B::cfg_nrounds() + C::cfg_nrounds()
    }

    #[inline]
    fn get_bucket(&self, round: usize) -> usize {
        if round < C::cfg_nrounds() {
            self.2.get_bucket(round)
        } else if round < B::cfg_nrounds() + C::cfg_nrounds() {
            self.1.get_bucket(round - C::cfg_nrounds())
        } else {
            self.0.get_bucket(round - B::cfg_nrounds() - C::cfg_nrounds())
        }
    }

    #[inline]
    fn reverse(round: usize, bucket: usize) -> bool {
        if round < C::cfg_nrounds() {
            C::reverse(round, bucket)
        } else if round < B::cfg_nrounds() + C::cfg_nrounds() {
            B::reverse(round - C::cfg_nrounds(), bucket)
        } else {
            A::reverse(round - B::cfg_nrounds() - C::cfg_nrounds(), bucket)
        }
    }
}