%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% %A operatio.tex GAP documentation Martin Schoenert %% %A @(#)$Id: operatio.tex,v 3.20 1994/06/30 15:33:58 vfelsch Rel $ %% %Y Copyright 1990-1992, Lehrstuhl D fuer Mathematik, RWTH Aachen, Germany %% %% This file describes the functions that deal with operations of groups. %% %H $Log: operatio.tex,v $ %H Revision 3.20 1994/06/30 15:33:58 vfelsch %H included final changes for version 3.4 %H %H Revision 3.19 1994/06/17 19:03:10 vfelsch %H examples adjusted to version 3.4 %H %H Revision 3.18 1994/06/10 02:51:09 vfelsch %H updated examples %H %H Revision 3.17 1993/07/27 08:18:35 fceller %H added missing index entries %H %H Revision 3.16 1993/02/19 10:48:42 gap %H adjustments in line length and spelling %H %H Revision 3.15 1993/02/18 14:52:37 felsch %H another example fixed %H %H Revision 3.14 1993/02/15 14:16:36 felsch %H DefineName eliminated %H %H Revision 3.13 1993/02/12 14:44:22 felsch %H examples adjusted to line length 72 %H %H Revision 3.12 1993/02/02 12:48:20 felsch %H long lines fixed %H %H Revision 3.11 1993/02/01 13:50:06 felsch %H example fixed %H %H Revision 3.10 1993/01/27 10:09:53 fceller %H various fixes of examples %H %H Revision 3.9 1992/06/01 19:18:28 martin %H changed 'Blocks', must operate transitively on %H %H Revision 3.8 1992/04/30 11:59:19 martin %H changed a few sentences to avoid bold non-roman fonts %H %H Revision 3.7 1992/04/06 16:18:22 martin %H fixed some more typos %H %H Revision 3.6 1992/03/25 15:37:32 martin %H added new sections for group homomorphisms %H %H Revision 3.5 1992/03/11 15:59:44 martin %H added the examples %H %H Revision 3.4 1992/03/10 12:17:54 martin %H added 'IsRegular' and 'IsSemiRegular' %H %H Revision 3.3 1992/01/31 16:33:08 martin %H added 'IsEquivalentOperation' %H %H Revision 3.2 1992/01/23 13:06:56 martin %H changed a reference %H %H Revision 3.1 1992/01/15 13:22:44 martin %H changed a reference %H %H Revision 3.0 1991/12/27 16:10:27 martin %H initial revision under RCS %H %% \Chapter{Operations of Groups} One of the most important tools in group theory is the *operation* or *action* of a group on a certain set. We say that a group $G$ operates on a set $D$ if we have a function that takes each $d \in D$ and each $g \in G$ to another element $d^g \in D$, which we call the image of $d$ under $g$, such that $d^{identity} = d$ and $(d^g)^h = d^{gh}$ for each $d \in D$ and $g,h \in G$. This is equivalent to saying that an operation is a homomorphism of the group $G$ into the full symmetric group on $D$. We usually call $D$ the *domain* of the operation and its elements *points*. An example of the usage of the functions in this package can be found in the introduction to {\GAP} (see "About Operations of Groups"). In {\GAP} group elements usually operate through the power operator, which is denoted by the caret '\^'. It is possible however to specify other operations (see "Other Operations"). First this chapter describes the functions that take a single element of the group and compute cycles of this group element and related information (see "Cycle", "CycleLength", "Cycles", and "CycleLengths"), and the function that describes how a group element operates by a permutation that operates the same way on '[1..]' (see "Permutation"). Next come the functions that test whether an orbit has minimal or maximal length and related functions (see "IsFixpoint", "IsFixpointFree", "DegreeOperation", "IsTransitive", and "Transitivity"). Next this chapter describes the functions that take a group and compute orbits of this group and related information (see "Orbit", "OrbitLength", "Orbits", and "OrbitLengths"). Next are the functions that compute the permutation group

that operates on '[ 1 .. Length() ]' in the same way that operates on , and the corresponding homomorphism from to

(see "Operation", "OperationHomomorphism"). Next is the functions that compute block systems, i.e., partitions of such that operates on the sets of the partition (see "Blocks"), and the function that tests whether has such a nontrivial partitioning under the operation of (see "IsPrimitive"). Finally come the functions that relate an orbit of on with the subgroup of that fixes the first point in the orbit (see "Stabilizer"), and the cosets of this subgroup in (see "RepresentativeOperation" and "RepresentativesOperation"). All functions described in this chapter are in 'LIBNAME/\"operatio.g\"'. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{Other Operations} \index{OnPoints}% \index{OnPairs}% \index{OnTuples}% \index{OnSets}% \index{OnRight}% \index{OnLeft}% \index{OnRightCosets}% \index{OnLeftCosets}% \index{OnLines}% The functions in the operation package generally compute with the operation of group elements defined by the canonical operation that is denoted with the caret ('\^') in {\GAP}. However they also allow you to specify other operations. Such operations are specified by functions, which are accepted as optional argument by all the operations package functions. This function must accept two arguments. The first argument will be the point and the second will be the group element. The function must return the image of the point under the group element. As an example, the function 'OnPairs' that specifies the operation on pairs could be defined as follows\\ | OnPairs := function ( pair, g ) return [ pair[1] ^ g, pair[2] ^ g ]; end; | The following operations are predefined. 'OnPoints':\\ specifies the canonical default operation. Passing this function is equivalent to specifying no operation. This function exists because there are places where the operation in not an option. 'OnPairs':\\ specifies the componentwise operation of group elements on pairs of points, which are represented by lists of length 2. 'OnTuples':\\ specifies the componentwise operation of group elements on tuples of points, which are represented by lists. 'OnPairs' is the special case of 'OnTuples' for tuples with two elements. 'OnSets':\\ specifies the operation of group elements on sets of points, which are represented by sorted lists of points without duplicates (see "Sets"). 'OnRight':\\ specifies that group elements operate by multiplication from the right. 'OnLeft':\\ specifies that group elements operate by multiplication from the left. 'OnRightCosets':\\ specifies that group elements operate by multiplication from the right on sets of points, which are represented by sorted lists of points without duplicates (see "Sets"). 'OnLeftCosets':\\ specifies that group elements operate by multiplication from the left on sets of points, which are represented by sorted lists of points without duplicates (see "Sets"). 'OnLines':\\ specifies that group elements, which must be matrices, operate on lines, which are represented by vectors with first nonzero coefficient one. That is, 'OnLines' multiplies the vector by the group element and then divides the vector by the first nonzero coefficient. Note that it is your responsibility to make sure that the elements of the domain on which you are operating are already in normal form. The reason is that all functions will compare points using the '=' operation. For example, if you are operating on sets with 'OnSets', you will get an error message it not all elements of the domain are sets. | gap> Cycle( (1,2), [2,1], OnSets ); Error, OnSets: must be a set | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{Cycle} 'Cycle( , )' \\ 'Cycle( , , )' 'Cycle' returns the orbit of the point , which may be an object of arbitrary type, under the group element as a list of points. The points in the cycle of under the group element are those for which a power $g^i$ exists such that $d^{g^i} = e$. The first point in the list returned by 'Cycle' is the point itself, the ordering of the other points is such that each point is the image of the previous point. 'Cycle' accepts a function of two arguments and as optional third argument, which specifies how the element operates (see "Other Operations"). | gap> Cycle( (1,5,3,8)(4,6,7), 3 ); [ 3, 8, 1, 5 ] gap> Cycle( (1,5,3,8)(4,6,7), [3,4], OnPairs ); [ [ 3, 4 ], [ 8, 6 ], [ 1, 7 ], [ 5, 4 ], [ 3, 6 ], [ 8, 7 ], [ 1, 4 ], [ 5, 6 ], [ 3, 7 ], [ 8, 4 ], [ 1, 6 ], [ 5, 7 ] ] | 'Cycle' calls \\ 'Domain([]).operations.Cycle( , , )' \\ and returns the value. Note that the third argument is not optional for the functions called this way. The default function called this way is 'GroupElementsOps.Cycle', which starts with and applies to the last point repeatedly until is reached again. Special categories of group elements overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{CycleLength} 'CycleLength( , )' \\ 'CycleLength( , , )' 'CycleLength' returns the length of the orbit of the point , which may be an object of arbitrary type, under the group elements . See "Cycle" for the definition of cycles. 'CycleLength' accepts a function of two arguments and as optional third argument, which specifies how the group element operates (see "Other Operations"). | gap> CycleLength( (1,5,3,8)(4,6,7), 3 ); 4 gap> CycleLength( (1,5,3,8)(4,6,7), [3,4], OnPairs ); 12 | 'CycleLength' calls \\ 'Domain([]).operations.CycleLength( , , )' \\ and returns the value. Note that the third argument is not optional for the functions called this way. The default function called this way is 'GroupElementsOps.CycleLength', which starts with and applies to the last point repeatedly until is reached again. Special categories of group elements overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{Cycles} 'Cycles( , )' \\ 'Cycles( , , )' 'Cycles' returns the set of cycles of the group element on the domain , which must be a list of points of arbitrary type, as a set of lists of points. See "Cycle" for the definition of cycles. It is allowed that is a proper subset of a domain, i.e., that is not invariant under the operation of . In this case is silently replaced by the smallest superset of which is invariant. The first point in each cycle is the smallest point of in this cycle. The ordering of the other points is such that each point is the image of the previous point. If is invariant under , then because 'Cycles' returns a set of cycles, i.e., a sorted list, and because cycles are compared lexicographically, and because the first point in each cycle is the smallest point in that cycle, the list returned by 'Cycles' is in fact sorted with respect to the smallest point in the cycles. 'Cycles' accepts a function of two arguments and as optional third argument, which specifies how the element operates (see "Other Operations"). | gap> Cycles( (1,5,3,8)(4,6,7), [3,5,7] ); [ [ 3, 8, 1, 5 ], [ 7, 4, 6 ] ] gap> Cycles( (1,5,3,8)(4,6,7), [[1,3],[4,6]], OnPairs ); [ [ [ 1, 3 ], [ 5, 8 ], [ 3, 1 ], [ 8, 5 ] ], [ [ 4, 6 ], [ 6, 7 ], [ 7, 4 ] ] ] | 'Cycles' calls \\ 'Domain([]).operations.Cycles( , , )' \\ and returns the value. Note that the third argument is not optional for the functions called this way. The default function called this way is 'GroupElementsOps.Cycles', which takes elements from , computes their orbit, removes all points in the orbit from , and repeats this until has been emptied. Special categories of group elements overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{CycleLengths} 'CycleLengths( , )' \\ 'CycleLengths( , , )' 'CycleLengths' returns a list of the lengths of the cycles of the group element on the domain , which must be a list of points of arbitrary type. See "Cycle" for the definition of cycles. It is allowed that is a proper subset of a domain, i.e., that is not invariant under the operation of . In this case is silently replaced by the smallest superset of which is invariant. The ordering of the lengths of cycles in the list returned by 'CycleLengths' corresponds to the list of cycles returned by 'Cycles', which is ordered with respect to the smallest point in each cycle. 'CycleLengths' accepts a function of two arguments and as optional third argument, which specifies how the element operates (see "Other Operations"). | gap> CycleLengths( (1,5,3,8)(4,6,7), [3,5,7] ); [ 4, 3 ] gap> CycleLengths( (1,5,3,8)(4,6,7), [[1,3],[4,6]], OnPairs ); [ 4, 3 ] | 'CycleLengths' calls \\ 'Domain([]).operations.CycleLengths( , , )' \\ and returns the value. Note that the third argument is not optional for the functions called this way. The default function called this way is 'GroupElementsOps.CycleLengths', which takes elements from , computes their orbit, removes all points in the orbit from , and repeats this until has been emptied. Special categories of group elements overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{Permutation} 'Permutation( , )' \\ 'Permutation( , , )' 'Permutation' returns a permutation that operates on the points '[1..Length(D)]' in the same way that the group element operates on the domain , which may be a list of arbitrary type. It is not allowed that is a proper subset of a domain, i.e., must be invariant under the element . 'Permutation' accepts a function of two arguments and as optional third argument, which specifies how the element operates (see "Other Operations"). | gap> Permutation( (1,5,3,8)(4,6,7), [4,7,6] ); (1,3,2) gap> D := [ [1,4], [1,6], [1,7], [3,4], [3,6], [3,7], > [4,5], [5,6], [5,7], [4,8], [6,8], [7,8] ];; gap> Permutation( (1,5,3,8)(4,6,7), D, OnSets ); ( 1, 8, 6,10, 2, 9, 4,11, 3, 7, 5,12) | 'Permutation' calls \\ 'Domain([]).operations.Permutation( , , )' \\ and returns the value. Note that the third argument is not optional for the functions called this way. The default function called this way is 'GroupElementsOps.Permutation', which simply applies to all the points of , finds the position of the image in , and finally applies 'PermList' (see "PermList") to the list of those positions. Actually this is not quite true. Because finding the position of an image in a sorted list is so much faster than finding it in , 'GroupElementsOps.Permutation' first sorts a copy of and remembers how it had to rearrange the elements of to achieve this. Special categories of group elements overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{IsFixpoint} 'IsFixpoint( , )' \\ 'IsFixpoint( , , )' 'IsFixpoint' returns 'true' if the point is a fixpoint under the operation of the group . We say that is a *fixpoint* under the operation of if every element of maps to itself. This is equivalent to saying that each generator of maps to itself. As a special case it is allowed that the first argument is a single group element, though this does not make a lot of sense, since in this case 'IsFixpoint' simply has to test '\^ = '. 'IsFixpoint' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> IsFixpoint( g, 1 ); false gap> IsFixpoint( g, [6,7,8], OnSets ); true | 'IsFixpoint' is so simple that it does all the work by itself, and, unlike the other functions described in this chapter, does not dispatch to another function. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{IsFixpointFree} 'IsFixpointFree( , )' \\ 'IsFixpointFree( , , )' 'IsFixpointFree' returns 'true' if the group operates without a fixpoint (see "IsFixpoint") on the domain , which must be a list of points of arbitrary type. We say that operates *fixpoint free* on the domain if each point of is moved by at least one element of . This is equivalent to saying that each point of is moved by at least one generator of . This definition also applies in the case that is a proper subset of a domain, i.e., that is not invariant under the operation of . As a special case it is allowed that the first argument is a single group element. 'IsFixpointFree' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> IsFixpointFree( g, [1..8] ); true gap> sets := Combinations( [1..8], 3 );; Length( sets ); 56 # a list of all three element subsets of '[1..8]' gap> IsFixpointFree( g, sets, OnSets ); false | 'IsFixpointFree' calls \\ '.operations.IsFixpointFree( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.IsFixpointFree', which simply loops over the elements of and applies to each all generators of , and tests whether each is moved by at least one generator. This function is seldom overlaid, because it is very difficult to improve it. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{DegreeOperation} 'DegreeOperation( , )' \\ 'DegreeOperation( , , )' 'DegreeOperation' returns the degree of the operation of the group on the domain , which must be a list of points of arbitrary type. The *degree* of the operation of on is defined as the number of points of that are properly moved by at least one element of . This definition also applies in the case that is a proper subset of a domain, i.e., that is not invariant under the operation of . 'DegreeOperation' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> DegreeOperation( g, [1..10] ); 8 gap> sets := Combinations( [1..8], 3 );; Length( sets ); 56 # a list of all three element subsets of '[1..8]' gap> DegreeOperation( g, sets, OnSets ); 55 | 'DegreeOperation' calls \\ '.operations.DegreeOperation( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.DegreeOperation', which simply loops over the elements of and applies to each all generators of , and counts those that are moved by at least one generator. This function is seldom overlaid, because it is very difficult to improve it. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{IsTransitive} 'IsTransitive( , )' \\ 'IsTransitive( , , )' 'IsTransitive' returns 'true' if the group operates transitively on the domain , which must be a list of points of arbitrary type. We say that a group acts *transitively* on a domain if and only if for every pair of points and there is an element of such that $d^g = e$. An alternative characterization of this property is to say that as a set is equal to the orbit of every single point. It is allowed that is a proper subset of a domain, i.e., that is not invariant under the operation of . In this case 'IsTransitive' checks whether for every pair of points , of there is an element of , such that $d^g = e$. This can also be characterized by saying that is a subset of the orbit of every single point. 'IsTransitive' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> IsTransitive( g, [1..8] ); false gap> IsTransitive( g, [1,6] ); false # note that the domain need not be invariant gap> sets := Combinations( [1..5], 3 );; Length( sets ); 10 # a list of all three element subsets of '[1..5]' gap> IsTransitive( g, sets, OnSets ); true | 'IsTransitive' calls \\ '.operations.IsTransitive( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.IsTransitive', which tests whether is a subset of the orbit of the first point in . This function is seldom overlaid, because it is difficult to improve it. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{Transitivity} 'Transitivity( , )' \\ 'Transitivity( , , )' 'Transitivity' returns the degree of transitivity of the group on the domain , which must be a list of points of arbitrary type. If does not operate transitively on then 'Transitivity' returns 0. The *degree of transitivity* of the operation of on is the largest such that operates -fold transitively on . We say that operates -*fold transitively* on if it operates transitively on (see "IsTransitive") and the stabilizer of one point of operates '-1'-fold transitively on 'Difference(,[])'. Because the stabilizers of the points of are conjugate this is equivalent to saying that the stabilizer of *each* point of operates '-1'-fold transitively on 'Difference(,[])'. It is not allowed that is a proper subset of a domain, i.e., must be invariant under the operation of . 'Transitivity' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> Transitivity( g, [1..8] ); 0 gap> Transitivity( g, [1..5] ); 3 gap> sets := Combinations( [1..5], 3 );; Length( sets ); 10 # a list of all three element subsets of '[1..5]' gap> Transitivity( g, sets, OnSets ); 1 | 'Transitivity' calls \\ '.operations.Transitivity( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.Transitivity', which first tests whether operates transitively on . If so, it returns \\ 'Transitivity(Stabilizer(,Difference(,[[1]]),)+1'; \\ if not, it simply returns 0. Special categories of groups overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{IsRegular} 'IsRegular( , )' 'IsRegular( , , )' 'IsRegular' returns 'true' if the group operates regularly on the domain , which must be a list of points of arbitrary type, and 'false' otherwise. A group operates *regularly* on a domain if it operates transitively and no element of other than the idenity leaves a point of fixed. An equal characterisation is that operates transitively on and the stabilizer of any point of is trivial. Yet another characterisation is that the operation of on is equivalent to the operation of on its elements by multiplication from the right. It is not allowed that is a proper subset of a domain, i.e., must be invariant under the operation of . 'IsRegular' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> IsRegular( g, [1..5] ); false gap> IsRegular( g, Elements(g), OnRight ); true gap> g := Group( (1,2,3), (3,4,5) );; gap> IsRegular( g, Orbit( g, [1,2,3], OnTuples ), OnTuples ); true | 'IsRegular' calls \\ '.operations.IsRegular( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.IsRegular', which tests if operates transitively and semiregularly on (see "IsTransitive" and "IsSemiRegular"). %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{IsSemiRegular} 'IsSemiRegular( , )'\\ 'IsSemiRegular( , , )' 'IsSemiRegular' returns 'true' if the group operates semiregularly on the domain , which must be a list of points of arbitrary type, and 'false' otherwise. A group operates *semiregularly* on a domain if no element of other than the idenity leaves a point of fixed. An equal characterisation is that the stabilizer of any point of is trivial. Yet another characterisation is that the operation of on is equivalent to multiple copies of the operation of on its elements by multiplication from the right. It is not allowed that is a proper subset of a domain, i.e., must be invariant under the operation of . 'IsSemiRegular' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2)(3,4)(5,7)(6,8), (1,3)(2,4)(5,6)(7,8) );; gap> IsSemiRegular( g, [1..8] ); true gap> g := Group( (1,2)(3,4)(5,7)(6,8), (1,3)(2,4)(5,6,7,8) );; gap> IsSemiRegular( g, [1..8] ); false gap> IsSemiRegular( g, Orbit( g, [1,5], OnSets ), OnSets ); true | 'IsSemiRegular' calls \\ '.operations.IsSemiRegular( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.IsSemiRegular', which computes a permutation group

that operates on '[1..Length()]' in the same way that operates on (see "Operation") and then checks if this permutation group operations semiregularly. This of course only works because this default function is overlaid for permutation groups (see "Operations of Permutation Groups"). %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{Orbit} 'Orbit( , )' \\ 'Orbit( , , )' 'Orbit' returns the orbit of the point , which may be an object of arbitrary type, under the group as a list of points. The points in the orbit of under the group are those points for which a group element of exists such that $d^g = e$. Suppose has generators. First we order the words of the free monoid with abstract generators according to length and for words with equal length lexicographically. So if has two generators called and the ordering is $identity, a, b, a^2, ab, ba, b^2, a^3, ...$. Next we order the elements of that can be written as a product of the generators, i.e., without inverses of the generators, according to the first occurence of a word representing the element in the above ordering. Then the ordering of points in the orbit returned by 'Orbit' is according to the order of the first representative of each point , i.e., the smallest such that $d^g = e$. Note that because the orbit is finite there is for every point in the orbit at least one representative that can be written as a product in the generators of . 'Orbit' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> Orbit( g, 1 ); [ 1, 2, 3, 4, 5 ] gap> Orbit( g, 2 ); [ 2, 3, 1, 4, 5 ] gap> Orbit( g, [1,6], OnPairs ); [ [ 1, 6 ], [ 2, 7 ], [ 3, 6 ], [ 2, 8 ], [ 1, 7 ], [ 4, 6 ], [ 3, 8 ], [ 2, 6 ], [ 1, 8 ], [ 4, 7 ], [ 5, 6 ], [ 3, 7 ], [ 5, 8 ], [ 5, 7 ], [ 4, 8 ] ] | 'Orbit' calls \\ '.operations.Orbit( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.Orbit', which performs an ordinary orbit algorithm. Special categories of groups overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{OrbitLength} 'OrbitLength( , )' \\ 'OrbitLength( , , )' 'OrbitLength' returns the length of the orbit of the point , which may be an object of arbitrary type, under the group . See "Orbit" for the definition of orbits. 'OrbitLength' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> OrbitLength( g, 1 ); 5 gap> OrbitLength( g, 10 ); 1 gap> OrbitLength( g, [1,6], OnPairs ); 15 | 'OrbitLength' calls \\ '.operations.OrbitLength( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.OrbitLength', which performs an ordinary orbit algorithm. Special categories of groups overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{Orbits} 'Orbits( , )' \\ 'Orbits( , , )' 'Orbits' returns the orbits of the group on the domain , which must be a list of points of arbitrary type, as a set of lists of points. See "Orbit" for the definition of orbits. It is allowed that is a proper subset of a domain, i.e., that is not invariant under the operation of . In this case is silently replaced by the smallest superset of which is invariant. The first point in each orbit is the smallest point, the other points of each orbit are ordered in the standard order defined for orbits (see "Orbit"). Because 'Orbits' returns a set of orbits, i.e., a sorted list, and because those orbits are compared lexicographically, and because the first point in each orbit is the smallest point in that orbit, the list returned by 'Orbits' is in fact sorted with respect to the smallest points the orbits. 'Orbits' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> Orbits( g, [1..8] ); [ [ 1, 2, 3, 4, 5 ], [ 6, 7, 8 ] ] gap> Orbits( g, [1,6] ); [ [ 1, 2, 3, 4, 5 ], [ 6, 7, 8 ] ] # the domain is not invariant gap> sets := Combinations( [1..8], 3 );; Length( sets ); 56 # a list of all three element subsets of '[1..8]' gap> Orbits( g, sets, OnSets ); [ [ [ 1, 2, 3 ], [ 1, 2, 4 ], [ 2, 3, 4 ], [ 1, 2, 5 ], [ 1, 3, 4 ], [ 2, 4, 5 ], [ 2, 3, 5 ], [ 1, 4, 5 ], [ 3, 4, 5 ], [ 1, 3, 5 ] ], [ [ 1, 2, 6 ], [ 2, 3, 7 ], [ 1, 3, 6 ], [ 2, 4, 8 ], [ 1, 2, 7 ], [ 1, 4, 6 ], [ 3, 4, 8 ], [ 2, 5, 7 ], [ 2, 3, 6 ], [ 1, 2, 8 ], [ 2, 4, 7 ], [ 1, 5, 6 ], [ 1, 4, 8 ], [ 4, 5, 7 ], [ 3, 5, 6 ], [ 2, 3, 8 ], [ 1, 3, 7 ], [ 2, 4, 6 ], [ 3, 4, 6 ], [ 2, 5, 8 ], [ 1, 5, 7 ], [ 4, 5, 6 ], [ 3, 5, 8 ], [ 1, 3, 8 ], [ 3, 4, 7 ], [ 2, 5, 6 ], [ 1, 4, 7 ], [ 1, 5, 8 ], [ 4, 5, 8 ], [ 3, 5, 7 ] ], [ [ 1, 6, 7 ], [ 2, 6, 7 ], [ 1, 6, 8 ], [ 3, 6, 7 ], [ 2, 6, 8 ], [ 2, 7, 8 ], [ 4, 6, 8 ], [ 3, 7, 8 ], [ 3, 6, 8 ], [ 4, 7, 8 ], [ 5, 6, 7 ], [ 1, 7, 8 ], [ 4, 6, 7 ], [ 5, 7, 8 ], [ 5, 6, 8 ] ], [ [ 6, 7, 8 ] ] ] | 'Orbits' calls \\ '.operations.Orbits( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.Orbits', which takes an element from , computes its orbit, removes all points in the orbit from , and repeats this until has been emptied. Special categories of groups overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{OrbitLengths} 'OrbitLengths( , )' \\ 'OrbitLengths( , , )' 'OrbitLengths' returns a list of the lengths of the orbits of the group on the domain , which may be a list of points of arbitrary type. See "Orbit" for the definition of orbits. It is allowed that is proper subset of a domain, i.e., that is not invariant under the operation of . In this case is silently replaced by the smallest superset of which is invariant. The ordering of the lengths of orbits in the list returned by 'OrbitLengths' corresponds to the list of cycles returned by 'Orbits', which is ordered with respect to the smallest point in each orbit. 'OrbitLengths' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> OrbitLengths( g, [1..8] ); [ 5, 3 ] gap> sets := Combinations( [1..8], 3 );; Length( sets ); 56 # a list of all three element subsets of '[1..8]' gap> OrbitLengths( g, sets, OnSets ); [ 10, 30, 15, 1 ] | 'OrbitLengths' calls \\ '.operations.OrbitLenghts( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.OrbitLengths', which takes an element from , computes its orbit, removes all points in the orbit from , and repeats this until has been emptied. Special categories of groups overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{Operation} 'Operation( , )' \\ 'Operation( , , )' 'Operation' returns a permutation group with the same number of generators as , such that each generator of the permutation group operates on the set '[1..Length(D)]' in the same way that the corresponding generator of the group operates on the domain , which may be a list of arbitrary type. It is not allowed that is a proper subset of a domain, i.e., must be invariant under the element . 'Operation' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). The function 'OperationHomomorphism' (see "OperationHomomorphism") can be used to compute the homomorphism that maps onto the new permutation group. Of course if you are only interested in mapping single elements of into the new permutation group you may also use 'Permutation' (see "Permutation"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> Operation( g, [1..5] ); Group( (1,2,3), (3,4,5) ) gap> Operation( g, Orbit( g, [1,6], OnPairs ), OnPairs ); Group( ( 1, 2, 3, 5, 8,12)( 4, 7, 9)( 6,10)(11,14), ( 2, 4)( 3, 6,11) ( 5, 9)( 7,10,13,12,15,14) ) | 'Operation' calls \\ '.operations.Operation( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.Operation', which simply applies each generator of to all the points of , finds the position of the image in , and finally applies 'PermList' (see "PermList") to the list of those positions. Actually this is not quite true. Because finding the position on an image in a sorted list is so much faster than finding it in , 'GroupElementsOps.Operation' first sorts a copy of and remembers how it had to rearrange the elements of to achieve this. Special categories of groups overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{OperationHomomorphism} \index{group homomorphisms!operation} \index{homomorphisms!operation, group} \index{Image!for OperationHomomorphism} \index{PreImage!for OperationHomomorphism} \index{Kernel!for OperationHomomorphism} 'OperationHomomorphism( ,

)' 'OperationHomomorphism' returns the group homomorphism (see "Group Homomorphisms") from the group to the permutation group

, which must be the result of a prior call to 'Operation' (see "Operation") with or a group of which is a subgroup (see "IsSubgroup") as first argument. | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> h := Operation( g, [1..5] ); Group( (1,2,3), (3,4,5) ) gap> p := OperationHomomorphism( g, h ); OperationHomomorphism( Group( (1,2,3)(6,7), (3,4,5)(7,8) ), Group( (1,2,3), (3,4,5) ) ) gap> (1,4,2,5,3)(6,7,8) ^ p; (1,4,2,5,3) gap> h := Operation( g, Orbit( g, [1,6], OnPairs ), OnPairs ); Group( ( 1, 2, 3, 5, 8,12)( 4, 7, 9)( 6,10)(11,14), ( 2, 4)( 3, 6,11) ( 5, 9)( 7,10,13,12,15,14) ) gap> p := OperationHomomorphism( g, h );; gap> s := SylowSubgroup( g, 2 ); Subgroup( Group( (1,2,3)(6,7), (3,4,5)(7,8) ), [ (7,8), (7,8), (2,5)(3,4), (2,3)(4,5) ] ) gap> Images( p, s ); Subgroup( Group( ( 1, 2, 3, 5, 8,12)( 4, 7, 9)( 6,10)(11,14), ( 2, 4) ( 3, 6,11)( 5, 9)( 7,10,13,12,15,14) ), [ ( 2, 4)( 5, 9)( 7,12)(10,15)(13,14), ( 2, 4)( 5, 9)( 7,12)(10,15)(13,14), ( 2,14)( 3, 6)( 4,13)( 7,15)( 8,11)(10,12), ( 2,12)( 3, 8)( 4, 7)( 6,11)(10,14)(13,15) ] ) gap> OperationHomomorphism( g, Group( (1,2,3), (3,4,5) ) ); Error, Record: element 'operation' must have an assigned value | 'OperationHomomorphism' calls \\ '

.operations.OperationHomomorphism( ,

)' \\ and returns the value. The default function called this way is 'GroupOps.OperationHomomorphism', which uses the fields '

.operationGroup', '

.operationDomain', and '

.operationOperation' (the arguments to the 'Operation' call that created

) to construct a generic homomorphism . This homomorphism uses \\ 'Permutation(,.range.operationDomain,.range.operationOperation)' \\ to compute the image of an element of under . It uses 'Representative' to compute the preimages of an element

of

under . And it computes the kernel by intersecting the cores (see "Core") of the stabilizers (see "Stabilizer") of representatives of the orbits of . Look under *OperationHomomorphism* in the index to see for which groups and operations this function is overlaid. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{Blocks} 'Blocks( , , )' \\ 'Blocks( , , , )' In this form 'Blocks' returns a block system of the domain , which may be a list of points of arbitrary type, under the group , such that the points in the list all lie in the same block. If no such nontrivial block system exists, 'Blocks' returns '[ ]'. must operate transitively on , otherwise an error is signalled. 'Blocks( , )' \\ 'Blocks( , , )' In this form 'Blocks' returns a minimal block system of the domain , which may be a list of points of arbitrary type, under the group . If no nontrivial block system exists, 'Blocks' returns '[ ]'. must operate transitively on , otherwise an error is signalled. A *block system* is a list of blocks with the following properties. Each block of is a subset of . The blocks are pairwise disjoint. The union of blocks is . The image of each block under each element of is as a set equal to some block of the block system. Note that this implies that all blocks contain the same number of elements as operates transitive on . Put differently a block system of is a partition of such that operates with 'OnSets' (see "Other Operations") on . The block system that consists of only singleton sets and the block system consisting only of are called *trivial*. A block system is called *minimal* if there is no nontrivial block system whose blocks are all subsets of the blocks of and whose number of blocks is larger than the number of blocks of . 'Blocks' accepts a function of two arguments and as optional third, resp. fourth, argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> Blocks( g, [1..5] ); [ [ 1 .. 5 ] ] gap> Blocks( g, Orbit( g, [1,2], OnPairs ), OnPairs ); [ [ [ 1, 2 ], [ 3, 2 ], [ 4, 2 ], [ 5, 2 ] ], [ [ 1, 3 ], [ 2, 3 ], [ 4, 3 ], [ 5, 3 ] ], [ [ 1, 4 ], [ 2, 4 ], [ 3, 4 ], [ 5, 4 ] ], [ [ 1, 5 ], [ 2, 5 ], [ 3, 5 ], [ 4, 5 ] ], [ [ 2, 1 ], [ 3, 1 ], [ 4, 1 ], [ 5, 1 ] ] ] | 'Blocks' calls \\ '.operations.Blocks( , , , )' \\ and returns the value. If no seed was given as argument to 'Blocks' it passes the empty list. Note that the fourth argument is not optional for functions called this way. The default function called this way is 'GroupOps.Blocks', which computes a permutation group

that operates on '[1..Length()]' in the same way that operates on (see "Operation") and leaves it to this permutation group to find the blocks. This of course works only because this default function is overlaid for permutation groups (see "Operations of Permutation Groups"). %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{IsPrimitive} 'IsPrimitive( , )' \\ 'IsPrimitive( , , )' 'IsPrimitive' returns 'true' if the group operates primitively on the domain , which may be a list of points of arbitrary type, and 'false' otherwise. A group operates *primitively* on a domain if and only if operates transitively (see "IsTransitive") and has only the trivial block systems (see "Blocks"). 'IsPrimitive' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> IsPrimitive( g, [1..5] ); true gap> IsPrimitive( g, Orbit( g, [1,2], OnPairs ), OnPairs ); false | 'IsPrimitive' calls \\ '.operations.IsPrimitive( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.IsPrimitive', which simply calls 'Blocks( , , )' and tests whether the returned block system is '[ ]'. This function is seldom overlaid, because all the important work is done in 'Blocks'. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{Stabilizer} 'Stabilizer( , )' \\ 'Stabilizer( , , )' 'Stabilizer' returns the stabilizer of the point under the operation of the group . The *stabilizer* $S$ of $d$ in $G$ is the subgroup of those elements $g$ of $G$ that fix $d$, i.e., for which $d^g = d$. The right cosets of $S$ correspond in a canonical way to the points $p$ in the orbit $O$ of $d$ under $G$; namely all elements from a right coset $S g$ map $d$ to the same point $d^g \in O$, and elements from different right cosets $S g$ and $S h$ map $d$ to different points $d^g$ and $d^h$. Thus the index of the stabilizer $S$ in $G$ is equal to the length of the orbit $O$. 'RepresentativesOperation' (see "RepresentativesOperation") computes a system of representatives of the right cosets of $S$ in $G$. 'Stabilizer' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> g.name := "G";; gap> Stabilizer( g, 1 ); Subgroup( G, [ (3,4,5)(7,8), (2,5,3)(6,7) ] ) gap> Stabilizer( g, [1,2,3], OnSets ); Subgroup( G, [ (7,8), (6,8), (2,3)(4,5)(6,7,8), (1,2)(4,5)(6,7,8) ] )| 'Stabilizer' calls \\ '.operations.Stabilizer( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.Stabilizer', which computes the orbit of $d$ under $G$, remembers a representative $r_e$ for each point $e$ in the orbit, and uses Schreier\'s theorem, which says that the stabilizer is generated by the elements $r_e g r_{e^g}^{-1}$. Special categories of groups overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{RepresentativeOperation} 'RepresentativeOperation( , , )' \\ 'RepresentativeOperation( , , , )' 'RepresentativeOperation' returns a representative of the point in the orbit of the point under the group . If = then 'RepresentativeOperation' returns '.identity', otherwise it is not specified which group element 'RepresentativeOperation' will return if there are several that map to . If is not in the orbit of under , 'RepresentativeOperation' returns 'false'. An element $g$ of $G$ is called a *representative* for the point $e$ in the orbit of $d$ under $G$ if $g$ maps $d$ to $e$, i.e., $d^g = e$. Note that the set of such representatives that map $d$ to $e$ forms a right coset of the stabilizer of $d$ in $G$ (see "Stabilizer"). 'RepresentativeOperation' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> RepresentativeOperation( g, 1, 5 ); (1,5,4,3,2)(6,8,7) gap> RepresentativeOperation( g, 1, 6 ); false gap> RepresentativeOperation( g, [1,2,3], [3,4,5], OnSets ); (1,3,5,2,4) gap> RepresentativeOperation( g, [1,2,3,4], [3,4,5,2], OnTuples ); false | 'RepresentativeOperation' calls \\ '.operations.RepresentativeOperation( , , , )' \\ and returns the value. Note that the fourth argument is not optional for functions called this way. The default function called this way is 'GroupOper.RepresentativeOperation', which starts a normal orbit calculation to compute the orbit of under , and remembers for each point how it was obtained, i.e., which generator of took which orbit point to this new point. When the point appears this information can be traced back to write down the representative of as a word in the generators. Special categories of groups overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{RepresentativesOperation} 'RepresentativesOperation( , )' \\ 'RepresentativesOperation( , , )' 'RepresentativesOperation' returns a list of representatives of the points in the orbit of the point under the group . The ordering of the representatives corresponds to the ordering of the points in the orbit as returned by 'Orbit' (see "Orbit"). Therefore 'List( RepresentativesOperation(,), r->\^r ) = Orbit(,)'. An element $g$ of $G$ is called a *representative* for the point $e$ in the orbit of $d$ under $G$ if $g$ maps $d$ to $e$, i.e., $d^g = e$. Note that the set of such representatives that map $d$ to $e$ forms a right coset of the stabilizer of $d$ in $G$ (see "Stabilizer"). The set of all representatives of the orbit of $d$ under $G$ thus forms a system of representatives of the right cosets of the stabilizer of $d$ in $G$. 'RepresentativesOperation' accepts a function of two arguments and as optional third argument, which specifies how the elements of operate (see "Other Operations"). | gap> g := Group( (1,2,3)(6,7), (3,4,5)(7,8) );; gap> RepresentativesOperation( g, 1 ); [ (), (1,2,3)(6,7), (1,3,2), (1,4,5,3,2)(7,8), (1,5,4,3,2) ] gap> Orbit( g, [1,2], OnSets ); [ [ 1, 2 ], [ 2, 3 ], [ 1, 3 ], [ 2, 4 ], [ 1, 4 ], [ 3, 4 ], [ 2, 5 ], [ 1, 5 ], [ 4, 5 ], [ 3, 5 ] ] gap> RepresentativesOperation( g, [1,2], OnSets ); [ (), (1,2,3)(6,7), (1,3,2), (1,2,4,5,3)(6,8,7), (1,4,5,3,2)(7,8), (1,3,2,4,5)(6,8), (1,2,5,4,3)(6,7), (1,5,4,3,2), (1,4,3,2,5)(6,7,8), (1,3,2,5,4) ] | 'RepresentativesOperation' calls \\ '.operations.RepresentativesOperation( , , )' \\ and returns the value. Note that the third argument is not optional for functions called this way. The default function called this way is 'GroupOps.RepresentativesOperation', which computes the orbit of with the normal algorithm, but remembers for each point $e$ in the orbit a representative $r_e$. When a generator $g$ of $G$ takes an old point $e$ to a point $f$ not yet in the orbit, the representative $r_f$ for $f$ is computed as $r_e g$. Special categories of groups overlay this default function with more efficient functions. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \Section{IsEquivalentOperation} 'IsEquivalentOperation( , , , )' \\ 'IsEquivalentOperation( , , , , )' \\ 'IsEquivalentOperation( , , , , )' \\ 'IsEquivalentOperation( , , , , , )' 'IsEquivalentOperation' returns 'true' if operates on in like operates on , and 'false' otherwise. The operations of $G$ on $D$ and $H$ on $E$ are equivalent if they have the same number of generators and there is a permutation $F$ of the elements of $E$ such that for every generator $g$ of $G$ and the corresponding generator $h$ of $H$ we have $Position( D, D_i^g ) = Position( F, F_i^h )$. Note that this assumes that the mapping defined by mapping $G.generators$ to $H.generators$ is a homomorphism (actually an isomorphism of factor groups of $G$ and $H$ represented by the respective operation). 'IsEquivalentOperation' accepts functions and of two arguments and as optional third and sixth arguments, which specify how the elements of and operate (see "Other Operations"). | gap> g := Group( (1,2)(4,5), (1,2,3)(4,5,6) );; gap> h := Group( (2,3)(4,5), (1,2,3)(4,5,6) );; gap> IsEquivalentOperation( g, [1..6], h, [1..6] ); true gap> h := Group( (1,2), (1,2,3) );; gap> IsEquivalentOperation(g,[[1,4],[2,5],[3,6]],OnPairs,h,[1..3]); true gap> h := Group( (1,2), (1,2,3)(4,5,6) );; gap> IsEquivalentOperation( g, [1..6], h, [1..6] ); false gap> h := Group( (1,2,3)(4,5,6), (1,2)(4,5) );; gap> IsEquivalentOperation( g, [1..6], h, [1..6] ); false # the generators must correspond | 'IsEquivalentOperation' calls \\ '.operations.IsEquivalentOperation(,,,,,)' and returns the value. Note that the third and sixth argument are not optional for functions called this way. The default function called this way is 'GroupOps.IsEquivalentOperation', which tries to rearrange so that the above condition is satisfied. This is done one orbit of at a time, and for each such orbit all the orbits of of the same length are tried to see if there is one which can be rearranged as necessary. Special categories of groups overlay this function with more efficient ones. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% %E Emacs . . . . . . . . . . . . . . . . . . . . . local Emacs variables %% %% Local Variables: %% mode: outline %% outline-regexp: "\\\\Chapter\\|\\\\Section" %% fill-column: 73 %% eval: (hide-body) %% End: %%