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

another sparsify

Browse Source
This commit is contained in:
kalmar 2017-07-05 13:19:12 +02:00
parent c0e2577684
commit e648d2c370
1 changed files with 13 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
src/OrbitDecomposition.jl
View File

@ -140,6 +140,19 @@ function Cstar_repr(x::GroupRingElem, mreps)
return res return res
end end
dens(M::SparseMatrixCSC) = length(M.nzval)/length(M)
dens(M::AbstractArray) = sum(abs.(M) .!= 0)/length(M)
function sparsify2(M::AbstractArray)
println("Density before sparsification: \t$(dens(M))")
M = deepcopy(M)
M[M .< eps(eltype(M))] .= zero(eltype(M))
M = sparse(M)
dropzeros!(M)
println("Density after sparsification: \t$(dens(M))")
return M
end
function orthSVD(M::AbstractMatrix) function orthSVD(M::AbstractMatrix)
M = full(M) M = full(M)
fact = svdfact(M) fact = svdfact(M)