matrix 0.22.0

The package provides a matrix laboratory.
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

use random::{self, Source};
use test::Bencher;

use matrix::decomposition::SymmetricEigen;
use matrix::format::Conventional;

#[bench]
fn symmetric_eigen_0010(b: &mut Bencher) {
    symmetric_eigen(10, b);
}

#[bench]
fn symmetric_eigen_0100(b: &mut Bencher) {
    symmetric_eigen(100, b);
}

#[bench]
fn symmetric_eigen_1000(b: &mut Bencher) {
    symmetric_eigen(1000, b);
}

fn symmetric_eigen(size: usize, bencher: &mut Bencher) {
    let matrix = make_symmetric(size);
    bencher.iter(|| matrix.decompose());
}

fn make_symmetric(size: usize) -> Conventional<f64> {
    let mut source = random::default().seed([69, 42]);
    let values = source.iter().take(size * size).collect::<Vec<f64>>();
    Conventional::from_vec(size, values)
}