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