aeon-tensor 0.3.3

Helper package for working with tensor-based PDEs in aeon-tk binaries
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
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

pub trait TensorStorageRef {
    /// Reference to underlying buffer for tensor storage.
    /// This is only gaurenteed to be >= the most recent call to length.
    fn buffer(&self) -> &[f64];
}

pub trait TensorStorageMut: TensorStorageRef {
    fn buffer_mut(&mut self) -> &mut [f64];
}

pub trait TensorStorageOwned: TensorStorageMut {
    fn resize(&mut self, length: usize);
}

// *********************************
// Implementations

impl TensorStorageRef for f64 {
    fn buffer(&self) -> &[f64] {
        std::slice::from_ref(self)
    }
}

impl TensorStorageMut for f64 {
    fn buffer_mut(&mut self) -> &mut [f64] {
        std::slice::from_mut(self)
    }
}

impl TensorStorageOwned for f64 {
    fn resize(&mut self, length: usize) {
        assert!(
            length <= 1,
            "f64 storage does not have capacity for {} elements",
            length
        );
    }
}

impl<const L: usize> TensorStorageRef for [f64; L] {
    fn buffer(&self) -> &[f64] {
        self.as_slice()
    }
}

impl<const L: usize> TensorStorageMut for [f64; L] {
    fn buffer_mut(&mut self) -> &mut [f64] {
        self.as_mut_slice()
    }
}

impl<const L: usize> TensorStorageOwned for [f64; L] {
    fn resize(&mut self, length: usize) {
        assert!(
            length <= L,
            "static storage of length {} does not have capacity for {} elements",
            L,
            length
        );
    }
}

impl<'a> TensorStorageRef for &'a [f64] {
    fn buffer(&self) -> &[f64] {
        self
    }
}

impl<'a> TensorStorageRef for &'a mut [f64] {
    fn buffer(&self) -> &[f64] {
        self
    }
}

impl<'a> TensorStorageMut for &'a mut [f64] {
    fn buffer_mut(&mut self) -> &mut [f64] {
        self
    }
}

impl TensorStorageRef for Vec<f64> {
    fn buffer(&self) -> &[f64] {
        self
    }
}

impl<'a> TensorStorageMut for Vec<f64> {
    fn buffer_mut(&mut self) -> &mut [f64] {
        self
    }
}

impl<'a> TensorStorageOwned for Vec<f64> {
    fn resize(&mut self, length: usize) {
        assert!(
            length <= self.len(),
            "slice storage of length {} does not have capacity for {} elements",
            self.len(),
            length
        );
    }
}