make orbits into vectors

OrbitDecomposition.jl

@@ -62,18 +62,18 @@ function orbit_decomposition(G::Nemo.Group, E::Vector, rdict=GroupRings.reverse_
     elts = collect(elements(G))
 
     tovisit = trues(E);
-    orbits = Vector{Set{Int}}()
+    orbits = Vector{Vector{Int}}()
 
     for i in 1:endof(E)
         if tovisit[i]
-            orbit = Set{Int}()
+            orbit = Vector{Int}()
             a = E[i]
             for g in elts
                 idx = rdict[g(a)]
                 tovisit[idx] = false
                 push!(orbit,idx)
             end
-            push!(orbits, orbit)
+            push!(orbits, unique(orbit))
         end
     end
     return orbits