private void GenerateCliques(int[] variableSizes, IList <int> startSet, ICollection <int> alreadyRepresented, GraphicalModel model, SourceOfRandomness randomness)
            {
                if (alreadyRepresented.Count == variableSizes.Length)
                {
                    return;
                }
                // Generate the clique variable set
                IList <int> cliqueContents = new List <int>();

                Sharpen.Collections.AddAll(cliqueContents, startSet);
                Sharpen.Collections.AddAll(alreadyRepresented, startSet);
                while (true)
                {
                    if (alreadyRepresented.Count == variableSizes.Length)
                    {
                        break;
                    }
                    if (cliqueContents.Count == 0 || randomness.NextDouble(0, 1) < 0.7)
                    {
                        int gen;
                        do
                        {
                            gen = randomness.NextInt(variableSizes.Length);
                        }while (alreadyRepresented.Contains(gen));
                        alreadyRepresented.Add(gen);
                        cliqueContents.Add(gen);
                    }
                    else
                    {
                        break;
                    }
                }
                // Create the actual table
                int[] neighbors     = new int[cliqueContents.Count];
                int[] neighborSizes = new int[neighbors.Length];
                for (int j = 0; j < neighbors.Length; j++)
                {
                    neighbors[j]     = cliqueContents[j];
                    neighborSizes[j] = variableSizes[neighbors[j]];
                }
                ConcatVectorTable table = new ConcatVectorTable(neighborSizes);

                foreach (int[] assignment in table)
                {
                    // Generate a vector
                    ConcatVector v = new ConcatVector(ConcatVecComponents);
                    for (int x = 0; x < ConcatVecComponents; x++)
                    {
                        if (randomness.NextBoolean())
                        {
                            v.SetSparseComponent(x, randomness.NextInt(32), randomness.NextDouble());
                        }
                        else
                        {
                            double[] val = new double[randomness.NextInt(12)];
                            for (int y = 0; y < val.Length; y++)
                            {
                                val[y] = randomness.NextDouble();
                            }
                            v.SetDenseComponent(x, val);
                        }
                    }
                    // set vec in table
                    table.SetAssignmentValue(assignment, null);
                }
                model.AddFactor(table, neighbors);
                // Pick the number of children
                IList <int> availableVariables = new List <int>();

                Sharpen.Collections.AddAll(availableVariables, cliqueContents);
                availableVariables.RemoveAll(startSet);
                int numChildren = randomness.NextInt(0, availableVariables.Count);

                if (numChildren == 0)
                {
                    return;
                }
                IList <IList <int> > children = new List <IList <int> >();

                for (int i = 0; i < numChildren; i++)
                {
                    children.Add(new List <int>());
                }
                // divide up the shared variables across the children
                int cursor = 0;

                while (true)
                {
                    if (availableVariables.Count == 0)
                    {
                        break;
                    }
                    if (children[cursor].Count == 0 || randomness.NextBoolean())
                    {
                        int gen = randomness.NextInt(availableVariables.Count);
                        children[cursor].Add(availableVariables[gen]);
                        availableVariables.Remove(availableVariables[gen]);
                    }
                    else
                    {
                        break;
                    }
                    cursor = (cursor + 1) % numChildren;
                }
                foreach (IList <int> shared1 in children)
                {
                    foreach (int i_1 in shared1)
                    {
                        foreach (IList <int> shared2 in children)
                        {
                            System.Diagnostics.Debug.Assert((shared1 == shared2 || !shared2.Contains(i_1)));
                        }
                    }
                }
                foreach (IList <int> shared in children)
                {
                    if (shared.Count > 0)
                    {
                        GenerateCliques(variableSizes, shared, alreadyRepresented, model, randomness);
                    }
                }
            }