1
0
mirror of https://github.com/kalmarek/PropertyT.jl.git synced 2024-11-30 10:00:28 +01:00

overwrite Q to allow gc

This commit is contained in:
kalmarek 2017-10-09 18:06:29 +02:00
parent 0693455656
commit 344e11a974

View File

@ -116,20 +116,20 @@ end
function rationalize_and_project{T}(Q::AbstractArray{T}, δ::T, logger) function rationalize_and_project{T}(Q::AbstractArray{T}, δ::T, logger)
info(logger, "") info(logger, "")
info(logger, "Rationalizing with accuracy ") info(logger, "Rationalizing with accuracy ")
t = @timed Q_ = (Q, δ) t = @timed Q = (Q, δ)
info(logger, timed_msg(t)) info(logger, timed_msg(t))
info(logger, "Projecting columns of the rationalized Q to the augmentation ideal...") info(logger, "Projecting columns of the rationalized Q to the augmentation ideal...")
t = @timed Q_int = correct_to_augmentation_ideal(Q_) t = @timed Q = correct_to_augmentation_ideal(Q)
info(logger, timed_msg(t)) info(logger, timed_msg(t))
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_int, :, i)) for i in 1:size(Q_int, 2)]) check = all([0.0 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
return Q_int return Q
end end
function check_distance_to_positive_cone(Δ::GroupRingElem, λ, Q, wlen; function check_distance_to_positive_cone(Δ::GroupRingElem, λ, Q, wlen;