2018-02-02 19:26:35 +08:00
|
|
|
use na::{DMatrix, Hessenberg, Matrix4};
|
2018-11-07 01:32:20 +08:00
|
|
|
use test::{self, Bencher};
|
2017-08-03 01:37:44 +08:00
|
|
|
|
2018-02-02 19:26:35 +08:00
|
|
|
#[path = "../common/macros.rs"]
|
2017-08-03 01:37:44 +08:00
|
|
|
mod macros;
|
|
|
|
|
|
|
|
// Without unpack.
|
|
|
|
#[bench]
|
|
|
|
fn hessenberg_decompose_4x4(bh: &mut Bencher) {
|
|
|
|
let m = Matrix4::<f64>::new_random();
|
|
|
|
bh.iter(|| test::black_box(Hessenberg::new(m.clone())))
|
|
|
|
}
|
|
|
|
|
|
|
|
#[bench]
|
|
|
|
fn hessenberg_decompose_100x100(bh: &mut Bencher) {
|
|
|
|
let m = DMatrix::<f64>::new_random(100, 100);
|
|
|
|
bh.iter(|| test::black_box(Hessenberg::new(m.clone())))
|
|
|
|
}
|
|
|
|
|
|
|
|
#[bench]
|
|
|
|
fn hessenberg_decompose_200x200(bh: &mut Bencher) {
|
|
|
|
let m = DMatrix::<f64>::new_random(200, 200);
|
|
|
|
bh.iter(|| test::black_box(Hessenberg::new(m.clone())))
|
|
|
|
}
|
|
|
|
|
|
|
|
#[bench]
|
|
|
|
fn hessenberg_decompose_500x500(bh: &mut Bencher) {
|
|
|
|
let m = DMatrix::<f64>::new_random(500, 500);
|
|
|
|
bh.iter(|| test::black_box(Hessenberg::new(m.clone())))
|
|
|
|
}
|
|
|
|
|
|
|
|
// With unpack.
|
|
|
|
#[bench]
|
|
|
|
fn hessenberg_decompose_unpack_100x100(bh: &mut Bencher) {
|
|
|
|
let m = DMatrix::<f64>::new_random(100, 100);
|
|
|
|
bh.iter(|| {
|
|
|
|
let hess = Hessenberg::new(m.clone());
|
|
|
|
let _ = hess.unpack();
|
|
|
|
})
|
|
|
|
}
|
|
|
|
|
|
|
|
#[bench]
|
|
|
|
fn hessenberg_decompose_unpack_200x200(bh: &mut Bencher) {
|
|
|
|
let m = DMatrix::<f64>::new_random(200, 200);
|
|
|
|
bh.iter(|| {
|
|
|
|
let hess = Hessenberg::new(m.clone());
|
|
|
|
let _ = hess.unpack();
|
|
|
|
})
|
|
|
|
}
|
|
|
|
|
|
|
|
#[bench]
|
|
|
|
fn hessenberg_decompose_unpack_500x500(bh: &mut Bencher) {
|
|
|
|
let m = DMatrix::<f64>::new_random(500, 500);
|
|
|
|
bh.iter(|| {
|
|
|
|
let hess = Hessenberg::new(m.clone());
|
|
|
|
let _ = hess.unpack();
|
|
|
|
})
|
|
|
|
}
|