public void IndexOf2(int size)
        {
            var random = RandomRegistry.GetRandom();

            var props   = new double[size];
            var divisor = props.Length * (props.Length + 1) / 2.0;

            for (var i = 0; i < props.Length; ++i)
            {
                props[i] = (i + 1) / divisor;
            }
            Randomize(props, random);

            var incremental = ProbabilitySelector.Incremental((double[])props.Clone());

            const int samples = 100000;

            for (var i = 0; i < samples; ++i)
            {
                var value  = random.NextDouble();
                var index1 = ProbabilitySelector.IndexOf(incremental, value);
                var index2 = IndexOf(props, value);

                Assert.Equal(index2, index1);
            }
        }
        public void IndexOf1()
        {
            var random = RandomRegistry.GetRandom();

            var props   = new double[10];
            var divisor = props.Length * (props.Length + 1) / 2.0;

            for (var i = 0; i < props.Length; ++i)
            {
                props[i] = (i + 1) / divisor;
            }
            Randomize(props, random);

            var incremental = ProbabilitySelector.Incremental((double[])props.Clone());

            const double samples = 1000000;
            var          indices = new double[props.Length];

            Array.Fill(indices, 0);

            for (var i = 0; i < samples; ++i)
            {
                indices[ProbabilitySelector.IndexOf(incremental, random.NextDouble())] += 1;
            }

            for (var i = 0; i < props.Length; ++i)
            {
                indices[i] /= samples;
            }

            for (var i = 0; i < props.Length; ++i)
            {
                Assert.Equal(props[i], indices[i], 1);
            }
        }
Exemplo n.º 3
0
        public virtual Population <TGene, TAllele> Select(Population <TGene, TAllele> population, int count, Optimize opt)
        {
            if (count < 0)
            {
                throw new ArgumentOutOfRangeException(nameof(count),
                                                      $"Selection count must be greater or equal then zero, but was {count}.");
            }

            var selection = new Population <TGene, TAllele>(count);

            if (count > 0 && !population.IsEmpty)
            {
                var pop = Copy(population);

                var prob = Probabilities(pop, count, opt);

                CheckAndCorrect(prob);

                ProbabilitySelector.Incremental(prob);

                var random = RandomRegistry.GetRandom();
                selection.Fill(() => pop[ProbabilitySelector.IndexOf(prob, random.NextDouble())], count);
            }

            return(selection);
        }
Exemplo n.º 4
0
        public void IndexOfSerialEqualBinary(int size)
        {
            var probabilities = Array(size, new Random(12));

            Arithmetics.Normalize(probabilities);
            ProbabilitySelector.Incremental(probabilities);

            Assert.Equal(
                ProbabilitySelector.IndexOfSerial(probabilities, 0.5),
                ProbabilitySelector.IndexOfBinary(probabilities, 0.5)
                );
        }