kalmar/PropertyT.jl
1
0
Fork 0
mirror of https://github.com/kalmarek/PropertyT.jl.git synced 2024-11-23 08:15:29 +01:00
Code Issues Releases Wiki Activity

use the new definition of base-less mul! in GroupRings

Browse Source
This commit is contained in:
kalmar 2017-07-21 16:07:08 +02:00
parent b8a40c42c8
commit e74b4c49ef
2 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/CheckSolution.jl
View File

@ -27,10 +27,10 @@ function EOI{T<:Number}(Δ::GroupRingElem{T}, λ::T)
return Δ*Δ - λ*Δ return Δ*Δ - λ*Δ
end end
function groupring_square(vect::Vector, elt::GroupRingElem) function groupring_square(vect::AbstractVector, elt::GroupRingElem)
zzz = zeros(eltype(vect), length(elt.coeffs)) zzz = zeros(eltype(vect), length(elt.coeffs))
zzz[1:length(vect)] = vect zzz[1:length(vect)] = vect
return GroupRings.mul(zzz, zzz, parent(elt).pm) return GroupRings.mul!(similar(zzz), zzz, zzz, parent(elt).pm)
end end
function compute_SOS(sqrt_matrix, elt) function compute_SOS(sqrt_matrix, elt)

3
src/Orbit-wise.jl
View File

@ -29,7 +29,8 @@ function OrbitData(name::String)
splap = load(joinpath(name, "delta.jld"), "Δ"); splap = load(joinpath(name, "delta.jld"), "Δ");
pm = load(joinpath(name, "pm.jld"), "pm"); pm = load(joinpath(name, "pm.jld"), "pm");
cnstr = PropertyT.constraints_from_pm(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"); Uπs = load(joinpath(name, "U_pis.jld"), "Uπs");
#dimensions of the corresponding πs: #dimensions of the corresponding πs: