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
|
#![feature(test)]
extern crate test;
use ndarray::s;
use ndarray::IntoNdProducer;
use ndarray::{Array3, ShapeBuilder, Zip};
use test::{black_box, Bencher};
pub fn zip_copy<'a, A, P, Q>(data: P, out: Q)
where
P: IntoNdProducer<Item = &'a A>,
Q: IntoNdProducer<Item = &'a mut A, Dim = P::Dim>,
A: Copy + 'a,
{
Zip::from(data).and(out).for_each(|&i, o| {
*o = i;
});
}
pub fn zip_copy_split<'a, A, P, Q>(data: P, out: Q)
where
P: IntoNdProducer<Item = &'a A>,
Q: IntoNdProducer<Item = &'a mut A, Dim = P::Dim>,
A: Copy + 'a,
{
let z = Zip::from(data).and(out);
let (z1, z2) = z.split();
let (z11, z12) = z1.split();
let (z21, z22) = z2.split();
let f = |&i: &A, o: &mut A| *o = i;
z11.for_each(f);
z12.for_each(f);
z21.for_each(f);
z22.for_each(f);
}
pub fn zip_indexed(data: &Array3<f32>, out: &mut Array3<f32>)
{
Zip::indexed(data).and(out).for_each(|idx, &i, o| {
let _ = black_box(idx);
*o = i;
});
}
// array size in benchmarks
const SZ3: (usize, usize, usize) = (100, 110, 100);
#[bench]
fn zip_cc(b: &mut Bencher)
{
let data: Array3<f32> = Array3::zeros(SZ3);
let mut out = Array3::zeros(data.dim());
b.iter(|| zip_copy(&data, &mut out));
}
#[bench]
fn zip_cf(b: &mut Bencher)
{
let data: Array3<f32> = Array3::zeros(SZ3);
let mut out = Array3::zeros(data.dim().f());
b.iter(|| zip_copy(&data, &mut out));
}
#[bench]
fn zip_fc(b: &mut Bencher)
{
let data: Array3<f32> = Array3::zeros(SZ3.f());
let mut out = Array3::zeros(data.dim());
b.iter(|| zip_copy(&data, &mut out));
}
#[bench]
fn zip_ff(b: &mut Bencher)
{
let data: Array3<f32> = Array3::zeros(SZ3.f());
let mut out = Array3::zeros(data.dim().f());
b.iter(|| zip_copy(&data, &mut out));
}
#[bench]
fn zip_indexed_cc(b: &mut Bencher)
{
let data: Array3<f32> = Array3::zeros(SZ3);
let mut out = Array3::zeros(data.dim());
b.iter(|| zip_indexed(&data, &mut out));
}
#[bench]
fn zip_indexed_ff(b: &mut Bencher)
{
let data: Array3<f32> = Array3::zeros(SZ3.f());
let mut out = Array3::zeros(data.dim().f());
b.iter(|| zip_indexed(&data, &mut out));
}
#[bench]
fn slice_zip_cc(b: &mut Bencher)
{
let data: Array3<f32> = Array3::zeros(SZ3);
let mut out = Array3::zeros(data.dim());
let data = data.slice(s![1.., 1.., 1..]);
let mut out = out.slice_mut(s![1.., 1.., 1..]);
b.iter(|| zip_copy(&data, &mut out));
}
#[bench]
fn slice_zip_ff(b: &mut Bencher)
{
let data: Array3<f32> = Array3::zeros(SZ3.f());
let mut out = Array3::zeros(data.dim().f());
let data = data.slice(s![1.., 1.., 1..]);
let mut out = out.slice_mut(s![1.., 1.., 1..]);
b.iter(|| zip_copy(&data, &mut out));
}
#[bench]
fn slice_split_zip_cc(b: &mut Bencher)
{
let data: Array3<f32> = Array3::zeros(SZ3);
let mut out = Array3::zeros(data.dim());
let data = data.slice(s![1.., 1.., 1..]);
let mut out = out.slice_mut(s![1.., 1.., 1..]);
b.iter(|| zip_copy_split(&data, &mut out));
}
#[bench]
fn slice_split_zip_ff(b: &mut Bencher)
{
let data: Array3<f32> = Array3::zeros(SZ3.f());
let mut out = Array3::zeros(data.dim().f());
let data = data.slice(s![1.., 1.., 1..]);
let mut out = out.slice_mut(s![1.., 1.., 1..]);
b.iter(|| zip_copy_split(&data, &mut out));
}
|
