Commit Graph

578 Commits

Author SHA1 Message Date
Marek Kaluba a05a0798ec
add tests for levels and their symmetrization 2023-05-22 22:46:22 +02:00
Marek Kaluba 66f10612de
streamline constraints generation
and make it a bit less dependent on the ordering of elements in basis
2023-05-22 22:40:32 +02:00
Marek Kaluba d1a58d3890
remove unnecessary @time from certify
It's all blazing fast now
2023-05-22 22:39:20 +02:00
Marek Kaluba cd7901b455
add G₂ to classify_root_system 2023-05-10 15:20:07 +02:00
Marek Kaluba d40a0fe117
fix error with printing of roots 2023-05-09 01:02:39 +02:00
Marek Kaluba 58f0ccb141
use IntervalMatrices
IntervalMatrices use Rump algorithm to matrix multiplication
This brings time to Qint'*Qint down to ~40s which is
5-8 × slower than Q'*Q (for size n=2^13).

The naive version is ~100 × slower than Q'*Q even for n = 2^10.
2023-04-06 13:17:52 +02:00
Marek Kaluba f0986982ce
reorganize Roots module 2023-04-06 11:39:54 +02:00
Marek Kaluba 5b4a7f6804
reshuffle sos_sdps for clarity 2023-04-04 23:50:48 +02:00
Marek Kaluba 150b5c2cba
skip the identity constraint if augmented 2023-04-04 23:48:59 +02:00
Marek Kaluba 3f2be20152
make nzpairs(::ConstraintMatrix) type stable 2023-04-04 23:15:40 +02:00
Marek Kaluba 132802feeb
use sparse matrices for invariant constraint 2023-04-04 23:14:56 +02:00
Marek Kaluba 1a43a1b1be
in reconstruct: average the sum, not sum the averages! 2023-04-04 19:58:51 +02:00
Marek Kaluba fb3b51fd6e
fix sos_problem_dual 2023-03-20 01:37:19 +01:00
Marek Kaluba bacd170504
preserve trace when diagonalizing M_orb 2023-03-20 01:37:18 +01:00
Marek Kaluba a5f5a4ea35
lots of re-formatting 2023-03-20 01:37:18 +01:00
Marek Kaluba a1de0ecc85
improve __fast_recursive_dot a bit 2023-03-20 01:37:17 +01:00
Marek Kaluba b5fa1ac0ef
use Cartan matrix to classify root-subsystems 2023-03-20 01:37:16 +01:00
Marek Kaluba 1fb324b49a
update constraints to StarAlgebras-0.2 2023-03-20 01:36:40 +01:00
Marek Kaluba 4e43811ea3
fix a nasty bug with negatives in ConstraintMatrix 2023-03-20 01:36:39 +01:00
Marek Kaluba 2cc9444667
small fixes 2022-11-17 11:22:17 +01:00
Marek Kaluba d1d46d13ef
use spzeros constructor compatible with julia-1.6 2022-11-16 00:14:06 +01:00
Marek Kaluba f053bffefe
rewrite reconstruct! with a better architecture 2022-11-15 19:36:39 +01:00
Marek Kaluba 35c5110a37
move solve to a separate file 2022-11-15 18:51:43 +01:00
Marek Kaluba ca6a17acca
use ProgressMeter instead of poormans progressbar 2022-11-15 00:27:09 +01:00
Marek Kaluba e6bd862e7a
small tweaks 2022-11-14 19:50:50 +01:00
Marek Kaluba 703e69fc62
fix actions on SpNs 2022-11-14 19:50:27 +01:00
Marek Kaluba 971e07b819
move to using sparse matrices in symmetrized sdp
dense are faster for small sizes only
2022-11-14 19:50:09 +01:00
Marek Kaluba 2f89538eb0
add __fast_recursive_dot to speed up constraints 2022-11-14 19:49:42 +01:00
Marek Kaluba 227e82d551
make augmented compute_sos fast 2022-11-14 19:45:41 +01:00
Marek Kaluba 65286e09d2
fix _positive_direction so that Adj works for N=2 2022-11-13 14:27:58 +01:00
Marek Kaluba fee9c38537
remove old unsed definition for Positive 2022-11-08 12:00:20 +01:00
Marek Kaluba 4e93e78a9b
test: old sq,adj,op = the graded Adj for SL/SAut 2022-11-08 11:49:04 +01:00
Marek Kaluba 903c1683ff
remove unused (old) constraints function 2022-11-08 09:59:26 +01:00
Marek Kaluba 6e9cc8e0f6
action on FPGroupElement via action on the underlying word 2022-11-08 09:59:07 +01:00
Marek Kaluba 7907137fb5
define group_algebra 2022-11-07 18:44:42 +01:00
Marek Kaluba 633f065488
cleanup Project.toml and fix imports 2022-11-07 18:44:27 +01:00
Marek Kaluba f00bfb7ca9
remove old laplacians.jl, 1712.* 1812.* 2022-11-07 16:29:26 +01:00
Marek Kaluba 0e5799862b
add action by alphabet permutation 2022-11-07 16:21:58 +01:00
Marek Kaluba 147211ea7a
add graded-by-root-system Adj 2022-11-07 16:17:41 +01:00
Marek Kaluba 4b8efd2a40
adjust sq, adj, op to the new world 2022-11-07 16:12:05 +01:00
Marek Kaluba 085f6bce3c
replace checksoltuion by certify 2022-11-07 16:01:35 +01:00
Marek Kaluba ecea3dfbcb
add ConstraintMatrix 2022-11-07 15:45:18 +01:00
Marek Kaluba bb0354d3a0
rewrite sos_spds.jl
This includes changes related to:
* SymbolicWedderburn
* new formulation of constraints
* general warmstart for JuMP-^1.3
2022-11-07 15:42:47 +01:00
Marek Kaluba 9511e34de4
replace BlockDecomposition by SymbolicWedderburn 2022-11-07 15:34:30 +01:00
kalmarek 649894236c
make verbose kwarg in BlockDecomposition 2020-10-18 01:08:23 +02:00
kalmarek 1e14380942
rename OrbitData → BlockDecomposition 2020-10-17 11:39:19 +02:00
kalmarek 0d48ff698c
use ^ for group action everywhere 2020-10-17 02:13:59 +02:00
kalmarek 984161c91a
cleanup orbitdata 2020-10-17 02:13:05 +02:00
kalmarek eeecc35232
use @SDconstraint and SOS_problem_(primal|dual) 2020-10-17 02:11:34 +02:00
kalmarek 7ab971eb76
update to SCS-0.7.0 2020-10-17 02:08:51 +02:00