poulpy/poulpy-core/src/conversion/lwe_to_glwe.rs
poulpy-core/src/conversion/lwe_to_glwe.rsRaw
1use poulpy_hal::{
2 api::{ScratchTakeBasic, VecZnxNormalize, VecZnxNormalizeTmpBytes},
3 layouts::{Backend, DataMut, Module, Scratch, VecZnx, ZnxView, ZnxViewMut, ZnxZero},
4};
5
6use crate::{
7 GLWEKeyswitch, ScratchTakeCore,
8 layouts::{GGLWEInfos, GGLWEPreparedToRef, GLWE, GLWEInfos, GLWELayout, GLWEToMut, LWE, LWEInfos, LWEToRef},
9};
10
11impl<BE: Backend> GLWEFromLWE<BE> for Module<BE>
12where
13 Self: GLWEKeyswitch<BE> + VecZnxNormalizeTmpBytes + VecZnxNormalize<BE>,
14 Scratch<BE>: ScratchTakeCore<BE>,
15{
16 fn glwe_from_lwe_tmp_bytes<R, A, K>(&self, glwe_infos: &R, lwe_infos: &A, key_infos: &K) -> usize
17 where
18 R: GLWEInfos,
19 A: LWEInfos,
20 K: GGLWEInfos,
21 {
22 let ct: usize = GLWE::bytes_of(
23 self.n().into(),
24 key_infos.base2k(),
25 lwe_infos.k().max(glwe_infos.k()),
26 1u32.into(),
27 );
28
29 let ks: usize = self.glwe_keyswitch_tmp_bytes(glwe_infos, glwe_infos, key_infos);
30 if lwe_infos.base2k() == key_infos.base2k() {
31 ct + ks
32 } else {
33 let a_conv = VecZnx::bytes_of(self.n(), 1, lwe_infos.size()) + self.vec_znx_normalize_tmp_bytes();
34 ct + a_conv + ks
35 }
36 }
37
38 fn glwe_from_lwe<R, A, K>(&self, res: &mut R, lwe: &A, ksk: &K, scratch: &mut Scratch<BE>)
39 where
40 R: GLWEToMut,
41 A: LWEToRef,
42 K: GGLWEPreparedToRef<BE> + GGLWEInfos,
43 {
44 let res: &mut GLWE<&mut [u8]> = &mut res.to_mut();
45 let lwe: &LWE<&[u8]> = &lwe.to_ref();
46
47 assert_eq!(res.n(), self.n() as u32);
48 assert_eq!(ksk.n(), self.n() as u32);
49 assert!(lwe.n() <= self.n() as u32);
50
51 let (mut glwe, scratch_1) = scratch.take_glwe(&GLWELayout {
52 n: ksk.n(),
53 base2k: ksk.base2k(),
54 k: lwe.k(),
55 rank: 1u32.into(),
56 });
57 glwe.data.zero();
58
59 let n_lwe: usize = lwe.n().into();
60
61 if lwe.base2k() == ksk.base2k() {
62 for i in 0..lwe.size() {
63 let data_lwe: &[i64] = lwe.data.at(0, i);
64 glwe.data.at_mut(0, i)[0] = data_lwe[0];
65 glwe.data.at_mut(1, i)[..n_lwe].copy_from_slice(&data_lwe[1..]);
66 }
67 } else {
68 let (mut a_conv, scratch_2) = scratch_1.take_vec_znx(self.n(), 1, lwe.size());
69 a_conv.zero();
70 for j in 0..lwe.size() {
71 let data_lwe: &[i64] = lwe.data.at(0, j);
72 a_conv.at_mut(0, j)[0] = data_lwe[0]
73 }
74
75 self.vec_znx_normalize(
76 &mut glwe.data,
77 ksk.base2k().into(),
78 0,
79 0,
80 &a_conv,
81 lwe.base2k().into(),
82 0,
83 scratch_2,
84 );
85
86 a_conv.zero();
87 for j in 0..lwe.size() {
88 let data_lwe: &[i64] = lwe.data.at(0, j);
89 a_conv.at_mut(0, j)[..n_lwe].copy_from_slice(&data_lwe[1..]);
90 }
91
92 self.vec_znx_normalize(
93 &mut glwe.data,
94 ksk.base2k().into(),
95 0,
96 1,
97 &a_conv,
98 lwe.base2k().into(),
99 0,
100 scratch_2,
101 );
102 }
103
104 self.glwe_keyswitch(res, &glwe, ksk, scratch_1);
105 }
106}
107
108pub trait GLWEFromLWE<BE: Backend>
109where
110 Self: GLWEKeyswitch<BE>,
111{
112 fn glwe_from_lwe_tmp_bytes<R, A, K>(&self, glwe_infos: &R, lwe_infos: &A, key_infos: &K) -> usize
113 where
114 R: GLWEInfos,
115 A: LWEInfos,
116 K: GGLWEInfos;
117
118 fn glwe_from_lwe<R, A, K>(&self, res: &mut R, lwe: &A, ksk: &K, scratch: &mut Scratch<BE>)
119 where
120 R: GLWEToMut,
121 A: LWEToRef,
122 K: GGLWEPreparedToRef<BE> + GGLWEInfos;
123}
124
125impl GLWE<Vec<u8>> {
126 pub fn from_lwe_tmp_bytes<R, A, K, M, BE: Backend>(module: &M, glwe_infos: &R, lwe_infos: &A, key_infos: &K) -> usize
127 where
128 R: GLWEInfos,
129 A: LWEInfos,
130 K: GGLWEInfos,
131 M: GLWEFromLWE<BE>,
132 {
133 module.glwe_from_lwe_tmp_bytes(glwe_infos, lwe_infos, key_infos)
134 }
135}
136
137impl<D: DataMut> GLWE<D> {
138 pub fn from_lwe<A, K, M, BE: Backend>(&mut self, module: &M, lwe: &A, ksk: &K, scratch: &mut Scratch<BE>)
139 where
140 M: GLWEFromLWE<BE>,
141 A: LWEToRef,
142 K: GGLWEPreparedToRef<BE> + GGLWEInfos,
143 Scratch<BE>: ScratchTakeCore<BE>,
144 {
145 module.glwe_from_lwe(self, lwe, ksk, scratch);
146 }
147}