This commit is contained in:
kalmar 2017-06-06 16:23:36 +02:00
parent 39a6b39a2f
commit 3059142db0

View File

@ -271,7 +271,8 @@ function compute_orbit_data(name::String, N::Int)
save(joinpath(name, "delta.jld"), "Δ", Δ.coeffs)
save(joinpath(name, "pm.jld"), "pm", pm)
BN = WreathProduct(Nemo.FiniteField(2,1, "a")[1], PermutationGroup(N))
info("WreathProduct")
@time BN = WreathProduct(Nemo.FiniteField(2,1, "a")[1], PermutationGroup(N))
info("Decomposing E into orbits of B$(N)")
@time orbs = orbit_decomposition(BN, E4, E_dict)
@ -290,9 +291,12 @@ function compute_orbit_data(name::String, N::Int)
info("Uπs...")
@time Uπs = Uπ_matrices(π_E_projections);
@show multiplicities = [size(U,2) for U in Uπs];
@show dimensions = [Int(p[BN()]*Int(order(BN))) for p in BN_mps];
multiplicities = [size(U,2) for U in Uπs];
@show multiplicities
dimensions = [Int(p[BN()]*Int(order(BN))) for p in BN_mps];
@show dimensions
@assert dot(multiplicities, dimensions) == sizes[2]
save(joinpath(name, "U_pis.jld"), "Uπs", Uπs, "dims", dimensions)
return 0
end