numpy 0.17.2

PyO3-based Rust bindings of the NumPy C-API
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

#![feature(test)]
#![allow(deprecated)]

extern crate test;
use test::{black_box, Bencher};

use ndarray::Zip;
use numpy::{npyiter::NpyMultiIterBuilder, PyArray};
use pyo3::Python;

fn numpy_iter(bencher: &mut Bencher, size: usize) {
    Python::with_gil(|py| {
        let x = PyArray::<f64, _>::zeros(py, size, false);
        let y = PyArray::<f64, _>::zeros(py, size, false);
        let z = PyArray::<f64, _>::zeros(py, size, false);

        let x = x.readonly();
        let y = y.readonly();
        let mut z = z.readwrite();

        bencher.iter(|| {
            let iter = NpyMultiIterBuilder::new()
                .add_readonly(black_box(&x))
                .add_readonly(black_box(&y))
                .add_readwrite(black_box(&mut z))
                .build()
                .unwrap();

            for (x, y, z) in iter {
                *z = x + y;
            }
        });
    });
}

#[bench]
fn numpy_iter_small(bencher: &mut Bencher) {
    numpy_iter(bencher, 2_usize.pow(5));
}

#[bench]
fn numpy_iter_medium(bencher: &mut Bencher) {
    numpy_iter(bencher, 2_usize.pow(10));
}

#[bench]
fn numpy_iter_large(bencher: &mut Bencher) {
    numpy_iter(bencher, 2_usize.pow(15));
}

fn ndarray_iter(bencher: &mut Bencher, size: usize) {
    Python::with_gil(|py| {
        let x = PyArray::<f64, _>::zeros(py, size, false);
        let y = PyArray::<f64, _>::zeros(py, size, false);
        let z = PyArray::<f64, _>::zeros(py, size, false);

        let x = x.readonly();
        let y = y.readonly();
        let mut z = z.readwrite();

        bencher.iter(|| {
            Zip::from(black_box(x.as_array()))
                .and(black_box(y.as_array()))
                .and(black_box(z.as_array_mut()))
                .for_each(|x, y, z| {
                    *z = x + y;
                });
        });
    });
}

#[bench]
fn ndarray_iter_small(bencher: &mut Bencher) {
    ndarray_iter(bencher, 2_usize.pow(5));
}

#[bench]
fn ndarray_iter_medium(bencher: &mut Bencher) {
    ndarray_iter(bencher, 2_usize.pow(10));
}

#[bench]
fn ndarray_iter_large(bencher: &mut Bencher) {
    ndarray_iter(bencher, 2_usize.pow(15));
}