From 3059142db0703ca3f60eeff10572f64b8c3435de Mon Sep 17 00:00:00 2001 From: kalmar Date: Tue, 6 Jun 2017 16:23:36 +0200 Subject: [PATCH] cleaning --- OrbitDecomposition.jl | 10 +++++++--- 1 file changed, 7 insertions(+), 3 deletions(-) diff --git a/OrbitDecomposition.jl b/OrbitDecomposition.jl index 58aba99..cca93d1 100644 --- a/OrbitDecomposition.jl +++ b/OrbitDecomposition.jl @@ -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