From 578a75d807993720ff29a9985d05e82fdd148cc4 Mon Sep 17 00:00:00 2001
From: kalmar <kalmar@amu.edu.pl>
Date: Fri, 9 Jun 2017 22:03:35 +0200
Subject: [PATCH] Add action of permutation group on SLn

---
 OrbitDecomposition.jl | 14 --------------
 SL_orbit.jl           |  9 +++++++++
 2 files changed, 9 insertions(+), 14 deletions(-)

diff --git a/OrbitDecomposition.jl b/OrbitDecomposition.jl
index c45e777..fbeb6f2 100644
--- a/OrbitDecomposition.jl
+++ b/OrbitDecomposition.jl
@@ -2,7 +2,6 @@ push!(LOAD_PATH, "./")
 
 using Nemo
 using Groups
-using WreathProducts
 
 using GroupRings
 using PropertyT
@@ -52,19 +51,6 @@ end
 import Nemo.elements
 elements(F::Nemo.FinField) = FFEltsIter(F)
 
-###############################################################################
-#
-#  Action of Premutations on Nemo.MatElem
-#
-###############################################################################
-
-function (p::Nemo.perm)(A::Nemo.MatElem)
-    length(p.d) == A.r == A.c || throw("Can't act via $p on matrix of size ($(A.r), $(A.c))")
-    R = parent(A)
-    inv_p = inv(p)
-    return R(Nemo.matrix_repr(p))*A*R(Nemo.matrix_repr(inv_p))
-end
-
 ###############################################################################
 #
 #  Orbit stuff
diff --git a/SL_orbit.jl b/SL_orbit.jl
index 76c71d3..44694a1 100644
--- a/SL_orbit.jl
+++ b/SL_orbit.jl
@@ -19,6 +19,14 @@ function (g::WreathProducts.WreathProductElem)(A::Nemo.MatElem)
     return G*A*inv_G
 end
 
+function (p::Nemo.perm)(A::Nemo.MatElem)
+    length(p.d) == A.r == A.c || throw("Can't act via $p on matrix of size ($(A.r), $(A.c))")
+    R = parent(A)
+    inv_p = inv(p)
+    return R(Nemo.matrix_repr(p))*A*R(Nemo.matrix_repr(inv_p))
+end
+
+
 ###############################################################################
 #
 #  Generating set
@@ -143,6 +151,7 @@ function main()
    info(logger, "Symmetric generating set of size $(length(S))")
    info(logger, S)
    AutS = WreathProduct(Nemo.FiniteField(2,1, "a")[1], PermutationGroup(N))
+   # AutS = PermutationGroup(N)
 
    solver = SCS.SCSSolver(eps=tol, max_iters=iterations, verbose=true, linearsolver=SCS.Indirect)