use base2k::{ Backend, FFT64, MatZnxDft, MatZnxDftAlloc, MatZnxDftOps, MatZnxDftToMut, MatZnxDftToRef, Module, ScalarZnx, ScalarZnxDft, ScalarZnxDftToRef, ScalarZnxToRef, Scratch, VecZnxAlloc, VecZnxDft, VecZnxDftAlloc, VecZnxDftOps, VecZnxDftToMut, VecZnxDftToRef, VecZnxOps, ZnxInfos, ZnxZero, }; use sampling::source::Source; use crate::{ elem::{GetRow, Infos, SetRow}, glwe_ciphertext::GLWECiphertext, glwe_ciphertext_fourier::GLWECiphertextFourier, glwe_plaintext::GLWEPlaintext, keys::SecretKeyFourier, utils::derive_size, }; pub struct GGSWCiphertext { pub data: MatZnxDft, pub basek: usize, pub k: usize, } impl GGSWCiphertext, B> { pub fn new(module: &Module, basek: usize, k: usize, rows: usize, rank: usize) -> Self { Self { data: module.new_mat_znx_dft(rows, rank + 1, rank + 1, derive_size(basek, k)), basek: basek, k: k, } } } impl Infos for GGSWCiphertext { type Inner = MatZnxDft; fn inner(&self) -> &Self::Inner { &self.data } fn basek(&self) -> usize { self.basek } fn k(&self) -> usize { self.k } } impl GGSWCiphertext { pub fn rank(&self) -> usize { self.data.cols_out() - 1 } } impl MatZnxDftToMut for GGSWCiphertext where MatZnxDft: MatZnxDftToMut, { fn to_mut(&mut self) -> MatZnxDft<&mut [u8], B> { self.data.to_mut() } } impl MatZnxDftToRef for GGSWCiphertext where MatZnxDft: MatZnxDftToRef, { fn to_ref(&self) -> MatZnxDft<&[u8], B> { self.data.to_ref() } } impl GGSWCiphertext, FFT64> { pub fn encrypt_sk_scratch_space(module: &Module, rank: usize, size: usize) -> usize { GLWECiphertext::encrypt_sk_scratch_space(module, size) + module.bytes_of_vec_znx(rank + 1, size) + module.bytes_of_vec_znx(1, size) + module.bytes_of_vec_znx_dft(rank + 1, size) } pub fn external_product_scratch_space( module: &Module, out_size: usize, in_size: usize, ggsw_size: usize, rank: usize, ) -> usize { let tmp_in: usize = module.bytes_of_vec_znx_dft(rank + 1, in_size); let tmp_out: usize = module.bytes_of_vec_znx_dft(rank + 1, out_size); let ggsw: usize = GLWECiphertextFourier::external_product_scratch_space(module, out_size, in_size, ggsw_size, rank); tmp_in + tmp_out + ggsw } pub fn external_product_inplace_scratch_space( module: &Module, out_size: usize, ggsw_size: usize, rank: usize, ) -> usize { let tmp: usize = module.bytes_of_vec_znx_dft(rank + 1, out_size); let ggsw: usize = GLWECiphertextFourier::external_product_inplace_scratch_space(module, out_size, ggsw_size, rank); tmp + ggsw } } impl GGSWCiphertext where MatZnxDft: MatZnxDftToMut + MatZnxDftToRef, { pub fn encrypt_sk( &mut self, module: &Module, pt: &ScalarZnx, sk_dft: &SecretKeyFourier, source_xa: &mut Source, source_xe: &mut Source, sigma: f64, scratch: &mut Scratch, ) where ScalarZnx: ScalarZnxToRef, ScalarZnxDft: ScalarZnxDftToRef, { #[cfg(debug_assertions)] { assert_eq!(self.rank(), sk_dft.rank()); assert_eq!(self.n(), module.n()); assert_eq!(pt.n(), module.n()); assert_eq!(sk_dft.n(), module.n()); } let size: usize = self.size(); let basek: usize = self.basek(); let k: usize = self.k(); let cols: usize = self.rank() + 1; let (tmp_znx_pt, scratch_1) = scratch.tmp_vec_znx(module, 1, size); let (tmp_znx_ct, scrach_2) = scratch_1.tmp_vec_znx(module, cols, size); let mut vec_znx_pt: GLWEPlaintext<&mut [u8]> = GLWEPlaintext { data: tmp_znx_pt, basek: basek, k: k, }; let mut vec_znx_ct: GLWECiphertext<&mut [u8]> = GLWECiphertext { data: tmp_znx_ct, basek: basek, k, }; (0..self.rows()).for_each(|row_j| { // Adds the scalar_znx_pt to the i-th limb of the vec_znx_pt module.vec_znx_add_scalar_inplace(&mut vec_znx_pt, 0, row_j, pt, 0); module.vec_znx_normalize_inplace(basek, &mut vec_znx_pt, 0, scrach_2); (0..cols).for_each(|col_i| { // rlwe encrypt of vec_znx_pt into vec_znx_ct vec_znx_ct.encrypt_sk_private( module, Some((&vec_znx_pt, col_i)), sk_dft, source_xa, source_xe, sigma, scrach_2, ); // Switch vec_znx_ct into DFT domain { let (mut vec_znx_dft_ct, _) = scrach_2.tmp_vec_znx_dft(module, cols, size); (0..cols).for_each(|i| { module.vec_znx_dft(&mut vec_znx_dft_ct, i, &vec_znx_ct, i); }); module.vmp_prepare_row(self, row_j, col_i, &vec_znx_dft_ct); } }); vec_znx_pt.data.zero(); // zeroes for next iteration }); } pub fn external_product( &mut self, module: &Module, lhs: &GGSWCiphertext, rhs: &GGSWCiphertext, scratch: &mut Scratch, ) where MatZnxDft: MatZnxDftToRef, MatZnxDft: MatZnxDftToRef, { #[cfg(debug_assertions)] { assert_eq!( self.rank(), lhs.rank(), "ggsw_out rank: {} != ggsw_in rank: {}", self.rank(), lhs.rank() ); assert_eq!( self.rank(), rhs.rank(), "ggsw_in rank: {} != ggsw_apply rank: {}", self.rank(), rhs.rank() ); } let (tmp_in_data, scratch1) = scratch.tmp_vec_znx_dft(module, lhs.rank() + 1, lhs.size()); let mut tmp_in: GLWECiphertextFourier<&mut [u8], FFT64> = GLWECiphertextFourier::<&mut [u8], FFT64> { data: tmp_in_data, basek: lhs.basek(), k: lhs.k(), }; let (tmp_out_data, scratch2) = scratch1.tmp_vec_znx_dft(module, self.rank() + 1, self.size()); let mut tmp_out: GLWECiphertextFourier<&mut [u8], FFT64> = GLWECiphertextFourier::<&mut [u8], FFT64> { data: tmp_out_data, basek: self.basek(), k: self.k(), }; (0..self.rank() + 1).for_each(|col_i| { (0..self.rows()).for_each(|row_j| { lhs.get_row(module, row_j, col_i, &mut tmp_in); tmp_out.external_product(module, &tmp_in, rhs, scratch2); self.set_row(module, row_j, col_i, &tmp_out); }); }); tmp_out.data.zero(); (self.rows().min(lhs.rows())..self.rows()).for_each(|row_i| { (0..self.rank() + 1).for_each(|col_j| { self.set_row(module, row_i, col_j, &tmp_out); }); }); } pub fn external_product_inplace( &mut self, module: &Module, rhs: &GGSWCiphertext, scratch: &mut Scratch, ) where MatZnxDft: MatZnxDftToRef, { #[cfg(debug_assertions)] { assert_eq!( self.rank(), rhs.rank(), "ggsw_out rank: {} != ggsw_apply: {}", self.rank(), rhs.rank() ); } let (tmp_data, scratch1) = scratch.tmp_vec_znx_dft(module, self.rank() + 1, self.size()); let mut tmp: GLWECiphertextFourier<&mut [u8], FFT64> = GLWECiphertextFourier::<&mut [u8], FFT64> { data: tmp_data, basek: self.basek(), k: self.k(), }; (0..self.rank() + 1).for_each(|col_i| { (0..self.rows()).for_each(|row_j| { self.get_row(module, row_j, col_i, &mut tmp); tmp.external_product_inplace(module, rhs, scratch1); self.set_row(module, row_j, col_i, &tmp); }); }); } } impl GetRow for GGSWCiphertext where MatZnxDft: MatZnxDftToRef, { fn get_row( &self, module: &Module, row_i: usize, col_j: usize, res: &mut GLWECiphertextFourier, ) where VecZnxDft: VecZnxDftToMut, { module.vmp_extract_row(res, self, row_i, col_j); } } impl SetRow for GGSWCiphertext where MatZnxDft: MatZnxDftToMut, { fn set_row(&mut self, module: &Module, row_i: usize, col_j: usize, a: &GLWECiphertextFourier) where VecZnxDft: VecZnxDftToRef, { module.vmp_prepare_row(self, row_i, col_j, a); } }