1
0
mirror of https://github.com/kalmarek/PropertyT.jl.git synced 2024-12-25 18:25:30 +01:00

make the laplacian generation code more generic

This commit is contained in:
kalmar 2017-02-11 13:34:28 +01:00
parent 7f707e60a0
commit f9f96acb9a

View File

@ -37,24 +37,21 @@ function constraints_from_pm(pm, total_length=maximum(pm))
end
return constraints
end
function Laplacian_sparse(S::Array{Array{Float64,2},1},
basis::Array{Array{Float64,2},1})
function splaplacian_coeff(S, basis, n=length(basis))
squares = unique(vcat([basis[1]], S, products(S,S)))
result = spzeros(length(basis))
result = spzeros(n)
result[1] = length(S)
for s in S
ind = find(basis, s)
ind = findfirst(basis, s)
result[ind] += -1
end
return result
end
function Laplacian(S::Array{Array{Float64,2},1},
basis:: Array{Array{Float64,2},1})
return full(Laplacian_sparse(S,basis))
function laplacian_coeff(S, basis)
return full(splaplacian_coeff(S,basis))
end