unitforge 0.4.0

A library for unit and quantity consistent computations in Rust
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

#[cfg(test)]
mod matrix2_tests {
    use ndarray::arr2;
    use unitforge::small_linalg::Matrix2;
    use unitforge::PhysicsQuantity;

    fn approx_eq(a: f64, b: f64, tol: f64) -> bool {
        (a - b).abs() < tol
    }

    #[test]
    fn test_matrix2_add() {
        let a: Matrix2<f64> = Matrix2::from_f64([[1.0, 2.0], [3.0, 4.0]]);
        let b = Matrix2::from_f64([[5.0, 6.0], [7.0, 8.0]]);
        let expected = Matrix2::from_f64([[6.0, 8.0], [10.0, 12.0]]);
        assert_eq!(a + b, expected);
    }

    #[test]
    fn test_matrix2_sub() {
        let a: Matrix2<f64> = Matrix2::from_f64([[5.0, 7.0], [9.0, 11.0]]);
        let b = Matrix2::from_f64([[1.0, 2.0], [3.0, 4.0]]);
        let expected = Matrix2::from_f64([[4.0, 5.0], [6.0, 7.0]]);
        assert_eq!(a - b, expected);
    }

    #[test]
    fn test_matrix2_add_assign() {
        let mut a: Matrix2<f64> = Matrix2::from_f64([[1.0, 2.0], [3.0, 4.0]]);
        let b = Matrix2::from_f64([[1.0, 1.0], [1.0, 1.0]]);
        a += b;
        assert_eq!(a, Matrix2::from_f64([[2.0, 3.0], [4.0, 5.0]]));
    }

    #[test]
    fn test_matrix2_sub_assign() {
        let mut a: Matrix2<f64> = Matrix2::from_f64([[5.0, 5.0], [5.0, 5.0]]);
        let b = Matrix2::from_f64([[1.0, 2.0], [3.0, 4.0]]);
        a -= b;
        assert_eq!(a, Matrix2::from_f64([[4.0, 3.0], [2.0, 1.0]]));
    }

    #[test]
    fn test_matrix2_scalar_mul() {
        let m: Matrix2<f64> = Matrix2::from_f64([[1.0, 2.0], [3.0, 4.0]]);
        let result = m * 2.0;
        assert_eq!(result, Matrix2::from_f64([[2.0, 4.0], [6.0, 8.0]]));
    }

    #[test]
    fn test_matrix2_neg() {
        let m: Matrix2<f64> = Matrix2::from_f64([[1.0, -2.0], [-3.0, 4.0]]);
        let result = -m;
        assert_eq!(result, Matrix2::from_f64([[-1.0, 2.0], [3.0, -4.0]]));
    }

    #[test]
    fn test_matrix2_det() {
        let m: Matrix2<f64> = Matrix2::from_f64([[4.0, 3.0], [6.0, 3.0]]);
        let det = m.det();
        assert!(approx_eq(det.as_f64(), -6.0, 1e-10));
    }

    #[test]
    fn test_matrix2_inverse_exists() {
        let m: Matrix2<f64> = Matrix2::from_f64([[4.0, 7.0], [2.0, 6.0]]);
        let inv = m.inverse().unwrap();
        let expected = Matrix2::from_f64([[0.6, -0.7], [-0.2, 0.4]]);
        for i in 0..2 {
            for j in 0..2 {
                assert!(approx_eq(inv[(i, j)].as_f64(), expected[(i, j)], 1e-10));
            }
        }
    }

    #[test]
    fn test_matrix2_inverse_singular() {
        let m: Matrix2<f64> = Matrix2::from_f64([[1.0, 2.0], [2.0, 4.0]]);
        assert!(m.inverse().is_none());
    }

    #[test]
    fn test_matrix2_transpose() {
        let m: Matrix2<f64> = Matrix2::from_f64([[1.0, 2.0], [3.0, 4.0]]);
        let t = m.transpose();
        assert_eq!(t, Matrix2::from_f64([[1.0, 3.0], [2.0, 4.0]]));
    }

    #[test]
    fn test_matrix2_frobenius_norm() {
        let m: Matrix2<f64> = Matrix2::from_f64([[1.0, 2.0], [3.0, 4.0]]);
        let norm = m.frobenius_norm().as_f64();
        assert!(approx_eq(norm, (1.0f64 + 4.0 + 9.0 + 16.0).sqrt(), 1e-10));
    }

    #[test]
    fn test_matrix2_abs() {
        let m: Matrix2<f64> = Matrix2::from_f64([[-1.0, 2.0], [-3.0, 4.0]]);
        let abs = m.abs();
        assert_eq!(abs, Matrix2::from_f64([[1.0, 2.0], [3.0, 4.0]]));
    }

    #[test]
    fn test_identity() {
        let eye = Matrix2::identity();
        assert_eq!(eye, Matrix2::from_f64([[1.0, 0.0], [0.0, 1.0]]));
    }

    #[test]
    fn test_to_ndarray() {
        let mat = Matrix2::new([[1.0, 2.0], [3.0, 4.0]]);
        let arr = mat.to_ndarray().unwrap();
        assert_eq!(mat[(0, 0)], arr[[0, 0]]);
        assert_eq!(mat[(1, 0)], arr[[1, 0]]);
        assert_eq!(mat[(0, 1)], arr[[0, 1]]);
        assert_eq!(mat[(1, 1)], arr[[1, 1]]);
    }

    #[test]
    fn test_from_ndarray() {
        let arr = arr2(&[[1.0, 2.0], [3.0, 4.0]]);
        let mat = Matrix2::from_ndarray(arr.view()).unwrap();
        assert_eq!(mat[(0, 0)], arr[[0, 0]]);
        assert_eq!(mat[(1, 0)], arr[[1, 0]]);
        assert_eq!(mat[(0, 1)], arr[[0, 1]]);
        assert_eq!(mat[(1, 1)], arr[[1, 1]]);
    }
}