From 8f88035bb72ac21fb1c2bb0908a906c109955713 Mon Sep 17 00:00:00 2001 From: kalmar Date: Fri, 4 Aug 2017 20:50:06 +0200 Subject: [PATCH] in-place projection to augmentation ideal --- src/CheckSolution.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/CheckSolution.jl b/src/CheckSolution.jl index f261788..80b6d73 100644 --- a/src/CheckSolution.jl +++ b/src/CheckSolution.jl @@ -50,9 +50,9 @@ function compute_SOS(sqrt_matrix, elt::GroupRingElem) return GroupRingElem(result, parent(elt)) end -function correct_to_augmentation_ideal{T<:Rational}(sqrt_matrix::Array{T,2}) +function correct_to_augmentation_ideal!{T<:Rational}(sqrt_matrix::Array{T,2}) l = size(sqrt_matrix, 2) - sqrt_corrected = copy(sqrt_matrix) + sqrt_corrected = sqrt_matrix Threads.@threads for j in 1:l col = sum(view(sqrt_matrix, :,j))//l for i in 1:l