slight tweaks

This commit is contained in:
kalmar 2017-04-08 13:48:00 +02:00 committed by kalmarek
parent da236c32ce
commit 57f92113ac

9
SL.jl
View File

@ -55,21 +55,18 @@ function products{T}(U::AbstractVector{T}, V::AbstractVector{T})
end end
function ΔandSDPconstraints(Id, S, radius) function ΔandSDPconstraints(Id, S, radius)
radius *=2
sizes = Vector{Int}() sizes = Vector{Int}()
S = vcat([Id], S) S = vcat([Id], S)
B = S B = S
push!(sizes,length(B)) push!(sizes,length(B))
for i in 2:radius for i in 2:2*radius
B = products(S, B); B = products(S, B);
push!(sizes, length(B)) push!(sizes, length(B))
end end
println("Generated balls of sizes $sizes") println("Generated balls of sizes $sizes")
k = div(radius,2) basis = B[1:sizes[radius]]
basis = B[1:sizes[k]]
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)) sdp_constraints = PropertyT.constraints_from_pm(product_matrix, length(B))
L_coeff = PropertyT.splaplacian_coeff(S, basis, length(B)); L_coeff = PropertyT.splaplacian_coeff(S, basis, length(B));
Δ = GroupAlgebraElement(L_coeff, product_matrix) Δ = GroupAlgebraElement(L_coeff, product_matrix)