From 677ae32ee2233e1562098e806f530cef08d04ad1 Mon Sep 17 00:00:00 2001 From: kalmarek Date: Wed, 15 Aug 2018 19:13:26 +0200 Subject: [PATCH] rankOne_projections take a GroupRing now --- src/OrbitDecomposition.jl | 2 +- src/Projections.jl | 14 ++++++-------- 2 files changed, 7 insertions(+), 9 deletions(-) diff --git a/src/OrbitDecomposition.jl b/src/OrbitDecomposition.jl index 41d766e..b0c5ad4 100644 --- a/src/OrbitDecomposition.jl +++ b/src/OrbitDecomposition.jl @@ -162,7 +162,7 @@ function compute_orbit_data{T<:GroupElem}(logger, name::String, S::Vector{T}, au reps = matrix_reps(reps) info(logger, "Projections") - @logtime logger autS_mps = Projections.rankOne_projections(autS); + @logtime logger autS_mps = Projections.rankOne_projections(GroupRing(autS)); @logtime logger π_E_projections = [Cstar_repr(p, reps) for p in autS_mps] diff --git a/src/Projections.jl b/src/Projections.jl index d2b5064..cfdba18 100644 --- a/src/Projections.jl +++ b/src/Projections.jl @@ -142,11 +142,9 @@ function rankOne_projection(chi::PermCharacter, idems::Vector{T}) where {T<:Grou throw("Couldn't find rank-one projection for $chi") end -function rankOne_projections(G::Generic.PermGroup, T::Type=Rational{Int}) - if G.n == 1 - return [one(GroupRing(G), T)] - else - RG = GroupRing(G, fastm=false) +function rankOne_projections(RG::GroupRing{G}, T::Type=Rational{Int}) where G<:Generic.PermGroup + if RG.group.n == 1 + return [GroupRingElem([one(T)], RG)] end RGidems = idempotents(RG, T) @@ -165,14 +163,14 @@ function orbit_selector(n::Integer, k::Integer, return Projections.DirectProdCharacter(ntuple(i -> (i <= k ? chi : psi), n)) end -function rankOne_projections(Bn::WreathProduct, T::Type=Rational{Int}) +function rankOne_projections(RBn::GroupRing{G}, T::Type=Rational{Int}) where {G<:WreathProduct} + Bn = RBn.group N = Bn.P.n # projections as elements of the group rings RSₙ - Sn_rankOnePr = [rankOne_projections(PermutationGroup(Int8(i))) for i in 1:N] + Sn_rankOnePr = [rankOne_projections(GroupRing(PermutationGroup(i))) for i in typeof(N)(1):N] # embedding into group ring of BN - RBn = GroupRing(Bn) RN = GroupRing(Bn.N) sign, id = collect(characters(Bn.N.group))