kalmar/PropertyT.jl
1
0
Fork 0
mirror of https://github.com/kalmarek/PropertyT.jl.git synced 2024-12-26 02:30:29 +01:00
Code Issues Releases Wiki Activity

add verbose=true to OrbitData constructor

Browse Source
This commit is contained in:
kalmarek 2018-09-16 18:03:25 +02:00
parent 3df3967fad
commit 9d6a6fa0f8
1 changed files with 12 additions and 33 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

45
src/orbitdata.jl
View File

@ -11,47 +11,26 @@ struct OrbitData{T<:AbstractArray{Float64, 2}, GEl<:GroupElem, P<:perm}
dims::Vector{Int} dims::Vector{Int}
end end
function OrbitData(RG::GroupRing, autS::Group) function OrbitData(RG::GroupRing, autS::Group, verbose=true)
orbs = orbit_decomposition(autS, RG.basis, RG.basis_dict) verbose && info("Decomposing basis of RG into orbits of $(autS)")
@assert sum(length(o) for o in orbs) == length(RG.basis)
autS_mps = Projections.rankOne_projections(GroupRing(autS))
preps = perm_reps(autS, RG.basis[1:size(RG.pm,1)], RG.basis_dict)
mreps = matrix_reps(preps)
Uπs = [orthSVD(matrix_repr(p, mreps)) for p in autS_mps]
multiplicities = size.(Uπs,2)
dimensions = [Int(p[autS()]*Int(order(autS))) for p in autS_mps]
@assert dot(multiplicities, dimensions) == size(RG.pm,1)
nzros = [i for i in 1:length(Uπs) if size(Uπs[i],2) !=0]
return OrbitData(orbs, preps, Uπs[nzros], dims[nzros])
end
function compute_OrbitData(RG::GroupRing, autS::Group)
info("Decomposing E into orbits of $(autS)")
@time orbs = orbit_decomposition(autS, RG.basis, RG.basis_dict) @time orbs = orbit_decomposition(autS, RG.basis, RG.basis_dict)
@assert sum(length(o) for o in orbs) == length(RG.basis) @assert sum(length(o) for o in orbs) == length(RG.basis)
info("E consists of $(length(orbs)) orbits!") verbose && info("The action has $(length(orbs)) orbits")
info("Action matrices") verbose && info("Projections in the Group Ring of AutS")
@time autS_mps = Projections.rankOne_projections(GroupRing(autS))
verbose && info("AutS-action matrix representatives")
@time preps = perm_reps(autS, RG.basis[1:size(RG.pm,1)], RG.basis_dict) @time preps = perm_reps(autS, RG.basis[1:size(RG.pm,1)], RG.basis_dict)
mreps = matrix_reps(preps) @time mreps = matrix_reps(preps)
info("Projections") verbose && info("Projection matrices Uπs")
@time autS_mps = Projections.rankOne_projections(GroupRing(autS));
info("Uπs...")
@time Uπs = [orthSVD(matrix_repr(p, mreps)) for p in autS_mps] @time Uπs = [orthSVD(matrix_repr(p, mreps)) for p in autS_mps]
multiplicities = size.(Uπs,2) multiplicities = size.(Uπs,2)
info("multiplicities = $multiplicities") verbose && info("multiplicities = $multiplicities")
dimensions = [Int(p[autS()]*Int(order(autS))) for p in autS_mps]; dimensions = [Int(p[autS()]*Int(order(autS))) for p in autS_mps]
info("dimensions = $dimensions") verbose && info("dimensions = $dimensions")
@assert dot(multiplicities, dimensions) == size(RG.pm,1) @assert dot(multiplicities, dimensions) == size(RG.pm,1)
return OrbitData(orbs, preps, Uπs, dimensions) return OrbitData(orbs, preps, Uπs, dimensions)