ndarray-linalg 0.4.1

Linear algebra package for rust-ndarray using LAPACK
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

include!("header.rs");

macro_rules! impl_test {
    ($funcname:ident, $random:path, $n:expr, $m:expr, $t:expr) => {
#[test]
fn $funcname() {
    use ndarray_linalg::prelude::*;
    let a = $random($n, $m, $t);
    let ans = a.clone();
    let (p, l, u) = a.lu().unwrap();
    println!("P = \n{:?}", &p);
    println!("L = \n{:?}", &l);
    println!("U = \n{:?}", &u);
    println!("LU = \n{:?}", l.dot(&u));
    assert_close_l2!(&l.dot(&u).permutated(&p), &ans, 1e-7);
}
}} // impl_test

macro_rules! impl_test_lu {
    ($modname:ident, $random:path) => {
mod $modname {
    impl_test!(lu_square, $random, 3, 3, false);
    impl_test!(lu_square_t, $random, 3, 3, true);
    impl_test!(lu_3x4, $random, 3, 4, false);
    impl_test!(lu_3x4_t, $random, 3, 4, true);
    impl_test!(lu_4x3, $random, 4, 3, false);
    impl_test!(lu_4x3_t, $random, 4, 3, true);
}
}} // impl_test_lu

impl_test_lu!(owned, super::random_owned);
impl_test_lu!(shared, super::random_shared);