poulpy/core/src/layouts/compressed/gglwe_tsk.rs
core/src/layouts/compressed/gglwe_tsk.rsRaw
1use backend::hal::{
2 api::{FillUniform, Reset, VecZnxCopy, VecZnxFillUniform},
3 layouts::{Backend, Data, DataMut, DataRef, MatZnx, Module, ReaderFrom, WriterTo},
4};
5
6use crate::layouts::{GGLWETensorKey, Infos, compressed::GGLWESwitchingKeyCompressed};
7use byteorder::{LittleEndian, ReadBytesExt, WriteBytesExt};
8use std::fmt;
9
10use crate::trait_families::Decompress;
11
12#[derive(PartialEq, Eq, Clone)]
13pub struct GGLWETensorKeyCompressed<D: Data> {
14 pub(crate) keys: Vec<GGLWESwitchingKeyCompressed<D>>,
15}
16
17impl<D: DataRef> fmt::Debug for GGLWETensorKeyCompressed<D> {
18 fn fmt(&self, f: &mut fmt::Formatter) -> fmt::Result {
19 write!(f, "{}", self)
20 }
21}
22
23impl<D: DataMut> FillUniform for GGLWETensorKeyCompressed<D> {
24 fn fill_uniform(&mut self, source: &mut sampling::source::Source) {
25 self.keys
26 .iter_mut()
27 .for_each(|key: &mut GGLWESwitchingKeyCompressed<D>| key.fill_uniform(source))
28 }
29}
30
31impl<D: DataMut> Reset for GGLWETensorKeyCompressed<D>
32where
33 MatZnx<D>: Reset,
34{
35 fn reset(&mut self) {
36 self.keys
37 .iter_mut()
38 .for_each(|key: &mut GGLWESwitchingKeyCompressed<D>| key.reset())
39 }
40}
41
42impl<D: DataRef> fmt::Display for GGLWETensorKeyCompressed<D> {
43 fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result {
44 writeln!(f, "(GLWETensorKeyCompressed)",)?;
45 for (i, key) in self.keys.iter().enumerate() {
46 write!(f, "{}: {}", i, key)?;
47 }
48 Ok(())
49 }
50}
51
52impl GGLWETensorKeyCompressed<Vec<u8>> {
53 pub fn alloc(n: usize, basek: usize, k: usize, rows: usize, digits: usize, rank: usize) -> Self {
54 let mut keys: Vec<GGLWESwitchingKeyCompressed<Vec<u8>>> = Vec::new();
55 let pairs: usize = (((rank + 1) * rank) >> 1).max(1);
56 (0..pairs).for_each(|_| {
57 keys.push(GGLWESwitchingKeyCompressed::alloc(
58 n, basek, k, rows, digits, 1, rank,
59 ));
60 });
61 Self { keys: keys }
62 }
63
64 pub fn bytes_of(n: usize, basek: usize, k: usize, rows: usize, digits: usize, rank: usize) -> usize {
65 let pairs: usize = (((rank + 1) * rank) >> 1).max(1);
66 pairs * GGLWESwitchingKeyCompressed::bytes_of(n, basek, k, rows, digits, 1)
67 }
68}
69
70impl<D: Data> Infos for GGLWETensorKeyCompressed<D> {
71 type Inner = MatZnx<D>;
72
73 fn inner(&self) -> &Self::Inner {
74 &self.keys[0].inner()
75 }
76
77 fn basek(&self) -> usize {
78 self.keys[0].basek()
79 }
80
81 fn k(&self) -> usize {
82 self.keys[0].k()
83 }
84}
85
86impl<D: Data> GGLWETensorKeyCompressed<D> {
87 pub fn rank(&self) -> usize {
88 self.keys[0].rank()
89 }
90
91 pub fn digits(&self) -> usize {
92 self.keys[0].digits()
93 }
94
95 pub fn rank_in(&self) -> usize {
96 self.keys[0].rank_in()
97 }
98
99 pub fn rank_out(&self) -> usize {
100 self.keys[0].rank_out()
101 }
102}
103
104impl<D: DataMut> ReaderFrom for GGLWETensorKeyCompressed<D> {
105 fn read_from<R: std::io::Read>(&mut self, reader: &mut R) -> std::io::Result<()> {
106 let len: usize = reader.read_u64::<LittleEndian>()? as usize;
107 if self.keys.len() != len {
108 return Err(std::io::Error::new(
109 std::io::ErrorKind::InvalidData,
110 format!("self.keys.len()={} != read len={}", self.keys.len(), len),
111 ));
112 }
113 for key in &mut self.keys {
114 key.read_from(reader)?;
115 }
116 Ok(())
117 }
118}
119
120impl<D: DataRef> WriterTo for GGLWETensorKeyCompressed<D> {
121 fn write_to<W: std::io::Write>(&self, writer: &mut W) -> std::io::Result<()> {
122 writer.write_u64::<LittleEndian>(self.keys.len() as u64)?;
123 for key in &self.keys {
124 key.write_to(writer)?;
125 }
126 Ok(())
127 }
128}
129
130impl<D: DataMut> GGLWETensorKeyCompressed<D> {
131 pub(crate) fn at_mut(&mut self, mut i: usize, mut j: usize) -> &mut GGLWESwitchingKeyCompressed<D> {
132 if i > j {
133 std::mem::swap(&mut i, &mut j);
134 };
135 let rank: usize = self.rank();
136 &mut self.keys[i * rank + j - (i * (i + 1) / 2)]
137 }
138}
139
140impl<D: DataMut, DR: DataRef, B: Backend> Decompress<B, GGLWETensorKeyCompressed<DR>> for GGLWETensorKey<D> {
141 fn decompress(&mut self, module: &Module<B>, other: &GGLWETensorKeyCompressed<DR>)
142 where
143 Module<B>: VecZnxFillUniform + VecZnxCopy,
144 {
145 #[cfg(debug_assertions)]
146 {
147 assert_eq!(
148 self.keys.len(),
149 other.keys.len(),
150 "invalid receiver: self.keys.len()={} != other.keys.len()={}",
151 self.keys.len(),
152 other.keys.len()
153 );
154 }
155
156 self.keys
157 .iter_mut()
158 .zip(other.keys.iter())
159 .for_each(|(a, b)| {
160 a.decompress(module, b);
161 });
162 }
163}