|
#############################################################################
##
## isomorph.gi
## Copyright (C) 2014-19 James D. Mitchell
## Wilf A. Wilson
##
## Licensing information can be found in the README file of this package.
##
#############################################################################
##
# This file contains methods using bliss or nauty, for computing isomorphisms,
# automorphisms, and canonical labellings of digraphs.
InstallGlobalFunction(DigraphsUseBliss,
function()
# Just do nothing if NautyTracesInterface is not available.
if DIGRAPHS_NautyAvailable then
Info(InfoWarning,
1,
"Using bliss by default for AutomorphismGroup . . .");
if not DIGRAPHS_UsingBliss then
InstallMethod(AutomorphismGroup, "for a digraph", [IsDigraph],
BlissAutomorphismGroup);
InstallMethod(AutomorphismGroup, "for a digraph and vertex coloring",
[IsDigraph, IsHomogeneousList], BlissAutomorphismGroup);
MakeReadWriteGlobal("DIGRAPHS_UsingBliss");
DIGRAPHS_UsingBliss := true;
MakeReadOnlyGlobal("DIGRAPHS_UsingBliss");
fi;
fi;
end);
InstallGlobalFunction(DigraphsUseNauty,
function()
if DIGRAPHS_NautyAvailable then
if DIGRAPHS_UsingBliss then
InstallMethod(AutomorphismGroup, "for a digraph", [IsDigraph],
NautyAutomorphismGroup);
InstallMethod(AutomorphismGroup, "for a digraph and vertex coloring",
[IsDigraph, IsHomogeneousList], NautyAutomorphismGroup);
MakeReadWriteGlobal("DIGRAPHS_UsingBliss");
DIGRAPHS_UsingBliss := false;
MakeReadOnlyGlobal("DIGRAPHS_UsingBliss");
fi;
Info(InfoWarning,
1,
"Using NautyTracesInterface by default for AutomorphismGroup");
Info(InfoWarning,
1,
"bliss will be used for edge coloured automorphisms");
else
Info(InfoWarning,
1,
"NautyTracesInterface is not available!");
Info(InfoWarning,
1,
"Using bliss by default for AutomorphismGroup . . .");
fi;
end);
# Wrappers for the C-level functions
BindGlobal("BLISS_DATA_NC",
function(digraph, vert_colours, edge_colours)
local collapsed, mults, data, edge_gp;
if IsMultiDigraph(digraph) then
if edge_colours = fail then
collapsed := DIGRAPHS_CollapseMultipleEdges(digraph);
else
collapsed := DIGRAPHS_CollapseMultiColouredEdges(digraph, edge_colours);
fi;
digraph := collapsed[1];
edge_colours := collapsed[2];
mults := collapsed[3];
data := DIGRAPH_AUTOMORPHISMS(digraph,
vert_colours,
edge_colours);
if IsEmpty(data[1]) then
data[1] := [()];
fi;
data[1] := Group(data[1]);
if IsBound(data[3]) then
SetSize(data[1], data[3]);
fi;
if Length(mults) > 0 then
edge_gp := Group(Flat(List(mults,
x -> GeneratorsOfGroup(SymmetricGroup(x)))));
data[1] := DirectProduct(data[1], edge_gp);
else
data[1] := DirectProduct(data[1], Group(()));
fi;
data[2] := [data[2], ()];
return data;
else
data := DIGRAPH_AUTOMORPHISMS(digraph,
vert_colours,
edge_colours);
if IsEmpty(data[1]) then
data[1] := [()];
fi;
data[1] := Group(data[1]);
if IsBound(data[3]) then
SetSize(data[1], data[3]);
fi;
return data;
fi;
end);
## The argument <vert_colours> should be a list of colours of the vertices
## of <digraph>, and the argument <edge_colours> should be a list of
## lists of the same shape as OutNeighbours(<digraph>).
## See ValidateVertexColouring and ValidateEdgeColouring.
##
## Returns a list where the first position is the automorphism group, and the
## second is the canonical labelling.
BindGlobal("BLISS_DATA",
function(D, vert_colours, edge_colours)
if vert_colours <> fail then
vert_colours := DIGRAPHS_ValidateVertexColouring(DigraphNrVertices(D),
vert_colours);
fi;
if edge_colours <> fail then
DIGRAPHS_ValidateEdgeColouring(D, edge_colours);
fi;
return BLISS_DATA_NC(D, vert_colours, edge_colours);
end);
BindGlobal("BLISS_DATA_NO_COLORS", D -> BLISS_DATA(D, fail, fail));
if DIGRAPHS_NautyAvailable then
BindGlobal("NAUTY_DATA",
function(D, colors)
local data;
if colors <> false then
colors := DIGRAPHS_ValidateVertexColouring(DigraphNrVertices(D),
colors);
colors := NautyColorData(colors);
fi;
if DigraphHasNoVertices(D) then
# This circumvents Issue #17 in NautyTracesInterface, whereby a graph
# with 0 vertices causes a seg fault.
return [Group(()), ()];
fi;
data := NautyDense(DigraphSource(D),
DigraphRange(D),
DigraphNrVertices(D),
not IsSymmetricDigraph(D),
colors);
if IsEmpty(data[1]) then
data[1] := [()];
fi;
data[1] := Group(data[1]);
data[2] := data[2] ^ -1;
return data;
end);
BindGlobal("NAUTY_DATA_NO_COLORS", D -> NAUTY_DATA(D, false));
else
BindGlobal("NAUTY_DATA", ReturnFail);
BindGlobal("NAUTY_DATA_NO_COLORS", ReturnFail);
fi;
# Canonical labellings
InstallMethod(BlissCanonicalLabelling, "for a digraph",
[IsDigraph],
function(D)
local data;
data := BLISS_DATA_NO_COLORS(D);
if IsImmutableDigraph(D) then
SetBlissAutomorphismGroup(D, data[1]);
fi;
return data[2];
end);
InstallMethod(BlissCanonicalLabelling, "for a digraph and vertex coloring",
[IsDigraph, IsHomogeneousList],
{D, colors} -> BLISS_DATA(D, colors, fail)[2]);
InstallMethod(NautyCanonicalLabelling, "for a digraph",
[IsDigraph],
function(D)
local data;
if not DIGRAPHS_NautyAvailable or IsMultiDigraph(D) then
Info(InfoWarning, 1, "NautyTracesInterface is not available");
return fail;
fi;
data := NAUTY_DATA_NO_COLORS(D);
if IsImmutableDigraph(D) then
SetNautyAutomorphismGroup(D, data[1]);
fi;
return data[2];
end);
InstallMethod(NautyCanonicalLabelling,
"for a digraph and vertex coloring",
[IsDigraph, IsHomogeneousList],
function(D, colors)
if not DIGRAPHS_NautyAvailable or IsMultiDigraph(D) then
Info(InfoWarning, 1, "NautyTracesInterface is not available");
return fail;
fi;
return NAUTY_DATA(D, colors)[2];
end);
# Canonical digraphs
InstallMethodThatReturnsDigraph(BlissCanonicalDigraph, "for a digraph",
[IsDigraph],
function(D)
if IsMultiDigraph(D) then
return OnMultiDigraphs(D, BlissCanonicalLabelling(D));
fi;
return OnDigraphsNC(D, BlissCanonicalLabelling(D));
end);
InstallMethod(BlissCanonicalDigraph, "for a digraph and vertex coloring",
[IsDigraph, IsHomogeneousList],
function(D, colors)
if IsMultiDigraph(D) then
return OnMultiDigraphs(D, BlissCanonicalLabelling(D, colors));
fi;
return OnDigraphsNC(D, BlissCanonicalLabelling(D, colors));
end);
InstallMethodThatReturnsDigraph(NautyCanonicalDigraph, "for a digraph",
[IsDigraph],
function(D)
if not DIGRAPHS_NautyAvailable or IsMultiDigraph(D) then
Info(InfoWarning, 1, "NautyTracesInterface is not available");
return fail;
fi;
return OnDigraphsNC(D, NautyCanonicalLabelling(D));
end);
InstallMethod(NautyCanonicalDigraph, "for a digraph and vertex coloring",
[IsDigraph, IsHomogeneousList],
function(D, colors)
if not DIGRAPHS_NautyAvailable or IsMultiDigraph(D) then
Info(InfoWarning, 1, "NautyTracesInterface is not available");
return fail;
fi;
return OnDigraphsNC(D, NautyCanonicalLabelling(D, colors));
end);
# Automorphism group
InstallMethod(BlissAutomorphismGroup, "for a digraph", [IsDigraph],
function(D)
local data;
data := BLISS_DATA_NO_COLORS(D);
if IsImmutableDigraph(D) then
SetBlissCanonicalLabelling(D, data[2]);
if not HasDigraphGroup(D) then
if IsMultiDigraph(D) then
SetDigraphGroup(D, Range(Projection(data[1], 1)));
else
SetDigraphGroup(D, data[1]);
fi;
fi;
fi;
return data[1];
end);
InstallMethod(NautyAutomorphismGroup, "for a digraph", [IsDigraph],
function(D)
local data;
if not DIGRAPHS_NautyAvailable or IsMultiDigraph(D) then
Info(InfoWarning, 1, "NautyTracesInterface is not available");
return fail;
fi;
data := NAUTY_DATA_NO_COLORS(D);
if IsImmutableDigraph(D) then
SetNautyCanonicalLabelling(D, data[2]);
if not HasDigraphGroup(D) then
# Multidigraphs not allowed
SetDigraphGroup(D, data[1]);
fi;
fi;
return data[1];
end);
InstallMethod(BlissAutomorphismGroup, "for a digraph and vertex coloring",
[IsDigraph, IsHomogeneousList], {D, colors} -> BLISS_DATA(D, colors, fail)[1]);
InstallMethod(BlissAutomorphismGroup,
"for a digraph, vertex colouring, and edge colouring",
[IsDigraph, IsHomogeneousList, IsList],
{digraph, vert_colours, edge_colours}
-> BLISS_DATA(digraph, vert_colours, edge_colours)[1]);
InstallMethod(BlissAutomorphismGroup,
"for a digraph, fail, and edge colouring",
[IsDigraph, IsBool, IsList],
function(digraph, vert_colours, edge_colours)
if vert_colours <> fail then
TryNextMethod();
fi;
return BLISS_DATA(digraph, vert_colours, edge_colours)[1];
end);
InstallMethod(NautyAutomorphismGroup, "for a digraph and vertex coloring",
[IsDigraph, IsHomogeneousList],
function(D, colors)
if IsMultiDigraph(D) then
Info(
InfoWarning,
1,
"NautyTracesInterfaces is not compatible with MultiDigraphs. Please ",
"consider calling either DigraphsUseBliss followed by ",
"AutomorphismGroup, or BlissAutomorphismGroup.");
return fail;
elif not DIGRAPHS_NautyAvailable then
Info(InfoWarning, 1, "NautyTracesInterface is not available");
return fail;
fi;
return NAUTY_DATA(D, colors)[1];
end);
InstallMethod(AutomorphismGroup, "for a digraph", [IsDigraph],
BlissAutomorphismGroup);
InstallMethod(AutomorphismGroup, "for a digraph and vertex coloring",
[IsDigraph, IsHomogeneousList], BlissAutomorphismGroup);
InstallMethod(AutomorphismGroup, "for a digraph, vertex and edge coloring",
[IsDigraph, IsHomogeneousList, IsList], BlissAutomorphismGroup);
InstallMethod(AutomorphismGroup, "for a digraph, vertex and edge coloring",
[IsDigraph, IsBool, IsList], BlissAutomorphismGroup);
# Check if two digraphs are isomorphic
InstallMethod(IsIsomorphicDigraph, "for digraphs", [IsDigraph, IsDigraph],
function(C, D)
local act;
if C = D then
return true;
elif DigraphNrVertices(C) <> DigraphNrVertices(D)
or DigraphNrEdges(C) <> DigraphNrEdges(D)
or IsMultiDigraph(C) <> IsMultiDigraph(D) then
return false;
fi; # JDM more!
if IsMultiDigraph(C) then
act := OnMultiDigraphs;
else
act := OnDigraphsNC;
fi;
if HasBlissCanonicalLabelling(C) and HasBlissCanonicalLabelling(D)
or not ((HasNautyCanonicalLabelling(C)
and NautyCanonicalLabelling(C) <> fail)
or (HasNautyCanonicalLabelling(D)
and NautyCanonicalLabelling(D) <> fail)) then
# Both digraphs either know their bliss canonical labelling or
# neither know their Nauty canonical labelling.
return act(C, BlissCanonicalLabelling(C))
= act(D, BlissCanonicalLabelling(D));
fi;
return act(C, NautyCanonicalLabelling(C))
= act(D, NautyCanonicalLabelling(D));
end);
InstallMethod(IsIsomorphicDigraph, "for digraphs and homogeneous lists",
[IsDigraph, IsDigraph, IsHomogeneousList, IsHomogeneousList],
function(C, D, c1, c2)
local m, colour1, n, colour2, max, class_sizes, act, i;
m := DigraphNrVertices(C);
n := DigraphNrVertices(D);
colour1 := DIGRAPHS_ValidateVertexColouring(m, c1);
colour2 := DIGRAPHS_ValidateVertexColouring(n, c2);
max := Maximum(colour1);
if max <> Maximum(colour2) then
return false;
fi;
# check some invariants
if m <> n
or DigraphNrEdges(C) <> DigraphNrEdges(D)
or IsMultiDigraph(C) <> IsMultiDigraph(D) then
# JDM more!
return false;
elif C = D and colour1 = colour2 then
return true;
fi;
class_sizes := ListWithIdenticalEntries(max, 0);
for i in DigraphVertices(C) do
class_sizes[colour1[i]] := class_sizes[colour1[i]] + 1;
class_sizes[colour2[i]] := class_sizes[colour2[i]] - 1;
od;
if not ForAll(class_sizes, x -> x = 0) then
return false;
elif IsMultiDigraph(C) then
act := OnMultiDigraphs;
else
act := OnDigraphsNC;
fi;
if DIGRAPHS_UsingBliss or IsMultiDigraph(C) then
return act(C, BlissCanonicalLabelling(C, colour1))
= act(D, BlissCanonicalLabelling(D, colour2));
fi;
return act(C, NautyCanonicalLabelling(C, colour1))
= act(D, NautyCanonicalLabelling(D, colour2));
end);
# Isomorphisms between digraphs
InstallMethod(IsomorphismDigraphs, "for digraphs", [IsDigraph, IsDigraph],
function(C, D)
local label1, label2;
if not IsIsomorphicDigraph(C, D) then
return fail;
elif IsMultiDigraph(C) then
if C = D then
return [(), ()];
fi;
label1 := BlissCanonicalLabelling(C);
label2 := BlissCanonicalLabelling(D);
return [label1[1] / label2[1], label1[2] / label2[2]];
elif C = D then
return ();
elif HasBlissCanonicalLabelling(C) and HasBlissCanonicalLabelling(D)
or not ((HasNautyCanonicalLabelling(C)
and NautyCanonicalLabelling(C) <> fail)
or (HasNautyCanonicalLabelling(D)
and NautyCanonicalLabelling(D) <> fail)) then
# Both digraphs either know their bliss canonical labelling or
# neither know their Nauty canonical labelling.
return BlissCanonicalLabelling(C) / BlissCanonicalLabelling(D);
fi;
return NautyCanonicalLabelling(C) / NautyCanonicalLabelling(D);
end);
InstallMethod(IsomorphismDigraphs, "for digraphs and homogeneous lists",
[IsDigraph, IsDigraph, IsHomogeneousList, IsHomogeneousList],
function(C, D, c1, c2)
local m, colour1, n, colour2, max, class_sizes, label1, label2, i;
m := DigraphNrVertices(C);
colour1 := DIGRAPHS_ValidateVertexColouring(m, c1);
n := DigraphNrVertices(D);
colour2 := DIGRAPHS_ValidateVertexColouring(n, c2);
max := Maximum(colour1);
if max <> Maximum(colour2) then
return fail;
fi;
# check some invariants
if m <> n
or DigraphNrEdges(C) <> DigraphNrEdges(D)
or IsMultiDigraph(C) <> IsMultiDigraph(D) then
return fail;
elif C = D and colour1 = colour2 then
if IsMultiDigraph(C) then
return [(), ()];
fi;
return ();
fi;
class_sizes := ListWithIdenticalEntries(max, 0);
for i in DigraphVertices(C) do
class_sizes[colour1[i]] := class_sizes[colour1[i]] + 1;
class_sizes[colour2[i]] := class_sizes[colour2[i]] - 1;
od;
if not ForAll(class_sizes, x -> x = 0) then
return fail;
elif DIGRAPHS_UsingBliss or IsMultiDigraph(C) then
label1 := BlissCanonicalLabelling(C, colour1);
label2 := BlissCanonicalLabelling(D, colour2);
else
label1 := NautyCanonicalLabelling(C, colour1);
label2 := NautyCanonicalLabelling(D, colour2);
fi;
if IsMultiDigraph(C) then
if OnMultiDigraphs(C, label1) <> OnMultiDigraphs(D, label2) then
return fail;
fi;
return [label1[1] / label2[1], label1[2] / label2[2]];
fi;
if OnDigraphsNC(C, label1) <> OnDigraphsNC(D, label2) then
return fail;
fi;
return label1 / label2;
end);
# Given a non-negative integer <n> and a homogeneous list <partition>,
# this global function tests whether <partition> is a valid partition
# of the list [1 .. n]. A valid partition of [1 .. n] is either:
#
# 1. A list of length <n> consisting of a numbers, such that the set of these
# numbers is [1 .. m] for some m <= n.
# 2. A list of non-empty disjoint lists whose union is [1 .. n].
#
# If <partition> is a valid partition of [1 .. n] then this global function
# returns the partition, in form 1 (converting it to this form if necessary).
# If <partition> is invalid, then the function returns <fail>.
InstallGlobalFunction(DIGRAPHS_ValidateVertexColouring,
function(n, partition)
local colours, i, missing, seen, x;
if not IsInt(n) or n < 0 then
ErrorNoReturn("the 1st argument <n> must be a non-negative integer,");
elif not IsHomogeneousList(partition) then
ErrorNoReturn("the 2nd argument <partition> must be a homogeneous list,");
elif n = 0 then
if IsEmpty(partition) then
return partition;
fi;
ErrorNoReturn("the only valid partition of the vertices of the digraph ",
"with 0 vertices is the empty list,");
elif not IsEmpty(partition) then
if IsPosInt(partition[1]) and Length(partition) = n then
# <partition> seems to be a list of colours
colours := [];
for i in partition do
if not IsPosInt(i) then
ErrorNoReturn("the 2nd argument <partition> does not define a ",
"colouring of the vertices [1 .. ", n, "], since it ",
"contains the element ", i, ", which is not a ",
"positive integer,");
elif i > n then
ErrorNoReturn("the 2nd argument <partition> does not define ",
"a colouring of the vertices [1 .. ", n, "], since ",
"it contains the integer ", i,
", which is greater than ", n, ",");
fi;
AddSet(colours, i);
od;
i := Length(colours);
missing := Difference([1 .. i], colours);
if not IsEmpty(missing) then
ErrorNoReturn("the 2nd argument <partition> does not define a ",
"colouring ",
"of the vertices [1 .. ", n, "], since it contains the ",
"colour ", colours[i], ", but it lacks the colour ",
missing[1], ". A colouring must use precisely the ",
"colours [1 .. m], for some positive integer m <= ", n,
",");
fi;
return partition;
elif IsList(partition[1]) then
seen := BlistList([1 .. n], []);
colours := EmptyPlist(n);
for i in [1 .. Length(partition)] do
# guaranteed to be non-empty since <partition> is homogeneous
for x in partition[i] do
if not IsPosInt(x) or x > n then
ErrorNoReturn("the 2nd argument <partition> does not define a ",
"colouring of the vertices [1 .. ", n, "], since ",
"the entry in position ", i, " contains ", x,
" which is not an integer in the range [1 .. ", n,
"],");
elif seen[x] then
ErrorNoReturn("the 2nd argument <partition> does not define a ",
"colouring of the vertices [1 .. ", n, "], since ",
"it contains the vertex ", x, " more than once,");
fi;
seen[x] := true;
colours[x] := i;
od;
od;
i := First([1 .. n], x -> not seen[x]);
if i <> fail then
ErrorNoReturn("the 2nd argument <partition> does not define a ",
"colouring of the vertices [1 .. ", n, "], since ",
"it does not assign a colour to the vertex ", i, ",");
fi;
return colours;
fi;
fi;
ErrorNoReturn("the 2nd argument <partition> does not define a ",
"colouring of the vertices [1 .. ", n, "]. The 2nd ",
"argument must have one of the following forms: ",
"1. a list of length ", n, " consisting of ",
"every integer in the range [1 .. m], for some m <= ", n,
"; or 2. a list of non-empty disjoint lists ",
"whose union is [1 .. ", n, "].");
end);
InstallGlobalFunction(DIGRAPHS_ValidateEdgeColouring,
function(graph, edge_colouring)
local n, colours, m, adj_colours, i;
if not IsDigraph(graph) then
ErrorNoReturn("the 1st argument must be a digraph,");
fi;
# Check: shapes and values from [1 .. something]
if edge_colouring = fail then
return true;
fi;
n := DigraphNrVertices(graph);
if not IsList(edge_colouring) or Length(edge_colouring) <> n then
ErrorNoReturn("the 2nd argument must be a list of the same shape as ",
"OutNeighbours(graph), where graph is the 1st argument,");
fi;
if ForAny(DigraphVertices(graph), x -> not IsList(edge_colouring[x]) or
(Length(edge_colouring[x]) <>
Length(OutNeighbours(graph)[x]))) then
ErrorNoReturn("the 2nd argument must be a list of the same shape as ",
"OutNeighbours(graph), where graph is the 1st argument,");
fi;
colours := [];
for adj_colours in edge_colouring do
for i in adj_colours do
if not IsPosInt(i) then
ErrorNoReturn("the 2nd argument should be a list of lists of ",
"positive integers,");
fi;
AddSet(colours, i);
od;
od;
m := Length(colours);
if ForAny([1 .. m], i -> i <> colours[i]) then
ErrorNoReturn("the 2nd argument should be a list of lists whose union ",
"is [1 .. number of colours],");
fi;
return true;
end);
InstallGlobalFunction(DIGRAPHS_CollapseMultiColouredEdges,
function(D, edge_colours)
local n, mults, out, new_cols, new_cols_set, idx, nbs, adjv, indices, colsv,
p, run, cur, range, cols, C, v, i;
n := DigraphNrVertices(D);
mults := [];
out := List([1 .. n], x -> []);
new_cols := List([1 .. n], x -> []);
new_cols_set := [];
idx := 1;
for v in [1 .. n] do
nbs := OutNeighbours(D)[v];
if not IsEmpty(nbs) then
adjv := ShallowCopy(nbs);
indices := [idx .. idx + Length(adjv) - 1];
colsv := ShallowCopy(edge_colours[v]);
p := Sortex(adjv);
colsv := Permuted(colsv, p);
indices := Permuted(indices, p);
run := 1;
cur := 1;
while cur <= Length(adjv) do
if cur < Length(adjv) and adjv[cur + 1] = adjv[cur] then
run := run + 1;
else
Add(out[v], adjv[cur]);
range := [cur - run + 1 .. cur];
cols := colsv{range};
C := List([1 .. Maximum(cols)], x -> []);
for i in range do
Add(C[colsv[i]], indices[i]);
od;
Append(mults, Filtered(C, x -> Length(x) > 1));
Sort(cols);
AddSet(new_cols_set, cols);
Add(new_cols[v], cols);
run := 1;
fi;
cur := cur + 1;
od;
idx := idx + Length(nbs);
fi;
od;
for cols in new_cols do
for i in [1 .. Length(cols)] do
cols[i] := Position(new_cols_set, cols[i]);
od;
od;
return [Digraph(out), new_cols, mults];
end);
InstallGlobalFunction(DIGRAPHS_CollapseMultipleEdges,
function(D)
local n, mults, out, new_cols, idx, nbs, adjv, indices, run, cur, v;
n := DigraphNrVertices(D);
mults := [];
out := List([1 .. n], x -> []);
new_cols := List([1 .. n], x -> []);
idx := 1;
for v in [1 .. n] do
nbs := OutNeighbours(D)[v];
if not IsEmpty(nbs) then
adjv := ShallowCopy(nbs);
indices := [idx .. idx + Length(adjv) - 1];
SortParallel(adjv, indices);
run := 1;
cur := 1;
while cur <= Length(adjv) do
if cur < Length(adjv) and adjv[cur + 1] = adjv[cur] then
run := run + 1;
else
Add(out[v], adjv[cur]);
Add(new_cols[v], run);
if run > 1 then
Add(mults, indices{[cur - run + 1 .. cur]});
fi;
run := 1;
fi;
cur := cur + 1;
od;
fi;
idx := idx + Length(nbs);
od;
return [Digraph(out), new_cols, mults];
end);
InstallMethod(IsDigraphIsomorphism, "for digraph, digraph, and permutation",
[IsDigraph, IsDigraph, IsPerm],
function(src, ran, x)
if IsMultiDigraph(src) or IsMultiDigraph(ran) then
ErrorNoReturn("the 1st and 2nd arguments <src> and <ran> must not have ",
"multiple edges,");
fi;
return IsDigraphHomomorphism(src, ran, x)
and IsDigraphHomomorphism(ran, src, x ^ -1);
end);
InstallMethod(IsDigraphIsomorphism,
"for digraph, digraph, permutation, list, and list",
[IsDigraph, IsDigraph, IsPerm, IsList, IsList],
function(src, ran, x, cols1, cols2)
return IsDigraphIsomorphism(src, ran, x)
and DigraphsRespectsColouring(src, ran, x, cols1, cols2);
end);
InstallMethod(IsDigraphAutomorphism, "for a digraph and a permutation",
[IsDigraph, IsPerm],
{D, x} -> IsDigraphIsomorphism(D, D, x));
InstallMethod(IsDigraphAutomorphism,
"for a digraph, a permutation, and a list",
[IsDigraph, IsPerm, IsList],
{D, x, c} -> IsDigraphIsomorphism(D, D, x, c, c));
InstallMethod(IsDigraphIsomorphism, "for digraph, digraph, and transformation",
[IsDigraph, IsDigraph, IsTransformation],
function(src, ran, x)
local y;
y := AsPermutation(RestrictedTransformation(x, DigraphVertices(src)));
if y = fail then
return false;
fi;
return IsDigraphIsomorphism(src, ran, y);
end);
InstallMethod(IsDigraphIsomorphism,
"for digraph, digraph, transformation, list, and list",
[IsDigraph, IsDigraph, IsTransformation, IsList, IsList],
function(src, ran, x, cols1, cols2)
return IsDigraphIsomorphism(src, ran, x)
and DigraphsRespectsColouring(src, ran, x, cols1, cols2);
end);
InstallMethod(IsDigraphAutomorphism, "for a digraph and a transformation",
[IsDigraph, IsTransformation],
{D, x} -> IsDigraphIsomorphism(D, D, x));
InstallMethod(IsDigraphAutomorphism,
"for a digraph, a transformation, and a list",
[IsDigraph, IsTransformation, IsList],
{D, x, c} -> IsDigraphIsomorphism(D, D, x, c, c));
[ Dauer der Verarbeitung: 0.43 Sekunden
(vorverarbeitet)
]
|
