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

threaded matrix_reps

Browse Source
This commit is contained in:
kalmar 2017-08-04 18:25:10 +02:00
parent 666fac1775
commit b97a0e30dd
1 changed files with 9 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
src/OrbitDecomposition.jl
View File

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