Magma

MAGMA Computational Algebra System

Magma
 •  How to get it
 •  Download
 •  Online Demo
 
Resources
 •  Online Help
 •  Discovering Mathematics with Magma
 •  Citations
 •  How to cite Magma
 •  Links
 •  Contact us
 
[Next][Prev] [Right] [Left] [Up] [Index] [Root]

Aschbacher Reduction

Subsections

Introduction

.

An on-going international research project seeks to develop algorithms to explore the structure of groups having either large order or large degree. The approach relies on the following theorem of Aschbacher [Asc84]:

A matrix group G acting on the finite dimensional K[G]-module V over a finite field K satisfies at least one of the following conditions (which we have simplified slightly for brevity):

(i)
G acts reducibly on V;
(ii)
G acts semilinearly over an extension field of K;
(iii)
G acts imprimitively on V;
(iv)
G preserves a nontrivial tensor-product decomposition of V;
(v)
G has a normal subgroup N, acting absolutely irreducibly on V, which is an extraspecial p-group or 2-group of symplectic type;
(vi)
G preserves a tensor-induced decomposition of V;
(vii)
G acts (modulo scalars) linearly over a proper subfield of K;
(viii)
G contains a classical group in its natural action over K;
(ix)
G is almost simple modulo scalars.

The philosophy underpinning the research program is to attempt to decide that G lies in at least one of the above categories, and to calculate the associated isomorphism or decomposition explicitly.

Groups in Category (i) can be recognised easily by means of the Meataxe functions described in the chapter on R-modules.

Groups which act irreducibly but not absolutely irreducibly on V fall theoretically into Category (ii), and furthermore act linearly over an extension field of K. In fact, absolute irreducibility can be tested using the built-in Magma functions and, by redefining their field to be an extension field L of K and reducing, they can be rewritten as absolutely irreducible groups of smaller dimension, but over L instead of K. We can therefore concentrate on absolutely irreducible matrix groups.

The Magma Aschbacher package currently includes functions which seek to decide membership of categories (ii)-(viii). It was prepared by E.A. O'Brien.

Primitivity

Let G be a subgroup of GL(d, q) and assume that G acts irreducibly on the underlying vector space V. Then G acts imprimitively on V if there is a non-trivial direct sum decomposition V = V1 direct-sum V2 direct-sum ... direct-sum Vr where V1, ..., Vr are permuted by G. In such a case, each block Vi has the same dimension or size, and we have the block system {V1, ..., Vr}. If no such system exists, then G is primitive.

Theoretical details of the algorithm used may be found in Holt, Leedham-Green, O'Brien, & Rees [DFH96b].

SetVerbose ("Smash", 1) will provide information on the progress of the algorithm.

IsPrimitive(G: parameters) : GrpMat -> BoolElt
    BlockSizes: [RngIntElt]             Default: []
Return true if the matrix group G is primitive, false if G is not primitive, or "unknown" if no decision can be reached.

If BlockSizes is supplied, then we search for systems of imprimitivity whose block sizes are in BlockSizes only. Otherwise we consider all valid sizes.

ImprimitiveBasis (G) : GrpMat -> SeqEnum
If the matrix group G is imprimitive, return the change-of-basis matrix which exhibits block structure for G.
Blocks(G) : GrpMat -> SeqEnum
If the matrix group G is imprimitive, return the blocks of imprimitivity of G.
BlocksImage(G) : GrpMat -> GrpPerm
Return the group induced by the action of the matrix group G on the system of imprimitivity.
ImprimitiveAction(G, g) : GrpMat, GrpMatElt -> GrpPermElt
Return action of g on blocks of imprimitivity of the matrix group G.

Example GrpMatFF_IsPrimitive (H58E2)

> MG := GL (4, 7);
> PG := Sym (3);
> G := WreathProduct (MG, PG);
> 
> IsPrimitive (G);
false
We investigate the block system for G.

> B := Blocks (G);
> B;
> #B;
4
> B[1];
Vector space of degree 12, dimension 4 over GF(7)
Generators:
(0 0 0 0 1 0 0 0 0 0 0 0)
(0 0 0 0 0 1 0 0 0 0 0 0)
(0 0 0 0 0 0 1 0 0 0 0 0)
(0 0 0 0 0 0 0 1 0 0 0 0)
Echelonized basis:
(0 0 0 0 1 0 0 0 0 0 0 0)
(0 0 0 0 0 1 0 0 0 0 0 0)
(0 0 0 0 0 0 1 0 0 0 0 0)
(0 0 0 0 0 0 0 1 0 0 0 0)
Now we obtain a permutation representation of G in its action on the blocks.

> P := BlocksImage (G);
> P;
Permutation group P acting on a set of cardinality 3
    (1, 2, 3)
    (2, 3)
> g := G.4 * G.3;
> ImprimitiveAction (G, g);
(1, 2)

Semilinearity

Let G be a subgroup of GL(d, q) and assume that G acts absolutely irreducibly on the underlying vector space V. Assume that a normal subgroup N of G embeds in GL(d/e, qe), for e>1, and a d x d matrix C acts as multiplication by a scalar λ(a field generator of GF(qe)) for that embedding.

We say that G acts as a semilinear group of automorphisms on the d/e-dimensional space if and only if, for each generator g of G, there is an integer i = i(g) such that Cg = gCi, that is, g corresponds to the field automorphism λ-> λi. If so, we have a map from G to the (cyclic) group ( Aut)(GF(qe)), and C centralises the kernel of this map, which thus lies in GL(d, qe).

Theoretical details of the algorithm used may be found in Holt, Leedham-Green, O'Brien, & Rees [DFH96a]

SetVerbose ("SemiLinear", 1) will provide information on the progress of the algorithm.

IsSemiLinear(G) : GrpMat -> BoolElt
Return true if the matrix group G is semilinear, false if G is not semilinear, or "unknown" if no decision can be reached.
DegreeOfFieldExtension(G) : GrpMat -> RngIntElt
The matrix group G is defined over K = GL(d, q). Return the degree e of the extension field of K over which G is semilinear.
CentralisingMatrix(G) : GrpMat -> AlgMatElt
Return the matrix C which centralises the normal subgroup of the matrix group G which acts linearly over the extension field.
FrobeniusAutomorphisms(G) : GrpMat -> SeqEnum
Return a sequence S of positive integers, one for each generator of the matrix group G. The element S[i] is the least positive integer such that G.i - 1 C G.i = CS[i].
WriteOverLargerField(G) : GrpMat -> GrpMat, GrpAb, SeqEnum
Return
(i)
the normal subgroup N of the matrix group G which is the kernel of the map from G to Ce; this subgroup acts linearly over the extension field of K and is precisely the centraliser of C in G;
(ii)
a cyclic group E of order e which is isomorphic to G/N; and
(iii)
a sequence of images of the generators of G in E.

Example GrpMatFF_Semilinearity (H58E3)

We analyze a semilinear group.

> P := GL(6,3);
> g1 := P![0,1,0,0,0,0,-1,0,0,0,0,0,
>          0,0,1,0,0,0,0,0,0,1,0,0,0,0,0,0,1,0,0,0,0,0,0,1];
> g2 := P![-1,0,0,0,1,0,0,-1,0,0,0,1,
>          -1,0,0,0,0,0,0,-1,0,0,0,0,0,0,-1,0,0,0,0,0,0,-1,0,0];
> g3 := P![1,0,0,0,0,0,0,-1,0,0,0,0,
>          0,0,1,0,0,0,0,0,0,-1,0,0,0,0,0,0,1,0,0,0,0,0,0,-1];
> G := sub <P | g1, g2, g3 >;
> 
> IsSemiLinear (G);
true
> DegreeOfFieldExtension (G);
2
> CentralisingMatrix (G);
[2 2 0 0 0 0]
[1 2 0 0 0 0]
[0 0 2 2 0 0]
[0 0 1 2 0 0]
[0 0 0 0 2 2]
[0 0 0 0 1 2]
> FrobeniusAutomorphisms (G);
[ 1, 1, 3 ]
> K, E, phi := WriteOverLargerField (G);

The group K is the kernel of the homomorphism from G into E.

> K.1;
[0 1 0 0 0 0]
[2 0 0 0 0 0]
[0 0 1 0 0 0]
[0 0 0 1 0 0]
[0 0 0 0 1 0]
[0 0 0 0 0 1]

E is the cyclic group of order e while phi gives the sequence of images of G.i in E.

> E;
Abelian Group isomorphic to Z/2
Defined on 1 generator
Relations:
    2*E.1 = 0
> 
> phi;
[ 0, 0, E.1 ]

Tensor Products

Let G be a subgroup of GL(d, K), where K = GF(q), and let V be the natural K[G]-module. We say that G preserves a tensor decomposition of V as U tensor W if there is an isomorphism of V onto U tensor W such that the induced image of G in GL(U tensor W) lies in GL(U) GL(W).

Theoretical details of the algorithm used may be found in Leedham-Green & O'Brien [O'B97a][O'B97b].

SetVerbose ("Tensor", 1) will provide information on the progress of the algorithm.

IsTensor(G: parameters) : GrpMat -> BoolElt
    Factors: [SeqEnum]                  Default: []
Return true if the matrix group G preserves a non-trivial tensor decomposition, false if G is does not preserve a tensor decomposition, or "unknown" if no decision can be reached.

A sequence of valid dimensions Factors for potential factors may be supplied; for all elements [u, w] of Factors, we search for decompositions of V as U tensor W, where U has dimension u and W has dimension w only. Otherwise we consider all valid factorisations.

TensorBasis(G) : GrpMat -> GrpMatElt
Return the change-of-basis matrix which exhibits the tensor decomposition of the matrix group G.
TensorFactors(G) : GrpMat -> GrpMat, GrpMat
Return the tensor factors of the matrix group G.
IsProportional(X, k) : Mtrx, RngIntElt -> BoolElt, Tup
Return true iff the matrix X is composed of k x k blocks which differ only by scalars; if so, return also the tensor decomposition of X.

Example GrpMatFF_Tensor (H58E4)

> P := GL(6, 3);
> 
> g := P![ 0, 1, 1, 2, 1, 0, 2, 2, 1, 2, 1, 1, 1, 0, 2, 1, 2, 2, 1, 2, 2,
>          2, 2, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 2, 2 ];
> 
> h := P![ 1, 0, 2, 1, 1, 2, 0, 0, 2, 0, 0, 2, 2, 0, 1, 0, 2, 1, 2, 1, 2,
>          2, 1, 1, 0, 2, 0, 1, 0, 0, 0, 0, 0, 2, 1, 2 ];
> 
> G := sub< P | g, h >;
> IsTensor(G);
true
> C := TensorBasis(G);

So C is the change-of-basis matrix. If we conjugate G.1 by C, we obtain a visible Kronecker product.

> G.1^C;
[0 0 2 0 2 0]
[0 0 2 2 2 2]
[2 0 0 0 2 0]
[2 2 0 0 2 2]
[0 0 0 0 1 0]
[0 0 0 0 1 1]
> 

We use the function IsProportional to verify that G.1C is a Kronecker product.

> IsProportional(G.1^C, 2);
true 
<
 [2 0]
 [2 2],

 [0 1 1]
 [1 0 1]
 [0 0 2]
>

Finally, we display the tensor factors.

> A := TensorFactors(G);
> A[1];
MatrixGroup(2, GF(3))
Generators:
    [1 2]
    [2 2]


    [2 0]
    [2 2]
> A[2];
MatrixGroup(3, GF(3))
Generators:
    [0 1 0]
    [1 2 1]
    [1 2 0]


    [0 1 1]
    [1 0 1]
    [0 0 2]

Tensor-induced Groups

Let G be a subgroup of GL(d, K), where K = GF(q) and q = pe for some prime p, and let V be the natural K[G]-module. Assume that d has a proper factorisation as ur. We say that G is tensor-induced if G preserves a decomposition of V as U1 tensor U2 tensor ... tensor Ur where each Ui has dimension u > 1 and r > 1, and the set of Ui is permuted by G. If G is tensor-induced, then there is a homomorphism of G into the symmetric group Sr.

Theoretical details of the algorithm used may be found in Leedham-Green & O'Brien [LGO02].

SetVerbose ("TensorInduced", 1) will provide information on the progress of the algorithm.

IsTensorInduced(G : parameters) : GrpMat -> BoolElt
    InducedDegree: RngIntElt            Default: "All"
Return true if the matrix group G is tensor-induced, false if G is not tensor-induced, or "unknown" if no decision can be reached.

If the value of InducedDegree is r, then we search for homomorphisms into the symmetric group of degree r only. Otherwise we consider all valid degrees.

TensorInducedBasis(G) : GrpMat -> GrpMatElt
Return the change-of-basis matrix which exhibits that the matrix group G is tensor-induced.
TensorInducedPermutations(G) : GrpMat -> SeqEnum
Return a sequence whose i-th entry is the homomorphic image of G.i in Sr.
TensorInducedAction(G, g) : GrpMat, GrpMatElt -> GrpPermElt
Return tensor induced action of the matrix group element g.

Example GrpMatFF_TensorInduced (H58E5)

We illustrate the use of the functions for determining if a matrix group is tensor induced.

> G := GL(2, 3);
> S := Sym(3);
> G := TensorWreathProduct(G, S);
> IsTensorInduced(G);
true

We recover the permutations.

> TensorInducedPermutations(G);
[
    Id(S),
    Id(S),
    (1, 2, 3),
    (1, 2)
]

Hence G.1 and G.2 are in the kernel of the homomorphism from G to S. We extract the change-of-basis matrix C and then conjugate G.1 by C, thereby obtaining a visible Kronecker product.

> C := TensorInducedBasis(G);
> x := G.1^C;
> x;
[2 0 0 0 0 0 0 0]
[0 2 0 0 0 0 0 0]
[0 0 2 0 0 0 0 0]
[0 0 0 2 0 0 0 0]
[1 0 0 0 1 0 0 0]
[0 1 0 0 0 1 0 0]
[0 0 1 0 0 0 1 0]
[0 0 0 1 0 0 0 1]

Finally, we verify that x = G.1C is a Kronecker product for each of 2 and 4.

> IsProportional(x, 2);
true
<[2 0]
[0 2], [1 0 0 0]
[0 1 0 0]
[2 0 2 0]
[0 2 0 2]>        
> IsProportional(x, 4);
true
<[2 0 0 0]
[0 2 0 0]
[0 0 2 0]
[0 0 0 2], [1 0]
[2 2]>                 

Normalisers of Extraspecial r-groups and Symplectic 2-groups

Let G ≤GL(d, q), where d=rm for some prime r. If G is contained in the normaliser of an r-group R, of order either r2m + 1 or 22m + 2, then either R is extraspecial (in the first case), or R is a 2-group of symplectic type (that is, a central product of an extraspecial 2-group with a cyclic group of order 4).

If d = r an odd prime, we use the Monte-Carlo algorithm of Niemeyer [Nie04] to decide whether or not G normalises such a subgroup. Otherwise, IsExtraSpecialNormaliser searches for elements of the normal subgroup, and can only reach negative conclusions in certain limited cases. If it cannot reach a conclusion it returns "unknown".

IsExtraSpecialNormaliser(G) : GrpMat -> BoolElt
Return true if the matrix group G normalises an extraspecial r-group or 2-group of symplectic type, false if G is known not to normalise an extraspecial r-group or a 2-group of symplectic type, or "unknown" if it cannot reach a conclusion.

ExtraSpecialParameters(G) : GrpMat -> [RngIntElt, RngIntElt]
Return sequence of integers, r and n, where the extraspecial or symplectic subgroup R normalised by the matrix group G has order rn.

ExtraSpecialGroup(G) : GrpMat -> GrpMat
Return the extraspecial or symplectic subgroup normalised by the matrix group G.
ExtraSpecialNormaliser(G) : GrpMat -> SeqEnum
Return the action of generators of the matrix group G on its normal extraspecial or symplectic subgroup as a sequence of matrices, each of degree 2r, one for each generator of G.
ExtraSpecialAction(G, g) : GrpMat, GrpMatElt -> GrpMatElt
Matrix of degree 2r describing action of element g on extraspecial or symplectic group normalised by the matrix group G.
ExtraSpecialBasis(G) : GrpMat -> GrpMatElt
In the odd prime degree case, return change-of-basis matrix which conjugates normal extraspecial subgroup into a "nice" representation, generated by a diagonal and a permutation matrix.

Example GrpMatFF_ExtraSpecialNormaliser (H58E6)

> F:=GF(8);
> P:=GL(7,F);
> w := PrimitiveElement(F);
> g1:=P![
> w,0,w^2,w^5,0,w^3,w,w,1,w^6,w^3,0,w^4,w,w^2,w^6,w^4,1,w^3,w^3,w^5,
> w^6,w,w^3,1,w^5,0,w^4,1,w^6,w^3,w^6,w^3,w^2,w^2,w^3,w^6,w^6,w^4,1,w^2,w^4,
> w^5,w^4,w^2,w^6,1,w^5,w ];
> g2:=P![w^3,w^4,w^2,w^6,w,w,w^3,w^3,w^4,w,w,w^2,w^3,w^3,w,w^3,w^5,w,1,w^3,w,
> 0,w^2,w^6,w,w^5,1,w,w^6,0,w^3,0,w^4,w,w^5,w^3,w^3,1,w^3,w^5,w^5,w^3,
> w^4,w^6,w,w^6,w^4,w^4,0 ];
> g3:=P![w^5,w^6,w^2,w,w,w^4,w^6,w^6,w^6,w,w^6,w,1,w^3,w,w^6,w^2,w,w^6,w^3,w^6,
> w^2,w^6,w^6,w^3,w,w^6,w^5,0,w^4,w^6,w^6,w,w^2,0,w,w^3,w^5,w^2,w^3,w^4,w^6,
> 0,w^3,w,w^3,w^4,w^3,1];
> gens := [g1,g2,g3];
> G := sub< P | gens >;
> IsExtraSpecialNormaliser(G);
true
> ExtraSpecialParameters (G);
[ 7, 3 ]
> N:=ExtraSpecialNormaliser(G);
> N;
[
    [3 4]
    [1 4],

    [4 3]
    [0 2],

    [1 0]
    [0 1]
]

Writing Representations over Subfields

The algorithm implemented by these functions is due to Glasby, Leedham-Green, and O'Brien [GLGO05]. We also provide access to an earlier algorithm for the non-scalar case developed by Glasby and Howlett [GH97].

IsOverSmallerField (G: parameters) : GrpMat -> BoolElt, GrpMat
    Scalars: BoolElt                    Default: false
    Algorithm: MonStgElt                Default: "GLO"
Decide whether or not an absolutely irreducible group G ≤GL(d, K) has an equivalent representation over a subfield of K. If so, it returns true and the representation over the smallest possible subfield, else it returns false. If the optional argument Scalars is true then decide whether or G modulo scalars has an equivalent representation over a subfield of K. If the optional argument Algorithm is set to "GH", then the non-scalar case uses the original Glasby and Howlett algorithm.
IsOverSmallerField (G, k: parameters) : GrpMat -> BoolElt, GrpMat
    Scalars: BoolElt                    Default: false
    Algorithm: MonStgElt                Default: "GLO"
Decide whether or not an absolutely irreducible group G ≤GL(d, K) has an equivalent representation over a proper subfield of K having degree k. If so, it returns true and the representation over this subfield, else it returns false. If the optional argument Scalars is true then decide whether or G modulo scalars has an equivalent representation over a subfield of K. If the optional argument Algorithm is set to "GH", then the non-scalar case uses the original Glasby and Howlett algorithm.
SmallerField(G) : GrpMat -> FLdFin
If the matrix group G (possibly mod scalars) has a representation over a subfield, return the subfield.
SmallerFieldBasis (G) : GrpMat -> GrpMatElt
Return change of basis matrix for the matrix group G so that G (possibly mod scalars) can be written over a smaller field.
SmallerFieldImage (G, g) : GrpMat, GrpMatElt -> GrpMatElt
The matrix group G can be rewritten (possibly mod scalars) over smaller field; return image of g ∈G in group defined over smaller field.

Example GrpMatFF_IsOverSmallerField (H58E7)

> G := GL (2, GF (3, 2));
> H := GL (2, GF (3, 8));
> K := sub < H | G.1, G.2 >;
> K;
MatrixGroup(2, GF(3^8))
Generators:
    [ $.1^820        0]
    [       0        1]


    [       2        1]
    [       2        0]
> flag, M := IsOverSmallerField (K);    
> flag;
true
> M;
MatrixGroup(2, GF(3^2))
Generators:
    [$.1^7 $.1^2]
    [    1     2]


    [$.1^7 $.1^6]
    [$.1^2 $.1^5]
> F := GF(3, 4);
> G := MatrixGroup<2, F | [ F.1^52, F.1^72, F.1^32, 0 ],
>                                [ 1, 0, F.1^20, 2 ] >;
> flag, X := IsOverSmallerField (G);
> flag;
false
> /* decide if G has an equivalent representation mod scalars */
> flag, X := IsOverSmallerField (G: Scalars := true);
> flag;
true
> X;
MatrixGroup(2, GF(3))
Generators:
    [2 1]
    [1 0]

    [2 1]
    [1 1]
> SmallerField (G);
Finite field of size 3
> SmallerFieldBasis (G);
[F.1^33 F.1^23]
[F.1^43 F.1^63]
> g := G.1 * G.2^2; g;
[F.1^52 F.1^72]
[F.1^32      0]
> SmallerFieldImage (G, g);
[1 2]
[2 0]

WriteOverSmallerField(G, F) : GrpMat, FldFin -> GrpMat, Map
Given a group G of d x d matrices over a finite field E having degree e and a subfield F of E having degree f, write the matrices of G as de/f by de/f matrices over F and return the group and the isomorphism.

Example GrpMatFF_WriteOverSmallerField (H58E8)

> G := GL(2, 4);
> H := WriteOverSmallerField(G, GF(2));
> H;
MatrixGroup(4, GF(2))
Generators:
    [0 1 0 0]
    [1 1 0 0]
    [0 0 1 0]
    [0 0 0 1]


    [1 0 1 0]
    [0 1 0 1]
    [1 0 0 0]
    [0 1 0 0]

Decompositions with Respect to a Normal Subgroup

SearchForDecomposition(G, S) : GrpMat, [GrpMatElt] -> BoolElt
The normal closure N of S in G where S is a sequence of elements in G is constructed by this function which seeks to decide whether or not G, with respect to N, has a decomposition corresponding to one of the categories (ii)--(vi) in the theorem of Aschbacher stated at the beginning of this section. Theoretical details of the algorithms used may be found in Holt, Leedham-Green, O'Brien, & Rees [DFH96a].

In summary, it tests for one of the following possibilities:

(ii)
G acts semilinearly over an extension field L of K, and N acts linearly over L;
(iii)
G acts imprimitively on V and N fixes each block of imprimitivity;
(iv)
G preserves a tensor product decomposition U tensor W of V, where N acts as scalar matrices on U;
(v)
N acts absolutely irreducibly on V and is an extraspecial p-group for some prime p, or 2-group of symplectic type;
(vi)
G preserves a tensor-induced decomposition V = tensor mU of V for some m>1, where N acts absolutely irreducibly on V and fixes each of the m factors.

If one of the listed decompositions is found, then the function reports the type found and returns true; if no decomposition is found with respect to N, then the function returns false. The answer provided by the function is conclusive for decompositions of types (ii)--(v), but a negative answer for (vi) is not necessarily conclusive.

Each involves a decomposition of G with respect to the normal subgroup N (which may sometimes be trivial or scalar). In (ii), N is the subgroup of G acting linearly over the extension field irreducibly on V. In (iii), N is the subgroup which fixes each of the subspaces in the imprimitive decomposition of V. In (iv), it is the subgroup acting as scalar matrices on one of the factors in the tensor-product decomposition. In (v), N is already described, and in (vi), it is the subgroup fixing each of the factors in the tensor-induced decomposition (so N itself falls in Category (iv)).

If any one of these decompositions can be found, then we may be able to obtain an explicit representation of G/N and hence reduce the study of G to a smaller problem. For example, in Category (iii), G/N acts as a permutation group on the subspaces in the imprimitive decomposition of V. Currently only limited facilities are provided to construct G/N.

Progress of the algorithm is provided by SetVerbose ("Smash", 1).

Accessing the Decomposition Information

The access functions described in the sections on Primitivity Testing, Semilinearity, Tensor Products, Tensor Induction, and Normalisers of Extraspecial groups may be used to extract information about decompositions of type (ii), (iii), (iv), (v) and (vi). We illustrate such decompositions below.


Example GrpMatFF_Decompose (H58E9)

We begin with an example where no decomposition exists.

> G := GL(4, 5);
> SearchForDecomposition (G, [G.1]);
Smash: No decomposition found
false

The second example is of an imprimitive decomposition.

> M := GL (4, 7);
> P := Sym (3);
> G := WreathProduct (M, P);
> SearchForDecomposition (G, [G.1, G.2]);
Smash: G is imprimitive
true
> IsPrimitive (G);
false
> BlocksImage (G);
Permutation group acting on a set of cardinality 3
    Id($)
    Id($)
    (1, 2, 3)
    (1, 2)

The third example admits a semilinear decomposition.

> P := GL(6,3);
> g1 := P![0,1,0,0,0,0,-1,0,0,0,0,0,
>          0,0,1,0,0,0,0,0,0,1,0,0,0,0,0,0,1,0,0,0,0,0,0,1];
> g2 := P![-1,0,0,0,1,0,0,-1,0,0,0,1,
>          -1,0,0,0,0,0,0,-1,0,0,0,0,0,0,-1,0,0,0,0,0,0,-1,0,0];
> g3 := P![1,0,0,0,0,0,0,-1,0,0,0,0,
>          0,0,1,0,0,0,0,0,0,-1,0,0,0,0,0,0,1,0,0,0,0,0,0,-1];
> G := sub <P | g1, g2, g3 >;
>
> SearchForDecomposition (G, [g1]);
Smash: G is semilinear
true
> IsSemiLinear (G);
true
> DegreeOfFieldExtension (G);
2
> CentralisingMatrix (G);
[2 2 0 0 0 0]
[1 2 0 0 0 0]
[0 0 2 2 0 0]
[0 0 1 2 0 0]
[0 0 0 0 2 2]
[0 0 0 0 1 2]
> FrobeniusAutomorphisms (G);
[ 1, 1, 3 ]

The fourth example admits a tensor product decomposition.

> F := GF(5);
> G := GL(5, F);
> H := GL(3, F);
> P := GL(15, F);
> A := MatrixAlgebra (F, 5);
> B := MatrixAlgebra (F, 3);
> g1 := A!G.1; g2 := A!G.2;  g3 := A!Identity(G);
> h1 := B!H.1; h2 := B!H.2; h3 := B!Identity(H);
> w := TensorProduct (g1, h3);
> x := TensorProduct (g2, h3);
> y := TensorProduct (g3, h1);
> z := TensorProduct (g3, h2);
> G := sub < P | w, x, y, z>;
> SearchForDecomposition (G, [G.1, G.2]);
Smash: G is a tensor product
true
> IsTensor (G);
true
> TensorBasis (G);
[1 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
[0 0 0 1 0 0 0 0 0 0 0 0 0 0 0]
[0 0 0 0 0 0 1 0 0 0 0 0 0 0 0]
[0 0 0 0 0 0 0 0 0 0 0 0 1 0 0]
[0 0 0 0 0 0 0 0 0 1 0 0 0 0 0]
[4 0 1 0 0 0 0 0 0 0 0 0 0 0 0]
[0 0 0 4 0 1 0 0 0 0 0 0 0 0 0]
[0 0 0 0 0 0 4 0 1 0 0 0 0 0 0]
[0 0 0 0 0 0 0 0 0 0 0 0 4 0 1]
[0 0 0 0 0 0 0 0 0 4 0 1 0 0 0]
[1 4 4 0 0 0 0 0 0 0 0 0 0 0 0]
[0 0 0 1 4 4 0 0 0 0 0 0 0 0 0]
[0 0 0 0 0 0 1 4 4 0 0 0 0 0 0]
[0 0 0 0 0 0 0 0 0 0 0 0 1 4 4]
[0 0 0 0 0 0 0 0 0 1 4 4 0 0 0]

Our fifth example is of a tensor-induced decomposition.

> M := GL (3, GF(2));
> P := Sym (3);
> G := TensorWreathProduct (M, P);
> SearchForDecomposition (G, [G.1]);
Smash: G is tensor induced
true
>
> IsTensorInduced (G);
true
> TensorInducedPermutations (G);
[ Id(P), Id(P), (1, 3, 2), (1, 3) ]

Our final example is of a normaliser of a symplectic group.

> F := GF(5);
> P := GL(4,F);
> g1 := P![ 1,0,0,0,0,4,0,0,2,0,2,3,3,0,4,3];
> g2 := P![ 4,0,0,1,2,4,4,0,1,0,1,2,0,0,0,1];
> g3 := P![ 4,0,1,1,0,1,0,0,0,1,3,4,0,4,3,2];
> g4 := P![ 2,0,4,3,4,4,2,4,0,1,3,4,4,2,0,1];
> g5 := P![ 1,1,3,4,0,0,3,4,2,0,0,4,3,1,3,4];
> g6 := P![ 2,0,0,0,0,2,0,0,0,0,2,0,0,0,0,2];
> G := sub < P | g1, g2, g3, g4, g5, g6 >;
> SearchForDecomposition (G, [G.4]);
Smash: G is normaliser of symplectic 2-group
true
> IsExtraSpecialNormaliser (G);
true
> ExtraSpecialParameters (G);
[2, 6]
> g := G.1 * G.2;
> ExtraSpecialAction(G, g);
[0 1 0 0]
[1 1 0 0]
[0 1 1 1]
[1 1 1 0]
 [Next][Prev] [Right] [Left] [Up] [Index] [Root]
                       

Version: V2.16 of Tue Jan 5 12:58:41 EST 2010

Valid HTML 4.01! Valid CSS!