kalmar/GroupsWithPropertyT
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

remove U\pi_matrices func

Browse Source
This commit is contained in:
kalmar 2017-06-21 16:44:55 +02:00
parent d2e2862ae6
commit b5c069b2a8
1 changed files with 1 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
OrbitDecomposition.jl
View File

@ -133,16 +133,6 @@ function orthSVD(M::AbstractMatrix)
return fact[:U][:,1:M_rank] return fact[:U][:,1:M_rank]
end end
function Uπ_matrices(P_matrices; orth=orthSVD)
U_p_matrices = Vector{Array{Float64,2}}()
for (i,p_mat) in enumerate(P_matrices)
U_p = orth(p_mat)
push!(U_p_matrices, U_p)
end
return U_p_matrices
end
function compute_orbit_data{T<:GroupElem}(logger, name::String, G::Group, S::Vector{T}, AutS; radius=2) function compute_orbit_data{T<:GroupElem}(logger, name::String, G::Group, S::Vector{T}, AutS; radius=2)
isdir(name) || mkdir(name) isdir(name) || mkdir(name)
@ -180,7 +170,7 @@ function compute_orbit_data{T<:GroupElem}(logger, name::String, G::Group, S::Vec
@time π_E_projections = [Cstar_repr(p, AutS_matrixreps) for p in AutS_mps] @time π_E_projections = [Cstar_repr(p, AutS_matrixreps) for p in AutS_mps]
info(logger, "Uπs...") info(logger, "Uπs...")
@time Uπs = Uπ_matrices(π_E_projections); @time Uπs = orthSVD.(π_E_projections)
multiplicities = [size(U,2) for U in Uπs]; multiplicities = [size(U,2) for U in Uπs];
info(logger, "multiplicities = $multiplicities") info(logger, "multiplicities = $multiplicities")