doctest failure in ~/work/AlgebraicSolving.jl/AlgebraicSolving.jl/src/interp/cuyt_lee.jl:229-240 ```jldoctest julia> using AlgebraicSolving julia> using AlgebraicSolving.Interpolation julia> R, _ = polynomial_ring(QQ, ["a", "b"]); julia> f = t -> (t[1] + t[2])^2 / (t[1] - t[2])^2; julia> cuyt_lee(R, f) (a^2 + 2*a*b + b^2)//(a^2 - 2*a*b + b^2) ``` Subexpression: cuyt_lee(R, f) Evaluated output: ERROR: TaskFailedException nested task error: DivideError: integer division error Stacktrace: [1] divexact(a::QQFieldElem, b::QQFieldElem; check::Bool) @ Nemo ~/.julia/packages/Nemo/XmujP/src/flint/fmpq.jl:553 [2] divexact @ ~/.julia/packages/Nemo/XmujP/src/flint/fmpq.jl:552 [inlined] [3] // @ ~/.julia/packages/Nemo/XmujP/src/flint/fmpq.jl:593 [inlined] [4] evaluate(f::AbstractAlgebra.Generic.FracFieldElem{Nemo.QQPolyRingElem}, v::QQFieldElem) @ AbstractAlgebra ~/.julia/packages/AbstractAlgebra/a2EqM/src/Fraction.jl:613 [5] thiele(R::Nemo.QQPolyRing, bb::AlgebraicSolving.Interpolation.var"#35#36"{Vector{ZZRingElem}, AlgebraicSolving.Interpolation.var"#49#50"{var"#2#3", Vector{Int64}}}; retry::Int64, show_progress::Bool, desc::String, offset::Int64) @ AlgebraicSolving.Interpolation ~/work/AlgebraicSolving.jl/AlgebraicSolving.jl/src/interp/thiele.jl:136 [6] thiele @ ~/work/AlgebraicSolving.jl/AlgebraicSolving.jl/src/interp/thiele.jl:36 [inlined] [7] (::AlgebraicSolving.Interpolation.var"#populate#33"{AlgebraicSolving.Interpolation.var"#populate#27#34"{Int64, Bool, AlgebraicSolving.Interpolation.var"#49#50"{var"#2#3", Vector{Int64}}, AlgebraicSolving.Progress.ProgressBar, ReentrantLock, Vector{Int64}, Int64, Int64, Nemo.QQPolyRing, Int64}})(cur::Vector{ZZRingElem}, dim::Int64; num_threads::Int64, offset::Int64) @ AlgebraicSolving.Interpolation ~/work/AlgebraicSolving.jl/AlgebraicSolving.jl/src/interp/cuyt_lee.jl:119 [8] macro expansion @ ~/work/AlgebraicSolving.jl/AlgebraicSolving.jl/src/interp/cuyt_lee.jl:148 [inlined] [9] (::AlgebraicSolving.Interpolation.var"#39#40"{AlgebraicSolving.Interpolation.var"#41#42"{Int64, Vector{ZZRingElem}, Int64, AlgebraicSolving.Interpolation.Atomic{Int64}, AlgebraicSolving.Interpolation.Atomic{Int64}, UnitRange{Int64}}})(tid::Int64; onethread::Bool) @ AlgebraicSolving.Interpolation ./threadingconstructs.jl:279 [10] #39 @ ./threadingconstructs.jl:246 [inlined] [11] (::Base.Threads.var"#threading_run##0#threading_run##1"{AlgebraicSolving.Interpolation.var"#39#40"{AlgebraicSolving.Interpolation.var"#41#42"{Int64, Vector{ZZRingElem}, Int64, AlgebraicSolving.Interpolation.Atomic{Int64}, AlgebraicSolving.Interpolation.Atomic{Int64}, UnitRange{Int64}}}, Int64})() @ Base.Threads ./threadingconstructs.jl:178 Stacktrace: [1] threading_run(fun::AlgebraicSolving.Interpolation.var"#39#40"{AlgebraicSolving.Interpolation.var"#41#42"{Int64, Vector{ZZRingElem}, Int64, AlgebraicSolving.Interpolation.Atomic{Int64}, AlgebraicSolving.Interpolation.Atomic{Int64}, UnitRange{Int64}}}, static::Bool) @ Base.Threads ./threadingconstructs.jl:199 [2] macro expansion @ ./threadingconstructs.jl:216 [inlined] [3] (::AlgebraicSolving.Interpolation.var"#populate#33"{AlgebraicSolving.Interpolation.var"#populate#27#34"{Int64, Bool, AlgebraicSolving.Interpolation.var"#49#50"{var"#2#3", Vector{Int64}}, AlgebraicSolving.Progress.ProgressBar, ReentrantLock, Vector{Int64}, Int64, Int64, Nemo.QQPolyRing, Int64}})(cur::Vector{ZZRingElem}, dim::Int64; num_threads::Int64, offset::Int64) @ AlgebraicSolving.Interpolation ~/work/AlgebraicSolving.jl/AlgebraicSolving.jl/src/interp/cuyt_lee.jl:147 [4] cuyt_lee_shifted(R::Nemo.QQMPolyRing, bb::AlgebraicSolving.Interpolation.var"#49#50"{var"#2#3", Vector{Int64}}; retry::Int64, nr_thrds::Int64, show_progress::Bool, desc::String) @ AlgebraicSolving.Interpolation ~/work/AlgebraicSolving.jl/AlgebraicSolving.jl/src/interp/cuyt_lee.jl:162 [5] cuyt_lee_shifted @ ~/work/AlgebraicSolving.jl/AlgebraicSolving.jl/src/i