diff --git a/src/gradings.jl b/src/gradings.jl index 1f20d6b..8f34efa 100644 --- a/src/gradings.jl +++ b/src/gradings.jl @@ -6,7 +6,8 @@ Roots.Root(e::MatrixGroups.ElementaryMatrix{N}) where {N} = function Roots.Root(s::MatrixGroups.ElementarySymplectic{N}) where {N} if s.symbol === :A return Roots.๐•–(N รท 2, s.i) - Roots.๐•–(N รท 2, s.j) - else#if s.symbol === :B + else + @assert s.symbol === :B n = N รท 2 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)) diff --git a/test/1703.09680.jl b/test/1703.09680.jl index 074d263..edf0306 100644 --- a/test/1703.09680.jl +++ b/test/1703.09680.jl @@ -79,15 +79,15 @@ end @test ฮป > 1 m = PropertyT.sos_problem_dual(elt, unit) - PropertyT.solve(m, scs_optimizer( - eps=1e-10, + PropertyT.solve(m, cosmo_optimizer( + eps=1e-6, max_iters=5_000, accel=50, alpha=1.9, )) @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 @testset "SAut(Fโ‚‚)" begin diff --git a/test/graded_adj.jl b/test/graded_adj.jl index 2c20506..bf5a804 100644 --- a/test/graded_adj.jl +++ b/test/graded_adj.jl @@ -85,7 +85,7 @@ ฮ”, ฮ”s end - @testset "Adj correctness: genus=$genus" begin + @testset "Adj numerics for genus=$genus" begin all_subtypes = ( :Aโ‚, :Cโ‚, Symbol("Aโ‚ร—Aโ‚"), Symbol("Cโ‚ร—Cโ‚"), Symbol("Aโ‚ร—Cโ‚"), :Aโ‚‚, :Cโ‚‚