\Chapter{An Example Program}

Here is a similar example to that in chapter "RDS:A basic example". But
now we do a little more handwork to see how things work.  Now we will
calculate the relative difference sets of ``Dembowski-Piper type d''
and order $16$. These difference sets represent projective planes
which admit a quasiregular collineation group such that the fixed
structure is an anti-flag. See \cite{DembowskiPiper}, \cite{Dembowski}
or \cite{RoederDiss} for details.

To have a little more output, you may want to increase "InfoRDS":
\begintt
gap> SetInfoLevel(InfoRDS,3);
\endtt

First, define some parameters and calculate the data needed:
\beginexample
gap> k:=16;;lambda:=1;;groupOrder:=255;; #Diffset parameters
gap> forbiddenGroupOrder:=15;;
gap> maxtest:=10^6;;                     #Bound for sig testing
gap> G:=CyclicGroup(groupOrder);
<pc group of size 255 with 3 generators>
gap> Gdata:=PermutationRepForDiffsetCalculations(G);;
gap> MakeImmutable(Gdata);; 
\endexample

Now the forbidden group is calculated in a very ineffective way. Then
we calculate admissible signatures. As there are only few normal
subgroups in $G$, we calculate them all. For other groups, one should
choose more wisely.

\beginexample
gap> N:=First(NormalSubgroups(Gdata.G),i->Size(i)=forbiddenGroupOrder);
Group([ f1*f3^9, f2*f3^10 ])
gap> globalSigData:=[];;
gap> normals:=Filtered(NormalSubgroups(Gdata.G),n->Size(n) in [2..groupOrder-1]);; 
gap> sigdat:=SignatureDataForNormalSubgroups(normals,globalSigData,
>                             N,Gdata,[k,lambda,maxtest,true]);;
\endexample

The last step gives better results, if a larger <maxtest> is chosen.
But it also takes more time. To find a suitable <maxtest>, the output
of "SignatureDataForNormalSubgroups" can be used, if
`InfoLevel(InfoRDS)' is at least $2$. Note that for recalculating
signatures, you will have to reset <globalSigData> to `[]'. Try experimenting
with <maxtest> to see the effect of signatures for the generation of
startsets.

Now examine the signatures found. Look if there is a normal subgroup
which has only one admissible signature (of course, you can also use
"NormalSgsHavingAtMostNSigs" here):

\beginexample
gap> Set(Filtered(sigdat,s->Size(s.sigs)=1 and Size(s.sigs[1])<=5),i->i.sigs);
[ [ [ 0, 4, 4, 4, 4 ] ], [ [ 4, 4, 8 ] ] ]
\endexample

Great! we'll take the subgroup of index $3$. The cosets modulo this
group will be used to generate startsets and we assume that the
trivial coset contains $8$ elements of the difference set.
So we generate startsets in $U$ and make a first reduction. For this,
we need $U$ and $N$ as lists of integers (recall that difference sets are
asumed to be lists of integers). We will call these lists $Up$ and
$Np$.
Furthermore, we will have to choose a suitable group of automorphisms
for reduction. As $G$ is cyclic, we may take $Gdata.Aac$ here. A good
choice is the stabilizer of all cosets modulo $U$. Yet sometimes
larger groups may be possible. For example if we want to generate
start sets in $U$ and then do a brute force search. In this case, we
may take the stabilizer of $U$ for reduction.
%notest
\beginexample
gap> U:=First(sigdat,s->s.sigs=[ [ 4, 4, 8 ] ]).subgroup; 
Group([ f2, f3 ])
gap> cosets:=RightCosets(G,U);
gap> U1:=cosets[2];;U2:=cosets[3];;
gap> Up:=GroupList2PermList(Set(U),Gdata);;  
gap> Np:=GroupList2PermList(Set(N),Gdata); 
[ 1, 12, 25, 43, 78, 97, 115, 116, 134, 151, 169, 188, 207, 238, 249 ]
gap> comps:=Difference(Up,Np);; 
gap> ssets:=List(comps,i->[i]);;
gap> ssets:=ReducedStartsets(ssets,[Gdata.Aac],sigdat,Gdata.diffTable); 
#I  Size 80
#I  2/ 0 @ 0:00:00.061
[ [ 3 ], [ 4 ] ]
\endexample

Observe that $1$ is assumed to be element of every difference set and
is not recorded explicitly. So the partial difference sets represented
by $ssets$ at this point are {\tt[ [ 1, 3 ], [ 1, 4 ] ]}. Now the
startsets are extended to size $7$ using elements of $Up$. The runtime
varies depending on the output of "SignatureDataForNormalSubgroups"
and hence on $maxtest$.
%notest
\beginexample
gap> repeat 
>     ssets:=ExtendedStartsets(ssets,comps,Np,7,Gdata); 
>     ssets:=ReducedStartsets(ssets,[Gdata.Aac],sigdat,Gdata.diffTable);; 
> until ssets=[] or Size(ssets[1])=7; 
#I  Size 133
#I  3/ 0 @ 0:00:00.133
#I  Size 847
#I  4/ 0 @ 0:00:00.949
#I  Size 6309
#I  4/ 0 @ 0:00:07.692
#I  Size 21527
#I  5/ 0 @ 0:00:28.337
#I  Size 15884
#I  4/ 0 @ 0:00:21.837
#I  Size 1216
#I  4/ 0 @ 0:00:01.758
gap> Size(ssets);
192
\endexample
At a higher level of "InfoRDS", the number of start sets which are
discarded because of wrong signatures are also shown.
Now the cosets are changed. Here we use the `NoSort' version of
"ExtendedStartsets".  This leads to a lot of start sets in the first
step. If the number of start sets in $U$ is very large, this could be
too much for a reduction. Then the only option is using the brute
force method. But also for the brute force search,
"ExtendedStartSetsNoSort" must be called first (remember that we chos
an enumeration of $G$ and assume the last element from each startset
to be the largeset ``interesting'' one for further completions).

%notest
\beginexample
gap> comps:=Difference(GroupList2PermList(Set(U1),Gdata),Np);;
gap> ssets:=ExtendedStartsetsNoSort(ssets,comps,Np,11,Gdata);;
gap> ssets:=ReducedStartsets(ssets,[Gdata.Aac],sigdat,Gdata.diffTable);; 
#I  Size 8640
#I  9/ 0 @ 0:00:14.159
gap> Size(ssets);  
6899
\endexample
%
And as above, we continue:
%
\begintt
repeat 
    ssets:=ExtendedStartsets(ssets,comps,Np,11,Gdata); 
    ssets:=ReducedStartsets(ssets,[Gdata.Aac],sigdat,Gdata.diffTable);; 
until ssets=[] or Size(ssets[1])=11; 
comps:=Difference(GroupList2PermList(Set(U2),Gdata),Np); 
RaiseStartSetLevelNoSort(ssets,comps,Np,15,Gdata); 
repeat 
    ssets:=ExtendedStartsets(ssets,comps,Np,15,Gdata); 
    ssets:=ReducedStartsets(ssets,[Gdata.Aac],sigdat,Gdata.diffTable);; 
until ssets=[] or Size(ssets[1])=15; 
\endtt

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
%E