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

small tweaks

Browse Source
This commit is contained in:
Marek Kaluba 2022-11-14 19:50:50 +01:00
parent 703e69fc62
commit e6bd862e7a
No known key found for this signature in database
GPG Key ID: 8BF1A3855328FC15
3 changed files with 6 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/gradings.jl
View File

@ -6,7 +6,8 @@ Roots.Root(e::MatrixGroups.ElementaryMatrix{N}) where {N} =
function Roots.Root(s::MatrixGroups.ElementarySymplectic{N}) where {N} function Roots.Root(s::MatrixGroups.ElementarySymplectic{N}) where {N}
if s.symbol === :A if s.symbol === :A
return Roots.𝕖(N ÷ 2, s.i) - Roots.𝕖(N ÷ 2, s.j) return Roots.𝕖(N ÷ 2, s.i) - Roots.𝕖(N ÷ 2, s.j)
else#if s.symbol === :B else
@assert s.symbol === :B
n = N ÷ 2 n = N ÷ 2
i, j = ifelse(s.i <= n, s.i, s.i - n), ifelse(s.j <= n, s.j, s.j - n) i, j = ifelse(s.i <= n, s.i, s.i - n), ifelse(s.j <= n, s.j, s.j - n)
return (-1)^(s.i > s.j) * (Roots.𝕖(n, i) + Roots.𝕖(n, j)) return (-1)^(s.i > s.j) * (Roots.𝕖(n, i) + Roots.𝕖(n, j))

6
test/1703.09680.jl
View File

@ -79,15 +79,15 @@ end
@test λ > 1 @test λ > 1
m = PropertyT.sos_problem_dual(elt, unit) m = PropertyT.sos_problem_dual(elt, unit)
PropertyT.solve(m, scs_optimizer( PropertyT.solve(m, cosmo_optimizer(
eps=1e-10, eps=1e-6,
max_iters=5_000, max_iters=5_000,
accel=50, accel=50,
alpha=1.9, alpha=1.9,
)) ))
@test JuMP.termination_status(m) in (JuMP.ALMOST_OPTIMAL, JuMP.OPTIMAL) @test JuMP.termination_status(m) in (JuMP.ALMOST_OPTIMAL, JuMP.OPTIMAL)
@test JuMP.objective_value(m) ≈ 1.5 atol = 1e-3 @test JuMP.objective_value(m) ≈ 1.5 atol = 1e-2
end end
@testset "SAut(F₂)" begin @testset "SAut(F₂)" begin

2
test/graded_adj.jl
View File

@ -85,7 +85,7 @@
Δ, Δs Δ, Δs
end end
@testset "Adj correctness: genus=$genus" begin @testset "Adj numerics for genus=$genus" begin
all_subtypes = ( all_subtypes = (
:A₁, :C₁, Symbol("A₁×A₁"), Symbol("C₁×C₁"), Symbol("A₁×C₁"), :A₂, :C₂ :A₁, :C₁, Symbol("A₁×A₁"), Symbol("C₁×C₁"), Symbol("A₁×C₁"), :A₂, :C₂