0

Is there an easier way to write the LCM part of this code? Gap only allows two arguments for LCM_INT.

SigmaSet := function(GROUP, g1, g2)
local S, g3;
S := [ ];
g3:= g2^-1 * g1^-1;
g_elts := [g1, g2, g3];
for g in GROUP do
    for j in [0..LCM_INT( Order(g1), LCM_INT( Order(g2), Order(g3) ) ) ] do
        for i in [1..3] do
            S := UnionSet(S, [g*g_elts[i]^j*g^-1 ]);
        od; 
    od; 
od;
return S;
end;
Olexandr Konovalov
  • 7,002
  • 2
  • 34
  • 72
ChrisG
  • 281

0 Answers0