use anyhow::Result; use rand::Rng; use rand_distr::{Normal, Uniform}; use std::ops::{Add, Mul}; use arith::{Ring, Rq, TR}; use crate::glwe::{PublicKey, SecretKey, GLWE}; const ERR_SIGMA: f64 = 3.2; // l GLWEs #[derive(Clone, Debug)] pub struct GLev(pub(crate) Vec>); impl GLev { pub fn encode(m: &Rq) -> Rq { m.remodule::() } pub fn decode(p: &Rq) -> Rq { p.remodule::() } pub fn encrypt( mut rng: impl Rng, beta: u32, l: u32, pk: &PublicKey, m: &Rq, ) -> Result { let glev: Vec> = (1..l + 1) .map(|i| { GLWE::::encrypt(&mut rng, pk, &(*m * (Q / beta.pow(i as u32) as u64))) }) .collect::>>()?; Ok(Self(glev)) } pub fn encrypt_s( mut rng: impl Rng, beta: u32, l: u32, sk: &SecretKey, m: &Rq, ) -> Result { let glev: Vec> = (1..l + 1) .map(|i| { GLWE::::encrypt_s(&mut rng, sk, &(*m * (Q / beta.pow(i as u32) as u64))) }) .collect::>>()?; Ok(Self(glev)) } pub fn decrypt(&self, sk: &SecretKey, beta: u32) -> Rq { let pt = self.0[0].decrypt(sk); pt.mul_div_round(beta as u64, Q) } } #[cfg(test)] mod tests { use anyhow::Result; use rand::distributions::Uniform; use super::*; #[test] fn test_encrypt_decrypt() -> Result<()> { const Q: u64 = 2u64.pow(16) + 1; const N: usize = 128; const T: u64 = 2; // plaintext modulus const K: usize = 16; type S = GLev; let beta: u32 = 2; let l: u32 = 16; let mut rng = rand::thread_rng(); for _ in 0..200 { let (sk, pk) = GLWE::::new_key(&mut rng)?; let msg_dist = Uniform::new(0_u64, T); let m = Rq::::rand_u64(&mut rng, msg_dist)?; let p: Rq = S::encode::(&m); // plaintext let c = S::encrypt(&mut rng, beta, l, &pk, &p)?; let p_recovered = c.decrypt::(&sk, beta); let m_recovered = S::decode::(&p_recovered); assert_eq!(m, m_recovered); } Ok(()) } }