Added a divgen() function to puritygen.jl