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
use crate::faer_tests::test_helper::*;
use crate::ArgminInv;
use approx::assert_relative_eq;
use faer::mat::AsMatRef;
use paste::item;
macro_rules! make_test {
($t:ty) => {
item! {
#[test]
fn [<test_inv_ $t>]() {
let a = matrix2_new(
2 as $t, 5 as $t,
1 as $t, 3 as $t,
);
let target = matrix2_new(
3 as $t, -5 as $t,
-1 as $t, 2 as $t,
);
let res = <_ as ArgminInv<_>>::inv(&a).unwrap();
let res1 = <_ as ArgminInv<_>>::inv(&a.as_mat_ref()).unwrap();
assert_eq!(res,res1);
assert_eq!(res.nrows(),2);
assert_eq!(res.ncols(),2);
for i in 0..2 {
for j in 0..2 {
//@note(geo-ant) the 20 epsilon are a bit arbitrary,
// but it's to avoid spurious errors due to numerical effects
// while keeping good accuracy.
assert_relative_eq!(res[(i, j)], target[(i, j)], epsilon = 20.*$t::EPSILON);
}
}
}
}
item! {
#[test]
fn [<test_inv_error $t>]() {
let a = matrix2_new(
2 as $t, 5 as $t,
4 as $t, 10 as $t,
);
let inv= <_ as ArgminInv<_>>::inv(&a);
assert!(inv.is_err());
}
}
};
}
make_test!(f32);
make_test!(f64);