arnaucube
arnaucube/poulpyPublic Self-hosted repository
poulpy/poulpy-core/src/layouts/gglwe_tsk.rs
poulpy-core/src/layouts/gglwe_tsk.rsRaw
  1use 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    // Returns a mutable reference to GLWESwitchingKey_{s}(s[i] * s[j])
104    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    // Returns a reference to GLWESwitchingKey_{s}(s[i] * s[j])
115    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}