From 13452577323fa30b7f4468de9e1f6e9c98731ad5 Mon Sep 17 00:00:00 2001
From: kalmarek <kalmar@amu.edu.pl>
Date: Fri, 28 Jun 2019 01:12:32 +0200
Subject: [PATCH] further speedups for WreathProduct on Mats

---
 src/orbitdata.jl | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/src/orbitdata.jl b/src/orbitdata.jl
index 41e08dd..c133a74 100644
--- a/src/orbitdata.jl
+++ b/src/orbitdata.jl
@@ -258,7 +258,10 @@ function (g::WreathProductElem{N})(A::MatElem) where N
     result = similar(A)
     @inbounds for i = 1:size(A,1)
         for j = 1:size(A,2)
-            result[i, j] = A[g.p[i], g.p[j]]*(flips[i]*flips[j])
+            x = A[g.p[i], g.p[j]]
+            result[i, j] = x*(flips[i]*flips[j])
+            # result[i, j] = AbstractAlgebra.mul!(x, x, flips[i]*flips[j])
+            # this mul! needs to be separately defined, but is 2x faster
         end
     end
     return result