diff --git a/OrbitDecomposition.jl b/OrbitDecomposition.jl
index bb501eb..23b5d6f 100644
--- a/OrbitDecomposition.jl
+++ b/OrbitDecomposition.jl
@@ -52,20 +52,19 @@ end
 import Nemo.elements
 elements(F::Nemo.FinField) = FFEltsIter(F)
 
+###############################################################################
+#
+#  Action of Premutations on Nemo.MatElem
+#
+###############################################################################
 
-function AutFG_emb(A::AutGroup, g::WreathProductElem)
-    isa(A.objectGroup, FreeGroup) || throw("Not an Aut(Fₙ)")
-    parent(g).P.n == length(A.objectGroup.gens) || throw("No natural embedding of $(parent(g)) into $A")
-    powers = [(elt == parent(elt)() ? 0: 1) for elt in g.n.elts]
-    elt = reduce(*, [A(Groups.flip_autsymbol(i))^pow for (i,pow) in enumerate(powers)])
-    Groups.r_multiply!(elt, [Groups.perm_autsymbol(g.p)])
-    return elt
+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
 
-function (g::WreathProductElem)(a::AutGroupElem)
-    g = AutFG_emb(parent(a),g)
-    return g*a*g^-1
-end
 
 function orbit_decomposition(G::Nemo.Group, E::Vector, rdict=GroupRings.reverse_dict(E))