main.py

import sys
import time
import random
import numpy

from genetics import Individual, Chromosome, Population, HEURISTICS
from Lineage import Lineage


def main(pol_n, len_i, len_g, lin_n, regen, n_threads, show, config_tag):

    THREADED = n_threads > 0
    THREADS = n_threads

    def lineageName(pol_n, gen_n):
        return "test/BernBase{}gen{}.lineage".format(pol_n, gen_n)

    #********** Generate the asked number of generations *********
    #   check

    if(regen):
        lin_n = pol_n
        lins = []
        # l is an int here
        individuals = [ Individual([Chromosome(), Chromosome()], pol_n-1) for _ in range(len_i) ]
        for indiv in individuals:
            for c in indiv.chromosomes:
                # generate a random set of weights for the chromosome
                weights_compressed = [[0]*pol_n]* (len(HEURISTICS)-1)
                for W_h in weights_compressed:
                    tmp_w = [0]*pol_n
                    for i in range(pol_n):
                        W_h[i] = max(0,1 -random.random() -tmp_w[i])
                        tmp_w[i] = tmp_w[i] + W_h[i]
                c.set_genes(numpy.array(weights_compressed))
        pop = Population(individuals, 1)
        print("First generation created for lineage number ", lin_n)
        lins.append(Lineage([pop], "0"))
        fname = lineageName(pol_n, lins[0].generations)
        print("saving to file : ", fname)
        lins[0].toFile(fname)

        # l is a lineage here
        for l in lins:
            while l.generations < len_g:
                print("Creating generation ",l.generations+1)
                l.nextGeneration(THREADED, THREADS)
                fname = lineageName(pol_n, l.generations)
                print("saving to file : ", fname)
                l.toFile(fname)
        lin = lins[-1]
    else:
        for g in range(len_g, 0, -1):
            try:
                fname = lineageName(pol_n, g)
                fh = open(fname, 'r')
                # Store configuration file values
                lin = Lineage.fromFile(fname)
                break
            except FileNotFoundError:
                if(g <= 1):
                    print("No previous generations were found, please add the -r or --regenerate flag")
                    exit("no files found : "+lineageName("<Bernstein base>", "<generationNumber>"))
    # print(lin)
    lin.beatYourElders(show, True, True, True, n_threads, config_tag)
    # Lineage.plotScorePerGeneration(lin.populations)


if __name__ == "__main__":
    import argparse
    parser = argparse.ArgumentParser(description='Test the algorithm for a given number of generations.')
    parser.add_argument('-n', type=int, default=0,
                        help='order of the Bernstein base (default: 0)')
    parser.add_argument('-i', type=int, default=2,
                        help='number of individuals in the population (default: 2)')
    parser.add_argument('-g', type=int, default=2,
                        help='number of succeeding generations (default: 2)')
    parser.add_argument('-l', type=int, default=1,
                        help='number of different lineages (default: 1)')
    parser.add_argument('-t', '--threads', type=int, default=-1,
                        help='number of different lineages (default: -1, non threaded)')
    parser.add_argument('--show', type=int, default=5,
                        help='show the last X genertions performance (default: 5)')
    parser.add_argument('-r', "--regenerate", action="store_true",
                        help="Regenerates the lineage from start (default : tries to load from file)")
    parser.add_argument('-m' "--mutation", type=int, default=10,
                        help='number of mutations per 32-bit weight (default: 10)')
    parser.add_argument('-c' "--cut", type=float, default=0.9,
                        help='approximate quota of dying individuals per generation (default: 0.9)')

    args = parser.parse_args()


    """
    Hyperparameter optimization

    mutation_strengths = [2, 4, 6, 8, 10]
    median_cuts = [0.2, 0.4, 0.6, 0.8, 0.9]

    for mut_strength in mutation_strengths:
        for median_cut in median_cuts:
            print("ms=" + str(mut_strength) + ",cut=" + str(median_cut))
            Population.MEAN_CUT_SELECTION = median_cut
            Chromosome.MUTATION_STRENGTH = mut_strength
            main(args.n+1, args.i, args.g, args.l, args.regenerate, args.threads, "ms=" + str(mut_strength) + ", cut=" + str(median_cut))

    test = input("SAVE ALL FIGURES FIRST! Then press enter to close")
    """
    Population.MEAN_CUT_SELECTION = 0.9
    Chromosome.MUTATION_STRENGTH = 10
    main(args.n+1, args.i, args.g, args.l, args.regenerate, args.threads, args.show, "ms = " + str(10) + ", cut = " + str(0.9))