matrix-mc 0.1.0

some useful matrix and vector operations performed on 2D Vecs of f64s
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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93

use matrix_mc::{identity, rref, matrix, normalize, dup};

#[test]
fn init_identity_3x3() {
    let expected = 
        matrix![1, 0, 0;
                0, 1, 0;
                0, 0, 1];
    
    let actual = identity(3usize);

    assert_eq!(expected, actual);
}

#[test]
fn rref_independent_3x3() {

    let expected = identity(3);

    let mut actual = 
        matrix![2,  3,  4 ;
                3,  1,  5 ;
                7,  11, 22];
    rref(&mut actual);

    assert_eq!(expected, actual);
}

#[test]
fn rref_dependent_3x3() {
    let expected = 
        matrix![1, 0, 0;
                0, 1, 1;
                0, 0, 0];
    
    let mut actual = 
        matrix![1, 1, 1;
                0, 1, 1;
                1, 1, 1];
    rref(&mut actual);
    
    assert_eq!(expected, actual);
}


#[test]
fn rref_dependent_3x2() {
    let expected = 
        matrix![1, 0;
                0, 1;
                0, 0];

    let mut actual = 
        matrix![5,  4;
                0,  3;
                10, 8];
    rref(&mut actual);

    assert_eq!(expected, actual);    
}

#[test]
fn normalize_3x3() {
    let expected =
        matrix![1f64/f64::sqrt(2f64), 1f64/f64::sqrt(5f64), 1;
                1f64/f64::sqrt(2f64),                    0, 0;
                                   0, 2f64/f64::sqrt(5f64), 0];
    
    let mut actual = 
        matrix![1, 1, 1;
                1, 0, 0;
                0, 2, 0];
    normalize(&mut actual);

    assert_eq!(expected, actual);
}

#[test]
fn dup_rref_3x3() {
    let expected = matrix!
        [1, 0, 0;
         0, 1, 0;
         0, 0, 1];
    
    let input = matrix! 
        [5, 4, 3;
         2, 1, 9;
         8, 7, 6];

    let actual = dup(rref, &input);

    assert_eq!(expected, actual);
}