remove restriction on nr of variables

Author: Rafael Mohr

src/siggb/hashtable.jl

@@ -34,11 +34,6 @@ end
 
 # initialize and set fields for basis hashtable
 function initialize_basis_hash_table(::Val{N}) where N
-
-    # for now at most 64 variables
-    if N > 64
-        error("At most 64 variables currently supported.")
-    end
 
     # not necessary to create `initial_size` exponents
     exponents = Vector{Monomial{N}}(undef, init_ht_size)
@@ -50,8 +45,7 @@ function initialize_basis_hash_table(::Val{N}) where N
     size = init_ht_size
 
     # initialize fast divisibility params
-    int64bits = 64
-    ndivbits = div(64, N)
+    ndivbits = div(DivMaskSize, N)
     # division mask stores at least 1 bit
     # per each of first ndivvars variables
     ndivbits == 0 && (ndivbits += 1)

src/siggb/monomials.jl

@@ -121,6 +121,7 @@ function divmask(e::Monomial{N},
     ctr = one(DivMask)
     res = zero(DivMask)
     o = one(DivMask)
+    lb = N > DivMaskSize ? N - DivMaskSize : 1
     for i in N:-1:1
         for j in 1:ndivbits
             @inbounds if e.exps[i] >= divmap[ctr]

src/siggb/typedefs.jl

@@ -7,6 +7,7 @@ const Exp = Int16
 const MonIdx = Int32
 const MonHash = UInt32
 const DivMask = UInt64
+const DivMaskSize = 64
 # stuff for matrix
 const ColIdx = UInt32
 const Coeff = UInt32