0xalivecow/kauma
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

refactor: clen edf code

Browse source
This commit is contained in:
Alivecow 2024-12-03 17:00:23 +01:00
parent 4c9adb9fdc
commit b632e0c803
1 changed files with 0 additions and 16 deletions
Download patch file Download diff file Expand all files Collapse all files

16
src/utils/edf.rs
View file

@ -4,34 +4,21 @@ use rand::Rng;
use super::poly::{gcd, Polynomial}; use super::poly::{gcd, Polynomial};
pub fn edf(f: Polynomial, d: u32) -> Vec<Polynomial> { pub fn edf(f: Polynomial, d: u32) -> Vec<Polynomial> {
eprintln!("edf started: {:?}", f.clone().to_c_array());
let q = BigUint::pow(&BigUint::from_u8(2).unwrap(), 128); let q = BigUint::pow(&BigUint::from_u8(2).unwrap(), 128);
let n: u32 = (f.degree() as u32) / (d); let n: u32 = (f.degree() as u32) / (d);
let mut z: Vec<Polynomial> = vec![f.clone()]; let mut z: Vec<Polynomial> = vec![f.clone()];
let one_cmp = Polynomial::one(); let one_cmp = Polynomial::one();
while (z.len() as u32) < n { while (z.len() as u32) < n {
eprintln!("z len {}", z.len());
eprintln!("n len {}", n);
let h = Polynomial::rand(&rand::thread_rng().gen_range(1..=f.degree())); let h = Polynomial::rand(&rand::thread_rng().gen_range(1..=f.degree()));
eprintln!("h: {:02X?}", h);
let exponent = (q.pow(d) - BigUint::one()) / BigUint::from_u8(3).unwrap(); let exponent = (q.pow(d) - BigUint::one()) / BigUint::from_u8(3).unwrap();
let g = h.bpow_mod(exponent, &f) + Polynomial::one(); let g = h.bpow_mod(exponent, &f) + Polynomial::one();
eprintln!("g before for {:0X?}", g);
eprintln!("z before for {:0X?}", z);
for i in (0..z.len()).rev() { for i in (0..z.len()).rev() {
if z[i].degree() as u32 > d { if z[i].degree() as u32 > d {
eprintln!("Inside if");
let j = gcd(&z[i], &g); let j = gcd(&z[i], &g);
eprintln!("j != one_cmp {:?}", j != one_cmp);
eprintln!("j != z[i] {:?}", j != z[i]);
eprintln!("Inside if");
if j != one_cmp && j != z[i] { if j != one_cmp && j != z[i] {
let intemediate = z[i].div(&j).0; let intemediate = z[i].div(&j).0;
z.remove(i); z.remove(i);
@ -40,11 +27,8 @@ pub fn edf(f: Polynomial, d: u32) -> Vec<Polynomial> {
} }
} }
} }
eprintln!("z after for {:0X?}", z);
} }
eprintln!("edf finished");
z z
} }