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

Save memory: ΔandSDPconstraints takes a lot of it

Browse Source
This commit is contained in:
kalmarek 2017-08-27 19:12:21 +02:00
parent 830ef95181
commit f543d847ea
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/Orbit-wise.jl
View File

@ -215,8 +215,6 @@ function check_property_T(sett::Settings)
init_orbit_data(logger, sett, radius=sett.radius)
Δ = PropertyT.ΔandSDPconstraints(sett.name, sett.G)[1]
fnames = PropertyT.λSDPfilenames(sett.name)
if all(isfile.(fnames))
@ -235,6 +233,10 @@ function check_property_T(sett::Settings)
info(logger, "minimum(P) = $(minimum(P))")
if λ > 0
pm_fname = joinpath(sett.name, "pm.jld")
RG = GroupRing(sett.G, load(pm_fname, "pm"))
Δ_fname = joinpath(sett.name, "delta.jld")
Δ = GroupRingElem(load(Δ_fname, "Δ")[:, 1], RG)
isapprox(eigvals(P), abs.(eigvals(P)), atol=sett.tol) ||
warn("The solution matrix doesn't seem to be positive definite!")