Added basis recomputing

Author: MaximusMellberg

src/SmoothNondegenerate.jl

@@ -33,7 +33,7 @@ function find_Ssmooth_model(f, M, S_target, params)
             bool = false
             return f_transformed, pseudo_inverse_mat_new
         catch e
-            (0 < verbose) && println("This f is not S smooth, changing to one that is")
+            (0 < params.verbose) && println("This f is not S smooth, changing to one that is")
             if isa(e, ArgumentError) && e.msg == "f is not smooth"
                 throw(ArgumentError("f is not smooth"))
                 bool = false  

src/ZetaFunction.jl

@@ -376,6 +376,15 @@ function zeta_function(f; S=[-1], verbose=false, givefrobmat=false, algorithm=:c
         f, pseudo_inverse_mat_new = find_Ssmooth_model(f, M, S, params)
     end
 
+    basis = compute_monomial_bases(f, R, PR, params.termorder) # basis of cohomology 
+    Basis = []
+    for i in 1:n
+        for j in basis[i]
+            push!(Basis,[j,i])
+        end
+    end
+    println(basis)
+
     #=
     BasisTLift = []
     for i in basis