poulpy_core/layouts/
gglwe_tsk.rs

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

poulpy_core/layouts/
gglwe_tsk.rs
