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

infos

Browse Source
This commit is contained in:
kalmar 2017-01-13 18:39:02 +01:00
parent 80b2ad0987
commit ad6c0a2365
1 changed files with 2 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
property(T).jl
View File

@ -44,6 +44,7 @@ function create_product_matrix(basis::Array{Array{Float64,2},1}, limit::Int)
for i in 1:limit for i in 1:limit
x_inv = inv(basis[i]) x_inv = inv(basis[i])
# info("$i of $limit")
for j in 1:limit for j in 1:limit
w = x_inv*basis[j] w = x_inv*basis[j]
@ -137,6 +138,7 @@ function resulting_SOS{T<:Number}(sqrt_matrix::Array{T,2},
zzz = zeros(elt.coefficients) zzz = zeros(elt.coefficients)
L = size(sqrt_matrix,2) L = size(sqrt_matrix,2)
for i in 1:L for i in 1:L
info("$i of $L")
zzz[1:L] = view(sqrt_matrix, :,i) zzz[1:L] = view(sqrt_matrix, :,i)
new_base = GroupAlgebraElement(zzz, elt.product_matrix) new_base = GroupAlgebraElement(zzz, elt.product_matrix)
result += (new_base*new_base).coefficients result += (new_base*new_base).coefficients