passagemath-gap-pkg-semigroups 10.6.30__cp311-cp311-macosx_13_0_arm64.whl

gap/pkg/semigroups/gap/elements/star.gd

@@ -0,0 +1,21 @@
+#############################################################################
+##
+##  elements/star.gd
+##  Copyright (C) 2015-2022                              James D. Mitchell
+##
+##  Licensing information can be found in the README file of this package.
+##
+#############################################################################
+##
+
+DeclareCategory("IsAssociativeElementWithStar", IsAssociativeElement);
+DeclareCategoryCollections("IsAssociativeElementWithStar");
+DeclareOperation("StarOp", [IsAssociativeElementWithStar]);
+DeclareAttribute("Star", IsAssociativeElementWithStar);
+
+DeclareProperty("IsStarSemigroup",
+                IsSemigroup and IsAssociativeElementWithStarCollection);
+DeclareSynonym("IsRegularStarSemigroup",
+               IsRegularSemigroup and IsStarSemigroup);
+
+DeclareOperation("InverseOp", [IsAssociativeElementWithStar]);

gap/pkg/semigroups/gap/elements/star.gi

@@ -0,0 +1,21 @@
+#############################################################################
+##
+##  elements/star.gi
+##  Copyright (C) 2015-2022                              James D. Mitchell
+##
+##  Licensing information can be found in the README file of this package.
+##
+#############################################################################
+##
+
+InstallMethod(Star, "for an associative element with star",
+[IsAssociativeElementWithStar],
+function(elm)
+  elm := StarOp(elm);
+  MakeImmutable(elm);
+  return elm;
+end);
+
+InstallMethod(IsStarSemigroup, "for a semigroup of elements with star",
+[IsSemigroup and IsAssociativeElementWithStarCollection],
+S -> ForAll(GeneratorsOfSemigroup(S), x -> Star(x) in S));

gap/pkg/semigroups/gap/elements/trans.gd

@@ -0,0 +1,13 @@
+############################################################################
+##
+##  elements/trans.gd
+##  Copyright (C) 2015-2022                              James D. Mitchell
+##
+##  Licensing information can be found in the README file of this package.
+##
+#############################################################################
+
+DeclareAttribute("CanonicalTransformation", IsTransformation);
+DeclareOperation("CanonicalTransformation", [IsTransformation, IsInt]);
+DeclareOperation("TransformationByImageAndKernel",
+                 [IsHomogeneousList, IsCyclotomicCollColl]);

gap/pkg/semigroups/gap/elements/trans.gi

@@ -0,0 +1,50 @@
+############################################################################
+##
+##  elements/trans.gi
+##  Copyright (C) 2015-2022                              James D. Mitchell
+##
+##  Licensing information can be found in the README file of this package.
+##
+#############################################################################
+
+InstallMethod(CanonicalTransformation, "for a transformation",
+[IsTransformation], f -> CanonicalTransformation(f, DegreeOfTransformation(f)));
+
+InstallMethod(CanonicalTransformation, "for a transformation",
+[IsTransformation, IsInt],
+function(f, n)
+  local digraph;
+  if n < DegreeOfTransformation(f) then
+    ErrorNoReturn("the second argument (an integer) must be at least ",
+                  "the degree of the first argument (a transformation)");
+  fi;
+
+  digraph := AsDigraph(f, n);
+  return AsTransformation(OnDigraphs(digraph,
+                                     BlissCanonicalLabelling(digraph)));
+end);
+
+InstallMethod(TransformationByImageAndKernel, "for an image and partition",
+[IsHomogeneousList, IsCyclotomicCollColl],
+function(im, ker)
+  local flat, i;
+
+  if not ForAll(ker, class -> ForAll(class, IsPosInt)) then
+    ErrorNoReturn("the argument must be a list of lists of pos ints");
+  fi;
+
+  flat := Union(ker);
+
+  if flat <> [1 .. Length(flat)] then
+    ErrorNoReturn("the union of the second argument (a partition) must ",
+                  "be [1 .. ", Length(flat), "]");
+  fi;
+
+  for i in [1 .. Length(ker)] do
+    flat{ker[i]} := [1 .. Length(ker[i])] * 0 + i;
+  od;
+  return TransformationByImageAndKernel(im, flat);
+end);
+
+InstallMethod(IndexPeriodOfSemigroupElement, "for a transformation",
+[IsTransformation], IndexPeriodOfTransformation);

gap/pkg/semigroups/gap/fp/freeband.gd

@@ -0,0 +1,22 @@
+###############################################################################
+##
+##  fp/freeband.gd
+##  Copyright (C) 2013-2022                                 Julius Jonusas
+##
+##  Licensing information can be found in the README file of this package.
+##
+###############################################################################
+
+DeclareCategory("IsFreeBandElement", IsAssociativeElement);
+DeclareCategoryCollections("IsFreeBandElement");
+DeclareCategory("IsFreeBandCategory", IsSemigroup and IsBand);
+DeclareProperty("IsFreeBand", IsSemigroup);
+DeclareSynonym("IsFreeBandSubsemigroup",
+               IsSemigroup and IsFreeBandElementCollection);
+
+DeclareGlobalFunction("FreeBand");
+DeclareAttribute("ContentOfFreeBandElement", IsFreeBandElement);
+DeclareAttribute("ContentOfFreeBandElementCollection",
+                 IsFreeBandElementCollection);
+
+DeclareOperation("EqualInFreeBand", [IsList, IsList]);

gap/pkg/semigroups/gap/fp/freeband.gi

@@ -0,0 +1,502 @@
+###############################################################################
+##
+##  fp/freeband.gi
+##  Copyright (C) 2013-2022                                Julius Jonusas
+##
+##  Licensing information can be found in the README file of this package.
+##
+###############################################################################
+
+# TODO(later): this is not really finished.
+
+# TODO(later)
+# InstallMethod(FreeBandOfFreeBandElement,
+
+InstallMethod(ContentOfFreeBandElement, "for a free band element",
+[IsFreeBandElement],
+w -> ListBlist([1 .. Length(w!.cont)], w!.cont));
+
+InstallMethod(ContentOfFreeBandElementCollection,
+"for a free band element collection",
+[IsFreeBandElementCollection],
+function(coll)
+  local n, content, w;
+
+  n := Length(coll[1]!.cont);
+  content := BlistList([1 .. n], []);
+
+  for w in coll do
+    UniteBlist(content, w!.cont);
+    if SizeBlist(content) = n then
+      break;
+    fi;
+  od;
+
+  return ListBlist([1 .. n], content);
+end);
+
+###############################################################################
+## Internal
+###############################################################################
+
+SEMIGROUPS.FreeBandElmToWord := function(elem)
+  local tuple, out, first, pre_tuple, last, su_tuple, word1, word2;
+
+  tuple := elem!.tuple;
+
+  if elem!.word <> fail then
+    return elem!.word;
+  elif tuple![2] = 0 then  # tuple corresponds to one the generators
+    out := [tuple![1]];
+  else
+    first := tuple![1];
+    pre_tuple := tuple![2];  # tuple corresponding to the prefix
+    last := tuple![3];
+    su_tuple := tuple![4];  # tuple corresponding to the suffix
+
+    # if first = last we only need a single letter
+    if first = last then
+      out := Concatenation(SEMIGROUPS.FreeBandElmToWord(pre_tuple), [first],
+                           SEMIGROUPS.FreeBandElmToWord(su_tuple));
+      # otherwise we need distinct letters for both first and last
+    else
+      word1 := Concatenation(SEMIGROUPS.FreeBandElmToWord(pre_tuple), [first]);
+      word2 := Concatenation([last], SEMIGROUPS.FreeBandElmToWord(su_tuple));
+      if word1 = word2 then
+        out := word1;
+      else
+        out := Concatenation(word1, word2);
+      fi;
+    fi;
+  fi;
+  elem!.word := out;
+  return out;
+end;
+
+SEMIGROUPS.HashFunctionForFreeBandElements := function(x, data)
+  return ORB_HashFunctionForPlainFlatList(SEMIGROUPS.FreeBandElmToWord(x),
+                                          data);
+end;
+
+InstallMethod(IsGeneratorsOfInverseSemigroup, "for a free band element coll",
+[IsFreeBandElementCollection], ReturnFalse);
+
+InstallTrueMethod(IsFinite, IsFreeBandSubsemigroup);
+
+InstallGlobalFunction(FreeBand,
+function(arg...)
+  local names, F, type, ngens, gens, filts, S, m;
+
+  # Get and check the argument list, and construct names if necessary.
+  if Length(arg) = 1 and IsInt(arg[1]) and 0 < arg[1] then
+    names := List([1 .. arg[1]],
+                  i -> Concatenation("x", String(i)));
+  elif Length(arg) = 2 and IsInt(arg[1]) and 0 < arg[1]
+      and IsString(arg[2]) then
+    names := List([1 .. arg[1]],
+                  i -> Concatenation(arg[2], String(i)));
+  elif 1 <= Length(arg) and ForAll(arg, IsString) then
+    names := arg;
+  elif Length(arg) = 1 and IsList(arg[1])
+      and ForAll(arg[1], IsString) then
+    names := arg[1];
+  else
+    ErrorNoReturn("FreeBand(<name1>,<name2>..) or FreeBand(<rank> [, name])");
+  fi;
+
+  MakeImmutable(names);
+
+  F := NewFamily("FreeBandElementsFamily", IsFreeBandElement,
+                                           CanEasilySortElements);
+
+  type := NewType(F, IsFreeBandElement and IsPositionalObjectRep);
+
+  ngens := Length(names);
+  gens := EmptyPlist(ngens);
+
+  for m in [1 .. ngens] do
+    gens[m] := Objectify(type, rec(tuple := [m, 0, m, 0],
+                                   cont := BlistList([1 .. ngens], [m]),
+                                   word := [m]));
+  od;
+
+  StoreInfoFreeMagma(F, names, IsFreeBandElement);
+
+  filts := IsFreeBandCategory and IsAttributeStoringRep and IsWholeFamily and
+           IsFreeBand;
+
+  S := Objectify(NewType(FamilyObj(gens), filts),
+                 rec(opts := SEMIGROUPS.DefaultOptionsRec));
+
+  SetGeneratorsOfMagma(S, gens);
+
+  FamilyObj(S)!.semigroup := S;
+  F!.semigroup := S;
+
+  return S;
+end);
+
+InstallMethod(IsFreeBand, "for a semigroup", [IsSemigroup],
+function(S)
+  local used, occurred, gens, max_d, g;
+
+  if not IsBand(S) then
+    return false;
+  elif IsFreeBandSubsemigroup(S) then
+    gens := Generators(S);
+    used := BlistList([1 .. Length(gens[1]!.cont)], []);
+    occurred := BlistList([1 .. Length(gens[1]!.cont)], []);
+    for g in gens do
+      used := IntersectionBlist(used, g!.cont);
+      if g!.tuple[2] = 0 then
+        occurred[g!.tuple[1]] := true;
+      fi;
+    od;
+    if used = occurred then
+      return true;
+    fi;
+  fi;
+
+  if not IsActingSemigroup(S) then
+    S := AsSemigroup(IsTransformationSemigroup, S);
+  fi;
+
+  max_d := MaximalDClasses(S);
+  return Size(S) = Size(FreeBand(Length(max_d)));
+end);
+
+InstallMethod(Iterator, "for a Greens D-class of a free band",
+[IsFreeBandElementCollection and IsGreensDClass],
+function(dclass)
+  local NextIterator_FreeBandDClass, NewIterator_FreeBandDClass,
+  ShallowCopyLocal, record, s, content, rep,
+  NextIterator_FreeBandDClassWithPrint;
+
+  s := Parent(dclass);
+  rep := Representative(dclass);
+  content := rep!.cont;
+
+  NextIterator_FreeBandDClass := function(iter)
+    local output, i, content, tempcont, tuple;
+
+    if iter!.element <> fail then
+      # output := StructuralCopy(iter!.element); doesn't work
+      output := rec(tuple := StructuralCopy(iter!.element!.tuple),
+                    content := ShallowCopy(iter!.element!.content),
+                    word := ShallowCopy(iter!.element!.word));
+    else
+      return fail;
+    fi;
+
+    content := iter!.content;
+    tuple := iter!.element!.tuple;
+    iter!.element!.word := fail;
+
+    if tuple[2] = 0 then
+      iter!.element := fail;
+    elif iter!.iter1!.element <> fail then
+      # Prefix word is not done yet
+      tuple[2] := NextIterator_FreeBandDClass(iter!.iter1);
+      iter!.element!.tuple := tuple;
+    elif Position(content, true, tuple[1]) <> fail then
+      # Update the first component
+      i := Position(content, true, tuple[1]);
+      tuple[1] := i;
+      tempcont := ShallowCopy(content);
+      tempcont[i] := false;
+      iter!.iter1 := NewIterator_FreeBandDClass(iter!.semigroup, tempcont);
+      tuple[2] := NextIterator_FreeBandDClass(iter!.iter1);
+      iter!.element!.tuple := tuple;
+    elif iter!.iter2!.element <> fail then
+      # Suffix word is not done yet
+      tuple[4] := NextIterator_FreeBandDClass(iter!.iter2);
+      iter!.element!.tuple := tuple;
+      # Restart the prefix
+      i := Position(content, true);
+      tuple[1] := i;
+      tempcont := ShallowCopy(content);
+      tempcont[i] := false;
+      iter!.iter1 := NewIterator_FreeBandDClass(iter!.semigroup, tempcont);
+      tuple[2] := NextIterator_FreeBandDClass(iter!.iter1);
+      iter!.element!.tuple := tuple;
+    elif Position(content, true, tuple[3]) <> fail then
+      # Update the third component
+      i := Position(content, true, tuple[3]);
+      tuple[3] := i;
+      tempcont := ShallowCopy(content);
+      tempcont[i] := false;
+      iter!.iter2 := NewIterator_FreeBandDClass(iter!.semigroup, tempcont);
+      tuple[4] := NextIterator_FreeBandDClass(iter!.iter2);
+      # Restart the prefix
+      i := Position(content, true);
+      tuple[1] := i;
+      tempcont := ShallowCopy(content);
+      tempcont[i] := false;
+      iter!.iter1 := NewIterator_FreeBandDClass(iter!.semigroup, tempcont);
+      tuple[2] := NextIterator_FreeBandDClass(iter!.iter1);
+      iter!.element!.tuple := tuple;
+    else
+      iter!.element := fail;
+    fi;
+    return output;
+  end;
+
+  NextIterator_FreeBandDClassWithPrint := function(iter)
+    local next_value;
+
+    next_value := NextIterator_FreeBandDClass(iter);
+
+    if next_value = fail then
+      return fail;
+    else
+      return Product(List(SEMIGROUPS.FreeBandElmToWord(next_value),
+                     x -> GeneratorsOfSemigroup(iter!.semigroup)[x]));
+    fi;
+  end;
+
+  NewIterator_FreeBandDClass := function(s, content)
+    local record, first, tempcont, elem;
+
+    first := Position(content, true);
+    elem := Objectify(TypeObj(s.1), rec(tuple := [],
+                                        content := content,
+                                        word := fail));
+
+    # If the content is of size one
+    if Position(content, true, first) = fail then
+      elem!.tuple := [first, 0, first, 0];
+      elem!.word := [first];
+      record := rec(element := elem,
+                    iter1 := fail,
+                    iter2 := fail,
+                    semigroup := s,
+                    content := content) ;
+    else
+      tempcont := ShallowCopy(content);
+      tempcont[first] := false;
+      record := rec(element := elem,
+                    iter1 := NewIterator_FreeBandDClass(s, tempcont),
+                    iter2 := NewIterator_FreeBandDClass(s, tempcont),
+                    semigroup := s,
+                    content := content);
+      record!.element!.tuple := [first,
+                                 NextIterator_FreeBandDClass(record!.iter1),
+                                 first,
+                                 NextIterator_FreeBandDClass(record!.iter2)];
+    fi;
+    return record;
+  end;
+
+  ShallowCopyLocal := function(it)
+    local record;
+    record := NewIterator_FreeBandDClass(s, content);
+    record.NextIterator := it!.NextIterator;
+    record.ShallowCopy := it!.ShallowCopy;
+    return record;
+  end;
+
+  record := NewIterator_FreeBandDClass(s, content);
+  record.NextIterator := NextIterator_FreeBandDClassWithPrint;
+  record.ShallowCopy := ShallowCopyLocal;
+  return IteratorByNextIterator(record);
+end);
+
+InstallMethod(GreensDClassOfElement, "for a free band and element",
+[IsFreeBandCategory, IsFreeBandElement],
+function(S, x)
+  local filt, type, D;
+  if not x in S then
+    ErrorNoReturn("the 2nd argument (a free band element) does not ",
+                  "belong to 1st argument (a free band category)");
+  fi;
+
+  filt := IsEquivalenceClass and IsEquivalenceClassDefaultRep
+          and IsGreensDClass;
+
+  if Length(GeneratorsOfSemigroup(S)) < 4 then
+    filt := filt and IsGreensClassOfSemigroupThatCanUseFroidurePinRep;
+  fi;
+
+  type := NewType(FamilyObj(S), filt);
+  D := Objectify(type, rec());
+  SetParent(D, S);
+  SetRepresentative(D, x);
+  SetEquivalenceClassRelation(D, GreensDRelation(S));
+
+  return D;
+end);
+
+InstallMethod(ViewString, "for a free band element",
+[IsFreeBandElement], PrintString);
+
+InstallMethod(PrintString, "for a free band element",
+[IsFreeBandElement],
+function(elem)
+  return Concatenation(List(SEMIGROUPS.FreeBandElmToWord(elem),
+                            x -> FamilyObj(elem)!.names[x]));
+end);
+
+InstallMethod(ViewObj, "for a free band",
+[IsFreeBandCategory],
+function(S)
+  if GAPInfo.ViewLength * 10 < Length(GeneratorsOfMagma(S)) then
+    Print("<free band with ", Length(GeneratorsOfSemigroup(S)),
+          " generators>");
+  else
+    Print("<free band on the generators ",
+          GeneratorsOfSemigroup(S), ">");
+  fi;
+end);
+
+InstallMethod(EqualInFreeBand, "for two strings",
+[IsString, IsString],
+function(w1_in, w2_in)
+  if IsEmpty(w1_in) or IsEmpty(w2_in) then
+    return IsEmpty(w1_in) and IsEmpty(w2_in);
+    # Necessary since [] is recognised as a string yet the below
+    # doesn't work with it.
+  fi;
+  w1_in := List(w1_in, IntChar);
+  w2_in := List(w2_in, IntChar);
+  return EqualInFreeBand(w1_in, w2_in);
+end);
+
+# This is an implementation of the algorithm described in Jakub Radoszewski
+# and Wojciech Rytter's paper 'Efficient Testing of Equivalence of Words
+# in a Free Idempotent Semigroup⋆'
+
+InstallMethod(EqualInFreeBand, "for two lists of positive integers",
+[IsHomogeneousList, IsHomogeneousList],
+{w1_in, w2_in} -> libsemigroups.freeband_equal_to(w1_in, w2_in));
+
+# TODO(later) Is there a more efficient way to compare elements? JJ
+
+InstallMethod(\=, "for elements of a free band",
+IsIdenticalObj, [IsFreeBandElement, IsFreeBandElement],
+function(x, y)
+  local i;
+
+  if x!.cont <> y!.cont then
+    return false;
+  fi;
+  for i in [1 .. 4] do
+    if x!.tuple![i] <> y!.tuple![i] then
+      return false;
+    fi;
+  od;
+  return true;
+end);
+
+# TODO(later) Is it possible to find a non recursive way to compare elements?
+# JJ
+
+InstallMethod(\<, "for elements of a free band",
+IsIdenticalObj, [IsFreeBandElement, IsFreeBandElement],
+function(x, y)
+  local i, tuple1, tuple2;
+
+  tuple1 := x!.tuple;
+  tuple2 := y!.tuple;
+
+  i := 1;
+  while i <= 3 and tuple1![i] = tuple2![i] do
+    i := i + 1;
+  od;
+
+  if tuple2[i] = 0  then
+    return false;
+  elif tuple1[i] = 0 then
+    return true;
+  fi;
+  return tuple1[i] < tuple2[i];
+end);
+
+InstallMethod(\*, "for elements of a free band", IsIdenticalObj,
+[IsFreeBandElement, IsFreeBandElement],
+function(x, y)
+  local type, out, diff, new_first, new_prefix, new_last, new_sufix, copy;
+
+  type := TypeObj(x);
+  # if the content of two elements is the same we only need the prefix of the
+  # first and the suffix of the second one
+  # cont = blist
+  if IsSubsetBlist(x!.cont, y!.cont) then
+    out := [x!.tuple[1], x!.tuple[2]];
+  else
+    diff := DifferenceBlist(y!.cont, x!.cont);
+    # new_first is the last letter to occur first in the product
+    new_first := y!.tuple[1];
+    new_prefix := y!.tuple[2];
+    while true do
+      if diff[new_first] and new_prefix = 0 then
+        copy := ShallowCopy(x);  # are shallow copies necessary?
+        out := [new_first, copy];
+        break;
+      elif diff[new_first] then
+        copy := ShallowCopy(x * new_prefix);
+        out := [new_first, copy];
+        break;
+      else
+        new_first := new_prefix!.tuple[1];
+        new_prefix := new_prefix!.tuple[2];
+      fi;
+    od;
+  fi;
+
+  if IsSubsetBlist(y!.cont, x!.cont) then
+    out{[3, 4]} := [y!.tuple[3], y!.tuple[4]];
+  else
+    diff := DifferenceBlist(x!.cont, y!.cont);
+    # new_last is the first letter to occur for the last time in the product
+    new_last := x!.tuple[3];
+    new_sufix := x!.tuple[4];
+    repeat
+      if diff[new_last] and new_sufix = 0 then
+        copy := ShallowCopy(y);
+        out{[3 .. 4]} := [new_last, copy];
+        break;
+      elif diff[new_last] then
+        copy := ShallowCopy(new_sufix * y);
+        out{[3, 4]} := [new_last, copy];
+        break;
+      else
+        new_last := new_sufix!.tuple[3];
+        new_sufix := new_sufix!.tuple[4];
+      fi;
+    until IsBound(out[3]);
+  fi;
+
+  return Objectify(type, rec(tuple := out,
+                             cont := UnionBlist(x!.cont, y!.cont),
+                             word := fail));
+end);
+
+InstallMethod(Size, "for a free band",
+[IsFreeBandCategory and IsFinite and HasGeneratorsOfSemigroup],
+function(S)
+  local c, output, k, i, n;
+
+  n := Length(FamilyObj(S.1)!.names);
+
+  c := [];
+  for k in [1 .. n] do
+    c[k] := 1;
+    for i in [1 .. k - 1] do
+      c[k] := c[k] * (k - i + 1) ^ (2 ^ i);
+    od;
+  od;
+
+  output := 0;
+  for k in [1 .. n] do
+    output := output + Binomial(n, k) * c[k];
+  od;
+
+  return output;
+end);
+
+InstallMethod(ChooseHashFunction, "for a free band element and int",
+[IsFreeBandElement, IsInt],
+function(_, hashlen)
+  return rec(func := SEMIGROUPS.HashFunctionForFreeBandElements,
+             data := hashlen);
+end);

gap/pkg/semigroups/gap/fp/freeinverse.gd

@@ -0,0 +1,30 @@
+###############################################################################
+##
+##  fp/freeinverse.gd
+##  Copyright (C) 2013-2022                                Julius Jonusas
+##
+##  Licensing information can be found in the README file of this package.
+##
+###############################################################################
+
+DeclareUserPreference(rec(
+  name := "FreeInverseSemigroupElementDisplay",
+  description := ["options for the display of free inverse semigroup elements"],
+  default := "minimal",
+  values := ["minimal", "canonical"],
+  multi := false,
+  package := "semigroups"));
+
+DeclareCategory("IsFreeInverseSemigroupElement", IsAssociativeElement);
+DeclareCategoryCollections("IsFreeInverseSemigroupElement");
+
+DeclareCategory("IsFreeInverseSemigroupCategory",
+                IsInverseSemigroup
+                and IsGeneratorsOfInverseSemigroup
+                and IsFreeInverseSemigroupElementCollection);
+DeclareProperty("IsFreeInverseSemigroup", IsSemigroup);
+
+DeclareGlobalFunction("FreeInverseSemigroup");
+
+DeclareAttribute("MinimalWord", IsFreeInverseSemigroupElement);
+DeclareAttribute("CanonicalForm", IsFreeInverseSemigroupElement);