passagemath-gap-pkg-semigroups 10.6.29__cp312-abi3-manylinux_2_24_x86_64.manylinux_2_28_x86_64.whl

gap/pkg/semigroups/tst/standard/elements/elements.tst

@@ -0,0 +1,117 @@
+#############################################################################
+##
+#W  standard/elements/elements.tst
+#Y  Copyright (C) 2016-2022                                 Wilf A. Wilson
+##
+##  Licensing information can be found in the README file of this package.
+##
+#############################################################################
+##
+
+#@local S, x
+gap> START_TEST("Semigroups package: standard/elements/elements.tst");
+gap> LoadPackage("semigroups", false);;
+
+#
+gap> SEMIGROUPS.StartTest();
+
+# elements: SEMIGROUPS.IndexPeriodByRank
+gap> x := Transformation(
+> [71, 14, 60, 68, 67, 74, 61, 19, 81, 10, 17, 21, 9, 49, 78,
+>  73, 81, 22, 59, 50, 29, 31, 6, 27, 40, 21, 75, 44, 71, 62, 22, 56, 90, 53,
+>  64, 53, 78, 73, 18, 96, 58, 46, 11, 19, 72, 12, 33, 83, 21, 34, 94, 94, 51,
+>  6, 52, 24, 66, 46, 61, 95, 40, 81, 41, 38, 26, 39, 90, 93, 30, 37, 55, 51,
+>  99, 38, 38, 79, 33, 58, 90, 89, 98, 54, 18, 58, 28, 27, 61, 82, 96, 26, 89,
+>  93, 44, 4, 96, 11, 39, 14, 57, 34]);
+<transformation on 100 pts with rank 64>
+gap> SEMIGROUPS.IndexPeriodByRank(x, RankOfTransformation);
+[ 12, 42 ]
+gap> x := Transformation([1, 1, 2, 3, 4, 5, 6, 7, 8]);;
+gap> SEMIGROUPS.IndexPeriodByRank(x, RankOfTransformation);
+[ 8, 1 ]
+gap> x := Bipartition([[1], [2, -5], [3, -4], [4, -3], [5, -2], [-1]]);
+<bipartition: [ 1 ], [ 2, -5 ], [ 3, -4 ], [ 4, -3 ], [ 5, -2 ], [ -1 ]>
+gap> SEMIGROUPS.IndexPeriodByRank(x, RankOfBipartition);
+[ 1, 2 ]
+gap> x := Bipartition([[1, 2, -1], [3, -2], [4, -3], [5, -4], [6, -5],
+> [7, -6], [8, -7], [9, -8], [-9]]);;
+gap> SEMIGROUPS.IndexPeriodByRank(x, RankOfBipartition);
+[ 8, 1 ]
+gap> x := Bipartition([[1, -2], [2, -1]]);;
+gap> SEMIGROUPS.IndexPeriodByRank(x, RankOfBipartition);
+[ 1, 2 ]
+
+# elements: IndexPeriodOfSemigroupElement
+gap> x := PBR([[-3, -2, -1, 1, 2, 3], [-3, -2, -1, 1, 3], [-3, -2, -1, 1, 2]],
+> [[-3, -2, -1, 1, 2], [-3, 2, 3], [-2, -1, 2, 3]]);;
+gap> IndexPeriodOfSemigroupElement(x);
+[ 2, 1 ]
+gap> IsIdempotent(x);
+false
+gap> IndexPeriodOfSemigroupElement(GeneratorsOfMagma(FreeMagma(1))[1]);
+Error, the argument (a mult. elt.) is not the generator of a semigroup
+gap> x := Matrix(IsBooleanMat, [[0, 1], [1, 0]]);;
+gap> IndexPeriodOfSemigroupElement(x);
+[ 1, 2 ]
+gap> x;;
+gap> IndexPeriodOfSemigroupElement(x);
+[ 1, 2 ]
+gap> x := Matrix(IsBooleanMat,
+> [[0, 0, 0, 1, 0, 0, 0],
+>  [0, 0, 0, 1, 0, 0, 0],
+>  [0, 0, 0, 1, 0, 0, 0],
+>  [0, 0, 0, 0, 0, 1, 0],
+>  [0, 0, 0, 0, 1, 0, 0],
+>  [0, 1, 0, 0, 0, 0, 0],
+>  [0, 0, 0, 0, 0, 1, 0]]);;
+gap> IndexPeriodOfSemigroupElement(x);
+[ 1, 3 ]
+
+# elements: SmallestIdempotentPower
+gap> x := PBR([[-3, -2, -1, 1, 2, 3], [-3, -2, -1, 1, 3], [-3, -2, -1, 1, 2]],
+> [[-3, -2, -1, 1, 2], [-3, 2, 3], [-2, -1, 2, 3]]);;
+gap> SmallestIdempotentPower(x);
+2
+gap> First([1 .. 100], i -> IsIdempotent(x ^ i)) = last;
+true
+gap> x := PBR([[-2, 2], [-1]], [[1], [1, 2]]);;
+gap> SmallestIdempotentPower(x);
+4
+gap> First([1 .. 100], i -> IsIdempotent(x ^ i)) = last;
+true
+gap> x := PBR([[-2], [-1]], [[1], [1, 2]]);;
+gap> SmallestIdempotentPower(x);
+2
+gap> First([1 .. 100], i -> IsIdempotent(x ^ i)) = last;
+true
+gap> x := Transformation(
+> [71, 14, 60, 68, 67, 74, 61, 19, 81, 10, 17, 21, 9, 49, 78,
+>  73, 81, 22, 59, 50, 29, 31, 6, 27, 40, 21, 75, 44, 71, 62, 22, 56, 90, 53,
+>  64, 53, 78, 73, 18, 96, 58, 46, 11, 19, 72, 12, 33, 83, 21, 34, 94, 94, 51,
+>  6, 52, 24, 66, 46, 61, 95, 40, 81, 41, 38, 26, 39, 90, 93, 30, 37, 55, 51,
+>  99, 38, 38, 79, 33, 58, 90, 89, 98, 54, 18, 58, 28, 27, 61, 82, 96, 26, 89,
+>  93, 44, 4, 96, 11, 39, 14, 57, 34]);;
+gap> SmallestIdempotentPower(x);
+42
+gap> First([1 .. 100], i -> IsIdempotent(x ^ i)) = last;
+true
+gap> x := Transformation([2, 6, 4, 5, 3, 7, 8, 8]);;
+gap> SmallestIdempotentPower(x);
+6
+gap> First([1 .. 100], i -> IsIdempotent(x ^ i)) = last;
+true
+gap> SmallestIdempotentPower(GeneratorsOfMagma(FreeMagma(1))[1]);
+Error, the argument (a mult. elt.) is not the generator of a semigroup
+
+# elements: IsMultiplicativeZero
+gap> S := SymmetricInverseMonoid(3);;
+gap> x := MultiplicativeZero(S);
+<empty partial perm>
+gap> IsMultiplicativeZero(S, x);
+true
+gap> IsMultiplicativeZero(S, PartialPerm([1]));
+false
+
+# 
+gap> SEMIGROUPS.StopTest();
+gap> STOP_TEST("Semigroups package: standard/elements/elements.tst");

gap/pkg/semigroups/tst/standard/elements/ffmat.tst

@@ -0,0 +1,349 @@
+#############################################################################
+##
+#W  standard/elements/ffmat.tst
+#Y  Copyright (C) 2015-2022                                Markus Pfeiffer
+##
+##  Licensing information can be found in the README file of this package.
+##
+#############################################################################
+##
+
+#@local M, S, coll, ht, mat, mat2, r1, r2, rs, x, y, z
+gap> START_TEST("Semigroups package: standard/elements/ffmat.tst");
+gap> LoadPackage("semigroups", false);;
+
+# Set info levels and user preferences
+gap> SEMIGROUPS.StartTest();
+
+# MatrixTest1: Create
+gap> M := Matrix(GF(2),
+> Z(2) * [[1, 0, 1, 0], [0, 1, 0, 1], [0, 1, 0, 0], [0, 0, 0, 1]]);
+<a 4x4 matrix over GF2>
+
+# MatrixTest2: EvalPrintString
+gap> EvalString(String(M)) = M;
+true
+
+# MatrixTest3: IndexPeriodOfSemigroupElement
+gap> M := Matrix(GF(3),
+> [[Z(3) ^ 0, 0 * Z(3), 0 * Z(3), 0 * Z(3), 0 * Z(3)],
+>  [Z(3) ^ 0, 0 * Z(3), 0 * Z(3), 0 * Z(3), 0 * Z(3)],
+>  [0 * Z(3), Z(3), 0 * Z(3), 0 * Z(3), 0 * Z(3)],
+>  [0 * Z(3), 0 * Z(3), Z(3), 0 * Z(3), 0 * Z(3)],
+>  [0 * Z(3), 0 * Z(3), 0 * Z(3), 0 * Z(3), 0 * Z(3)]]);;
+gap> IndexPeriodOfSemigroupElement(M);
+[ 3, 1 ]
+gap> M := Matrix(GF(2 ^ 2),
+> [[Z(2 ^ 2), 0 * Z(2), 0 * Z(2), 0 * Z(2)],
+>  [Z(2 ^ 2), 0 * Z(2), 0 * Z(2), 0 * Z(2)],
+>  [0 * Z(2), Z(2) ^ 0, 0 * Z(2), 0 * Z(2)],
+>  [0 * Z(2), 0 * Z(2), Z(2) ^ 0, 0 * Z(2)]]);;
+gap> IndexPeriodOfSemigroupElement(M);
+[ 3, 3 ]
+
+# MatrixTest4: SmallestIdempotentPower
+gap> M := Matrix(GF(5),
+> [[Z(5) ^ 3, 0 * Z(5), Z(5) ^ 0],
+>  [Z(5) ^ 3, 0 * Z(5), 0 * Z(5)],
+>  [0 * Z(5), 0 * Z(5), 0 * Z(5)]]);;
+gap> SmallestIdempotentPower(M);
+4
+gap> M := Matrix(GF(2 ^ 2),
+> [[Z(2 ^ 2), 0 * Z(2), Z(2 ^ 2) ^ 2, Z(2 ^ 2) ^ 2],
+>  [Z(2 ^ 2) ^ 2, 0 * Z(2), 0 * Z(2), Z(2) ^ 0],
+>  [Z(2 ^ 2) ^ 2, 0 * Z(2), 0 * Z(2), 0 * Z(2)],
+>  [0 * Z(2), 0 * Z(2), 0 * Z(2), 0 * Z(2)]]);;
+gap> SmallestIdempotentPower(M);
+15
+
+# Test ELM_LIST and IsBound
+gap> mat := Matrix(GF(3), Z(3) * [[1, 0, 0], [0, 1, 0], [0, 0, 1]]);;
+gap> mat[1];
+[ Z(3), 0*Z(3), 0*Z(3) ]
+gap> mat[4];
+Error, row index 4 exceeds 3, the number of rows
+
+# Test \*
+gap> x := Matrix(GF(3), Z(3) * [[1, 2], [1, 2]]);;
+gap> y := Matrix(GF(3), Z(3) * [[0, 0], [0, 0]]);;
+gap> x * y;
+[ [ 0*Z(3), 0*Z(3) ], [ 0*Z(3), 0*Z(3) ] ]
+gap> y * x;
+[ [ 0*Z(3), 0*Z(3) ], [ 0*Z(3), 0*Z(3) ] ]
+
+# Test OneImmutable
+gap> x := Matrix(GF(3), Z(3) * [[1, 2], [1, 2]]);;
+gap> y := OneImmutable(x);
+[ [ Z(3)^0, 0*Z(3) ], [ 0*Z(3), Z(3)^0 ] ]
+gap> x * y = x;
+true
+gap> y * x = x;
+true
+
+# Test RandomMatrix
+gap> x := RandomMatrix(GF(8), 1);;
+gap> NrRows(x);
+1
+gap> x := RandomMatrix(GF(8), 10);;
+gap> NrRows(x);
+10
+gap> x := RandomMatrix(GF(3), 3, [3]);;
+gap> Rank(x);
+3
+gap> x := RandomMatrix(GF(3), 10, [3 .. 8]);;
+gap> Rank(x) in [3 .. 8];
+true
+gap> x := RandomMatrix(GF(3), 5, [0]);
+[ [ 0*Z(3), 0*Z(3), 0*Z(3), 0*Z(3), 0*Z(3) ], 
+  [ 0*Z(3), 0*Z(3), 0*Z(3), 0*Z(3), 0*Z(3) ], 
+  [ 0*Z(3), 0*Z(3), 0*Z(3), 0*Z(3), 0*Z(3) ], 
+  [ 0*Z(3), 0*Z(3), 0*Z(3), 0*Z(3), 0*Z(3) ], 
+  [ 0*Z(3), 0*Z(3), 0*Z(3), 0*Z(3), 0*Z(3) ] ]
+gap> x := RandomMatrix(GF(3), 5, [1, -2, 1314]);
+Error, the list of ranks has to consist of numbers > 0 and < n
+
+# Test \ = , \< for row basis
+gap> mat := Matrix(GF(3 ^ 2),
+> [[Z(3 ^ 2) ^ 7, Z(3 ^ 2) ^ 3, Z(3 ^ 2) ^ 5], [Z(3 ^ 2) ^ 3, Z(3),
+>   Z(3 ^ 2) ^ 2], [Z(3 ^ 2) ^ 7, Z(3 ^ 2) ^ 7, Z(3 ^ 2)]]);;
+gap> r1 := RowSpaceBasis(mat);
+<rowbasis of rank 3 over GF(3^2)>
+gap> r2 := RowSpaceBasis(TransposedMat(mat));
+<rowbasis of rank 3 over GF(3^2)>
+gap> r1 = r2;
+true
+gap> r2 < r1;
+false
+gap> mat := Matrix(GF(3 ^ 2),
+> [[Z(3 ^ 2) ^ 7, Z(3 ^ 2) ^ 7, Z(3 ^ 2) ^ 6], [Z(3 ^ 2), Z(3 ^ 2) ^ 2,
+>   Z(3) ^ 0], [Z(3 ^ 2) ^ 6, Z(3 ^ 2) ^ 7, Z(3 ^ 2) ^ 5]]);;
+gap> RowSpaceBasis(TransposedMat(mat)) = RowSpaceBasis(mat);
+false
+gap> RowSpaceBasis(TransposedMat(mat)) < RowSpaceBasis(mat);
+true
+gap> RowSpaceBasis(TransposedMat(mat)) > RowSpaceBasis(mat);
+false
+
+# Test printing viewing etc for row basis
+gap> mat := Matrix(GF(3 ^ 2),
+> [[Z(3 ^ 2) ^ 7, Z(3 ^ 2) ^ 3, Z(3 ^ 2) ^ 5], [Z(3 ^ 2) ^ 3, Z(3),
+>   Z(3 ^ 2) ^ 2], [Z(3 ^ 2) ^ 7, Z(3 ^ 2) ^ 7, Z(3 ^ 2)]]);;
+gap> rs := RowSpaceBasis(mat);;
+gap> EvalString(String(rs)) = rs;
+true
+
+# Test RightInverse and LeftInverse
+gap> mat := Matrix(GF(2 ^ 2),
+> [[Z(2 ^ 2), Z(2) ^ 0, Z(2 ^ 2) ^ 2, Z(2) ^ 0, Z(2 ^ 2) ^ 2],
+>  [Z(2 ^ 2) ^ 2, Z(2 ^ 2), Z(2 ^ 2) ^ 2, Z(2 ^ 2) ^ 2, Z(2) ^ 0],
+>  [Z(2) ^ 0, 0 * Z(2), Z(2) ^ 0, 0 * Z(2), 0 * Z(2)],
+>  [Z(2 ^ 2) ^ 2, Z(2 ^ 2) ^ 2, Z(2) ^ 0, Z(2) ^ 0, 0 * Z(2)],
+>  [Z(2) ^ 0, 0 * Z(2), Z(2 ^ 2), 0 * Z(2), Z(2 ^ 2) ^ 2]]);;
+gap> RightInverse(mat);
+[ [ 0*Z(2), 0*Z(2), Z(2^2)^2, 0*Z(2), 0*Z(2) ], 
+  [ Z(2^2)^2, Z(2^2), Z(2)^0, Z(2^2), 0*Z(2) ], 
+  [ 0*Z(2), 0*Z(2), Z(2^2), 0*Z(2), 0*Z(2) ], 
+  [ Z(2^2), Z(2)^0, Z(2^2)^2, 0*Z(2), 0*Z(2) ], 
+  [ 0*Z(2), Z(2)^0, Z(2^2)^2, Z(2^2)^2, 0*Z(2) ] ]
+gap> LeftInverse(mat);
+[ [ 0*Z(2), Z(2^2), Z(2^2), Z(2)^0, Z(2)^0 ], 
+  [ Z(2^2)^2, Z(2^2), Z(2)^0, Z(2^2), 0*Z(2) ], 
+  [ 0*Z(2), Z(2^2), Z(2^2)^2, Z(2)^0, Z(2)^0 ], 
+  [ Z(2^2), Z(2^2), Z(2)^0, Z(2^2), Z(2^2) ], 
+  [ 0*Z(2), 0*Z(2), 0*Z(2), 0*Z(2), 0*Z(2) ] ]
+gap> Inverse(mat);
+fail
+
+# Test Inverse
+gap> mat := Matrix(GF(2 ^ 2), [[Z(2 ^ 2) ^ 2]]);;
+gap> Inverse(mat);
+[ [ Z(2^2) ] ]
+gap> mat := Matrix(GF(2 ^ 2),
+> [[Z(2 ^ 2), Z(2 ^ 2) ^ 2, Z(2 ^ 2)], [Z(2) ^ 0, Z(2 ^ 2), Z(2 ^ 2)],
+>  [Z(2) ^ 0, Z(2 ^ 2) ^ 2, Z(2 ^ 2)]]);;
+gap> Inverse(mat);
+[ [ Z(2^2), 0*Z(2), Z(2^2) ], [ 0*Z(2), Z(2)^0, Z(2)^0 ], 
+  [ Z(2)^0, Z(2^2), 0*Z(2) ] ]
+
+# Test ComputeRowSpaceAndTransformation
+gap> ComputeRowSpaceAndTransformation(2);
+Error, Assertion failure
+gap> ComputeRowSpaceAndTransformation(Matrix(GF(9), [[Z(9)]]));
+[ <rowbasis of rank 1 over GF(3^2)>, [ [ Z(3^2)^7 ] ], [ [ Z(3^2) ] ] ]
+
+# Test Matrix, checker, 1/1
+gap> mat := Matrix(GF(5), [[0 * Z(5), Z(5) ^ 3], [Z(5) ^ 2, Z(5) ^ 0]]);
+[ [ 0*Z(5), Z(5)^3 ], [ Z(5)^2, Z(5)^0 ] ]
+
+# Test One, 1/1
+gap> mat := Matrix(GF(5), [[0 * Z(5), Z(5) ^ 3], [Z(5) ^ 2, Z(5) ^ 0]]);
+[ [ 0*Z(5), Z(5)^3 ], [ Z(5)^2, Z(5)^0 ] ]
+gap> One(mat);
+[ [ Z(5)^0, 0*Z(5) ], [ 0*Z(5), Z(5)^0 ] ]
+
+# Test RandomMatrix, 1/1
+gap> mat := RandomMatrix(GF(5), 10);
+< mutable compressed matrix 10x10 over GF(5) >
+
+# Test \*, works, 1/3
+gap> mat := Matrix(GF(5), [[0 * Z(5), Z(5) ^ 3], [Z(5), 0 * Z(5)]]);
+[ [ 0*Z(5), Z(5)^3 ], [ Z(5), 0*Z(5) ] ]
+gap> mat ^ 2;
+[ [ Z(5)^0, 0*Z(5) ], [ 0*Z(5), Z(5)^0 ] ]
+
+# Test \*, wrong base field, 3/3
+gap> mat := Matrix(GF(5), [[0 * Z(5), Z(5) ^ 3], [Z(5), 0 * Z(5)]]);
+[ [ 0*Z(5), Z(5)^3 ], [ Z(5), 0*Z(5) ] ]
+gap> mat2 := Matrix(GF(7), [[0 * Z(7), Z(7) ^ 3], [Z(7), 0 * Z(7)]]);
+[ [ 0*Z(7), Z(7)^3 ], [ Z(7), 0*Z(7) ] ]
+gap> mat * mat2;
+Error, <x> and <y> have different characteristic
+
+# Test Display, 1/1
+gap> mat := Matrix(GF(11), [[Z(11) ^ 9, 0 * Z(11), Z(11), Z(11) ^ 9, 0 * Z(11)],
+>   [Z(11) ^ 3, Z(11) ^ 4, 0 * Z(11), Z(11) ^ 2, Z(11) ^ 7],
+>   [Z(11) ^ 9, Z(11) ^ 3, Z(11) ^ 5, Z(11) ^ 4, Z(11) ^ 4],
+>   [Z(11) ^ 6, Z(11), Z(11) ^ 7, Z(11) ^ 3, Z(11) ^ 5],
+>   [0 * Z(11), Z(11) ^ 8, Z(11) ^ 3, Z(11) ^ 6, Z(11) ^ 0]]);
+[ [ Z(11)^9, 0*Z(11), Z(11), Z(11)^9, 0*Z(11) ], 
+  [ Z(11)^3, Z(11)^4, 0*Z(11), Z(11)^2, Z(11)^7 ], 
+  [ Z(11)^9, Z(11)^3, Z(11)^5, Z(11)^4, Z(11)^4 ], 
+  [ Z(11)^6, Z(11), Z(11)^7, Z(11)^3, Z(11)^5 ], 
+  [ 0*Z(11), Z(11)^8, Z(11)^3, Z(11)^6, Z(11)^0 ] ]
+gap> Display(mat);
+  6  .  2  6  .
+  8  5  .  4  7
+  6  8 10  5  5
+  9  2  7  8 10
+  .  3  8  9  1
+
+# Test RandomMatrix for finite field, dim, rank
+gap> x := RandomMatrix(GF(3), 3, 1);;
+gap> Rank(x);
+1
+gap> Rank(TransposedMat(x));
+1
+
+# Test MatrixNC
+gap> x := Matrix(GF(3), [[Z(3) ^ 0, 0 * Z(3), Z(3)], [Z(3), 0 * Z(3), Z(3)],
+> [Z(3) ^ 0, Z(3) ^ 0, Z(3)]]);;
+gap> y := Matrix(GF(3), [[0 * Z(3), Z(3), Z(3) ^ 0], [0 * Z(3), Z(3) ^ 0, 0 * Z(3)],
+>  [Z(3), 0 * Z(3), Z(3)]]);;
+gap> z := Matrix(Unpack(x), y);
+[ [ Z(3)^0, 0*Z(3), Z(3) ], [ Z(3), 0*Z(3), Z(3) ], [ Z(3)^0, Z(3)^0, Z(3) ] ]
+gap> z = x;
+true
+
+# Test RandomMatrix
+gap> RandomMatrix(GF(11), 0, []);
+Error, no method found! For debugging hints type ?Recovery from NoMethodFound
+Error, no 1st choice method found for `RandomMatrixOp' on 3 arguments
+gap> RandomMatrix(GF(11), 0, [0]);
+Error, no method found! For debugging hints type ?Recovery from NoMethodFound
+Error, no 1st choice method found for `RandomMatrixOp' on 3 arguments
+gap> RandomMatrix(GF(11), 0, [1]);
+Error, no method found! For debugging hints type ?Recovery from NoMethodFound
+Error, no 1st choice method found for `RandomMatrixOp' on 3 arguments
+
+# Test AsList for IsPlistMatrixRep 
+gap> S := Semigroup(Transformation([2, 3, 1]));;
+gap> IsGroupAsSemigroup(S);
+true
+gap> IsGroupAsSemigroup(AsSemigroup(IsMatrixOverFiniteFieldSemigroup, S));
+true
+
+# Test BaseDomain (error)
+gap> coll := [Matrix(GF(3), [[Z(3), Z(3), 0 * Z(3)],
+>                            [Z(3), Z(3) ^ 0, Z(3)],
+>                            [0 * Z(3), 0 * Z(3), Z(3)]]),
+>             Matrix(GF(2 ^ 2), [[Z(2 ^ 2) ^ 2, Z(2 ^ 2) ^ 2, Z(2) ^ 0],
+>                              [Z(2) ^ 0, 0 * Z(2), Z(2) ^ 0],
+>                              [0 * Z(2), 0 * Z(2), Z(2 ^ 2) ^ 2]])];;
+gap> BaseDomain(coll);
+Error, no method found! For debugging hints type ?Recovery from NoMethodFound
+Error, no 1st choice method found for `BaseDomain' on 1 arguments
+
+# Test OneMutable
+gap> x := Matrix(GF(3), [[Z(3) ^ 0, 0 * Z(3), Z(3)], [Z(3), 0 * Z(3), Z(3)],
+> [Z(3) ^ 0, Z(3) ^ 0, Z(3)]]);;
+gap> OneMutable(x);
+[ [ Z(3)^0, 0*Z(3), 0*Z(3) ], [ 0*Z(3), Z(3)^0, 0*Z(3) ], 
+  [ 0*Z(3), 0*Z(3), Z(3)^0 ] ]
+
+# Test \< 
+gap> x := Matrix(GF(3), [[Z(3) ^ 0, 0 * Z(3), Z(3)], [Z(3), 0 * Z(3), Z(3)],
+> [Z(3) ^ 0, Z(3) ^ 0, Z(3)]]);;
+gap> y := Matrix(GF(3), [[0 * Z(3), Z(3), Z(3) ^ 0], [0 * Z(3), Z(3) ^ 0, 0 * Z(3)],
+>  [Z(3), 0 * Z(3), Z(3)]]);;
+gap> x < y;
+false
+gap> y < x;
+true
+gap> x < x;
+false
+
+# Test TransposedMat
+gap> x := Matrix(GF(2 ^ 2),
+> [[0 * Z(2), Z(2) ^ 0, 0 * Z(2)], [Z(2) ^ 0, Z(2 ^ 2) ^ 2, Z(2 ^ 2)],
+>  [0 * Z(2), 0 * Z(2), 0 * Z(2)]]);;
+gap> TransposedMat(x);
+[ [ 0*Z(2), Z(2)^0, 0*Z(2) ], [ Z(2)^0, Z(2^2)^2, 0*Z(2) ], 
+  [ 0*Z(2), Z(2^2), 0*Z(2) ] ]
+gap> TransposedMat(TransposedMat(x)) = x;
+true
+
+# Test hashing (for finite field row basis)
+gap> x := Matrix(GF(2 ^ 2),
+> [[0 * Z(2), Z(2) ^ 0, 0 * Z(2)], [Z(2) ^ 0, Z(2 ^ 2) ^ 2, Z(2 ^ 2)],
+>  [0 * Z(2), 0 * Z(2), 0 * Z(2)]]);;
+gap> ht := HTCreate(RowSpaceBasis(x));;
+gap> HTAdd(ht, RowSpaceBasis(x), true);;
+gap> HTAdd(ht, RowSpaceBasis(Matrix(GF(11), [[Z(11)]])), true);;
+gap> ht := HTCreate(RowSpaceBasis(Matrix(GF(11), [[Z(11)]])));;
+gap> HTAdd(ht, RowSpaceBasis(Matrix(GF(11), [[Z(11)]])), true);;
+gap> HTAdd(ht, RowSpaceBasis(Matrix(GF(7), [[Z(7)]])), true);;
+gap> HTAdd(ht, RowSpaceBasis(x), true);;
+gap> HTAdd(ht, RowSpaceBasis(ZeroMatrix(GF(4), 2, 2)), true);;
+gap> ht := HTCreate(RowSpaceBasis(ZeroMatrix(GF(4), 2, 2)));;
+gap> HTAdd(ht, RowSpaceBasis(ZeroMatrix(GF(4), 2, 2)), true);;
+
+# TraceMat
+gap> x := Matrix(GF(2 ^ 2),
+> [[0 * Z(2), Z(2) ^ 0, 0 * Z(2)], 
+>  [Z(2) ^ 0, Z(2 ^ 2) ^ 2, Z(2 ^ 2)],
+>  [0 * Z(2), 0 * Z(2), 0 * Z(2)]]);;
+gap> TraceMat(x);
+Z(2^2)^2
+
+# Finite field matrix attributes TryNextMethod
+gap> x := Matrix(Integers, [[-1, -2, 0], [0, 3, -1], [1, 0, -3]]);
+<3x3-matrix over Integers>
+gap> RowSpaceBasis(x);
+Error, no method found! For debugging hints type ?Recovery from NoMethodFound
+Error, no 2nd choice method found for `RowSpaceBasis' on 1 arguments
+gap> RowSpaceTransformation(x);
+Error, no method found! For debugging hints type ?Recovery from NoMethodFound
+Error, no 2nd choice method found for `RowSpaceTransformation' on 1 arguments
+gap> RightInverse(x);
+Error, no method found! For debugging hints type ?Recovery from NoMethodFound
+Error, no 2nd choice method found for `RightInverse' on 1 arguments
+gap> LeftInverse(x);
+Error, no method found! For debugging hints type ?Recovery from NoMethodFound
+Error, no 2nd choice method found for `LeftInverse' on 1 arguments
+gap> RowSpaceTransformationInv(x);
+Error, no method found! For debugging hints type ?Recovery from NoMethodFound
+Error, no 2nd choice method found for `RowSpaceTransformationInv' on 1 argumen\
+ts
+gap> x := Matrix(GF(2 ^ 2),
+> [[0 * Z(2), Z(2) ^ 0, 0 * Z(2)], [Z(2) ^ 0, Z(2 ^ 2) ^ 2, Z(2 ^ 2)],
+>  [0 * Z(2), 0 * Z(2), 0 * Z(2)]]);;
+gap> RowSpaceTransformation(x);
+[ [ Z(2^2)^2, Z(2)^0, 0*Z(2) ], [ Z(2)^0, 0*Z(2), 0*Z(2) ], 
+  [ 0*Z(2), 0*Z(2), Z(2)^0 ] ]
+gap> RowSpaceTransformationInv(x);
+[ [ 0*Z(2), Z(2)^0, 0*Z(2) ], [ Z(2)^0, Z(2^2)^2, 0*Z(2) ], 
+  [ 0*Z(2), 0*Z(2), Z(2)^0 ] ]
+
+#
+gap> SEMIGROUPS.StopTest();
+gap> STOP_TEST("Semigroups package: standard/elements/ffmat.tst");