qmat 0.5.1

A simple library for 2-dimensional matrices.
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

#[cfg(test)]
mod mul {
    use qmat::matrix;

    #[test]
    fn basic_scalar_2x3() {
        let a = matrix!(2, 3, [1, 5, 3, 2, 6, 7]);
        let b = 4;
        let output = a.mul_scalar(b);
        assert_eq!(*output.as_flat_array(), [4, 20, 12, 8, 24, 28]);
    }

    #[test]
    fn basic_mat_2x2() {
        let a = matrix!(2, 2, [3, 7, 4, 9]);
        let b = matrix!(2, 2, [6, 2, 5, 8]);
        let output = a.multiply::<2, 4, 4>(&b);
        assert_eq!(*output.as_flat_array(), [53, 62, 69, 80]);
    }

    #[test]
    fn basic_mat_3x3() {
        let a = matrix!(3, 3, [12, 8, 4, 3, 17, 14, 9, 8, 10]);
        let b = matrix!(3, 3, [5, 19, 3, 6, 15, 9, 7, 8, 16]);
        let output = a.multiply::<3, 9, 9>(&b);
        assert_eq!(
            *output.as_flat_array(),
            [136, 380, 172, 215, 424, 386, 163, 371, 259]
        );
    }
}

#[cfg(test)]
mod add {
    use qmat::matrix;

    #[test]
    fn basic_mat_2x2() {
        let a = matrix!(3, 3, [12, 8, 4, 3, 17, 14, 9, 8, 10]);
        let b = matrix!(3, 3, [5, 19, 3, 6, 15, 9, 7, 8, 16]);
        let output = a + b;
        assert_eq!(*output.as_flat_array(), [17, 27, 7, 9, 32, 23, 16, 16, 26]);
    }
}

#[cfg(test)]
mod sub {
    use qmat::matrix;

    #[test]
    fn basic_mat_2x2() {
        let a = matrix!(3, 3, [12, 8, 4, 3, 17, 14, 9, 8, 10]);
        let b = matrix!(3, 3, [5, 19, 3, 6, 15, 9, 7, 8, 16]);
        let output = a - b;
        assert_eq!(*output.as_flat_array(), [7, -11, 1, -3, 2, 5, 2, 0, -6]);
    }
}

#[cfg(test)]
mod inverse {
    use qmat::{errors::MatrixOperationError, matrix};

    #[test]
    fn nonsingular_2x2() {
        let mat = matrix!(2, 2, [3.0, 7.0, 1.0, -4.0]);
        let inv = mat.inverse().unwrap();
        assert_eq!(
            *inv.as_flat_array(),
            [4.0 / 19.0, 7.0 / 19.0, 1.0 / 19.0, -3.0 / 19.0]
        )
    }

    #[test]
    fn singular_2x2() {
        let mat = matrix!(2, 2, [3, 6, 2, 4]);
        let res = mat.inverse();
        assert!(res.is_err());
        if let Err(e) = res {
            assert_eq!(e, MatrixOperationError::InvalidDeterminant);
        }
    }
}