kalmar
/
PropertyT.jl
mirror of https://github.com/kalmarek/PropertyT.jl.git
1
0
Fork 0
Code Issues Releases Wiki Activity

remove l (length) from args of

This commit is contained in:
kalmarek 2018-08-14 18:17:38 +02:00
parent cdb3c39d2a
commit 6b0cbd0114
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/CheckSolution.jl
View File

@ -25,14 +25,14 @@ end
EOI{T<:Number}(Δ::GroupRingElem{T}, λ::T) = Δ*Δ - λ*Δ EOI{T<:Number}(Δ::GroupRingElem{T}, λ::T) = Δ*Δ - λ*Δ
function groupring_square(vect::AbstractVector, l, pm) function groupring_square(vect::AbstractVector, pm)
zzz = zeros(eltype(vect), l) zzz = zeros(eltype(vect), maximum(pm))
return GroupRings.mul!(zzz, vect, vect, pm) return GroupRings.mul!(zzz, vect, vect, pm)
end end
function compute_SOS(Q::AbstractArray, pm::Array{Int,2}, l::Int) function compute_SOS(Q::AbstractArray, pm::Array{Int,2})
# result = zeros(eltype(Q), l) # result = zeros(eltype(Q), maximum(pm))
# r = similar(result) # r = similar(result)
# for i in 1:size(Q,2) # for i in 1:size(Q,2)
# print(" $i") # print(" $i")
@ -42,15 +42,15 @@ function compute_SOS(Q::AbstractArray, pm::Array{Int,2}, l::Int)
@everywhere groupring_square = PropertyT.groupring_square @everywhere groupring_square = PropertyT.groupring_square
result = @parallel (+) for i in 1:size(Q,2) result = @parallel (+) for i in 1:size(Q,2)
groupring_square(Q[:,i], l, pm) groupring_square(Q[:,i], pm)
end end
return result return result
end end
function compute_SOS(Q::AbstractArray, RG::GroupRing, l::Int) function compute_SOS(Q::AbstractArray, RG::GroupRing)
result = compute_SOS(Q, RG.pm, l) result = compute_SOS(Q, RG.pm)
return GroupRingElem(result, RG) return GroupRingElem(result, RG)
end end