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

update autIdproj

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

13
src/CheckSolution.jl
View File

@ -64,23 +64,24 @@ function distances_to_cone(elt::GroupRingElem, wlen::Int)
end end
function augIdproj{T, I<:AbstractInterval}(S::Type{I}, Q::AbstractArray{T,2}) function augIdproj{T, I<:AbstractInterval}(S::Type{I}, Q::AbstractArray{T,2})
function augIdproj(Q::AbstractArray{T,2}) where {T<:Real}
R = zeros(Interval{T}, size(Q))
l = size(Q, 2) l = size(Q, 2)
R = zeros(S, (l,l))
Threads.@threads for j in 1:l Threads.@threads for j in 1:l
col = sum(view(Q, :,j))/l col = sum(view(Q, :,j))/l
for i in 1:l for i in 1:size(Q, 1)
R[i,j] = Q[i,j] - col ± eps(0.0) R[i,j] = @interval(Q[i,j] - col)
end end
end end
return R return R
end end
function augIdproj{T}(Q::AbstractArray{T,2}, logger) function augIdproj(Q::AbstractArray{T,2}, logger) where {T<:Real}
info(logger, "Projecting columns of Q to the augmentation ideal...") info(logger, "Projecting columns of Q to the augmentation ideal...")
@logtime logger Q = augIdproj(Interval{T}, Q) @logtime logger Q = augIdproj(Q)
info(logger, "Checking that sum of every column contains 0.0... ") info(logger, "Checking that sum of every column contains 0.0... ")
check = all([0.0 in sum(view(Q, :, i)) for i in 1:size(Q, 2)]) check = all([zero(T) in sum(view(Q, :, i)) for i in 1:size(Q, 2)])
info(logger, (check? "They do." : "FAILED!")) info(logger, (check? "They do." : "FAILED!"))
@assert check @assert check