poulpy/base2k/src/module.rs
base2k/src/module.rsRaw
1use crate::ffi::module::{delete_module_info, module_info_t, new_module_info, MODULE};
2use crate::GALOISGENERATOR;
3
4#[derive(Copy, Clone)]
5#[repr(u8)]
6pub enum BACKEND {
7 FFT64,
8 NTT120,
9}
10
11pub struct Module {
12 pub ptr: *mut MODULE,
13 pub n: usize,
14 pub backend: BACKEND,
15}
16
17impl Module {
18 // Instantiates a new module.
19 pub fn new(n: usize, module_type: BACKEND) -> Self {
20 unsafe {
21 let module_type_u32: u32;
22 match module_type {
23 BACKEND::FFT64 => module_type_u32 = 0,
24 BACKEND::NTT120 => module_type_u32 = 1,
25 }
26 let m: *mut module_info_t = new_module_info(n as u64, module_type_u32);
27 if m.is_null() {
28 panic!("Failed to create module.");
29 }
30 Self {
31 ptr: m,
32 n: n,
33 backend: module_type,
34 }
35 }
36 }
37
38 pub fn backend(&self) -> BACKEND {
39 self.backend
40 }
41
42 pub fn n(&self) -> usize {
43 self.n
44 }
45
46 pub fn log_n(&self) -> usize {
47 (usize::BITS - (self.n() - 1).leading_zeros()) as _
48 }
49
50 pub fn cyclotomic_order(&self) -> u64 {
51 (self.n() << 1) as _
52 }
53
54 // GALOISGENERATOR^|gen| * sign(gen)
55 pub fn galois_element(&self, gen: i64) -> i64 {
56 if gen == 0 {
57 return 1;
58 }
59
60 let mut gal_el: u64 = 1;
61 let mut gen_1_pow: u64 = GALOISGENERATOR;
62 let mut e: usize = gen.abs() as usize;
63 while e > 0 {
64 if e & 1 == 1 {
65 gal_el = gal_el.wrapping_mul(gen_1_pow);
66 }
67
68 gen_1_pow = gen_1_pow.wrapping_mul(gen_1_pow);
69 e >>= 1;
70 }
71
72 gal_el &= self.cyclotomic_order() - 1;
73
74 (gal_el as i64) * gen.signum()
75 }
76
77 pub fn free(self) {
78 unsafe { delete_module_info(self.ptr) }
79 drop(self);
80 }
81}