public void Cube2eneWithoutBonds()
{
AtomRefinable ar = Make(Cubene(), true);
EquitablePartitionRefiner refiner = new EquitablePartitionRefiner(ar);
Partition finer = refiner.Refine(Partition.Unit(8));
Partition expected = Partition.Unit(8);
Assert.AreEqual(expected, finer);
}
public void Cube2eneWithBonds()
{
AtomRefinable ar = Make(Cubene(), false);
EquitablePartitionRefiner refiner = new EquitablePartitionRefiner(ar);
Partition finer = refiner.Refine(Partition.Unit(8));
Partition expected = Partition.FromString("0,2,5,7|1,3,4,6");
Assert.AreEqual(expected, finer);
}
public void RefineTest()
{
EquitablePartitionRefiner refiner = new EquitablePartitionRefiner(MakeExampleTable());
Partition coarser = Partition.FromString("[0|1,2,3]");
Partition finer = refiner.Refine(coarser);
Partition expected = Partition.FromString("[0|1,2|3]");
Assert.AreEqual(expected, finer);
}
/// <summary>
/// Does the work of the class, that refines a coarse partition into a finer
/// one using the supplied automorphism group to prune the search.
/// </summary>
/// <param name="group">the automorphism group of the graph</param>
/// <param name="coarser">the partition to refine</param>
private void Refine(PermutationGroup group, Partition coarser)
{
int vertexCount = GetVertexCount();
Partition finer = equitableRefiner.Refine(coarser);
int firstNonDiscreteCell = finer.GetIndexOfFirstNonDiscreteCell();
if (firstNonDiscreteCell == -1)
{
firstNonDiscreteCell = vertexCount;
}
Permutation pi1 = new Permutation(firstNonDiscreteCell);
Result result = Result.Better;
if (bestExist)
{
pi1 = finer.SetAsPermutation(firstNonDiscreteCell);
result = CompareRowwise(pi1);
}
// partition is discrete
if (finer.Count == vertexCount)
{
if (!bestExist)
{
best = finer.ToPermutation();
first = finer.ToPermutation();
bestExist = true;
}
else
{
if (result == Result.Better)
{
best = new Permutation(pi1);
}
else if (result == Result.Equal)
{
group.Enter(pi1.Multiply(best.Invert()));
}
}
}
else
{
if (result != Result.Worse)
{
var blockCopy = finer.CopyBlock(firstNonDiscreteCell);
for (int vertexInBlock = 0; vertexInBlock < vertexCount; vertexInBlock++)
{
if (blockCopy.Contains(vertexInBlock))
{
Partition nextPartition = finer.SplitBefore(firstNonDiscreteCell, vertexInBlock);
this.Refine(group, nextPartition);
int[] permF = new int[vertexCount];
int[] invF = new int[vertexCount];
for (int i = 0; i < vertexCount; i++)
{
permF[i] = i;
invF[i] = i;
}
for (int j = 0; j <= firstNonDiscreteCell; j++)
{
int x = nextPartition.GetFirstInCell(j);
int i = invF[x];
int h = permF[j];
permF[j] = x;
permF[i] = h;
invF[h] = i;
invF[x] = j;
}
Permutation pPermF = new Permutation(permF);
group.ChangeBase(pPermF);
for (int j = 0; j < vertexCount; j++)
{
Permutation g = group[firstNonDiscreteCell, j];
if (g != null)
{
blockCopy.Remove(g[vertexInBlock]);
}
}
}
}
}
}
}