matrix_algebra: put basis in the usual (row,column) order.

diff --git a/mjo/hurwitz.py b/mjo/hurwitz.py
index 1f7c9dc3781a844d6d77a1e63168c52340bf6770..e2800345cf550b3f6970fa62592c0d43a1e85531 100644 (file)
--- a/mjo/hurwitz.py
+++ b/mjo/hurwitz.py
@@ -319,6 +319,7 @@ class HurwitzMatrixAlgebraElement(MatrixAlgebraElement):
             sage: A = HurwitzMatrixAlgebra(2, QQbar, ZZ)
             sage: M = A([ [ I,   2*I],
             ....:         [ 3*I, 4*I] ])
+            sage: M.conjugate_transpose()
             +------+------+
             | -1*I | -3*I |
             +------+------+

diff --git a/mjo/matrix_algebra.py b/mjo/matrix_algebra.py
index d67347b3aa1252acaf165448d1ea35ebb4a77cd7..84aa8d28bb3de9e6797b9b0a209eb96c3a14dc22 100644 (file)
--- a/mjo/matrix_algebra.py
+++ b/mjo/matrix_algebra.py
@@ -198,8 +198,6 @@ class MatrixAlgebra(CombinatorialFreeModule):
         # sticking a "1" in each position doesn't give us a basis for
         # the space. We actually need to stick each of e0, e1, ...  (a
         # basis for the entry algebra itself) into each position.
-        I = range(n)
-        J = range(n)
         self._entry_algebra = entry_algebra
 
         # Needs to make the (overridden) method call when, for example,
@@ -207,8 +205,8 @@ class MatrixAlgebra(CombinatorialFreeModule):
         # lies to us.
         entry_basis = self.entry_algebra_gens()
 
-        basis_indices = [(i,j,e) for i in range(n)
-                                 for j in range(n)
+        basis_indices = [(i,j,e) for j in range(n)
+                                 for i in range(n)
                                  for e in entry_basis]
 
         super().__init__(scalars,