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

allow to have 0s in completed pm

Browse Source
This commit is contained in:
kalmarek 2017-09-10 21:42:45 +02:00
parent e223763b75
commit 61bed0f3d0
1 changed files with 11 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
src/GroupRings.jl
View File

@ -556,12 +556,21 @@ function complete!(RG::GroupRing)
RG.basis = [elements(RG.group)...]
end
fastm!(RG, fill=true)
fastm!(RG, fill=false)
warning = false
for linidx in find(RG.pm .== 0)
i,j = ind2sub(size(RG.pm), linidx)
RG.pm[i,j] = RG.basis_dict[RG.basis[i]*RG.basis[j]]
g = RG.basis[i]*RG.basis[j]
if haskey(RG.basis_dict, g)
RG.pm[i,j] = RG.basis_dict[g]
else
if !warning
warning = true
end
end
end
warning && warn("Some products were not supported on basis")
return RG
end