feat(ntt): impl product plan

co-authored-by: sarah el kazdadi <[email protected]>

Cargo.toml

@@ -13,6 +13,7 @@ rust-version = "1.67"
 
 [dependencies]
 aligned-vec = { version = "0.5", default-features = false }
+bytemuck = "1"
 pulp = { version = "0.18.22", default-features = false }
 
 [features]

src/lib.rs

@@ -109,6 +109,8 @@ pub mod native_binary32;
 /// polynomial.
 pub mod native_binary64;
 
+pub mod product;
+
 // Fn arguments are (simd, z0, z1, w, w_shoup, p, neg_p, two_p)
 trait Butterfly<S: Copy, V: Copy>: Copy + Fn(S, V, V, V, V, V, V, V) -> (V, V) {}
 impl<F: Copy + Fn(S, V, V, V, V, V, V, V) -> (V, V), S: Copy, V: Copy> Butterfly<S, V> for F {}

src/prime32.rs

@@ -685,6 +685,10 @@ impl Plan {
         }
     }
 
+    pub(crate) fn p_div(&self) -> Div32 {
+        self.p_div
+    }
+
     /// Returns the polynomial size of the negacyclic NTT plan.
     #[inline]
     pub fn ntt_size(&self) -> usize {

src/prime64.rs

@@ -770,6 +770,10 @@ impl Plan {
         }
     }
 
+    pub(crate) fn p_div(&self) -> Div64 {
+        self.p_div
+    }
+
     /// Returns the polynomial size of the negacyclic NTT plan.
     #[inline]
     pub fn ntt_size(&self) -> usize {