From 2ef8de7d42fdc59a30d27b0e1a9fcc5d9c22653f Mon Sep 17 00:00:00 2001
From: kalmarek <kalmar@amu.edu.pl>
Date: Mon, 1 Jul 2019 01:46:32 +0200
Subject: [PATCH] =?UTF-8?q?add=20commented=20test=20on=20Adj=E2=82=84=20+?=
 =?UTF-8?q?=20100=20Op=E2=82=84=20in=20SAut(F=E2=82=84)?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 test/1812.03456.jl | 21 +++++++++++++++++++++
 1 file changed, 21 insertions(+)

diff --git a/test/1812.03456.jl b/test/1812.03456.jl
index f674227..77c9191 100644
--- a/test/1812.03456.jl
+++ b/test/1812.03456.jl
@@ -221,4 +221,25 @@ end
         end
     end
 
+    # @testset "Adj₄ + 100 Op₄ ∈ ISAut(F₄) is SOS" begin
+    #     N = 4
+    #     halfradius = 2
+    #     M = SAut(FreeGroup(N))
+    #     S = PropertyT.generating_set(M)
+    #     Δ = PropertyT.Laplacian(S, halfradius)
+    #     RG = parent(Δ)
+    #     orbit_data = PropertyT.OrbitData(RG, WreathProduct(PermGroup(2), PermGroup(N)))
+    #     orbit_data = PropertyT.decimate(orbit_data);
+    #
+    #     @time elt = PropertyT.Adj(RG) + 100PropertyT.Op(RG)
+    #     UB = 0.05
+    #
+    #     warm = nothing
+    #
+    #     residual, λ, warm = check_positivity(elt, Δ, orbit_data, UB, warm, with_solver=with_SCS(20_000, accel=10))
+    #     @info "obtained λ and residual" λ norm(residual, 1)
+    #
+    #     @test 2^2*norm(residual, 1) < λ # i.e. we can certify positivity
+    #     @test 2^2*norm(residual, 1) < λ/100
+    # end
 end