set alpha=1.95, acceleration_lookback=1 for numerical stability

MCG.jl

@@ -72,8 +72,6 @@ function main()
    radius = parsed_args["radius"]
    N = parsed_args["N"]
 
-   solver = SCSSolver(eps=tol, max_iters=iterations, linearsolver=SCS.Direct)
-
    prefix = "MCG($N)"
    name = "$(prefix)"
 
@@ -101,12 +99,13 @@ function main()
    ];
 
    relations = [relations; [inv(rel) for rel in relations]]
-
    Groups.add_rels!(MCGroup, Dict(rel => MCGroup() for rel in relations))
 
    S = gens(MCGroup)
    S = unique([S; [inv(s) for s in S]])
    Id = MCGroup()
+   
+   solver = SCSSolver(eps=tol, max_iters=iterations, linearsolver=SCS.Direct, alpha=1.9, acceleration_lookback=1)
 
    @time PropertyT.check_property_T(name, S, Id, solver, upper_bound, tol, radius)
    return 0