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
use core::cmp::min;
#[derive(Copy, Clone)]
pub struct RangeChunk { i: usize, n: usize, chunk: usize }
pub fn range_chunk(n: usize, chunk: usize) -> RangeChunk {
RangeChunk {
i: 0,
n: n,
chunk: chunk,
}
}
impl Iterator for RangeChunk {
type Item = (usize, usize);
#[inline]
fn next(&mut self) -> Option<Self::Item> {
if self.n == 0 {
None
} else {
let i = self.i;
let rem = min(self.n, self.chunk);
self.i += 1;
self.n -= rem;
Some((i, rem))
}
}
}
#[inline]
pub fn round_up_to(x: usize, multiple_of: usize) -> usize {
let (mut d, r) = (x / multiple_of, x % multiple_of);
if r > 0 { d += 1; }
d * multiple_of
}
impl RangeChunk {
#[cfg(feature="threading")]
pub(crate) fn part(self, index: usize, total: usize) -> Self {
debug_assert_eq!(self.i, 0, "range must be uniterated");
debug_assert_ne!(total, 0);
let (n, chunk) = (self.n, self.chunk);
let mut nchunks = n / chunk;
nchunks += (n % chunk != 0) as usize;
let mut chunks_per = nchunks / total;
chunks_per += (nchunks % total != 0) as usize;
let i = chunks_per * index;
let nn = min(n, (i + chunks_per) * chunk).saturating_sub(i * chunk);
RangeChunk { i, n: nn, chunk }
}
}