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
pub trait MaxDifference {
fn max_diff(self, other: Self) -> f64;
}
impl MaxDifference for f32 {
fn max_diff(self, other: Self) -> f64 {
f64::from((self - other).abs())
}
}
impl MaxDifference for f64 {
fn max_diff(self, other: Self) -> f64 {
(self - other).abs()
}
}
impl MaxDifference for Vec<f32> {
fn max_diff(self, other: Self) -> f64 {
let mut max: f64 = 0.0;
for (a, b) in self.iter().zip(other.iter()) {
let diff = f64::from((*a - *b).abs());
if diff > max {
max = diff;
}
}
max
}
}
impl MaxDifference for Vec<f64> {
fn max_diff(self, other: Self) -> f64 {
let mut max: f64 = 0.0;
for (a, b) in self.iter().zip(other.iter()) {
let diff = (*a - *b).abs();
if diff > max {
max = diff;
}
}
max
}
}
macro_rules! array_impls {
($($N:literal)+) => {
$(
impl MaxDifference for [f64; $N] {
fn max_diff(self, other: Self) -> f64 {
let mut max: f64 = 0.0;
for (a, b) in self.iter().zip(other.iter()) {
let diff = (*a - *b).abs();
if diff > max {
max = diff;
}
}
max
}
}
impl MaxDifference for [f32; $N] {
fn max_diff(self, other: Self) -> f64 {
let mut max: f64 = 0.0;
for (a, b) in self.iter().zip(other.iter()) {
let diff = f64::from((*a - *b).abs());
if diff > max {
max = diff;
}
}
max
}
}
)+
}
}
array_impls! {
0 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
}