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

reorganize autIdproj

Browse Source
This commit is contained in:
kalmarek 2019-01-08 04:59:41 +01:00
parent 38d80e63ed
commit d4613168c3
1 changed files with 15 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
src/checksolution.jl
View File

@ -16,7 +16,19 @@ function compute_SOS(RG::GroupRing, Q::AbstractArray)
return GroupRingElem(result, RG) return GroupRingElem(result, RG)
end end
function augIdproj(Q::AbstractArray{T,2}) where {T<:Real} function augIdproj(Q::AbstractMatrix{T}) where {T<:Real}
result = zeros(size(Q))
l = size(Q, 2)
Threads.@threads for j in 1:l
col = sum(view(Q, :,j))/l
for i in 1:size(Q, 1)
result[i,j] = Q[i,j] - col
end
end
return result
end
function augIdproj(::Interval, Q::AbstractMatrix{T}) where {T<:Real}
result = zeros(Interval{T}, size(Q)) result = zeros(Interval{T}, size(Q))
l = size(Q, 2) l = size(Q, 2)
Threads.@threads for j in 1:l Threads.@threads for j in 1:l
@ -25,8 +37,6 @@ function augIdproj(Q::AbstractArray{T,2}) where {T<:Real}
result[i,j] = @interval(Q[i,j] - col) result[i,j] = @interval(Q[i,j] - col)
end end
end end
return result check = all([zero(T) in sum(view(result, :, i)) for i in 1:size(result, 2)])
end return result, check
end
end end