diff --git a/SL.jl b/SL.jl
index 76aabc3..3855944 100644
--- a/SL.jl
+++ b/SL.jl
@@ -55,21 +55,18 @@ function products{T}(U::AbstractVector{T}, V::AbstractVector{T})
 end
 
 function ΔandSDPconstraints(Id, S, radius)
-    radius *=2
-
     sizes = Vector{Int}()
     S = vcat([Id], S)
     B = S
     push!(sizes,length(B))
-    for i in 2:radius
+    for i in 2:2*radius
         B = products(S, B);
         push!(sizes, length(B))
     end
     println("Generated balls of sizes $sizes")
-    k = div(radius,2)
-    basis = B[1:sizes[k]]
+    basis = B[1:sizes[radius]]
 
-    product_matrix = PropertyT.create_product_matrix(B, sizes[k]);
+    product_matrix = PropertyT.create_product_matrix(B, sizes[radius]);
     sdp_constraints = PropertyT.constraints_from_pm(product_matrix, length(B))
     L_coeff = PropertyT.splaplacian_coeff(S, basis, length(B));
     Δ = GroupAlgebraElement(L_coeff, product_matrix)