arnaucube
arnaucube/poulpyPublic Self-hosted repository
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}