Attempt to fix bug caused by LINS_MaxPGenerators being too low.

Unfortunately the method `PQuotient` ungracefully throws an error, instead of returning fail.

Author: FriedrichRober

gap/findPQuotients.gi

@@ -176,13 +176,6 @@ BindGlobal("LINS_PullBackH", function(GenM, p, Gens, O, Mu, Psi)
 end);
 
 
-#############################################################################
-## maximal generators of a p-quotient.
-#############################################################################
-
-BindGlobal("LINS_MaxPGenerators", 1000);
-
-
 #############################################################################
 ##  LINS_FindPModules
 #############################################################################
@@ -225,6 +218,8 @@ InstallGlobalFunction(LINS_FindPModules, function(gr, rH, p, opts)
     n,      # pos-int:      index bound of LINS graph `gr`.
     Iso,    # isomorphism:  from `H` into fp-group
     IH,     # fp-group:     image under `Iso`, isomorphic to `H`
+    logord, # int:          option for PQuotient, the maximal order
+            #               for P we consider is p^logord
     P,      # quotsys:      p-quotient of `IH` of class 1
     Mu,     # epimorphism:  from `IH` into `P`
     M,      # pc-group:     image under `Mu`
@@ -272,7 +267,8 @@ InstallGlobalFunction(LINS_FindPModules, function(gr, rH, p, opts)
     IH := Image(Iso);
 
     # Create the Isomorphism to the group structure of the `p`-Module `M`
-    P := PQuotient(IH, p, 1, LINS_MaxPGenerators);
+    logord := Int(Log(Float(n / rH!.Index))/Log(Float(p)));
+    P := PQuotient(IH, p, 1, logord);
     Mu := EpimorphismQuotientSystem(P);
     M := Image(Mu);
     GenM := GeneratorsOfGroup(M);