replace rug by num_bigint

Author: mwien

.gitignore

@@ -1 +1,2 @@
 target/
+cliquepicking_python/.env/

cliquepicking_python/.env/lib64

@@ -12,4 +12,3 @@ crate-type = ["cdylib"]
 cliquepicking_rs = { path = "../cliquepicking_rs" }
 num-bigint = "0.4.6"
 pyo3 = { version = "0.21.1", features = ["num-bigint"] }
-rug = "1.24.1"

cliquepicking_python/.env/pyvenv.cfg

@@ -6,7 +6,6 @@ use cliquepicking_rs::sample::sample_cpdag;
 use cliquepicking_rs::sample::sample_cpdag_orders;
 
 use num_bigint::BigUint;
-use rug::integer::Order;
 
 /// A Python module for counting and sampling Markov equivalent DAGs.
 #[pymodule]
@@ -23,8 +22,7 @@ fn mec_size(cpdag: Vec<(usize, usize)>) -> PyResult<BigUint> {
     let mx = max_element(&cpdag);
     let g = PartiallyDirectedGraph::from_edge_list(cpdag, mx + 1);
     let res = count_cpdag(&g);
-    let digits = res.to_digits::<u32>(Order::Lsf);
-    Ok(BigUint::new(digits))
+    Ok(res)
 }
 
 /// Sample k DAGs uniformly from the Markov equivalence class represented by CPDAG cpdag.

cliquepicking_python/Cargo.lock

@@ -7,4 +7,5 @@ edition = "2021"
 
 [dependencies]
 rand = "0.8.5"
-rug = "1.24.1"
+num-bigint = { version = "0.4.6", features = ["rand"] }
+num-traits = "0.2.19"

cliquepicking_python/Cargo.toml

@@ -1,12 +1,13 @@
-use rug::Integer;
+use num_bigint::BigUint;
+use num_traits::One;
 
 use crate::memoization::Memoization;
 
-pub fn factorial(n: usize, memoization: &mut [Integer]) -> Integer {
-    if memoization[n] != 0 {
+pub fn factorial(n: usize, memoization: &mut [BigUint]) -> BigUint {
+    if memoization[n] != BigUint::ZERO {
         return memoization[n].clone();
     }
-    let mut result = Integer::from(1);
+    let mut result = BigUint::one();
     for i in 1..n + 1 {
         result *= i;
     }
@@ -15,7 +16,7 @@ pub fn factorial(n: usize, memoization: &mut [Integer]) -> Integer {
 }
 
 // TODO: maybe implement without recursion?
-pub fn rho(x: &[usize], memoization: &mut Memoization) -> Integer {
+pub fn rho(x: &[usize], memoization: &mut Memoization) -> BigUint {
     let x_vec = x.to_vec();
     if let Some(res) = memoization.rho.get(&x_vec) {
         return res.clone();