matrix-42 0.1.0

A linear algebra library written in Rust with basic matrix and vector operations
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

use matrix::{Complex, Matrix, C, M};

#[test]
fn test_trace() {
    let u = M!([[1., 0.], [0., 1.]]);
    let a = u.trace();
    assert_eq!(a.unwrap(), 2.);

    let u = M!([[2., -5., 0.], [4., 3., 7.], [-2., 3., 4.]]);
    let a = u.trace();
    assert_eq!(a.unwrap(), 9.);

    let u = M!([[-2., -8., 4.], [1., -23., 4.], [0., 6., 4.]]);
    let a = u.trace();
    assert_eq!(a.unwrap(), -21.);
}

#[test]
#[should_panic(expected = "matrix must be squared")]
fn test_trace_fail() {
    let u = M!([[1., 0.], [0., 1.], [1., 0.]]);
    let _a = u.trace(); // This should panic
}

#[test]
fn test_trace_complex() {
    let u = M!([[C!(1., 0.), C!(0., 1.)], [C!(0., 0.), C!(1., 0.)]]);
    let a = u.trace();
    assert_eq!(a.unwrap(), C!(2., 0.));

    let u = M!([[C!(2., -5.), C!(0., 0.)], [C!(4., 3.), C!(7., 0.)],]);
    let a = u.trace();
    assert_eq!(a.unwrap(), C!(9., -5.));

    let u = M!([
        [C!(-2., 0.), C!(-8., 0.), C!(4., 0.)],
        [C!(1., 0.), C!(-23., 0.), C!(4., 0.)],
        [C!(0., 0.), C!(6., 0.), C!(4., 0.)]
    ]);
    let a = u.trace();
    assert_eq!(a.unwrap(), C!(-21., 0.));
}