1use poulpy_hal::{
  2    api::{VmpPMatAlloc, VmpPMatBytesOf, VmpPrepare, VmpPrepareTmpBytes},
  3    layouts::{Backend, Data, DataMut, DataRef, Module, Scratch, VmpPMat, VmpPMatToMut, VmpPMatToRef, ZnxInfos},
  4};
  5
  6use crate::layouts::{
  7    Base2K, Degree, Dnum, Dsize, GGLWE, GGLWEInfos, GGLWEToRef, GLWEInfos, GetDegree, LWEInfos, Rank, TorusPrecision,
  8};
  9
 10#[derive(PartialEq, Eq)]
 11pub struct GGLWEPrepared<D: Data, B: Backend> {
 12    pub(crate) data: VmpPMat<D, B>,
 13    pub(crate) k: TorusPrecision,
 14    pub(crate) base2k: Base2K,
 15    pub(crate) dsize: Dsize,
 16}
 17
 18impl<D: Data, B: Backend> LWEInfos for GGLWEPrepared<D, B> {
 19    fn n(&self) -> Degree {
 20        Degree(self.data.n() as u32)
 21    }
 22
 23    fn base2k(&self) -> Base2K {
 24        self.base2k
 25    }
 26
 27    fn k(&self) -> TorusPrecision {
 28        self.k
 29    }
 30
 31    fn size(&self) -> usize {
 32        self.data.size()
 33    }
 34}
 35
 36impl<D: Data, B: Backend> GLWEInfos for GGLWEPrepared<D, B> {
 37    fn rank(&self) -> Rank {
 38        self.rank_out()
 39    }
 40}
 41
 42impl<D: Data, B: Backend> GGLWEInfos for GGLWEPrepared<D, B> {
 43    fn rank_in(&self) -> Rank {
 44        Rank(self.data.cols_in() as u32)
 45    }
 46
 47    fn rank_out(&self) -> Rank {
 48        Rank(self.data.cols_out() as u32 - 1)
 49    }
 50
 51    fn dsize(&self) -> Dsize {
 52        self.dsize
 53    }
 54
 55    fn dnum(&self) -> Dnum {
 56        Dnum(self.data.rows() as u32)
 57    }
 58}
 59
 60pub trait GGLWEPreparedAlloc<B: Backend>
 61where
 62    Self: GetDegree + VmpPMatAlloc<B> + VmpPMatBytesOf,
 63{
 64    fn alloc_gglwe_prepared(
 65        &self,
 66        base2k: Base2K,
 67        k: TorusPrecision,
 68        rank_in: Rank,
 69        rank_out: Rank,
 70        dnum: Dnum,
 71        dsize: Dsize,
 72    ) -> GGLWEPrepared<Vec<u8>, B> {
 73        let size: usize = k.0.div_ceil(base2k.0) as usize;
 74        debug_assert!(
 75            size as u32 > dsize.0,
 76            "invalid gglwe: ceil(k/base2k): {size} <= dsize: {}",
 77            dsize.0
 78        );
 79
 80        assert!(
 81            dnum.0 * dsize.0 <= size as u32,
 82            "invalid gglwe: dnum: {} * dsize:{} > ceil(k/base2k): {size}",
 83            dnum.0,
 84            dsize.0,
 85        );
 86
 87        GGLWEPrepared {
 88            data: self.vmp_pmat_alloc(dnum.into(), rank_in.into(), (rank_out + 1).into(), size),
 89            k,
 90            base2k,
 91            dsize,
 92        }
 93    }
 94
 95    fn alloc_gglwe_prepared_from_infos<A>(&self, infos: &A) -> GGLWEPrepared<Vec<u8>, B>
 96    where
 97        A: GGLWEInfos,
 98    {
 99        assert_eq!(self.n(), infos.n());
100        self.alloc_gglwe_prepared(
101            infos.base2k(),
102            infos.k(),
103            infos.rank_in(),
104            infos.rank_out(),
105            infos.dnum(),
106            infos.dsize(),
107        )
108    }
109
110    fn bytes_of_gglwe_prepared(
111        &self,
112        base2k: Base2K,
113        k: TorusPrecision,
114        rank_in: Rank,
115        rank_out: Rank,
116        dnum: Dnum,
117        dsize: Dsize,
118    ) -> usize {
119        let size: usize = k.0.div_ceil(base2k.0) as usize;
120        debug_assert!(
121            size as u32 > dsize.0,
122            "invalid gglwe: ceil(k/base2k): {size} <= dsize: {}",
123            dsize.0
124        );
125
126        assert!(
127            dnum.0 * dsize.0 <= size as u32,
128            "invalid gglwe: dnum: {} * dsize:{} > ceil(k/base2k): {size}",
129            dnum.0,
130            dsize.0,
131        );
132
133        self.bytes_of_vmp_pmat(dnum.into(), rank_in.into(), (rank_out + 1).into(), size)
134    }
135
136    fn bytes_of_gglwe_prepared_from_infos<A>(&self, infos: &A) -> usize
137    where
138        A: GGLWEInfos,
139    {
140        assert_eq!(self.n(), infos.n());
141        self.bytes_of_gglwe_prepared(
142            infos.base2k(),
143            infos.k(),
144            infos.rank_in(),
145            infos.rank_out(),
146            infos.dnum(),
147            infos.dsize(),
148        )
149    }
150}
151
152impl<B: Backend> GGLWEPreparedAlloc<B> for Module<B> where Module<B>: GetDegree + VmpPMatAlloc<B> + VmpPMatBytesOf {}
153
154impl<B: Backend> GGLWEPrepared<Vec<u8>, B>
155where
156    Module<B>: GGLWEPreparedAlloc<B>,
157{
158    pub fn alloc_from_infos<A>(module: &Module<B>, infos: &A) -> Self
159    where
160        A: GGLWEInfos,
161    {
162        module.alloc_gglwe_prepared_from_infos(infos)
163    }
164
165    pub fn alloc(
166        module: &Module<B>,
167        base2k: Base2K,
168        k: TorusPrecision,
169        rank_in: Rank,
170        rank_out: Rank,
171        dnum: Dnum,
172        dsize: Dsize,
173    ) -> Self {
174        module.alloc_gglwe_prepared(base2k, k, rank_in, rank_out, dnum, dsize)
175    }
176
177    pub fn bytes_of_from_infos<A>(module: &Module<B>, infos: &A) -> usize
178    where
179        A: GGLWEInfos,
180    {
181        module.bytes_of_gglwe_prepared_from_infos(infos)
182    }
183
184    pub fn bytes_of(
185        module: &Module<B>,
186        base2k: Base2K,
187        k: TorusPrecision,
188        rank_in: Rank,
189        rank_out: Rank,
190        dnum: Dnum,
191        dsize: Dsize,
192    ) -> usize {
193        module.bytes_of_gglwe_prepared(base2k, k, rank_in, rank_out, dnum, dsize)
194    }
195}
196
197pub trait GGLWEPrepare<B: Backend>
198where
199    Self: GetDegree + VmpPrepareTmpBytes + VmpPrepare<B>,
200{
201    fn prepare_gglwe_tmp_bytes<A>(&self, infos: &A) -> usize
202    where
203        A: GGLWEInfos,
204    {
205        self.vmp_prepare_tmp_bytes(
206            infos.dnum().into(),
207            infos.rank_in().into(),
208            (infos.rank() + 1).into(),
209            infos.size(),
210        )
211    }
212
213    fn prepare_gglwe<R, O>(&self, res: &mut R, other: &O, scratch: &mut Scratch<B>)
214    where
215        R: GGLWEPreparedToMut<B>,
216        O: GGLWEToRef,
217    {
218        let mut res: GGLWEPrepared<&mut [u8], B> = res.to_mut();
219        let other: GGLWE<&[u8]> = other.to_ref();
220
221        assert_eq!(res.n(), self.n());
222        assert_eq!(other.n(), self.n());
223        assert_eq!(res.base2k, other.base2k);
224        assert_eq!(res.k, other.k);
225        assert_eq!(res.dsize, other.dsize);
226
227        self.vmp_prepare(&mut res.data, &other.data, scratch);
228    }
229}
230
231impl<B: Backend> GGLWEPrepare<B> for Module<B> where Self: GetDegree + VmpPrepareTmpBytes + VmpPrepare<B> {}
232
233impl<D: DataMut, B: Backend> GGLWEPrepared<D, B>
234where
235    Module<B>: GGLWEPrepare<B>,
236{
237    pub fn prepare<O>(&mut self, module: &Module<B>, other: &O, scratch: &mut Scratch<B>)
238    where
239        O: GGLWEToRef,
240    {
241        module.prepare_gglwe(self, other, scratch);
242    }
243}
244
245impl<B: Backend> GGLWEPrepared<Vec<u8>, B> {
246    pub fn prepare_tmp_bytes(&self, module: &Module<B>) -> usize
247    where
248        Module<B>: GGLWEPrepare<B>,
249    {
250        module.prepare_gglwe_tmp_bytes(self)
251    }
252}
253
254pub trait GGLWEPreparedToMut<B: Backend> {
255    fn to_mut(&mut self) -> GGLWEPrepared<&mut [u8], B>;
256}
257
258impl<D: DataMut, B: Backend> GGLWEPreparedToMut<B> for GGLWEPrepared<D, B> {
259    fn to_mut(&mut self) -> GGLWEPrepared<&mut [u8], B> {
260        GGLWEPrepared {
261            k: self.k,
262            base2k: self.base2k,
263            dsize: self.dsize,
264            data: self.data.to_mut(),
265        }
266    }
267}
268
269pub trait GGLWEPreparedToRef<B: Backend> {
270    fn to_ref(&self) -> GGLWEPrepared<&[u8], B>;
271}
272
273impl<D: DataRef, B: Backend> GGLWEPreparedToRef<B> for GGLWEPrepared<D, B> {
274    fn to_ref(&self) -> GGLWEPrepared<&[u8], B> {
275        GGLWEPrepared {
276            k: self.k,
277            base2k: self.base2k,
278            dsize: self.dsize,
279            data: self.data.to_ref(),
280        }
281    }
282}