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
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

use matrix::{complex::Complex, Matrix, Vector, C, M, V};

#[test]
fn test_vector_add() {
    let mut u = V!([2., 3.]);
    let v = V!([5., 7.]);
    u.add(&v);
    assert!(u._d.iter().eq(&[7., 10.]));
}

#[test]
#[should_panic(expected = "vectors must be the same size")]
fn test_vector_add_fail() {
    let mut u = V!([2., 3.]);
    let v = V!([5., 7., 8.]);
    u.add(&v);
}

#[test]
#[should_panic(expected = "matrices must be the same size")]
fn test_matrix_add_fail() {
    let mut u = M!([[1., 2.], [3., 4.]]);
    let v = M!([[7., 4.], [-2., 2.], [1., 1.]]);
    u.add(&v);
}

#[test]
fn test_matrix_add() {
    let mut u = M!([[1., 2.], [3., 4.]]);
    let v = M!([[7., 4.], [-2., 2.]]);
    u.add(&v);
    assert!(u._d.iter().eq(&[8., 6., 1., 6.]));
}

#[test]
fn test_vector_add_2() {
    for i in 0..20 {
        let size = i;
        let mut v = vec![];
        for i in 0..size {
            v.push(i as f32);
        }
        let mut a = V!(v);
        let b = a.clone();
        a.add(&b);

        for i in 0..size {
            assert_eq!(a[i] / 2., b[i], "a[i] and b[i] must equal")
        }
    }
}

#[test]
fn test_add_complex_vector() {
    let mut u = V!([C!(1., 2.), C!(-1., -3.)]);
    let v = V!([C!(1., 3.), C!(2., 3.)]);
    u.add(&v);
    assert!(u._d.iter().eq(&[C!(2., 5.), C!(1., 0.)]));
}

#[test]
fn test_add_complex_matrix() {
    let mut u = M!([[C!(1., 0.), C!(2., 0.)], [C!(3., 0.), C!(5., 0.)],]);
    let v = M!([[C!(7., 0.), C!(3., 0.)], [C!(-2., 0.), C!(2., 0.)],]);
    u.add(&v);
    assert!(u
        ._d
        .iter()
        .eq(&[C!(8., 0.), C!(5., 0.), C!(1., 0.), C!(7., 0.)]));
}

#[test]
fn test_add_complex() {
    for i in 0..20 {
        let size = i;
        let mut v = vec![];
        for i in 0..size {
            v.push(i as f32);
        }
        let mut a: Vector<f32> = Vector::from(v);
        let b = a.clone();
        a.add(&b);

        for i in 0..size {
            assert_eq!(a[i] / 2., b[i], "a[i] and b[i] must equal")
        }
    }
}