Skip to content
Open
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
24 changes: 20 additions & 4 deletions lib/InduceReduce.gi
Original file line number Diff line number Diff line change
Expand Up @@ -509,7 +509,7 @@ InstallValue( IndRed , rec(
##
InstallGlobalFunction( InduceReduce,
function(GR,Opt)
local TR, RedTR, H, ccsizesH, temp, it;
local TR, RedTR, H, ccsizesH, temp, it, t, tElem, tInduce, tReduce;

TR:=IndRed.GroupTools; # get group tools and reduce tools
RedTR:=IndRed.ReduceTools;
Expand All @@ -529,27 +529,43 @@ local TR, RedTR, H, ccsizesH, temp, it;
return GR.Ir;
fi;


tElem := 0;
tInduce := 0;
tReduce := 0;
while Size(GR.Ir)<GR.k and not (GR.NumberOfPrimes=0 and GR.NumberOfCyclic=0) do
# if number of irr. characters=number of conjugacy classes,
# all irr. characters have been found.

Opt.LLLOffset:=Opt.LLLOffset-1;
# Opt.LLLOffset postpones the first LLL lattice reduction

t := NanosecondsSinceEpoch();

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Just for curiosity: Why not Runtime()?

IndRed.FindElementary(GR,Opt); # find elementary subgroup

IndRed.InitElementary(GR,TR); # determine information needed about elementary subgroup

t := NanosecondsSinceEpoch() - t;
Info(InfoCTUnger + InfoTiming, 2, "InduceReduce timing for Elementary: ", t/1000000., " ms");
tElem := tElem + t;

Info(InfoCTUnger, 1, "Induce/Restrict: Trying [|Z|, |P|, k(E)] = ",
[ GR.orders[GR.IndexCyc],
GR.Elementary.n/GR.orders[GR.IndexCyc], GR.Elementary.k ]);

t := NanosecondsSinceEpoch();
IndRed.Induce(GR); # append induced characters to GR.B
t := NanosecondsSinceEpoch() - t;
Info(InfoCTUnger + InfoTiming, 2, "InduceReduce timing for Induce: ", t/1000000., " ms");
tInduce := tInduce + t;

t := NanosecondsSinceEpoch();
IndRed.Reduce(GR,RedTR,Opt); # reduce GR.B by GR.Ir and do lattice reduction
t := NanosecondsSinceEpoch() - t;
Info(InfoCTUnger + InfoTiming, 2, "InduceReduce timing for Reduce: ", t/1000000., " ms");
tReduce := tReduce + t;
od;
Info(InfoCTUnger + InfoTiming, 2, "InduceReduce total: ", (tElem + tInduce + tReduce)/1000000., " ms");
Info(InfoCTUnger + InfoTiming, 2, " Elementary: ", tElem/1000000., " ms");
Info(InfoCTUnger + InfoTiming, 2, " Induce: ", tInduce/1000000., " ms");
Info(InfoCTUnger + InfoTiming, 2, " Reduce: ", tReduce/1000000., " ms");

if Size(GR.Ir)>=GR.k then
return GR.Ir;
Expand Down