use poulpy_hal::{ layouts::{Data, DataMut, DataRef, FillUniform, MatZnx, ReaderFrom, WriterTo, ZnxInfos}, source::Source, }; use std::fmt; use crate::layouts::{Base2K, BuildError, Degree, Dnum, Dsize, GLWECiphertext, GLWEInfos, LWEInfos, Rank, TorusPrecision}; pub trait GGSWInfos where Self: GLWEInfos, { fn dnum(&self) -> Dnum; fn dsize(&self) -> Dsize; fn layout(&self) -> GGSWCiphertextLayout { GGSWCiphertextLayout { n: self.n(), base2k: self.base2k(), k: self.k(), rank: self.rank(), dnum: self.dnum(), dsize: self.dsize(), } } } #[derive(PartialEq, Eq, Copy, Clone, Debug)] pub struct GGSWCiphertextLayout { pub n: Degree, pub base2k: Base2K, pub k: TorusPrecision, pub rank: Rank, pub dnum: Dnum, pub dsize: Dsize, } impl LWEInfos for GGSWCiphertextLayout { fn base2k(&self) -> Base2K { self.base2k } fn k(&self) -> TorusPrecision { self.k } fn n(&self) -> Degree { self.n } } impl GLWEInfos for GGSWCiphertextLayout { fn rank(&self) -> Rank { self.rank } } impl GGSWInfos for GGSWCiphertextLayout { fn dsize(&self) -> Dsize { self.dsize } fn dnum(&self) -> Dnum { self.dnum } } #[derive(PartialEq, Eq, Clone)] pub struct GGSWCiphertext { pub(crate) data: MatZnx, pub(crate) k: TorusPrecision, pub(crate) base2k: Base2K, pub(crate) dsize: Dsize, } impl LWEInfos for GGSWCiphertext { fn n(&self) -> Degree { Degree(self.data.n() as u32) } fn base2k(&self) -> Base2K { self.base2k } fn k(&self) -> TorusPrecision { self.k } fn size(&self) -> usize { self.data.size() } } impl GLWEInfos for GGSWCiphertext { fn rank(&self) -> Rank { Rank(self.data.cols_out() as u32 - 1) } } impl GGSWInfos for GGSWCiphertext { fn dsize(&self) -> Dsize { self.dsize } fn dnum(&self) -> Dnum { Dnum(self.data.rows() as u32) } } pub struct GGSWCiphertextBuilder { data: Option>, base2k: Option, k: Option, dsize: Option, } impl GGSWCiphertext { #[inline] pub fn builder() -> GGSWCiphertextBuilder { GGSWCiphertextBuilder { data: None, base2k: None, k: None, dsize: None, } } } impl GGSWCiphertextBuilder> { #[inline] pub fn layout(mut self, infos: &A) -> Self where A: GGSWInfos, { debug_assert!( infos.size() as u32 > infos.dsize().0, "invalid ggsw: ceil(k/base2k): {} <= dsize: {}", infos.size(), infos.dsize() ); assert!( infos.dnum().0 * infos.dsize().0 <= infos.size() as u32, "invalid ggsw: dnum: {} * dsize:{} > ceil(k/base2k): {}", infos.dnum(), infos.dsize(), infos.size(), ); self.data = Some(MatZnx::alloc( infos.n().into(), infos.dnum().into(), (infos.rank() + 1).into(), (infos.rank() + 1).into(), infos.size(), )); self.base2k = Some(infos.base2k()); self.k = Some(infos.k()); self.dsize = Some(infos.dsize()); self } } impl GGSWCiphertextBuilder { #[inline] pub fn data(mut self, data: MatZnx) -> Self { self.data = Some(data); self } #[inline] pub fn base2k(mut self, base2k: Base2K) -> Self { self.base2k = Some(base2k); self } #[inline] pub fn k(mut self, k: TorusPrecision) -> Self { self.k = Some(k); self } #[inline] pub fn dsize(mut self, dsize: Dsize) -> Self { self.dsize = Some(dsize); self } pub fn build(self) -> Result, BuildError> { let data: MatZnx = self.data.ok_or(BuildError::MissingData)?; let base2k: Base2K = self.base2k.ok_or(BuildError::MissingBase2K)?; let k: TorusPrecision = self.k.ok_or(BuildError::MissingK)?; let dsize: Dsize = self.dsize.ok_or(BuildError::MissingDigits)?; if base2k == 0_u32 { return Err(BuildError::ZeroBase2K); } if dsize == 0_u32 { return Err(BuildError::ZeroBase2K); } if k == 0_u32 { return Err(BuildError::ZeroTorusPrecision); } if data.n() == 0 { return Err(BuildError::ZeroDegree); } if data.cols() == 0 { return Err(BuildError::ZeroCols); } if data.size() == 0 { return Err(BuildError::ZeroLimbs); } Ok(GGSWCiphertext { data, base2k, k, dsize, }) } } impl fmt::Debug for GGSWCiphertext { fn fmt(&self, f: &mut fmt::Formatter) -> fmt::Result { write!(f, "{}", self.data) } } impl fmt::Display for GGSWCiphertext { fn fmt(&self, f: &mut fmt::Formatter<'_>) -> fmt::Result { write!( f, "(GGSWCiphertext: k: {} base2k: {} dsize: {}) {}", self.k().0, self.base2k().0, self.dsize().0, self.data ) } } impl FillUniform for GGSWCiphertext { fn fill_uniform(&mut self, log_bound: usize, source: &mut Source) { self.data.fill_uniform(log_bound, source); } } impl GGSWCiphertext { pub fn at(&self, row: usize, col: usize) -> GLWECiphertext<&[u8]> { GLWECiphertext::builder() .data(self.data.at(row, col)) .base2k(self.base2k()) .k(self.k()) .build() .unwrap() } } impl GGSWCiphertext { pub fn at_mut(&mut self, row: usize, col: usize) -> GLWECiphertext<&mut [u8]> { GLWECiphertext::builder() .base2k(self.base2k()) .k(self.k()) .data(self.data.at_mut(row, col)) .build() .unwrap() } } impl GGSWCiphertext> { pub fn alloc(infos: &A) -> Self where A: GGSWInfos, { Self::alloc_with( infos.n(), infos.base2k(), infos.k(), infos.rank(), infos.dnum(), infos.dsize(), ) } pub fn alloc_with(n: Degree, base2k: Base2K, k: TorusPrecision, rank: Rank, dnum: Dnum, dsize: Dsize) -> Self { let size: usize = k.0.div_ceil(base2k.0) as usize; debug_assert!( size as u32 > dsize.0, "invalid ggsw: ceil(k/base2k): {size} <= dsize: {}", dsize.0 ); assert!( dnum.0 * dsize.0 <= size as u32, "invalid ggsw: dnum: {} * dsize:{} > ceil(k/base2k): {size}", dnum.0, dsize.0, ); Self { data: MatZnx::alloc( n.into(), dnum.into(), (rank + 1).into(), (rank + 1).into(), k.0.div_ceil(base2k.0) as usize, ), k, base2k, dsize, } } pub fn alloc_bytes(infos: &A) -> usize where A: GGSWInfos, { Self::alloc_bytes_with( infos.n(), infos.base2k(), infos.k(), infos.rank(), infos.dnum(), infos.dsize(), ) } pub fn alloc_bytes_with(n: Degree, base2k: Base2K, k: TorusPrecision, rank: Rank, dnum: Dnum, dsize: Dsize) -> usize { let size: usize = k.0.div_ceil(base2k.0) as usize; debug_assert!( size as u32 > dsize.0, "invalid ggsw: ceil(k/base2k): {size} <= dsize: {}", dsize.0 ); assert!( dnum.0 * dsize.0 <= size as u32, "invalid ggsw: dnum: {} * dsize:{} > ceil(k/base2k): {size}", dnum.0, dsize.0, ); MatZnx::alloc_bytes( n.into(), dnum.into(), (rank + 1).into(), (rank + 1).into(), k.0.div_ceil(base2k.0) as usize, ) } } use byteorder::{LittleEndian, ReadBytesExt, WriteBytesExt}; impl ReaderFrom for GGSWCiphertext { fn read_from(&mut self, reader: &mut R) -> std::io::Result<()> { self.k = TorusPrecision(reader.read_u32::()?); self.base2k = Base2K(reader.read_u32::()?); self.dsize = Dsize(reader.read_u32::()?); self.data.read_from(reader) } } impl WriterTo for GGSWCiphertext { fn write_to(&self, writer: &mut W) -> std::io::Result<()> { writer.write_u32::(self.k.into())?; writer.write_u32::(self.base2k.into())?; writer.write_u32::(self.dsize.into())?; self.data.write_to(writer) } }