debug statements

Author: jacksonwalters

src/lib.rs

@@ -156,7 +156,9 @@ fn factorize(n: i64) -> HashMap<i64, u32> {
 
 /// Fast computation of a primitive root mod p^e
 pub fn primitive_root(p: i64, e: u32) -> i64 {
+    println!("primitive_root called");
     let g = primitive_root_mod_p(p);
+    println!("g: {}", g);
     let mut g_lifted = g; // Lift it to p^e
     for _ in 1..e {
         if g_lifted.pow((p - 1) as u32) % p.pow(e) == 1 {

src/main.rs

@@ -45,4 +45,8 @@ fn main() {
     println!("Standard polynomial mult. result: {:?}", c_std);
     println!("Polynomial multiplication method using NTT: {:?}", c_fast);
 
+    let q = 12289;
+    let n = 512;
+    println!("omega = {}", ntt::omega(q*q, n));
+
 }