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

slighlty more performant reconstruct_sol

Browse Source
This commit is contained in:
kalmarek 2018-08-20 03:55:40 +02:00
parent 5acc00a4d0
commit 66ac37a6a7
1 changed files with 6 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
src/OrbitDecomposition.jl
View File

@ -96,8 +96,7 @@ function perm_reps(S::Vector, autS::Group, radius::Int)
return perm_reps(autS, E) return perm_reps(autS, E)
end end
function reconstruct_sol{T<:GroupElem, S<:perm}(preps::Dict{T, S}, function reconstruct_sol(preps::Dict{T, S}, Us::Vector, Ps::Vector, dims::Vector) where {T<:GroupElem, S<:perm}
Us::Vector, Ps::Vector, dims::Vector)
l = length(Us) l = length(Us)
transfP = [dims[π].*Us[π]*Ps[π]*Us[π]' for π in 1:l] transfP = [dims[π].*Us[π]*Ps[π]*Us[π]' for π in 1:l]
@ -111,7 +110,11 @@ function reconstruct_sol{T<:GroupElem, S<:perm}(preps::Dict{T, S},
end end
recP = 1/length(perms) .* sum(tmp) recP = 1/length(perms) .* sum(tmp)
recP[abs.(recP) .< eps(eltype(recP))] = zero(eltype(recP)) for i in eachindex(recP)
if abs(recP[i]) .< eps(eltype(recP))*100
recP[i] = zero(eltype(recP))
end
end
return recP return recP
end end