poulpy/core/src/ggsw.rs
core/src/ggsw.rsRaw
1use base2k::{
2 Backend, FFT64, MatZnxDft, MatZnxDftAlloc, MatZnxDftOps, MatZnxDftScratch, MatZnxDftToMut, MatZnxDftToRef, Module, ScalarZnx,
3 ScalarZnxDft, ScalarZnxDftToRef, ScalarZnxToRef, Scratch, VecZnx, VecZnxAlloc, VecZnxBig, VecZnxBigOps, VecZnxBigScratch,
4 VecZnxDft, VecZnxDftAlloc, VecZnxDftOps, VecZnxDftToMut, VecZnxDftToRef, VecZnxOps, VecZnxToMut, VecZnxToRef, ZnxInfos,
5 ZnxZero,
6};
7use sampling::source::Source;
8
9use crate::{
10 elem::{GetRow, Infos, SetRow},
11 encryption::EncryptSkScratchSpace,
12 external_product::{
13 ExternalProduct, ExternalProductInplace, ExternalProductInplaceScratchSpace, ExternalProductScratchSpace,
14 },
15 glwe::{GLWECiphertext, GLWECiphertextFourier, GLWEPlaintext, encrypt_glwe_sk},
16 keys::SecretKeyFourier,
17 keyswitch::{KeySwitch, KeySwitchInplace, KeySwitchInplaceScratchSpace, KeySwitchScratchSpace},
18 keyswitch_key::GLWEKeySwitchKey,
19 utils::derive_size,
20 vec_glwe_product::{VecGLWEProduct, VecGLWEProductScratchSpace},
21};
22
23pub struct GGSWCiphertext<C, B: Backend> {
24 pub data: MatZnxDft<C, B>,
25 pub log_base2k: usize,
26 pub log_k: usize,
27}
28
29impl<B: Backend> GGSWCiphertext<Vec<u8>, B> {
30 pub fn new(module: &Module<B>, log_base2k: usize, log_k: usize, rows: usize) -> Self {
31 Self {
32 data: module.new_mat_znx_dft(rows, 2, 2, derive_size(log_base2k, log_k)),
33 log_base2k: log_base2k,
34 log_k: log_k,
35 }
36 }
37}
38
39impl<T, B: Backend> Infos for GGSWCiphertext<T, B> {
40 type Inner = MatZnxDft<T, B>;
41
42 fn inner(&self) -> &Self::Inner {
43 &self.data
44 }
45
46 fn basek(&self) -> usize {
47 self.log_base2k
48 }
49
50 fn k(&self) -> usize {
51 self.log_k
52 }
53}
54
55impl<C, B: Backend> MatZnxDftToMut<B> for GGSWCiphertext<C, B>
56where
57 MatZnxDft<C, B>: MatZnxDftToMut<B>,
58{
59 fn to_mut(&mut self) -> MatZnxDft<&mut [u8], B> {
60 self.data.to_mut()
61 }
62}
63
64impl<C, B: Backend> MatZnxDftToRef<B> for GGSWCiphertext<C, B>
65where
66 MatZnxDft<C, B>: MatZnxDftToRef<B>,
67{
68 fn to_ref(&self) -> MatZnxDft<&[u8], B> {
69 self.data.to_ref()
70 }
71}
72
73impl GGSWCiphertext<Vec<u8>, FFT64> {
74 pub fn encrypt_sk_scratch_space(module: &Module<FFT64>, size: usize) -> usize {
75 GLWECiphertext::encrypt_sk_scratch_space(module, size)
76 + module.bytes_of_vec_znx(2, size)
77 + module.bytes_of_vec_znx(1, size)
78 + module.bytes_of_vec_znx_dft(2, size)
79 }
80}
81
82pub fn encrypt_rgsw_sk<C, P, S>(
83 module: &Module<FFT64>,
84 ct: &mut GGSWCiphertext<C, FFT64>,
85 pt: &ScalarZnx<P>,
86 sk_dft: &SecretKeyFourier<S, FFT64>,
87 source_xa: &mut Source,
88 source_xe: &mut Source,
89 sigma: f64,
90 bound: f64,
91 scratch: &mut Scratch,
92) where
93 MatZnxDft<C, FFT64>: MatZnxDftToMut<FFT64>,
94 ScalarZnx<P>: ScalarZnxToRef,
95 ScalarZnxDft<S, FFT64>: ScalarZnxDftToRef<FFT64>,
96{
97 let size: usize = ct.size();
98 let log_base2k: usize = ct.basek();
99
100 let (tmp_znx_pt, scratch_1) = scratch.tmp_vec_znx(module, 1, size);
101 let (tmp_znx_ct, scrach_2) = scratch_1.tmp_vec_znx(module, 2, size);
102
103 let mut vec_znx_pt: GLWEPlaintext<&mut [u8]> = GLWEPlaintext {
104 data: tmp_znx_pt,
105 log_base2k: log_base2k,
106 log_k: ct.k(),
107 };
108
109 let mut vec_znx_ct: GLWECiphertext<&mut [u8]> = GLWECiphertext {
110 data: tmp_znx_ct,
111 log_base2k: log_base2k,
112 log_k: ct.k(),
113 };
114
115 (0..ct.rows()).for_each(|row_j| {
116 // Adds the scalar_znx_pt to the i-th limb of the vec_znx_pt
117 module.vec_znx_add_scalar_inplace(&mut vec_znx_pt, 0, row_j, pt, 0);
118 module.vec_znx_normalize_inplace(log_base2k, &mut vec_znx_pt, 0, scrach_2);
119
120 (0..ct.rank()).for_each(|col_i| {
121 // rlwe encrypt of vec_znx_pt into vec_znx_ct
122 encrypt_glwe_sk(
123 module,
124 &mut vec_znx_ct,
125 Some((&vec_znx_pt, col_i)),
126 sk_dft,
127 source_xa,
128 source_xe,
129 sigma,
130 bound,
131 scrach_2,
132 );
133
134 // Switch vec_znx_ct into DFT domain
135 {
136 let (mut vec_znx_dft_ct, _) = scrach_2.tmp_vec_znx_dft(module, 2, size);
137 module.vec_znx_dft(&mut vec_znx_dft_ct, 0, &vec_znx_ct, 0);
138 module.vec_znx_dft(&mut vec_znx_dft_ct, 1, &vec_znx_ct, 1);
139 module.vmp_prepare_row(ct, row_j, col_i, &vec_znx_dft_ct);
140 }
141 });
142
143 vec_znx_pt.data.zero(); // zeroes for next iteration
144 });
145}
146
147impl<C> GGSWCiphertext<C, FFT64> {
148 pub fn encrypt_sk<P, S>(
149 &mut self,
150 module: &Module<FFT64>,
151 pt: &ScalarZnx<P>,
152 sk_dft: &SecretKeyFourier<S, FFT64>,
153 source_xa: &mut Source,
154 source_xe: &mut Source,
155 sigma: f64,
156 bound: f64,
157 scratch: &mut Scratch,
158 ) where
159 MatZnxDft<C, FFT64>: MatZnxDftToMut<FFT64>,
160 ScalarZnx<P>: ScalarZnxToRef,
161 ScalarZnxDft<S, FFT64>: ScalarZnxDftToRef<FFT64>,
162 {
163 encrypt_rgsw_sk(
164 module, self, pt, sk_dft, source_xa, source_xe, sigma, bound, scratch,
165 )
166 }
167}
168
169impl<C> GetRow<FFT64> for GGSWCiphertext<C, FFT64>
170where
171 MatZnxDft<C, FFT64>: MatZnxDftToRef<FFT64>,
172{
173 fn get_row<R>(&self, module: &Module<FFT64>, row_i: usize, col_j: usize, res: &mut GLWECiphertextFourier<R, FFT64>)
174 where
175 VecZnxDft<R, FFT64>: VecZnxDftToMut<FFT64>,
176 {
177 module.vmp_extract_row(res, self, row_i, col_j);
178 }
179}
180
181impl<C> SetRow<FFT64> for GGSWCiphertext<C, FFT64>
182where
183 MatZnxDft<C, FFT64>: MatZnxDftToMut<FFT64>,
184{
185 fn set_row<R>(&mut self, module: &Module<FFT64>, row_i: usize, col_j: usize, a: &GLWECiphertextFourier<R, FFT64>)
186 where
187 VecZnxDft<R, FFT64>: VecZnxDftToRef<FFT64>,
188 {
189 module.vmp_prepare_row(self, row_i, col_j, a);
190 }
191}
192
193impl KeySwitchScratchSpace for GGSWCiphertext<Vec<u8>, FFT64> {
194 fn keyswitch_scratch_space(module: &Module<FFT64>, res_size: usize, lhs: usize, rhs: usize) -> usize {
195 <GLWEKeySwitchKey<Vec<u8>, FFT64> as VecGLWEProductScratchSpace>::prod_with_vec_glwe_scratch_space(
196 module, res_size, lhs, rhs,
197 )
198 }
199}
200
201impl<DataSelf, DataLhs, DataRhs> KeySwitch<DataLhs, DataRhs> for GGSWCiphertext<DataSelf, FFT64>
202where
203 MatZnxDft<DataSelf, FFT64>: MatZnxDftToMut<FFT64> + MatZnxDftToRef<FFT64>,
204 MatZnxDft<DataLhs, FFT64>: MatZnxDftToRef<FFT64>,
205 MatZnxDft<DataRhs, FFT64>: MatZnxDftToRef<FFT64>,
206{
207 type Lhs = GGSWCiphertext<DataLhs, FFT64>;
208 type Rhs = GLWEKeySwitchKey<DataRhs, FFT64>;
209
210 fn keyswitch(&mut self, module: &Module<FFT64>, lhs: &Self::Lhs, rhs: &Self::Rhs, scratch: &mut Scratch) {
211 rhs.prod_with_vec_glwe(module, self, lhs, scratch);
212 }
213}
214
215impl KeySwitchInplaceScratchSpace for GGSWCiphertext<Vec<u8>, FFT64> {
216 fn keyswitch_inplace_scratch_space(module: &Module<FFT64>, res_size: usize, rhs: usize) -> usize {
217 <GLWEKeySwitchKey<Vec<u8>, FFT64> as VecGLWEProductScratchSpace>::prod_with_vec_glwe_inplace_scratch_space(
218 module, res_size, rhs,
219 )
220 }
221}
222
223impl<DataSelf, DataRhs> KeySwitchInplace<DataRhs> for GGSWCiphertext<DataSelf, FFT64>
224where
225 MatZnxDft<DataSelf, FFT64>: MatZnxDftToMut<FFT64> + MatZnxDftToRef<FFT64>,
226 MatZnxDft<DataRhs, FFT64>: MatZnxDftToRef<FFT64>,
227{
228 type Rhs = GLWEKeySwitchKey<DataRhs, FFT64>;
229
230 fn keyswitch_inplace(&mut self, module: &Module<FFT64>, rhs: &Self::Rhs, scratch: &mut Scratch) {
231 rhs.prod_with_vec_glwe(module, self, rhs, scratch);
232 }
233}
234
235impl ExternalProductScratchSpace for GGSWCiphertext<Vec<u8>, FFT64> {
236 fn external_product_scratch_space(module: &Module<FFT64>, res_size: usize, lhs: usize, rhs: usize) -> usize {
237 <GGSWCiphertext<Vec<u8>, FFT64> as VecGLWEProductScratchSpace>::prod_with_vec_glwe_scratch_space(
238 module, res_size, lhs, rhs,
239 )
240 }
241}
242
243impl<DataSelf, DataLhs, DataRhs> ExternalProduct<DataLhs, DataRhs> for GGSWCiphertext<DataSelf, FFT64>
244where
245 MatZnxDft<DataSelf, FFT64>: MatZnxDftToMut<FFT64> + MatZnxDftToRef<FFT64>,
246 MatZnxDft<DataLhs, FFT64>: MatZnxDftToRef<FFT64>,
247 MatZnxDft<DataRhs, FFT64>: MatZnxDftToRef<FFT64>,
248{
249 type Lhs = GGSWCiphertext<DataLhs, FFT64>;
250 type Rhs = GGSWCiphertext<DataRhs, FFT64>;
251
252 fn external_product(&mut self, module: &Module<FFT64>, lhs: &Self::Lhs, rhs: &Self::Rhs, scratch: &mut Scratch) {
253 rhs.prod_with_vec_glwe(module, self, lhs, scratch);
254 }
255}
256
257impl ExternalProductInplaceScratchSpace for GGSWCiphertext<Vec<u8>, FFT64> {
258 fn external_product_inplace_scratch_space(module: &Module<FFT64>, res_size: usize, rhs: usize) -> usize {
259 <GGSWCiphertext<Vec<u8>, FFT64> as VecGLWEProductScratchSpace>::prod_with_glwe_inplace_scratch_space(
260 module, res_size, rhs,
261 )
262 }
263}
264
265impl<DataSelf, DataRhs> ExternalProductInplace<DataRhs> for GGSWCiphertext<DataSelf, FFT64>
266where
267 MatZnxDft<DataSelf, FFT64>: MatZnxDftToMut<FFT64> + MatZnxDftToRef<FFT64>,
268 MatZnxDft<DataRhs, FFT64>: MatZnxDftToRef<FFT64>,
269{
270 type Rhs = GGSWCiphertext<DataRhs, FFT64>;
271
272 fn external_product_inplace(&mut self, module: &Module<FFT64>, rhs: &Self::Rhs, scratch: &mut Scratch) {
273 rhs.prod_with_vec_glwe_inplace(module, self, scratch);
274 }
275}
276
277impl VecGLWEProductScratchSpace for GGSWCiphertext<Vec<u8>, FFT64> {
278 fn prod_with_glwe_scratch_space(module: &Module<FFT64>, res_size: usize, a_size: usize, rgsw_size: usize) -> usize {
279 module.bytes_of_vec_znx_dft(2, rgsw_size)
280 + ((module.bytes_of_vec_znx_dft(2, a_size) + module.vmp_apply_tmp_bytes(res_size, a_size, a_size, 2, 2, rgsw_size))
281 | module.vec_znx_big_normalize_tmp_bytes())
282 }
283}
284
285impl<C> VecGLWEProduct for GGSWCiphertext<C, FFT64>
286where
287 MatZnxDft<C, FFT64>: MatZnxDftToRef<FFT64> + ZnxInfos,
288{
289 fn prod_with_glwe<R, A>(
290 &self,
291 module: &Module<FFT64>,
292 res: &mut GLWECiphertext<R>,
293 a: &GLWECiphertext<A>,
294 scratch: &mut Scratch,
295 ) where
296 VecZnx<R>: VecZnxToMut,
297 VecZnx<A>: VecZnxToRef,
298 {
299 let log_base2k: usize = self.basek();
300
301 #[cfg(debug_assertions)]
302 {
303 assert_eq!(res.basek(), log_base2k);
304 assert_eq!(a.basek(), log_base2k);
305 assert_eq!(self.n(), module.n());
306 assert_eq!(res.n(), module.n());
307 assert_eq!(a.n(), module.n());
308 }
309
310 let (mut res_dft, scratch1) = scratch.tmp_vec_znx_dft(module, 2, self.size()); // Todo optimise
311
312 {
313 let (mut a_dft, scratch2) = scratch1.tmp_vec_znx_dft(module, 2, a.size());
314 module.vec_znx_dft(&mut a_dft, 0, a, 0);
315 module.vec_znx_dft(&mut a_dft, 1, a, 1);
316 module.vmp_apply(&mut res_dft, &a_dft, self, scratch2);
317 }
318
319 let res_big: VecZnxBig<&mut [u8], FFT64> = module.vec_znx_idft_consume(res_dft);
320
321 module.vec_znx_big_normalize(log_base2k, res, 0, &res_big, 0, scratch1);
322 module.vec_znx_big_normalize(log_base2k, res, 1, &res_big, 1, scratch1);
323 }
324}