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
pub mod evd;
use crate::number::Number;
#[derive(Clone, Debug, Default, Hash)]
pub struct CirculantMatrix<T = f64>
where
T: Number,
{
col_elems: Vec<T>,
}
impl<T> CirculantMatrix<T>
where
T: Number,
{
pub fn new(col_elems: Vec<T>) -> Self {
Self { col_elems }
}
pub fn from(row_elems: &[T]) -> Self {
let col_elems = if row_elems.len() <= 1 {
row_elems.to_vec()
} else {
row_elems[0..1]
.iter()
.chain(row_elems[1..].iter().rev())
.map(|&e| e)
.collect::<Vec<_>>()
};
Self::new(col_elems)
}
pub fn col_elems(&self) -> &[T] {
&self.col_elems
}
pub fn row_elems(&self) -> Vec<T> {
if self.col_elems.len() <= 1 {
return self.col_elems.clone();
}
self.col_elems[0..1]
.iter()
.chain(self.col_elems[1..].iter().rev())
.map(|&e| e)
.collect::<Vec<_>>()
}
}