sprs 0.6.2

#[macro_use]
extern crate bencher;
extern crate sprs;
extern crate ndarray;

use bencher::Bencher;
use sprs::{CsVec, CsMat};
use ndarray::{Array, Array2, ShapeBuilder};

fn sparse_dense_dotprod_default(bench: &mut Bencher) {
    let w = Array::range(0., 10., 0.00001);
    let x = CsVec::new(1000000, vec![0, 200000, 800000], vec![1., 2., 3.]);
    bench.iter(|| {
        x.dot(&w);
    });
}

fn sparse_dense_dotprod_specialized(bench: &mut Bencher) {
    let w = Array::range(0., 10., 0.00001);
    let x = CsVec::new(1000000, vec![0, 200000, 800000], vec![1., 2., 3.]);
    bench.iter(|| {
        x.dot_dense(w.view());
    });
}

fn sparse_dense_vec_matprod_default(bench: &mut Bencher) {
    let w = Array::range(0., 10., 0.00001);
    let a = CsMat::new((3, 1000000),
                       vec![0, 2, 4, 5],
                       vec![0, 1, 0, 2, 2],
                       vec![1., 2., 3., 4., 5.]);
    bench.iter(|| {
        &a * &w;
    });
}

fn sparse_dense_vec_matprod_specialized(bench: &mut Bencher) {
    let w = Array::range(0., 10., 0.00001);
    let a = CsMat::new((3, 1000000),
                       vec![0, 2, 4, 5],
                       vec![0, 1, 0, 2, 2],
                       vec![1., 2., 3., 4., 5.]);
    let rows = a.rows();
    let cols = w.shape()[0];
    let w_reshape = w.view().into_shape((1, cols)).unwrap();
    let w_t = w_reshape.t();
    let mut res = Array2::zeros((rows, 1).f());
    bench.iter(|| {
        sprs::prod::csr_mulacc_dense_colmaj(a.view(),
                                            w_t.view(),
                                            res.view_mut());
    });
}

benchmark_group!(benches,
                 sparse_dense_dotprod_default,
                 sparse_dense_dotprod_specialized,
                 sparse_dense_vec_matprod_specialized,
                 sparse_dense_vec_matprod_default);
benchmark_main!(benches);