poulpy/backend/src/vec_znx_dft.rs
backend/src/vec_znx_dft.rsRaw
1use std::marker::PhantomData;
2
3use crate::ffi::vec_znx_dft;
4use crate::znx_base::ZnxInfos;
5use crate::{
6 Backend, DataView, DataViewMut, FFT64, Module, VecZnxBig, ZnxSliceSize, ZnxView, alloc_aligned,
7};
8use std::fmt;
9
10pub struct VecZnxDft<D, B: Backend> {
11 pub(crate) data: D,
12 pub(crate) n: usize,
13 pub(crate) cols: usize,
14 pub(crate) size: usize,
15 pub(crate) _phantom: PhantomData<B>,
16}
17
18impl<D, B: Backend> VecZnxDft<D, B> {
19 pub fn into_big(self) -> VecZnxBig<D, B> {
20 VecZnxBig::<D, B>::from_data(self.data, self.n, self.cols, self.size)
21 }
22}
23
24impl<D, B: Backend> ZnxInfos for VecZnxDft<D, B> {
25 fn cols(&self) -> usize {
26 self.cols
27 }
28
29 fn rows(&self) -> usize {
30 1
31 }
32
33 fn n(&self) -> usize {
34 self.n
35 }
36
37 fn size(&self) -> usize {
38 self.size
39 }
40}
41
42impl<D> ZnxSliceSize for VecZnxDft<D, FFT64> {
43 fn sl(&self) -> usize {
44 self.n() * self.cols()
45 }
46}
47
48impl<D, B: Backend> DataView for VecZnxDft<D, B> {
49 type D = D;
50 fn data(&self) -> &Self::D {
51 &self.data
52 }
53}
54
55impl<D, B: Backend> DataViewMut for VecZnxDft<D, B> {
56 fn data_mut(&mut self) -> &mut Self::D {
57 &mut self.data
58 }
59}
60
61impl<D: AsRef<[u8]>> ZnxView for VecZnxDft<D, FFT64> {
62 type Scalar = f64;
63}
64
65pub(crate) fn bytes_of_vec_znx_dft<B: Backend>(module: &Module<B>, cols: usize, size: usize) -> usize {
66 unsafe { vec_znx_dft::bytes_of_vec_znx_dft(module.ptr, size as u64) as usize * cols }
67}
68
69impl<D: From<Vec<u8>>, B: Backend> VecZnxDft<D, B> {
70 pub(crate) fn new(module: &Module<B>, cols: usize, size: usize) -> Self {
71 let data = alloc_aligned::<u8>(bytes_of_vec_znx_dft(module, cols, size));
72 Self {
73 data: data.into(),
74 n: module.n(),
75 cols,
76 size,
77 _phantom: PhantomData,
78 }
79 }
80
81 pub(crate) fn new_from_bytes(module: &Module<B>, cols: usize, size: usize, bytes: impl Into<Vec<u8>>) -> Self {
82 let data: Vec<u8> = bytes.into();
83 assert!(data.len() == bytes_of_vec_znx_dft(module, cols, size));
84 Self {
85 data: data.into(),
86 n: module.n(),
87 cols,
88 size,
89 _phantom: PhantomData,
90 }
91 }
92}
93
94pub type VecZnxDftOwned<B> = VecZnxDft<Vec<u8>, B>;
95
96impl<D, B: Backend> VecZnxDft<D, B> {
97 pub(crate) fn from_data(data: D, n: usize, cols: usize, size: usize) -> Self {
98 Self {
99 data,
100 n,
101 cols,
102 size,
103 _phantom: PhantomData,
104 }
105 }
106}
107
108pub trait VecZnxDftToRef<B: Backend> {
109 fn to_ref(&self) -> VecZnxDft<&[u8], B>;
110}
111
112impl<D, B: Backend> VecZnxDftToRef<B> for VecZnxDft<D, B>
113where
114 D: AsRef<[u8]>,
115 B: Backend,
116{
117 fn to_ref(&self) -> VecZnxDft<&[u8], B> {
118 VecZnxDft {
119 data: self.data.as_ref(),
120 n: self.n,
121 cols: self.cols,
122 size: self.size,
123 _phantom: std::marker::PhantomData,
124 }
125 }
126}
127
128pub trait VecZnxDftToMut<B: Backend> {
129 fn to_mut(&mut self) -> VecZnxDft<&mut [u8], B>;
130}
131
132impl<D, B: Backend> VecZnxDftToMut<B> for VecZnxDft<D, B>
133where
134 D: AsRef<[u8]> + AsMut<[u8]>,
135 B: Backend,
136{
137 fn to_mut(&mut self) -> VecZnxDft<&mut [u8], B> {
138 VecZnxDft {
139 data: self.data.as_mut(),
140 n: self.n,
141 cols: self.cols,
142 size: self.size,
143 _phantom: std::marker::PhantomData,
144 }
145 }
146}
147
148impl<D: AsRef<[u8]>> fmt::Display for VecZnxDft<D, FFT64> {
149 fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result {
150 writeln!(
151 f,
152 "VecZnxDft(n={}, cols={}, size={})",
153 self.n, self.cols, self.size
154 )?;
155
156 for col in 0..self.cols {
157 writeln!(f, "Column {}:", col)?;
158 for size in 0..self.size {
159 let coeffs = self.at(col, size);
160 write!(f, " Size {}: [", size)?;
161
162 let max_show = 100;
163 let show_count = coeffs.len().min(max_show);
164
165 for (i, &coeff) in coeffs.iter().take(show_count).enumerate() {
166 if i > 0 {
167 write!(f, ", ")?;
168 }
169 write!(f, "{}", coeff)?;
170 }
171
172 if coeffs.len() > max_show {
173 write!(f, ", ... ({} more)", coeffs.len() - max_show)?;
174 }
175
176 writeln!(f, "]")?;
177 }
178 }
179 Ok(())
180 }
181}