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
#[macro_use]
pub mod approx_eq {
#[macro_export]
macro_rules! assert_approx_eq {
($x: expr, $y: expr) => {
let eps = 1e-6;
let (x, y): (f64, f64) = ($x, $y);
if x == 0. {
assert!(y.abs() < eps);
} else if y == 0. {
assert!(x.abs() < eps);
} else {
assert!(&x.signum() == &y.signum());
let (x, y): (f64, f64) = (x.abs(), y.abs());
assert!((&x - &y).abs() / [x, y].iter().cloned().fold(f64::NAN, f64::min) < eps);
}
};
($x: expr, $y: expr, $e: expr) => {
let (x, y): (f64, f64) = ($x, $y);
if x == 0. {
assert!(y.abs() < $e);
} else if y == 0. {
assert!(x.abs() < $e);
} else {
assert!(&x.signum() == &y.signum());
let (x, y): (f64, f64) = (x.abs(), y.abs());
assert!((&x - &y).abs() / [x, y].iter().cloned().fold(f64::NAN, f64::min) < $e);
}
};
}
}
#[test]
fn test_noeps() {
assert_approx_eq!(1., 1.);
assert_approx_eq!(1., 1.000001);
}
#[test]
fn test_witheps() {
assert_approx_eq!(1.0000000001, 1., 1e-5);
}
#[test]
fn test_nearzero() {
assert_approx_eq!(0., -0.0000000000245);
}
#[test]
fn test_zero() {
assert_approx_eq!(0., 0.);
}
#[test]
#[should_panic(expected = "assertion failed")]
fn test_invalid() {
assert_approx_eq!(1.0000000001, 1., 1e-10);
}
#[test]
#[should_panic(expected = "assertion failed")]
fn test_sign() {
assert_approx_eq!(1., -1.);
}