FIX: Ensure all points are in domain if not single orbit

This fixes #990
gap-system · Dec 7, 2016 · dcd687d · dcd687d
1 parent 041efa0
commit dcd687d
Show file tree

Hide file tree

Showing 2 changed files with 14 additions and 0 deletions.
diff --git a/lib/grppcaut.gi b/lib/grppcaut.gi
@@ -1218,6 +1218,11 @@ local spaceincl,outvecs,l,sub,yet,i,j,k,s,t,new,incl,min,rans,sofar,done,
  for i in osporb do
  # assumption: orbit is not too long... (i.e. TODO: improve)
  m:=Orbit(a,i[1],OnSubspacesByCanonicalBasis);
+ for yet in i do
+ if not yet in m then
+ m:=Union(m,Orbit(a,yet,OnSubspacesByCanonicalBasis));
+ fi;
+ od;
  if Length(m)>Length(i) then
  yet:=ActionHomomorphism(a,m,OnSubspacesByCanonicalBasis,"surjective");
  sub:=Stabilizer(Image(yet),Set(List(i,x->Position(m,x))),OnSets);

diff --git a/tst/teststandard/grpauto.tst b/tst/teststandard/grpauto.tst
@@ -11,6 +11,15 @@
 
 gap> START_TEST("grpauto.tst");
 
+gap> G:=PcGroupCode(2664057732308087629446918084678644731357911420184737672,
+> 6272);;
+gap> H:=PcGroupCode(849169324138369034371257796152864931485540637536648, 
+> 6272);;
+gap> IsomorphismGroups(G,H);
+fail
+gap> IsomorphismGroups(G,PcGroupCode(CodePcGroup(G),Size(G)))=fail;
+false
+
 #
 # hard-iso
 #