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}, encryption::EncryptSkScratchSpace, external_product::{ ExternalProduct, ExternalProductInplace, ExternalProductInplaceScratchSpace, ExternalProductScratchSpace, }, glwe::{GLWECiphertext, GLWECiphertextFourier, GLWEPlaintext, encrypt_glwe_sk}, keys::SecretKeyFourier, keyswitch::{KeySwitch, KeySwitchInplace, KeySwitchInplaceScratchSpace, KeySwitchScratchSpace}, keyswitch_key::GLWEKeySwitchKey, utils::derive_size, vec_glwe_product::{VecGLWEProduct, VecGLWEProductScratchSpace}, }; pub struct GGSWCiphertext { pub data: MatZnxDft, pub log_base2k: usize, pub log_k: usize, } impl GGSWCiphertext, B> { pub fn new(module: &Module, log_base2k: usize, log_k: usize, rows: usize) -> Self { Self { data: module.new_mat_znx_dft(rows, 2, 2, derive_size(log_base2k, log_k)), log_base2k: log_base2k, log_k: log_k, } } } impl Infos for GGSWCiphertext { type Inner = MatZnxDft; fn inner(&self) -> &Self::Inner { &self.data } fn basek(&self) -> usize { self.log_base2k } fn k(&self) -> usize { self.log_k } } 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, size: usize) -> usize { GLWECiphertext::encrypt_sk_scratch_space(module, size) + module.bytes_of_vec_znx(2, size) + module.bytes_of_vec_znx(1, size) + module.bytes_of_vec_znx_dft(2, size) } } pub fn encrypt_rgsw_sk( module: &Module, ct: &mut GGSWCiphertext, pt: &ScalarZnx

, sk_dft: &SecretKeyFourier, source_xa: &mut Source, source_xe: &mut Source, sigma: f64, bound: f64, scratch: &mut Scratch, ) where MatZnxDft: MatZnxDftToMut, ScalarZnx

: ScalarZnxToRef, ScalarZnxDft: ScalarZnxDftToRef, { let size: usize = ct.size(); let log_base2k: usize = ct.basek(); 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, 2, size); let mut vec_znx_pt: GLWEPlaintext<&mut [u8]> = GLWEPlaintext { data: tmp_znx_pt, log_base2k: log_base2k, log_k: ct.k(), }; let mut vec_znx_ct: GLWECiphertext<&mut [u8]> = GLWECiphertext { data: tmp_znx_ct, log_base2k: log_base2k, log_k: ct.k(), }; (0..ct.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(log_base2k, &mut vec_znx_pt, 0, scrach_2); (0..ct.rank()).for_each(|col_i| { // rlwe encrypt of vec_znx_pt into vec_znx_ct encrypt_glwe_sk( module, &mut vec_znx_ct, Some((&vec_znx_pt, col_i)), sk_dft, source_xa, source_xe, sigma, bound, scrach_2, ); // Switch vec_znx_ct into DFT domain { let (mut vec_znx_dft_ct, _) = scrach_2.tmp_vec_znx_dft(module, 2, size); module.vec_znx_dft(&mut vec_znx_dft_ct, 0, &vec_znx_ct, 0); module.vec_znx_dft(&mut vec_znx_dft_ct, 1, &vec_znx_ct, 1); module.vmp_prepare_row(ct, row_j, col_i, &vec_znx_dft_ct); } }); vec_znx_pt.data.zero(); // zeroes for next iteration }); } impl GGSWCiphertext { 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 MatZnxDft: MatZnxDftToMut, ScalarZnx

: ScalarZnxToRef, ScalarZnxDft: ScalarZnxDftToRef, { encrypt_rgsw_sk( module, self, pt, sk_dft, source_xa, source_xe, sigma, bound, scratch, ) } } 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); } } impl KeySwitchScratchSpace for GGSWCiphertext, FFT64> { fn keyswitch_scratch_space(module: &Module, res_size: usize, lhs: usize, rhs: usize) -> usize { , FFT64> as VecGLWEProductScratchSpace>::prod_with_vec_glwe_scratch_space( module, res_size, lhs, rhs, ) } } impl KeySwitch for GGSWCiphertext where MatZnxDft: MatZnxDftToMut + MatZnxDftToRef, MatZnxDft: MatZnxDftToRef, MatZnxDft: MatZnxDftToRef, { type Lhs = GGSWCiphertext; type Rhs = GLWEKeySwitchKey; fn keyswitch(&mut self, module: &Module, lhs: &Self::Lhs, rhs: &Self::Rhs, scratch: &mut Scratch) { rhs.prod_with_vec_glwe(module, self, lhs, scratch); } } impl KeySwitchInplaceScratchSpace for GGSWCiphertext, FFT64> { fn keyswitch_inplace_scratch_space(module: &Module, res_size: usize, rhs: usize) -> usize { , FFT64> as VecGLWEProductScratchSpace>::prod_with_vec_glwe_inplace_scratch_space( module, res_size, rhs, ) } } impl KeySwitchInplace for GGSWCiphertext where MatZnxDft: MatZnxDftToMut + MatZnxDftToRef, MatZnxDft: MatZnxDftToRef, { type Rhs = GLWEKeySwitchKey; fn keyswitch_inplace(&mut self, module: &Module, rhs: &Self::Rhs, scratch: &mut Scratch) { rhs.prod_with_vec_glwe(module, self, rhs, scratch); } } impl ExternalProductScratchSpace for GGSWCiphertext, FFT64> { fn external_product_scratch_space(module: &Module, res_size: usize, lhs: usize, rhs: usize) -> usize { , FFT64> as VecGLWEProductScratchSpace>::prod_with_vec_glwe_scratch_space( module, res_size, lhs, rhs, ) } } impl ExternalProduct for GGSWCiphertext where MatZnxDft: MatZnxDftToMut + MatZnxDftToRef, MatZnxDft: MatZnxDftToRef, MatZnxDft: MatZnxDftToRef, { type Lhs = GGSWCiphertext; type Rhs = GGSWCiphertext; fn external_product(&mut self, module: &Module, lhs: &Self::Lhs, rhs: &Self::Rhs, scratch: &mut Scratch) { rhs.prod_with_vec_glwe(module, self, lhs, scratch); } } impl ExternalProductInplaceScratchSpace for GGSWCiphertext, FFT64> { fn external_product_inplace_scratch_space(module: &Module, res_size: usize, rhs: usize) -> usize { , FFT64> as VecGLWEProductScratchSpace>::prod_with_glwe_inplace_scratch_space( module, res_size, rhs, ) } } impl ExternalProductInplace for GGSWCiphertext where MatZnxDft: MatZnxDftToMut + MatZnxDftToRef, MatZnxDft: MatZnxDftToRef, { type Rhs = GGSWCiphertext; fn external_product_inplace(&mut self, module: &Module, rhs: &Self::Rhs, scratch: &mut Scratch) { rhs.prod_with_vec_glwe_inplace(module, self, scratch); } } impl VecGLWEProductScratchSpace for GGSWCiphertext, FFT64> { fn prod_with_glwe_scratch_space(module: &Module, res_size: usize, a_size: usize, rgsw_size: usize) -> usize { module.bytes_of_vec_znx_dft(2, rgsw_size) + ((module.bytes_of_vec_znx_dft(2, a_size) + module.vmp_apply_tmp_bytes(res_size, a_size, a_size, 2, 2, rgsw_size)) | module.vec_znx_big_normalize_tmp_bytes()) } } impl VecGLWEProduct for GGSWCiphertext where MatZnxDft: MatZnxDftToRef + ZnxInfos, { fn prod_with_glwe( &self, module: &Module, res: &mut GLWECiphertext, a: &GLWECiphertext, scratch: &mut Scratch, ) where 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 a_dft, scratch2) = scratch1.tmp_vec_znx_dft(module, 2, a.size()); module.vec_znx_dft(&mut a_dft, 0, a, 0); module.vec_znx_dft(&mut a_dft, 1, a, 1); module.vmp_apply(&mut res_dft, &a_dft, self, scratch2); } let res_big: VecZnxBig<&mut [u8], FFT64> = module.vec_znx_idft_consume(res_dft); 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); } }