From e74b4c49efed1d8e2e817bb063b657bf56366143 Mon Sep 17 00:00:00 2001 From: kalmar Date: Fri, 21 Jul 2017 16:07:08 +0200 Subject: [PATCH] use the new definition of base-less mul! in GroupRings --- src/CheckSolution.jl | 4 ++-- src/Orbit-wise.jl | 3 ++- 2 files changed, 4 insertions(+), 3 deletions(-) diff --git a/src/CheckSolution.jl b/src/CheckSolution.jl index 14f764d..42d0786 100644 --- a/src/CheckSolution.jl +++ b/src/CheckSolution.jl @@ -27,10 +27,10 @@ function EOI{T<:Number}(Δ::GroupRingElem{T}, λ::T) return Δ*Δ - λ*Δ end -function groupring_square(vect::Vector, elt::GroupRingElem) +function groupring_square(vect::AbstractVector, elt::GroupRingElem) zzz = zeros(eltype(vect), length(elt.coeffs)) zzz[1:length(vect)] = vect - return GroupRings.mul(zzz, zzz, parent(elt).pm) + return GroupRings.mul!(similar(zzz), zzz, zzz, parent(elt).pm) end function compute_SOS(sqrt_matrix, elt) diff --git a/src/Orbit-wise.jl b/src/Orbit-wise.jl index 9a16741..a305db4 100644 --- a/src/Orbit-wise.jl +++ b/src/Orbit-wise.jl @@ -29,7 +29,8 @@ function OrbitData(name::String) splap = load(joinpath(name, "delta.jld"), "Δ"); pm = load(joinpath(name, "pm.jld"), "pm"); cnstr = PropertyT.constraints_from_pm(pm); - splap² = GroupRings.mul(splap, splap, pm); + splap² = similar(splap) + splap² = GroupRings.mul!(splap², splap, splap, pm); Uπs = load(joinpath(name, "U_pis.jld"), "Uπs"); #dimensions of the corresponding πs: