mymatrix 0.9.0

My simple matrix library that can perform fraction 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
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
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165

use mymatrix::Vector;
use pyinrs::Fraction;
use rstest::{fixture, rstest};

struct Fixture {
    empty: Vector,
    one: Vector,
    some: Vector,
}

#[fixture]
fn setup() -> Fixture {
    Fixture {
        empty: Vector::new(),
        one: Vector::from([1]),
        some: Vector::from([1, 2, 3, 4, 5]),
    }
}

#[rstest]
fn basics(setup: Fixture) {
    assert_eq!(setup.empty.size(), 0);
    assert!(setup.empty.is_empty());

    assert_eq!(setup.one.size(), 1);
    assert!(!setup.one.is_empty());

    assert_eq!(setup.some.size(), 5);
    assert!(!setup.some.is_empty());
}

#[rstest]
fn compare(setup: Fixture) {
    assert!(setup.some == Vector::from([1, 2, 3, 4, 5]));
    assert!(setup.some != Vector::from([1, 3, 5]));
}

#[rstest]
fn access(mut setup: Fixture) {
    // for
    for i in 0..setup.some.size() {
        assert_eq!(setup.some[i], ((i + 1) as i32).into());
    }

    // assignment
    setup.some[0] = 0.into();
    assert_eq!(setup.some[0], 0.into());
}

#[rstest]
fn is_zero() {
    assert!(Vector::from([0]).is_zero());
    assert!(Vector::from([0, 0, 0]).is_zero());

    assert!(!Vector::from([1]).is_zero());
    assert!(!Vector::from([0, 0, 1]).is_zero());
}

#[rstest]
fn is_orthogonal() {
    let zero = Vector::from([0, 0]);
    assert!(zero.is_orthogonal(&Vector::from([0, 0])));
    assert!(zero.is_orthogonal(&Vector::from([1, 1])));
    assert!(zero.is_orthogonal(&Vector::from([2, 3])));

    let one = Vector::from([1, 1]);
    assert!(!one.is_orthogonal(&Vector::from([1, 1])));
    assert!(one.is_orthogonal(&Vector::from([1, -1])));
    assert!(one.is_orthogonal(&Vector::from([-1, 1])));
    assert!(one.is_orthogonal(&Vector::from([-2, 2])));
}

#[rstest]
fn is_parallel() {
    let zero = Vector::from([0, 0]);
    assert!(zero.is_parallel(&Vector::from([0, 0])));
    assert!(zero.is_parallel(&Vector::from([1, 1])));
    assert!(zero.is_parallel(&Vector::from([2, 3])));

    let some = Vector::from([3, 4]);
    assert!(!some.is_parallel(&Vector::from([1, 1])));
    assert!(some.is_parallel(&Vector::from([3, 4])));
    assert!(some.is_parallel(&Vector::from([-3, -4])));
    assert!(some.is_parallel(&Vector::from([6, 8])));

    assert!(Vector::from([1, 1]).is_parallel(&Vector::from([2, 2])));
}

#[rstest]
fn norm() {
    assert_eq!(Vector::from([0]).norm(), 0.0);
    assert_eq!(Vector::from([1]).norm(), 1.0);
    assert_eq!(Vector::from([3, 4]).norm(), 5.0);
}

#[rstest]
fn count_leading_zeros() {
    assert_eq!(Vector::from([0]).count_leading_zeros(), 1);
    assert_eq!(Vector::from([0, 1]).count_leading_zeros(), 1);
    assert_eq!(Vector::from([0, 0]).count_leading_zeros(), 2);
    assert_eq!(Vector::from([0, 0, 1]).count_leading_zeros(), 2);
    assert_eq!(Vector::from([0, 0, 0, 1, 2, 3]).count_leading_zeros(), 3);
    assert_eq!(Vector::from([0, 0, 0, 1, 2, 3, 0, 0, 0]).count_leading_zeros(), 3);
}

#[rstest]
fn cross() {
    assert_eq!(Vector::cross(&[1, 2].into(), &[3, 4].into()), [-2].into());
    assert_eq!(Vector::cross(&[1, 2, 3].into(), &[4, 5, 6].into()), [-3, 6, -3].into());
}

#[rstest]
#[should_panic(expected = "Error: Incompatible dimensions for cross product.")]
fn bad_cross() {
    Vector::cross(&[1].into(), &[2].into());
}

#[rstest]
fn add() {
    assert_eq!(Vector::from([1]) + Vector::from([1]), Vector::from([2]));
    assert_eq!(Vector::from([1, 2, 3]) + Vector::from([4, 5, 6]), Vector::from([5, 7, 9]));
    assert_eq!(Vector::from([1, 2]) + Vector::from([1, 2]) + Vector::from([1, 2]), Vector::from([3, 6]));
}

#[rstest]
#[should_panic(expected = "Error: The dimensions mismatch.")]
fn bad_add() {
    let _ = Vector::from([1]) + Vector::from([1, 2]);
}

#[rstest]
fn sub() {
    assert_eq!(Vector::from([1]) - Vector::from([1]), Vector::from([0]));
    assert_eq!(Vector::from([1, 2, 3]) - Vector::from([4, 5, 6]), Vector::from([-3, -3, -3]));
    assert_eq!(Vector::from([1, 2]) - Vector::from([1, 2]) - Vector::from([1, 2]), Vector::from([-1, -2]));
}

#[rstest]
fn scalar_mul() {
    assert_eq!(Fraction::from(1) * Vector::from([1]), Vector::from([1]));
    assert_eq!(Vector::from([1]) * Fraction::from(1), Vector::from([1]));
    assert_eq!(Vector::from([1, 2, 3]) * Fraction::from(2), Vector::from([2, 4, 6]));
    assert_eq!(
        Vector::from([1, 2]) * Fraction::from(2) * Fraction::from((4, 10)),
        Vector::from([Fraction::from((8, 10)), Fraction::from((16, 10))])
    );
}

#[rstest]
fn dot_mul() {
    assert_eq!(Vector::from([1]) * Vector::from([1]), 1.into());
    assert_eq!(Vector::from([1, 2, 3]) * Vector::from([4, 5, 6]), 32.into());
}

#[rstest]
fn format(setup: Fixture) {
    assert_eq!(format!("{}", setup.empty), "[]");
    assert_eq!(format!("{}", setup.one), "[1]");
    assert_eq!(format!("{}", setup.some), "[1 2 3 4 5]");

    assert_eq!(
        format!("{}", Vector::from([Fraction::from((-3, 4)), Fraction::new(), Fraction::from((5, 6))])),
        "[-3/4    0  5/6]"
    );
}