kalmar/PropertyT.jl
1
0
Fork 0
mirror of https://github.com/kalmarek/PropertyT.jl.git synced 2024-10-15 08:05:35 +02:00
Code Issues Releases Wiki Activity

threaded generation of matrix reps

Browse Source
This commit is contained in:
kalmarek 2017-08-27 18:32:50 +02:00
parent 0cc0ad94a6
commit 6061f731b3
1 changed files with 6 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/OrbitDecomposition.jl
View File

@ -108,14 +108,16 @@ function matrix_reps{T<:GroupElem}(G::Group, S::Vector{T}, AutS::Group, radius::
Id = (isa(G, Nemo.Ring) ? one(G) : G())
E2, _ = Groups.generate_balls(S, Id, radius=radius)
Edict = GroupRings.reverse_dict(E2)
elts = collect(elements(AutS))
mreps = Vector{SparseMatrixCSC{Int, Int}}(length(elts))
Threads.@threads for i in 1:length(elts)
elts = collect(elements(AutS))
l = length(elts)
mreps = Vector{SparseMatrixCSC{Int, Int}}(l)
Threads.@threads for i in 1:l
mreps[i] = PropertyT.matrix_repr(elts[i], E2, Edict)
end
mreps_dict = Dict(elts[i]=>mreps[i] for i in 1:length(elts))
mreps_dict = Dict(elts[i]=>mreps[i] for i in 1:l)
return mreps_dict
end