Commit Graph

583 Commits

Author SHA1 Message Date
Marek Kaluba 083e1b4735
@info "checking in (...) arithmetic" before the actual check 2024-02-19 18:48:28 +01:00
Marek Kaluba 526def226e
PG.degree is not guaranteed to be constant 2024-02-19 18:42:42 +01:00
Marek Kaluba 3b9c75604c
in reconstruct: accumulate to a single matrix 2024-02-19 18:39:01 +01:00
Marek Kaluba c7491f37f9
make sure Q'*Q hits the fast path for IntervalMatrices 2024-02-19 18:20:35 +01:00
Marek Kaluba cd90e072c4
back to IntervalArithmetic-0.20 to keep julia-1.6 2024-02-16 12:09:11 +01:00
Marek Kaluba c9f68eee8f
update IntervalArithmetic to 0.21 2024-02-15 22:44:23 +01:00
Marek Kaluba 27ce603872
import packages as ... instead of using
This includes:
* PermutationGroups as PG
* SymbolicWedderburn as SW
* StarAlgebras as SA
2024-02-15 22:42:17 +01:00
Marek Kaluba b23975a627
update and fix versions for dependencies
i.e. PermutationGroups, GroupsCore, SymbolicWedderburn
2024-02-14 11:00:05 +01: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