From 76b7a02dc6767f4a7937b7da5686926c7b08680a Mon Sep 17 00:00:00 2001
From: Alexander Hulpke <hulpke@math.colostate.edu>
Date: Wed, 7 Dec 2016 11:39:50 -0700
Subject: [PATCH] FIX: Ensure all points are in domain if not single orbit

This fixes #990
---
 lib/grppcaut.gi              | 5 +++++
 tst/teststandard/grpauto.tst | 9 +++++++++
 2 files changed, 14 insertions(+)

diff --git a/lib/grppcaut.gi b/lib/grppcaut.gi
index 78ec806f24d..6165eb66110 100644
--- a/lib/grppcaut.gi
+++ b/lib/grppcaut.gi
@@ -1220,6 +1220,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
index 677c90085cc..87dcdfe30ef 100644
--- 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
 #