diff --git a/SL.jl b/SL.jl
index 3855944..bc52037 100644
--- a/SL.jl
+++ b/SL.jl
@@ -44,25 +44,8 @@ function SL_generatingset(n::Int, p::Int)
     return unique(S), one(G)
 end
 
-function products{T}(U::AbstractVector{T}, V::AbstractVector{T})
-    result = Vector{T}()
-    for u in U
-        for v in V
-            push!(result, u*v)
-        end
-    end
-    return unique(result)
-end
-
 function ΔandSDPconstraints(Id, S, radius)
-    sizes = Vector{Int}()
-    S = vcat([Id], S)
-    B = S
-    push!(sizes,length(B))
-    for i in 2:2*radius
-        B = products(S, B);
-        push!(sizes, length(B))
-    end
+    B, sizes = PropertyT.generate_balls(Id, S, radius=2*radius)
     println("Generated balls of sizes $sizes")
     basis = B[1:sizes[radius]]