nalgebra 0.27.0

General-purpose linear algebra library with transformations and statically-sized or dynamically-sized 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

use na::Matrix3;

#[test]
#[rustfmt::skip]
fn udu_simple() {
    let m = Matrix3::new(
        2.0, -1.0,  0.0,
       -1.0,  2.0, -1.0,
        0.0, -1.0,  2.0);

    let udu = m.udu().unwrap();
    
    // Rebuild
    let p = udu.u * udu.d_matrix() * udu.u.transpose();

    assert!(relative_eq!(m, p, epsilon = 3.0e-16));
}

#[test]
#[should_panic]
#[rustfmt::skip]
fn udu_non_sym_panic() {
    let m = Matrix3::new(
        2.0, -1.0,  0.0,
        1.0, -2.0,  3.0,
       -2.0,  1.0,  0.3);

    let udu = m.udu().unwrap();
    // Rebuild
    let p = udu.u * udu.d_matrix() * udu.u.transpose();

    assert!(relative_eq!(m, p, epsilon = 3.0e-16));
}

#[cfg(feature = "proptest-support")]
mod proptest_tests {
    #[allow(unused_imports)]
    use crate::core::helper::{RandComplex, RandScalar};

    macro_rules! gen_tests(
        ($module: ident, $scalar: expr) => {
            mod $module {
                #[allow(unused_imports)]
                use crate::core::helper::{RandScalar, RandComplex};
                use crate::proptest::*;
                use proptest::{prop_assert, proptest};

                proptest! {
                    #[test]
                    fn udu(m in dmatrix_($scalar)) {
                        let m = &m * m.adjoint();

                        if let Some(udu) = m.clone().udu() {
                            let p = &udu.u * &udu.d_matrix() * &udu.u.transpose();
                            println!("m: {}, p: {}", m, p);

                            prop_assert!(relative_eq!(m, p, epsilon = 1.0e-7));
                        }
                    }

                    #[test]
                    fn udu_static(m in matrix4_($scalar)) {
                        let m = m.hermitian_part();

                        if let Some(udu) = m.udu() {
                            let p = udu.u * udu.d_matrix() * udu.u.transpose();
                            prop_assert!(relative_eq!(m, p, epsilon = 1.0e-7));
                        }
                    }
                }
            }
        }
    );

    gen_tests!(f64, PROPTEST_F64);
}