use base2k::{ Backend, FFT64, MatZnxDft, MatZnxDftAlloc, MatZnxDftOps, MatZnxDftScratch, MatZnxDftToMut, MatZnxDftToRef, Module, ScalarZnx, ScalarZnxDft, ScalarZnxDftToRef, ScalarZnxToRef, Scratch, VecZnx, VecZnxAlloc, VecZnxBig, VecZnxBigOps, VecZnxBigScratch, VecZnxDft, VecZnxDftAlloc, VecZnxDftOps, VecZnxDftToMut, VecZnxDftToRef, VecZnxOps, VecZnxToMut, VecZnxToRef, 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, vec_glwe_product::{VecGLWEProduct, VecGLWEProductScratchSpace}, }; pub struct GGLWECiphertext { pub(crate) data: MatZnxDft, pub(crate) basek: usize, pub(crate) k: usize, } impl GGLWECiphertext, B> { pub fn new(module: &Module, base2k: usize, k: usize, rows: usize, rank_in: usize, rank_out: usize) -> Self { Self { data: module.new_mat_znx_dft(rows, rank_in, rank_out + 1, derive_size(base2k, k)), basek: base2k, k, } } } impl Infos for GGLWECiphertext { type Inner = MatZnxDft; fn inner(&self) -> &Self::Inner { &self.data } fn basek(&self) -> usize { self.basek } fn k(&self) -> usize { self.k } } impl GGLWECiphertext { pub fn rank(&self) -> usize { self.data.cols_out() - 1 } } impl MatZnxDftToMut for GGLWECiphertext where MatZnxDft: MatZnxDftToMut, { fn to_mut(&mut self) -> MatZnxDft<&mut [u8], B> { self.data.to_mut() } } impl MatZnxDftToRef for GGLWECiphertext where MatZnxDft: MatZnxDftToRef, { fn to_ref(&self) -> MatZnxDft<&[u8], B> { self.data.to_ref() } } impl GGLWECiphertext, FFT64> { pub fn encrypt_sk_scratch_space(module: &Module, rank: usize, size: usize) -> usize { GLWECiphertext::encrypt_sk_scratch_space(module, rank, 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 encrypt_pk_scratch_space(_module: &Module, _rank: usize, _pk_size: usize) -> usize { unimplemented!() } } impl GGLWECiphertext where MatZnxDft: MatZnxDftToMut + ZnxInfos, { pub fn encrypt_sk( &mut self, module: &Module, pt: &ScalarZnx, sk_dft: &SecretKeyFourier, source_xa: &mut Source, source_xe: &mut Source, sigma: f64, bound: 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!(sk_dft.n(), module.n()); assert_eq!(pt.n(), module.n()); } let rows: usize = self.rows(); 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, scrach_1) = scratch.tmp_vec_znx(module, 1, size); let (tmp_znx_ct, scrach_2) = scrach_1.tmp_vec_znx(module, cols, size); let (tmp_znx_dft_ct, scratch_3) = scrach_2.tmp_vec_znx_dft(module, cols, size); let mut vec_znx_pt: GLWEPlaintext<&mut [u8]> = GLWEPlaintext { data: tmp_znx_pt, basek, k, }; let mut vec_znx_ct: GLWECiphertext<&mut [u8]> = GLWECiphertext { data: tmp_znx_ct, basek, k, }; let mut vec_znx_ct_dft: GLWECiphertextFourier<&mut [u8], FFT64> = GLWECiphertextFourier { data: tmp_znx_dft_ct, basek, k, }; (0..rows).for_each(|row_i| { // 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_i, pt, 0); module.vec_znx_normalize_inplace(basek, &mut vec_znx_pt, 0, scratch_3); // rlwe encrypt of vec_znx_pt into vec_znx_ct vec_znx_ct.encrypt_sk( module, &vec_znx_pt, sk_dft, source_xa, source_xe, sigma, bound, scratch_3, ); vec_znx_pt.data.zero(); // zeroes for next iteration // Switch vec_znx_ct into DFT domain vec_znx_ct.dft(module, &mut vec_znx_ct_dft); // Stores vec_znx_dft_ct into thw i-th row of the MatZnxDft module.vmp_prepare_row(self, row_i, 0, &vec_znx_ct_dft); }); } } impl GetRow for GGLWECiphertext where MatZnxDft: MatZnxDftToRef, { fn get_row(&self, module: &Module, row_i: usize, col_j: usize, res: &mut GLWECiphertextFourier) where VecZnxDft: VecZnxDftToMut, { #[cfg(debug_assertions)] { assert_eq!(col_j, 0); } module.vmp_extract_row(res, self, row_i, col_j); } } impl SetRow for GGLWECiphertext where MatZnxDft: MatZnxDftToMut, { fn set_row(&mut self, module: &Module, row_i: usize, col_j: usize, a: &GLWECiphertextFourier) where VecZnxDft: VecZnxDftToRef, { #[cfg(debug_assertions)] { assert_eq!(col_j, 0); } module.vmp_prepare_row(self, row_i, col_j, a); } } impl VecGLWEProductScratchSpace for GGLWECiphertext, FFT64> { fn prod_with_glwe_scratch_space(module: &Module, res_size: usize, a_size: usize, grlwe_size: usize) -> usize { module.bytes_of_vec_znx_dft(2, grlwe_size) + (module.vec_znx_big_normalize_tmp_bytes() | (module.vmp_apply_tmp_bytes(res_size, a_size, a_size, 1, 2, grlwe_size) + module.bytes_of_vec_znx_dft(1, a_size))) } } impl VecGLWEProduct for GGLWECiphertext where MatZnxDft: MatZnxDftToRef + ZnxInfos, { fn prod_with_glwe( &self, module: &Module, res: &mut GLWECiphertext, a: &GLWECiphertext, scratch: &mut Scratch, ) where MatZnxDft: MatZnxDftToRef, VecZnx: VecZnxToMut, VecZnx: VecZnxToRef, { let log_base2k: usize = self.basek(); #[cfg(debug_assertions)] { assert_eq!(res.basek(), log_base2k); assert_eq!(a.basek(), log_base2k); assert_eq!(self.n(), module.n()); assert_eq!(res.n(), module.n()); assert_eq!(a.n(), module.n()); } let (mut res_dft, scratch1) = scratch.tmp_vec_znx_dft(module, 2, self.size()); // Todo optimise { let (mut a1_dft, scratch2) = scratch1.tmp_vec_znx_dft(module, 1, a.size()); module.vec_znx_dft(&mut a1_dft, 0, a, 1); module.vmp_apply(&mut res_dft, &a1_dft, self, scratch2); } let mut res_big: VecZnxBig<&mut [u8], FFT64> = module.vec_znx_idft_consume(res_dft); module.vec_znx_big_add_small_inplace(&mut res_big, 0, a, 0); module.vec_znx_big_normalize(log_base2k, res, 0, &res_big, 0, scratch1); module.vec_znx_big_normalize(log_base2k, res, 1, &res_big, 1, scratch1); } }