lineas 0.1.1

A mathematical linear algebra library for 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

use lineas::Matrix;
use lineas::decompose::*;

#[cfg(test)]
mod lu_decompose {
	use super::*;
	
	#[test]
	fn correct() {
		let a = Matrix::new([[2, -8, 4], [1, -1, 11], [-3, 10, -8]]).dtype::<f64>();
		let l_check = Matrix::new([[2, 0, 0], [1, 3, 0], [-3, -2, 4]]).dtype::<f64>();
		let u_check = Matrix::new([[1, -4, 2], [0, 1, 3], [0, 0, 1]]).dtype::<f64>();
		let (l, u) = a.lu_decompose().unwrap();
		assert_eq!(l, l_check);
		assert_eq!(u, u_check);
		assert_eq!(l * u, a)
	}
	
	#[test]
	fn incorrect() {
		let a = Matrix::new([[0, 1], [1, 1]]);
		assert!(a.lu_decompose().is_none())
	}
}

#[cfg(test)]
mod plu_decompose {
	use super::*;
	
	#[test]
	fn correct() {
		let a = Matrix::new([[0, 1], [1, 1]]);
		let p_check = Matrix::new([[0, 1], [1, 0]]);
		let l_check = Matrix::new([[1, 0], [0, 1]]);
		let u_check = Matrix::new([[1, 1], [0, 1]]);
		let (p, l, u) = a.plu_decompose().unwrap();
		assert_eq!(p, p_check);
		assert_eq!(l, l_check);
		assert_eq!(u, u_check);
		assert_eq!(p * l * u, a)
	}
	
	// #[test]
	// fn incorrect() {
	// 	let a = Matrix::new([[0, 1], [1, 1]]);
	// 	assert!(a.lu_decompose().is_none())
	// }
}