diff --git a/Project.toml b/Project.toml
index e0d8a35c..87bb45fd 100644
--- a/Project.toml
+++ b/Project.toml
@@ -48,7 +48,7 @@ Pkg = "1.10, 1.11"
 PrecompileTools = "1.2"
 Primes = "0.5"
 Random = "1.10, 1.11"
-RationalFunctionFields = "0.2.2"
+RationalFunctionFields = "0.2.3"
 SpecialFunctions = "2"
 SymbolicUtils = "4"
 Symbolics = "7"
diff --git a/src/parametrizations.jl b/src/parametrizations.jl
index ccdd8286..5d8a679a 100644
--- a/src/parametrizations.jl
+++ b/src/parametrizations.jl
@@ -158,6 +158,7 @@ The function accepts the following optional arguments.
 
 - `seed`: A float in the range from 0 to 1, random seed (default is `seed = 42`). 
 - `prob_threshold`: The probability of correctness (default is `prob_threshold = 0.99`).
+- `loglevel`: the level of logs to be displayed (default is `Logging.Info`).
 
 ## Example
 
@@ -208,7 +209,7 @@ end
 
 function _reparametrize_global(ode::ODE{P}; prob_threshold = 0.99, seed = 42) where {P}
     Random.seed!(seed)
-    id_funcs = find_identifiable_functions(
+    id_funcs = _find_identifiable_functions(
         ode,
         with_states = true,
         simplify = :strong,
diff --git a/test/raw_generators.jl b/test/raw_generators.jl
index a06b17d6..26edb6e4 100644
--- a/test/raw_generators.jl
+++ b/test/raw_generators.jl
@@ -47,14 +47,14 @@ import Groebner
     gb = Groebner.groebner(eqs, ordering = Groebner.DegRevLex())
     # GB is linear
     @test length(gb) == length(gens(parent(eqs[1])))
-    expected = 9202476
+    expected = 9220151
     str = join(map(string, eqs), ",")
     @info "" length(str)
     @test abs(length(str) - expected) / expected * 100 < 5
 
     # Part 2: over Q
     eqs = fractionfree_generators_raw(rff.mqs)[1]
-    expected = 21486079
+    expected = 14652203
     str = join(map(string, eqs), ",")
     @info "" length(str)
     @test abs(length(str) - expected) / expected * 100 < 5