1
0
mirror of https://github.com/kalmarek/PropertyT.jl.git synced 2024-11-19 07:20:28 +01:00
This commit is contained in:
kalmarek 2017-11-08 15:40:38 +01:00
parent 38ae1e0656
commit a35122f8bb
2 changed files with 5 additions and 5 deletions

View File

@ -195,6 +195,7 @@ function compute_orbit_data{T<:GroupElem}(logger, name::String, G::Nemo.Group, S
info(logger, "Decomposing E into orbits of $(autS)")
@logtime logger orbs = orbit_decomposition(autS, E_2R, E_rdict)
@assert sum(length(o) for o in orbs) == length(E_2R)
info(logger, "E consists of $(length(orbs)) orbits!")
save(joinpath(name, "orbits.jld"), "orbits", orbs)
info(logger, "Action matrices")

View File

@ -90,19 +90,19 @@ function idempotents(RG::GroupRing{PermGroup}, T::Type=Rational{Int})
return unique(idems)
end
function rankOne_projection{S}(chi::PropertyT.PermCharacter, idems::Vector{GroupRingElem{S}})
function rankOne_projection{T}(chi::PropertyT.PermCharacter, idems::Vector{GroupRingElem{T}})
RG = parent(first(idems))
ids = [[one(RG, S)]; idems]
ids = [[one(RG, T)]; idems]
for (i,j,k) in Base.product(ids, ids, ids)
if chi(i) == zero(S) || chi(j) == zero(S) || chi(k) == zero(S)
if chi(i) == zero(T) || chi(j) == zero(T) || chi(k) == zero(T)
continue
end
elt = i*j*k
elt^2 == elt || continue
if chi(elt) == one(S)
if chi(elt) == one(T)
return elt
# return (i,j,k)
end
@ -134,7 +134,6 @@ function rankOne_projections(G::PermutationGroup, T::Type=Rational{Int})
return min_projs
end
function rankOne_projections(BN::WreathProduct, T::Type=Rational{Int})
N = BN.P.n