Ejemplo n.º 1
0
def specie_h(population, p, n, pset, direccion):  # heuristic applied by specie
    gpo_specie=specie_gpo(population)
    for specie in gpo_specie:
        eval_prob(specie[0])
    for ind in population:
        if ind.bestspecie_get()==1:
            if random.random()<=ind.get_LS_prob():
                strg=ind.__str__()
                args=[]
                if len(pset.arguments) > 0:
                    for arg in pset.arguments:
                        args.append(arg)
                l_strg=add_subt_cf(strg, args)
                c = tree2f()
                cd=c.convert(l_strg)
                xdata,ydata=get_address(p, n,direccion)
                beta_opt, beta_cov, success, nfev = curve_fit_2(eval_,cd , xdata, ydata, p0=ind.get_params(), method='trf', max_nfev=40)
                if not success:
                    ind.LS_applied_set(1)
                    ind.LS_story_set(1)
                else:
                    ind.LS_applied_set(1)
                    ind.LS_story_set(1)
                    ind.params_set(beta_opt)
                    funcEval.cont_evalp += nfev
Ejemplo n.º 2
0
def specie_h(population, p, n, pset, direccion):  # heuristic applied by specie
    gpo_specie=specie_gpo(population)
    for specie in gpo_specie:
        eval_prob(specie[0])
    for ind in population:
        if ind.bestspecie_get()==1:
            if random.random()<=ind.get_LS_prob():
                strg=ind.__str__()
                args=[]
                if len(pset.arguments) > 0:
                    for arg in pset.arguments:
                        args.append(arg)
                l_strg=add_subt_cf(strg, args)
                c = tree2f()
                cd=c.convert(l_strg)
                xdata,ydata=get_address(p, n,direccion)
                beta_opt, beta_cov, success, nfev = curve_fit_2(eval_,cd , xdata, ydata, p0=ind.get_params(), method='trf', max_nfev=40)
                if not success:
                    ind.LS_applied_set(1)
                    ind.LS_story_set(1)
                else:
                    ind.LS_applied_set(1)
                    ind.LS_story_set(1)
                    ind.params_set(beta_opt)
                    funcEval.cont_evalp += nfev
Ejemplo n.º 3
0
def ls_random(population, p, n, pset, direccion, problem,
              benchmark_flag):  # size heuristic
    random_set = random_set_pop(population)
    for ind in random_set:
        strg = ind.__str__()
        args = []
        if len(pset.arguments) > 1:
            for arg in pset.arguments:
                args.append(arg)
        l_strg = add_subt_cf(strg, args)
        c = tree2f()
        cd = c.convert(l_strg)
        xdata, ydata = get_address(p, n, problem, direccion, benchmark_flag)
        beta_opt, beta_cov, success, nfev = curve_fit_2(eval_,
                                                        cd,
                                                        xdata,
                                                        ydata,
                                                        p0=ind.get_params(),
                                                        method='trf',
                                                        max_nfev=40)
        if not success:
            ind.LS_applied_set(1)
            ind.LS_story_set(1)
        else:
            ind.LS_applied_set(1)
            ind.LS_story_set(1)
            funcEval.cont_evalp += nfev
            ind.params_set(beta_opt)
Ejemplo n.º 4
0
def best_specie(population, p, n, pset, direccion, problem,
                benchmark_flag):  # best of each specie
    for ind in population:
        if ind.bestspecie_get() == 1:
            strg = ind.__str__()
            args = []
            if len(pset.arguments) > 0:
                for arg in pset.arguments:
                    args.append(arg)
            l_strg = add_subt_cf(strg, args)
            c = tree2f()
            cd = c.convert(l_strg)
            xdata, ydata = get_address(p, n, problem, direccion,
                                       benchmark_flag)
            beta_opt, beta_cov, success, nfev = curve_fit_2(
                eval_,
                cd,
                xdata,
                ydata,
                p0=ind.get_params(),
                method='trf',
                max_nfev=40)
            if not success:
                ind.LS_applied_set(1)
                ind.LS_story_set(1)
            else:
                ind.LS_applied_set(1)
                ind.LS_story_set(1)
                ind.params_set(beta_opt)
                funcEval.cont_evalp += nfev
Ejemplo n.º 5
0
def best_pop_ls(population, p,  n, pset, direccion):  # best of the pop
    ind = best_pop(population)
    strg = ind.__str__()
    args = []
    if len(pset.arguments) > 1:
        for arg in pset.arguments:
            args.append(arg)
    l_strg = add_subt_cf(strg, args)
    c = tree2f()
    cd = c.convert(l_strg)
    xdata,ydata = get_address(p, n, direccion)
    beta_opt, beta_cov, success, nfev = curve_fit_2(eval_, cd, xdata, ydata, p0=ind.get_params(), method='trf', max_nfev=40)
    if not success:
        ind.LS_applied_set(1)
        ind.LS_story_set(1)
    else:
        ind.LS_applied_set(1)
        ind.LS_story_set(1)
        ind.params_set(beta_opt)
        funcEval.cont_evalp += nfev
Ejemplo n.º 6
0
def best_pop_ls(population, p,  n, pset, direccion):  # best of the pop
    ind = best_pop(population)
    strg = ind.__str__()
    args = []
    if len(pset.arguments) > 1:
        for arg in pset.arguments:
            args.append(arg)
    l_strg = add_subt_cf(strg, args)
    c = tree2f()
    cd = c.convert(l_strg)
    xdata,ydata = get_address(p, n, direccion)
    beta_opt, beta_cov, success, nfev = curve_fit_2(eval_, cd, xdata, ydata, p0=ind.get_params(), method='trf', max_nfev=40)
    if not success:
        ind.LS_applied_set(1)
        ind.LS_story_set(1)
    else:
        ind.LS_applied_set(1)
        ind.LS_story_set(1)
        ind.params_set(beta_opt)
        funcEval.cont_evalp += nfev
Ejemplo n.º 7
0
def ls_random(population, p, n, pset, direccion, problem):  # size heuristic
    random_set=random_set_pop(population)
    for ind in random_set:
        strg = ind.__str__()
        args = []
        if len(pset.arguments) > 1:
            for arg in pset.arguments:
                args.append(arg)
        l_strg = add_subt_cf(strg, args)
        c = tree2f()
        cd = c.convert(l_strg)
        xdata, ydata = get_address(p, n, problem,direccion)
        beta_opt, beta_cov, success, nfev = curve_fit_2(eval_,cd , xdata, ydata, p0=ind.get_params(), method='trf', max_nfev=40)
        if not success:
            ind.LS_applied_set(1)
            ind.LS_story_set(1)
        else:
            ind.LS_applied_set(1)
            ind.LS_story_set(1)
            funcEval.cont_evalp += nfev
            ind.params_set(beta_opt)
Ejemplo n.º 8
0
def neat_GP_LS(population,
               toolbox,
               cxpb,
               mutpb,
               ngen,
               neat_alg,
               neat_cx,
               neat_h,
               neat_pelit,
               LS_flag,
               LS_select,
               cont_evalf,
               num_salto,
               SaveMatrix,
               GenMatrix,
               pset,
               n_corr,
               num_p,
               params,
               direccion,
               problem,
               testing,
               version,
               benchmark_flag,
               beta,
               random_speciation,
               config,
               stats=None,
               halloffame=None,
               verbose=__debug__):
    """This algorithm reproduce the simplest evolutionary algorithm as
    presented in chapter 7 of [Back2000]_.

    :param population: A list of individuals.
    :param toolbox: A :class:`~deap.base.Toolbox` that contains the evolution
                    operators.
    :param cxpb: The probability of mating two individuals.
    :param mutpb: The probability of mutating an individual.
    :param ngen: The number of generation.
    :param neat_alg: wheter or not to use species stuff.
    :param neat_cx: wheter or not to use neatGP cx
    :param neat_h: indicate the distance allowed between each specie
    :param neat_pelit: probability of being elitist, it's used in the neat cx and mutation
    :param LS_flag: wheter or not to use LocalSearchGP
    :param LS_select: indicate the kind of selection to use the LSGP on the population.
    :param cont_evalf: contador maximo del numero de evaluaciones
    :param n_corr: run number just to wirte the txt file
    :param p: problem number just to wirte the txt file
    :param params:indicate the params for the fitness sharing, the diffetent
                    options are:
                    -DontPenalize(str): 'best_specie' or 'best_of_each_specie'
                    -Penalization_method(int):
                        1.without penalization
                        2.penalization fitness sharing
                        3.new penalization
                    -ShareFitness(str): 'yes' or 'no'
    :param stats: A :class:`~deap.tools.Statistics` object that is updated
                  inplace, optional.
    :param halloffame: A :class:`~deap.tools.HallOfFame` object that will
                       contain the best individuals, optional.
    :param verbose: Whether or not to log the statistics.
    :returns: The final population.

    It uses :math:`\lambda = \kappa = \mu` and goes as follow.
    It first initializes the population (:math:`P(0)`) by evaluating
    every individual presenting an invalid fitness. Then, it enters the
    evolution loop that begins by the selection of the :math:`P(g+1)`
    population. Then the crossover operator is applied on a proportion of
    :math:`P(g+1)` according to the *cxpb* probability, the resulting and the
    untouched individuals are placed in :math:`P'(g+1)`. Thereafter, a
    proportion of :math:`P'(g+1)`, determined by *mutpb*, is
    mutated and placed in :math:`P''(g+1)`, the untouched individuals are
    transferred :math:`P''(g+1)`. Finally, those new individuals are evaluated
    and the evolution loop continues until *ngen* generations are completed.
    Briefly, the operators are applied in the following order ::

        evaluate(population)
        for i in range(ngen):
            offspring = select(population)
            offspring = mate(offspring)
            offspring = mutate(offspring)
            evaluate(offspring)
            population = offspring

    This function expects :meth:`toolbox.mate`, :meth:`toolbox.mutate`,
    :meth:`toolbox.select` and :meth:`toolbox.evaluate` aliases to be
    registered in the toolbox.

    .. [Back2000] Back, Fogel and Michalewicz, "Evolutionary Computation 1 :
       Basic Algorithms and Operators", 2000.
    """

    logbook = tools.Logbook()
    logbook.header = ['gen', 'nevals'] + (stats.fields if stats else [])

    #creating files to save data.
    d = './Results/%s/pop_file_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    pop_file = open(d, 'a')

    d = './Results/%s/bestind_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    best = open(d, 'w')  # save data

    d = './Results/%s/bestind_str_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    out = open(d, 'w')

    d = './Timing/%s/pop_file_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    time_file = open(d, 'w')

    d = './Timing/%s/specie_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    time_specie = open(d, 'w')

    d = './Timing/%s/cruce_s_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    time_cx = open(d, 'w')

    d = './Specie/%s/specieind_%d_%d.csv' % (problem, num_p, n_corr)
    ensure_dir(d)
    specie_file = open(d, 'w')

    d = './Specie/%s/speciepop_%d_%d.csv' % (problem, num_p, n_corr)
    ensure_dir(d)
    specie_file_2 = open(d, 'w')

    d = './Specie/%s/specist_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    specie_statis = open(d, 'w')

    d = './Results/%s/bestind_LStr_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    bestind = open(d, 'w')

    d = './Matrix/%s/' % (problem)
    ensure_dir(d)

    begin = time.time()

    if SaveMatrix:  # Saving data in matrix
        num_r = 11
        num_r_sp = 2
        if GenMatrix:
            num_salto = 1
            num_c = ngen + 1
            Matrix = np.empty((
                num_c,
                num_r,
            ))
            Matrix_specie = np.empty((
                num_c,
                num_r_sp,
            ), dtype=object)
            vector = np.arange(0, num_c, num_salto)
        else:
            num_c = (cont_evalf / num_salto) + 1
            num_c_sp = (cont_evalf / num_salto) + 1
            Matrix = np.empty((
                num_c,
                num_r,
            ))
            Matrix_specie = np.empty((
                num_c_sp,
                num_r_sp,
            ), dtype=object)
            vector = np.arange(1, cont_evalf + num_salto, num_salto)
        for i in range(len(vector)):
            Matrix[i, 0] = vector[i]
            Matrix_specie[i, 0] = vector[i]
        Matrix[:, 6] = 0.

    if neat_alg:
        begin_sp = time.time()
        if version == 1:
            for ind in population:
                bit = p_bin(ind)
                ind.binary_rep_set(bit)
        elif version != 1:
            for ind in population:
                level_info = level_node(ind)
                ind.nodefeat_set(level_info)

        if random_speciation:
            species_random(population, neat_h, version, beta)
        else:
            species(population, neat_h, version, beta)

        end_sp = time.time()
        time_specie.write(
            '\n%s;%s;%s;%s' %
            (0, begin_sp, end_sp, str(round(end_sp - begin_sp, 2))))

        for ind in population:
            specie_file_2.write('\n%s,%s,%s' %
                                (0, ind.get_specie(), ind.get_intracluster()))

        list_s = list_species(population)
        for specie in list_s:
            list_ind = get_ind_specie(specie, population)
            specie_file.write('\n%s,%s,%s' %
                              (0, specie, list_ind[0].get_intracluster()))

    if funcEval.LS_flag:
        for ind in population:
            sizep = len(ind) + 2
            param_ones = np.ones(sizep)
            param_ones[0] = 0
            ind.params_set(param_ones)

    # Evaluate the individuals with an invalid fitness
    invalid_ind = [ind for ind in population if not ind.fitness.valid]
    fitnesses = toolbox.map(toolbox.evaluate, invalid_ind)
    for ind, fit in zip(invalid_ind, fitnesses):
        funcEval.cont_evalp += 1
        ind.fitness.values = fit

    best_ind = best_pop(population)  # best individual of the population
    if testing:
        fitnesst_best = toolbox.map(toolbox.evaluate_test, [best_ind])
        best_ind.fitness_test.values = fitnesst_best[0]
    if testing:
        best.write(
            '\n%s;%s;%s;%s;%s;%s' %
            (0, funcEval.cont_evalp, best_ind.fitness_test.values[0],
             best_ind.fitness.values[0], len(best_ind), avg_nodes(population)))
    else:
        best.write('\n%s;%s;%s;%s;%s;%s' %
                   (0, funcEval.cont_evalp, None, best_ind.fitness.values[0],
                    len(best_ind), avg_nodes(population)))

    data_pop = avg_nodes(population)

    if SaveMatrix:
        idx = 0
        Matrix[idx, 1] = best_ind.fitness.values[0]
        if testing:
            Matrix[idx, 2] = best_ind.fitness_test.values[0]
        else:
            Matrix[idx, 2] = None
        Matrix[idx, 3] = len(best_ind)
        Matrix[idx, 4] = data_pop[0]
        Matrix[idx, 5] = 0.
        Matrix[idx, 6] = 1  # just an id to know if the current row is full
        Matrix[idx, 7] = data_pop[1]  # max size
        Matrix[idx, 8] = data_pop[2]  # min size
        Matrix[idx, 9] = ind.get_specie(
        )  # max number of species in the current population
        Matrix[idx, 10] = max(ind_specie(population), key=lambda x: x[0])[
            0]  # max number of species in the current population

        np.savetxt('./Matrix/%s/idx_%d_%d.txt' % (problem, num_p, n_corr),
                   Matrix,
                   delimiter=",",
                   fmt="%s")

    if neat_alg:
        SpeciesPunishment(population, params, neat_h)

    if funcEval.LS_flag == 1:
        strg = best_ind.__str__()
        l_strg = add_subt_cf(strg, args=[])
        c = tree2f()
        cd = c.convert(l_strg)
        out.write('\n%s;%s;%s;%s;%s;%s' %
                  (0, len(best_ind), best_ind.LS_applied_get(),
                   best_ind.get_params(), cd, best_ind))
    else:
        out.write('\n%s;%s;%s' % (0, len(best_ind), best_ind))

    for ind in population:
        pop_file.write('\n%s;%s;%s' % (0, ind.fitness.values[0], ind))

    print '---- Generation %d -----' % (0)
    print 'Problem: ', problem
    print 'Problem No.: ', num_p
    print 'Run No.: ', n_corr
    print 'neat-GP:', neat_alg
    print 'neat-cx:', neat_cx
    print 'Local Search:', funcEval.LS_flag
    if funcEval.LS_flag:
        ls_type = ''
        if LS_select == 1:
            ls_type = 'LSHS'
        elif LS_select == 2:
            ls_type = 'Best-Sp'
        elif LS_select == 3:
            ls_type = 'LSHS-Sp'
        elif LS_select == 4:
            ls_type = 'Best-Pop'
        elif LS_select == 5:
            ls_type = 'All-Pop'
        elif LS_select == 6:
            ls_type = 'LSHS-test'
        elif LS_select == 7:
            ls_type = 'Best set'
        elif LS_select == 8:
            ls_type = 'Random set'
        elif LS_select == 9:
            ls_type = "Best-Random set"
        print 'Local Search Heuristic: %s (%s)' % (LS_select, ls_type)
    print 'Best Ind.:', best_ind
    print 'Best Fitness:', best_ind.fitness.values[0]
    if testing:
        print 'Test fitness:', best_ind.fitness_test.values[0]
    print 'Avg Nodes:', avg_nodes(population)
    print 'Evaluations: ', funcEval.cont_evalp
    end_t = time.time()

    if SaveMatrix:
        idx = 0
        Matrix_specie[idx, 1] = ind_specie(population)
        np.savetxt('./Specie/%s/specist_%d_%d.txt' % (problem, num_p, n_corr),
                   Matrix_specie,
                   delimiter=";",
                   fmt="%s")
    if testing:
        time_file.write(
            '\n%s;%s;%s;%s;%s;%s' %
            (0, begin, end_t, str(round(end_t - begin, 2)),
             best_ind.fitness.values[0], best_ind.fitness_test.values[0]))
    else:
        time_file.write('\n%s;%s;%s;%s;%s' %
                        (0, begin, end_t, str(round(
                            end_t - begin, 2)), best_ind.fitness.values[0]))

    # Begin the generational process
    for gen in range(1, ngen + 1):
        print "hola"
        if not GenMatrix:
            if funcEval.cont_evalp > cont_evalf:
                break

        begin = time.time()
        print '---- Generation %d -----' % (gen)
        print 'Problem: ', problem
        print 'Problem No.: ', num_p
        print 'Run No.: ', n_corr
        print 'neat-GP:', neat_alg
        print 'neat-cx:', neat_cx
        print 'Local Search:', funcEval.LS_flag
        if funcEval.LS_flag:
            print 'Local Search Heuristic: %s (%s)' % (LS_select, ls_type)

        best_ind = copy.deepcopy(best_pop(population))
        if neat_alg:
            # survival_p = p_worst
            parents = p_selection(population, survival_p=0.5)
        else:
            parents = toolbox.select(population, len(population))

        begin_cx = time.time()
        if neat_alg:  # neat-Crossover neat_cx and
            n = len(parents)
            mut = 1
            cx = 1
            offspring = neatGP(toolbox, parents, cxpb, mutpb, n, mut, cx,
                               neat_pelit, neat_cx)
        else:
            offspring = varOr(parents, toolbox, cxpb, mutpb)
        end_cx = time.time()
        time_specie.write(
            '\n%s;%s;%s;%s' %
            (0, begin_cx, end_cx, str(round(end_cx - begin_cx, 2))))

        if neat_alg:  # Speciation of the descendants
            begin_sp = time.time()

            if version == 1:
                for ind in offspring:
                    if ind.binary_rep_get() is None:
                        bit = p_bin(ind)
                        ind.binary_rep_set(bit)

            elif version == 2:
                for ind in offspring:
                    if ind.nodefeat_get() is None:
                        level_info = level_node(ind)
                        ind.nodefeat_set(level_info)

            if random_speciation:
                specie_offspring_random(parents, offspring, neat_h, version,
                                        beta)
            else:
                specie_parents_child(parents, offspring, neat_h, version, beta)

            #calc_intracluster(population)
            offspring[:] = parents + offspring
            calc_intracluster(offspring)

            for ind in offspring:
                specie_file_2.write(
                    '\n%s,%s,%s' %
                    (gen, ind.get_specie(), ind.get_intracluster()))

            list_s = list_species(offspring)
            for specie in list_s:
                list_ind = get_ind_specie(specie, offspring)
                specie_file.write(
                    '\n%s,%s,%s' %
                    (gen, specie, list_ind[0].get_intracluster()))

            invalid_ind = [ind for ind in offspring]
            if funcEval.LS_flag:
                new_invalid_ind = []
                for ind in invalid_ind:
                    strg = ind.__str__()
                    l_strg = add_subt(strg, ind)
                    c = tree2f()
                    cd = c.convert(l_strg)
                    new_invalid_ind.append(cd)
                fitness_ls = toolbox.map(toolbox.evaluate, new_invalid_ind)
                for ind, ls_fit in zip(invalid_ind, fitness_ls):
                    funcEval.cont_evalp += 1
                    ind.fitness.values = ls_fit
            else:
                fitnesses = toolbox.map(toolbox.evaluate, invalid_ind)
                for ind, fit in zip(invalid_ind, fitnesses):
                    funcEval.cont_evalp += 1
                    ind.fitness.values = fit

            end_sp = time.time()
            time_specie.write(
                '\n%s;%s;%s;%s' %
                (gen, begin_sp, end_sp, str(round(end_sp - begin_sp, 2))))
            print "neat-alg"
        else:
            invalid_ind = [ind for ind in offspring]
            if funcEval.LS_flag:
                new_invalid_ind = []
                for ind in invalid_ind:
                    strg = ind.__str__()
                    l_strg = add_subt(strg, ind)
                    c = tree2f()
                    cd = c.convert(l_strg)
                    new_invalid_ind.append(cd)
                fitness_ls = toolbox.map(toolbox.evaluate, new_invalid_ind)
                for ind, ls_fit in zip(invalid_ind, fitness_ls):
                    funcEval.cont_evalp += 1
                    ind.fitness.values = ls_fit
            else:
                fitnesses = toolbox.map(toolbox.evaluate, invalid_ind)
                for ind, fit in zip(invalid_ind, fitnesses):
                    funcEval.cont_evalp += 1
                    ind.fitness.values = fit

            orderbyfit = sorted(offspring, key=lambda ind: ind.fitness.values)

            if best_ind.fitness.values[0] <= orderbyfit[0].fitness.values[0]:
                offspring[:] = [best_ind] + orderbyfit[:len(population) - 1]

        if neat_alg:
            SpeciesPunishment(offspring, params, neat_h)

        if SaveMatrix:
            if GenMatrix:
                idx_aux = np.searchsorted(Matrix[:, 0], gen)
                Matrix_specie[idx_aux, 1] = ind_specie(offspring)
            else:
                idx_aux = np.searchsorted(Matrix_specie[:, 0],
                                          funcEval.cont_evalp)
                try:
                    Matrix_specie[idx_aux, 1] = ind_specie(offspring)
                except IndexError:
                    Matrix_specie[-1, 1] = ind_specie(offspring)
            np.savetxt('./Specie/%s/specist_%d_%d.txt' %
                       (problem, num_p, n_corr),
                       Matrix_specie,
                       delimiter=";",
                       fmt="%s")

        population[:] = offspring  # population update

        cond_ind = 0
        cont_better = 0
        if funcEval.LS_flag:
            for ind in population:
                ind.LS_applied_set(0)

            if LS_select == 1:
                trees_h(population, num_p, n_corr, pset, direccion, problem,
                        benchmark_flag)
            elif LS_select == 2:
                best_specie(population, num_p, n_corr, pset, direccion,
                            problem, benchmark_flag)
            elif LS_select == 3:
                specie_h(population, num_p, n_corr, pset, direccion, problem,
                         benchmark_flag)
            elif LS_select == 4:
                best_pop_ls(population, num_p, n_corr, pset, direccion,
                            problem, benchmark_flag)
            elif LS_select == 5:
                all_pop(population, num_p, n_corr, pset, direccion, problem,
                        benchmark_flag)
            elif LS_select == 6:
                trees_h_wo(population, num_p, n_corr, pset, direccion, problem,
                           benchmark_flag)
            elif LS_select == 7:
                ls_bestset(population, num_p, n_corr, pset, direccion, problem,
                           benchmark_flag)
            elif LS_select == 8:
                ls_random(population, num_p, n_corr, pset, direccion, problem,
                          benchmark_flag)
            elif LS_select == 9:
                ls_randbestset(population, num_p, n_corr, pset, direccion,
                               problem, benchmark_flag)
            #
            invalid_ind = [ind for ind in population]
            new_invalid_ind = []
            for ind in population:
                strg = ind.__str__()
                l_strg = add_subt(strg, ind)
                c = tree2f()
                cd = c.convert(l_strg)
                new_invalid_ind.append(cd)
            fitness_ls = toolbox.map(toolbox.evaluate, new_invalid_ind)
            print 'Fitness comp.:',
            for ind, ls_fit in zip(invalid_ind, fitness_ls):
                if ind.LS_applied_get() == 1:
                    cond_ind += 1
                    if ind.fitness.values[0] < ls_fit:
                        print '-',
                    elif ind.fitness.values[0] > ls_fit:
                        cont_better += 1
                        print '+',
                    elif ind.fitness.values[0] == ls_fit:
                        print '=',
                funcEval.cont_evalp += 1
                pop_file.write(
                    '\n%s;%s;%s;%s;%s;%s;%s;%s' %
                    (gen, ind.fitness.values[0], ls_fit, ind.LS_applied_get(),
                     ind, [x for x in ind.get_parent()] if ind.get_parent()
                     is not None else [], [x for x in ind.get_params()
                                           ], ind.get_id()))
                ind.fitness.values = ls_fit
            print ''

        #for ind in population:

        else:

            for ind in population:
                pop_file.write('\n%s;%s;%s;%s;%s;%s;%s' %
                               (gen, ind.fitness.values[0],
                                ind.LS_applied_get(), ind.LS_story_get(),
                                ind.off_cx_get(), ind.off_mut_get(), ind))

        best_ind = best_pop(population)
        if funcEval.LS_flag:
            strg = best_ind.__str__()
            l_strg = add_subt(strg, best_ind)
            c = tree2f()
            cd = c.convert(l_strg)
            new_invalid_ind.append(cd)
            bestind.write('\n%s;%s;%s' % (gen, best_ind.fitness.values[0], cd))
            toolbox.map(toolbox.evaluate_best_t, [cd])
            if testing:
                fit_best = toolbox.map(toolbox.evaluate_best, [cd])
                best_ind.fitness_test.values = fit_best[0]

                best.write(
                    '\n%s;%s;%s;%s;%s;%s;%s' %
                    (gen, funcEval.cont_evalp, best_ind.fitness.values[0],
                     best_ind.LS_fitness_get(),
                     best_ind.fitness_test.values[0], len(best_ind),
                     avg_nodes(population)))
            else:
                best.write(
                    '\n%s;%s;%s;%s;%s;%s;%s' %
                    (gen, funcEval.cont_evalp, best_ind.fitness.values[0],
                     best_ind.LS_fitness_get(), None, len(best_ind),
                     avg_nodes(population)))
            out.write('\n%s;%s;%s;%s;%s;%s' %
                      (gen, len(best_ind), best_ind.LS_applied_get(),
                       best_ind.get_params(), cd, best_ind))
        else:
            toolbox.map(toolbox.evaluate_best_t, [best_ind])
            if testing:
                fitnesses_test = toolbox.map(toolbox.evaluate_best, [best_ind])
                best_ind.fitness_test.values = fitnesses_test[0]
                best.write(
                    '\n%s;%s;%s;%s;%s;%s' %
                    (gen, funcEval.cont_evalp, best_ind.fitness_test.values[0],
                     best_ind.fitness.values[0], len(best_ind),
                     avg_nodes(population)))
            else:
                best.write('\n%s;%s;%s;%s;%s;%s' %
                           (gen, funcEval.cont_evalp, None,
                            best_ind.fitness.values[0], len(best_ind),
                            avg_nodes(population)))
            bestind.write('\n%s;%s;%s' %
                          (gen, best_ind.fitness.values[0], best_ind))
            out.write('\n%s;%s;%s' % (gen, len(best_ind), best_ind))

        if funcEval.LS_flag:
            print 'Num. LS:', cond_ind
            print 'Ind. Improvement:', cont_better
            print 'Best Ind. LS:', best_ind.LS_applied_get()

        print 'Best Ind.:', best_ind
        print 'Best Fitness:', best_ind.fitness.values[0]
        if testing:
            print 'Test fitness:', best_ind.fitness_test.values[0]
        print 'Avg Nodes:', avg_nodes(population)
        print 'Evaluations: ', funcEval.cont_evalp

        if SaveMatrix:
            data_pop = avg_nodes(population)
            if GenMatrix:
                idx_aux = np.searchsorted(Matrix[:, 0], gen)
                Matrix[idx_aux, 1] = best_ind.fitness.values[0]
                if testing:
                    Matrix[idx_aux, 2] = best_ind.fitness_test.values[0]
                else:
                    Matrix[idx_aux, 2] = None
                Matrix[idx_aux, 3] = len(best_ind)
                Matrix[idx_aux, 4] = data_pop[0]
                Matrix[idx_aux, 5] = gen
                Matrix[idx_aux, 6] = 1
                Matrix[idx_aux, 7] = data_pop[1]  # max nodes
                Matrix[idx_aux, 8] = data_pop[2]  # min nodes
                if neat_alg:
                    Matrix[idx_aux, 9] = best_ind.get_specie()
                    Matrix[idx_aux, 10] = max(ind_specie(population),
                                              key=lambda x: x[0])[0]
            else:
                if funcEval.cont_evalp >= cont_evalf:
                    num_c -= 1
                    idx_aux = num_c
                    Matrix[num_c, 1] = best_ind.fitness.values[0]
                    if testing:
                        Matrix[num_c, 2] = best_ind.fitness_test.values[0]
                    else:
                        Matrix[num_c, 2] = None
                    Matrix[num_c, 3] = len(best_ind)
                    Matrix[num_c, 4] = data_pop[0]
                    Matrix[num_c, 5] = gen
                    Matrix[num_c, 6] = 1
                    Matrix[num_c, 7] = data_pop[1]  #max_nodes
                    Matrix[num_c, 8] = data_pop[2]  #min nodes
                    if neat_alg:
                        Matrix[num_c, 9] = max(ind_specie(population),
                                               key=lambda x: x[0])[0]

                else:
                    idx_aux = np.searchsorted(Matrix[:, 0],
                                              funcEval.cont_evalp)
                    Matrix[idx_aux, 1] = best_ind.fitness.values[0]
                    if testing:
                        Matrix[idx_aux, 2] = best_ind.fitness_test.values[0]
                    else:
                        Matrix[idx_aux, 2] = None
                    Matrix[idx_aux, 3] = len(best_ind)
                    Matrix[idx_aux, 4] = data_pop[0]
                    Matrix[idx_aux, 5] = gen
                    Matrix[idx_aux, 6] = 1
                    Matrix[idx_aux, 7] = data_pop[1]  #max nodes
                    Matrix[idx_aux, 8] = data_pop[2]  #min nodes
                    if neat_alg:
                        Matrix[idx_aux, 9] = max(ind_specie(population),
                                                 key=lambda x: x[0])[0]

                id_it = idx_aux - 1
                id_beg = 0
                flag = True
                flag2 = False
                while flag:
                    if Matrix[id_it, 6] == 0:
                        id_it -= 1
                        flag2 = True
                    else:
                        id_beg = id_it
                        flag = False
                if flag2:
                    x = Matrix[id_beg, 1:8]
                    Matrix[id_beg:idx_aux, 1:] = Matrix[id_beg, 1:]

            np.savetxt('./Matrix/%s/idx_%d_%d.txt' % (problem, num_p, n_corr),
                       Matrix,
                       delimiter=",",
                       fmt="%s")

        end_t = time.time()
        if testing:
            time_file.write(
                '\n%s;%s;%s;%s;%s;%s' %
                (gen, begin, end_t, str(round(end_t - begin, 2)),
                 best_ind.fitness.values[0], best_ind.fitness_test.values[0]))
        else:
            time_file.write('\n%s;%s;%s;%s;%s' %
                            (gen, begin, end_t, str(round(end_t - begin, 2)),
                             best_ind.fitness.values[0]))


#        import convert_order
#        convert_order.convert_(population)
    return population, logbook, best_ind.fitness.values[
        0], best_ind.fitness_test.values[0]
Ejemplo n.º 9
0
def GP_LS(population,
          toolbox,
          cxpb,
          mutpb,
          ngen,
          LS_flag,
          LS_select,
          cont_evalf,
          numsalto,
          SaveMatrix,
          GenMatrix,
          pset,
          n_corr,
          num_p,
          params,
          direccion,
          problem,
          testing,
          version,
          benchmark_flag,
          config,
          stats=None,
          halloffame=None,
          verbose=__debug__):
    """This algorithm reproduce the simplest evolutionary algorithm in GP with a local search.

    :param population: A list of individuals.
    :param toolbox: A :class:`~deap.base.Toolbox` that contains the evolution
                    operators.
    :param cxpb: The probability of mating two individuals.
    :param mutpb: The probability of mutating an individual.
    :param ngen: The number of generation.
    :param LS_flag: wheter or not to use LocalSearchGP
    :param LS_select: indicate the kind of selection to use the LSGP on the population.
    :param cont_evalf: maximum number of evaluations
    :param n_corr: run number just to write the txt file
    :param p: problem number just to write the txt file
    :param stats: A :class:`~deap.tools.Statistics` object that is updated
                  inplace, optional.
    :param halloffame: A :class:`~deap.tools.HallOfFame` object that will
                       contain the best individuals, optional.
    :param verbose: Whether or not to log the statistics.
    :returns: The final population.

    It uses :math:`\lambda = \kappa = \mu` and goes as follow.
    It first initializes the population (:math:`P(0)`) by evaluating
    every individual presenting an invalid fitness. Then, it enters the
    evolution loop that begins by the selection of the :math:`P(g+1)`
    population. Then the crossover operator is applied on a proportion of
    :math:`P(g+1)` according to the *cxpb* probability, the resulting and the
    untouched individuals are placed in :math:`P'(g+1)`. Thereafter, a
    proportion of :math:`P'(g+1)`, determined by *mutpb*, is
    mutated and placed in :math:`P''(g+1)`, the untouched individuals are
    transferred :math:`P''(g+1)`. Finally, those new individuals are evaluated
    and the evolution loop continues until *ngen* generations are completed.
    Briefly, the operators are applied in the following order ::

        evaluate(population)
        for i in range(ngen):
            offspring = select(population)
            offspring = mate(offspring)
            offspring = mutate(offspring)
            evaluate(offspring)
            population = offspring

    This function expects :meth:`toolbox.mate`, :meth:`toolbox.mutate`,
    :meth:`toolbox.select` and :meth:`toolbox.evaluate` aliases to be
    registered in the toolbox.

    .. Adapted from NGP-LS
    https://github.com/saarahy/NGP-LS
    """
    logbook = tools.Logbook()
    logbook.header = ['gen', 'nevals'] + (stats.fields if stats else [])

    #creating files to save data.
    d = './Results/%s/pop_file_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    pop_file = open(d, 'a')

    d = './Results/%s/bestind_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    best = open(d, 'w')  # save data

    d = './Results/%s/bestind_str_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    out = open(d, 'w')

    d = './Timing/%s/pop_file_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    time_file = open(d, 'w')

    d = './Results/%s/bestind_LStr_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    bestind = open(d, 'w')

    d = './Matrix/%s/' % (problem)
    ensure_dir(d)

    begin = time.time()

    if funcEval.LS_flag:
        for ind in population:
            sizep = len(ind) + 2
            param_ones = np.ones(sizep)
            param_ones[0] = 0
            ind.params_set(param_ones)

    # Evaluate the individuals with an invalid fitness
    invalid_ind = [ind for ind in population if not ind.fitness.valid]
    fitnesses = list(toolbox.map(toolbox.evaluate, invalid_ind))
    for ind, fit in zip(invalid_ind, fitnesses):
        funcEval.cont_evalp += 1
        ind.fitness.values = fit

    best_ind = best_pop(population)  # best individual of the population
    if testing:
        fitnesst_best = list(toolbox.map(toolbox.evaluate_test, [best_ind]))
        best_ind.fitness_test.values = fitnesst_best[0]
    if testing:
        best.write(
            '\n%s;%s;%s;%s;%s;%s' %
            (0, funcEval.cont_evalp, best_ind.fitness_test.values[0],
             best_ind.fitness.values[0], len(best_ind), avg_nodes(population)))
    else:
        best.write('\n%s;%s;%s;%s;%s;%s' %
                   (0, funcEval.cont_evalp, None, best_ind.fitness.values[0],
                    len(best_ind), avg_nodes(population)))

    data_pop = avg_nodes(population)

    if funcEval.LS_flag == 1:
        strg = best_ind.__str__()
        l_strg = add_subt_cf(strg, args=[])
        c = tree2f()
        cd = c.convert(l_strg)
        out.write('\n%s;%s;%s;%s;%s;%s' %
                  (0, len(best_ind), best_ind.LS_applied_get(),
                   best_ind.get_params(), cd, best_ind))
    else:
        out.write('\n%s;%s;%s' % (0, len(best_ind), best_ind))

    for ind in population:
        pop_file.write('\n%s;%s;%s' % (0, ind.fitness.values[0], ind))

    print('---- Generation %d -----' % (0))
    print('Problem: ', problem)
    print('Problem No.: ', num_p)
    print('Run No.: ', n_corr)
    print('Local Search:', funcEval.LS_flag)
    if funcEval.LS_flag:
        ls_type = ''
        if LS_select == 1:
            ls_type = 'LSHS'
        elif LS_select == 2:
            ls_type = 'Best-Pop'
        elif LS_select == 3:
            ls_type = 'All-Pop'
        elif LS_select == 4:
            ls_type = 'LSHS-test'
        elif LS_select == 5:
            ls_type = 'Best set'
        elif LS_select == 6:
            ls_type = 'Random set'
        elif LS_select == 7:
            ls_type = "Best-Random set"
        print('Local Search Heuristic: %s (%s)' % (LS_select, ls_type))
    print('Best Ind.:', best_ind)
    print('Best Fitness:', best_ind.fitness.values[0])
    if testing:
        print('Test fitness:', best_ind.fitness_test.values[0])
    print('Avg Nodes:', avg_nodes(population))
    print('Evaluations: ', funcEval.cont_evalp)
    end_t = time.time()

    # Begin the generational process
    for gen in range(1, ngen + 1):

        if not GenMatrix:
            if funcEval.cont_evalp > cont_evalf:
                break

        begin = time.time()
        print('---- Generation %d -----' % (gen))
        print('Problem: ', problem)
        print('Problem No.: ', num_p)
        print('Run No.: ', n_corr)
        print('Local Search:', funcEval.LS_flag)
        if funcEval.LS_flag:
            print('Local Search Heuristic: %s (%s)' % (LS_select, ls_type))
        best_ind = copy.deepcopy(best_pop(population))

        parents = toolbox.select(population, len(population))

        begin_cx = time.time()
        offspring = varOr(parents, toolbox, cxpb, mutpb)
        end_cx = time.time()
        #time_specie.write('\n%s;%s;%s;%s' % (0, begin_cx, end_cx, str(round(end_cx - begin_cx, 2))))

        invalid_ind = [ind for ind in offspring]
        if funcEval.LS_flag:
            new_invalid_ind = []
            for ind in invalid_ind:
                strg = ind.__str__()
                l_strg = add_subt(strg, ind)
                c = tree2f()
                cd = c.convert(l_strg)
                new_invalid_ind.append(cd)
            fitness_ls = toolbox.map(toolbox.evaluate, new_invalid_ind)
            for ind, ls_fit in zip(invalid_ind, fitness_ls):
                funcEval.cont_evalp += 1
                ind.fitness.values = ls_fit
        else:
            fitnesses = toolbox.map(toolbox.evaluate, invalid_ind)
            for ind, fit in zip(invalid_ind, fitnesses):
                funcEval.cont_evalp += 1
                ind.fitness.values = fit

        orderbyfit = sorted(offspring, key=lambda ind: ind.fitness.values)

        if best_ind.fitness.values[0] <= orderbyfit[0].fitness.values[0]:
            offspring[:] = [best_ind] + orderbyfit[:len(population) - 1]

        population[:] = offspring  # population update

        cond_ind = 0
        cont_better = 0
        if funcEval.LS_flag:
            for ind in population:
                ind.LS_applied_set(0)

            if LS_select == 1:
                trees_h(population, num_p, n_corr, pset, direccion, problem,
                        benchmark_flag)
            elif LS_select == 2:
                best_pop_ls(population, num_p, n_corr, pset, direccion,
                            problem, benchmark_flag)
            elif LS_select == 3:
                all_pop(population, num_p, n_corr, pset, direccion, problem,
                        benchmark_flag)
            elif LS_select == 4:
                trees_h_wo(population, num_p, n_corr, pset, direccion, problem,
                           benchmark_flag)
            elif LS_select == 5:
                ls_bestset(population, num_p, n_corr, pset, direccion, problem,
                           benchmark_flag)
            elif LS_select == 6:
                ls_random(population, num_p, n_corr, pset, direccion, problem,
                          benchmark_flag)
            elif LS_select == 7:
                ls_randbestset(population, num_p, n_corr, pset, direccion,
                               problem, benchmark_flag)
            #
            invalid_ind = [ind for ind in population]
            new_invalid_ind = []
            for ind in population:
                strg = ind.__str__()
                l_strg = add_subt(strg, ind)
                c = tree2f()
                cd = c.convert(l_strg)
                new_invalid_ind.append(cd)
            fitness_ls = toolbox.map(toolbox.evaluate, new_invalid_ind)
            print('Fitness comp.:'),
            for ind, ls_fit in zip(invalid_ind, fitness_ls):
                if ind.LS_applied_get() == 1:
                    cond_ind += 1
                    if ind.fitness.values[0] < ls_fit:
                        print('-'),
                    elif ind.fitness.values[0] > ls_fit:
                        cont_better += 1
                        print('+'),
                    elif ind.fitness.values[0] == ls_fit:
                        print('='),
                funcEval.cont_evalp += 1
                pop_file.write(
                    '\n%s;%s;%s;%s;%s;%s;%s;%s' %
                    (gen, ind.fitness.values[0], ls_fit, ind.LS_applied_get(),
                     ind, [], [x for x in ind.get_params()], []))
                ind.fitness.values = ls_fit
            print('')

        #for ind in population:

        else:

            for ind in population:
                pop_file.write('\n%s;%s;%s;%s;%s;%s;%s' %
                               (gen, ind.fitness.values[0],
                                ind.LS_applied_get(), ind.LS_story_get(),
                                ind.off_cx_get(), ind.off_mut_get(), ind))

        best_ind = best_pop(population)
        if funcEval.LS_flag:
            strg = best_ind.__str__()
            l_strg = add_subt(strg, best_ind)
            c = tree2f()
            cd = c.convert(l_strg)
            new_invalid_ind.append(cd)
            bestind.write('\n%s;%s;%s' % (gen, best_ind.fitness.values[0], cd))
            toolbox.map(toolbox.evaluate_best_t, [cd])
            if testing:
                fit_best = list(toolbox.map(toolbox.evaluate_best, [cd]))
                best_ind.fitness_test.values = fit_best[0]

                best.write(
                    '\n%s;%s;%s;%s;%s;%s;%s' %
                    (gen, funcEval.cont_evalp, best_ind.fitness.values[0],
                     best_ind.LS_fitness_get(),
                     best_ind.fitness_test.values[0], len(best_ind),
                     avg_nodes(population)))
            else:
                best.write(
                    '\n%s;%s;%s;%s;%s;%s;%s' %
                    (gen, funcEval.cont_evalp, best_ind.fitness.values[0],
                     best_ind.LS_fitness_get(), None, len(best_ind),
                     avg_nodes(population)))
            out.write('\n%s;%s;%s;%s;%s;%s' %
                      (gen, len(best_ind), best_ind.LS_applied_get(),
                       best_ind.get_params(), cd, best_ind))
        else:
            toolbox.map(toolbox.evaluate_best_t, [best_ind])
            if testing:
                fitnesses_test = toolbox.map(toolbox.evaluate_best, [best_ind])
                best_ind.fitness_test.values = fitnesses_test[0]
                best.write(
                    '\n%s;%s;%s;%s;%s;%s' %
                    (gen, funcEval.cont_evalp, best_ind.fitness_test.values[0],
                     best_ind.fitness.values[0], len(best_ind),
                     avg_nodes(population)))
            else:
                best.write('\n%s;%s;%s;%s;%s;%s' %
                           (gen, funcEval.cont_evalp, None,
                            best_ind.fitness.values[0], len(best_ind),
                            avg_nodes(population)))
            bestind.write('\n%s;%s;%s' %
                          (gen, best_ind.fitness.values[0], best_ind))
            out.write('\n%s;%s;%s' % (gen, len(best_ind), best_ind))

        if funcEval.LS_flag:
            print('Num. LS:', cond_ind)
            print('Ind. Improvement:', cont_better)
            print('Best Ind. LS:', best_ind.LS_applied_get())

        print('Best Ind.:', best_ind)
        print('Best Fitness:', best_ind.fitness.values[0])
        if testing:
            print('Test fitness:', best_ind.fitness_test.values[0])
        print('Avg Nodes:', avg_nodes(population))
        print('Evaluations: ', funcEval.cont_evalp)

        end_t = time.time()
        if testing:
            time_file.write(
                '\n%s;%s;%s;%s;%s;%s' %
                (gen, begin, end_t, str(round(end_t - begin, 2)),
                 best_ind.fitness.values[0], best_ind.fitness_test.values[0]))
        else:
            time_file.write('\n%s;%s;%s;%s;%s' %
                            (gen, begin, end_t, str(round(end_t - begin, 2)),
                             best_ind.fitness.values[0]))
    return population, logbook
Ejemplo n.º 10
0
def neat_GP_LS(population, toolbox, cxpb, mutpb, ngen, neat_alg, neat_cx, neat_h,neat_pelit, LS_flag, LS_select, cont_evalf,
               num_salto, SaveMatrix, GenMatrix, pset,n_corr, num_p, params, direccion, problem, testing, version,
               stats=None, halloffame=None, verbose=__debug__):
    """This algorithm reproduce the simplest evolutionary algorithm as
    presented in chapter 7 of [Back2000]_.

    :param population: A list of individuals.
    :param toolbox: A :class:`~deap.base.Toolbox` that contains the evolution
                    operators.
    :param cxpb: The probability of mating two individuals.
    :param mutpb: The probability of mutating an individual.
    :param ngen: The number of generation.
    :param neat_alg: wheter or not to use species stuff.
    :param neat_cx: wheter or not to use neat-crossover
    :param neat_h: indicate the distance allowed between each specie
    :param neat_pelit: probability of being elitist, it's used in the neat cx and mutation
    :param LS_flag: wheter or not to use GP-LS
    :param LS_select: indicate the heuristic to select individuals for the LS method.
            LS_select == 1:'LSHS'            -- Heurisitic proposed by Z-Flores (2015)
            LS_select == 2:'Best-Sp'         -- Best individual of the specie
            LS_select == 3:'LSHS-Sp'         -- Heuristic 1 in each specie
            LS_select == 4:'Best-Pop'        -- Best ind of the population
            LS_select == 5:'All-Pop'         -- All population
            LS_select == 6:'LSHS-test'       -- Heuristic 1 only in the test data
            LS_select == 7:'Best set'        -- Best individuals set
            LS_select == 8:'Random set'      -- Random individuals set
            LS_select == 9:'Best-Random set' -- Best individual set and then random selection of this individuals
    :param cont_evalf: Counter of the function evaluations
    :param num_salto: This is a parameter to save data in the Matrix, indicates the range for each row
    :param SaveMatrix: wheter or not to save data on a Matrix
    :param GenMatrix: wheter or not to save data in the Matrix on a generations form or function evals form
    :param pset: this is the primitive set
    :param n_corr: run number just to write the txt file
    :param num_p: problem number just to write the txt file
    :param params:indicate the params for the fitness sharing, the different
                    options are:
                    -DontPenalize(str): 'best_specie' or 'best_of_each_specie'
                    -Penalization_method(int):
                        1.without penalization
                        2.penalization fitness sharing
                        3.new penalization
                    -ShareFitness(str): 'yes' or 'no'
    :param stats: A :class:`~deap.tools.Statistics` object that is updated
                  inplace, optional.
    :param direccion: Path to find the training data. This parameter works in the g_address.py file.
    :param problem: Name of the problem. Save and Get data.
    :param testing: Wheter or not to use testing data.
    :param version: Number of the version to make the speciation.
    :param halloffame: A :class:`~deap.tools.HallOfFame` object that will
                       contain the best individuals, optional.
    :param verbose: Whether or not to log the statistics.
    :returns: The final population.

    It uses :math:`\lambda = \kappa = \mu` and goes as follow.
    It first initializes the population (:math:`P(0)`) by evaluating
    every individual presenting an invalid fitness. Then, it enters the
    evolution loop that begins by the selection of the :math:`P(g+1)`
    population. Then the crossover operator is applied on a proportion of
    :math:`P(g+1)` according to the *cxpb* probability, the resulting and the
    untouched individuals are placed in :math:`P'(g+1)`. Thereafter, a
    proportion of :math:`P'(g+1)`, determined by *mutpb*, is
    mutated and placed in :math:`P''(g+1)`, the untouched individuals are
    transferred :math:`P''(g+1)`. Finally, those new individuals are evaluated
    and the evolution loop continues until *ngen* generations are completed.
    Briefly, the operators are applied in the following order ::

        evaluate(population)
        for i in range(ngen):
            offspring = select(population)
            offspring = mate(offspring)
            offspring = mutate(offspring)
            evaluate(offspring)
            population = offspring

    This function expects :meth:`toolbox.mate`, :meth:`toolbox.mutate`,
    :meth:`toolbox.select` and :meth:`toolbox.evaluate` aliases to be
    registered in the toolbox.

    .. [Back2000] Back, Fogel and Michalewicz, "Evolutionary Computation 1 :
       Basic Algorithms and Operators", 2000.
    """

    logbook = tools.Logbook()
    logbook.header = ['gen', 'nevals'] + (stats.fields if stats else [])

    # creating files to save data.
    d = './Results/%s/pop_file_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    pop_file = open(d, 'a')

    d = './Results/%s/bestind_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    best = open(d, 'w')

    d = './Results/%s/bestind_str_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    out = open(d, 'w')

    d='./Timing/%s/pop_file_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    time_file = open(d, 'w')

    d = './Timing/%s/specie_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    time_specie = open(d, 'w')

    d = './Timing/%s/cruce_s_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    time_cx = open(d, 'w')

    d = './Specie/%s/specieind_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    specie_file = open(d, 'w')

    d = './Specie/%s/specist_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    specie_statis = open(d, 'w')

    d = './Results/%s/bestind_LStr_%d_%d.txt' % (problem, num_p, n_corr)
    ensure_dir(d)
    bestind = open(d, 'w')

    begin = time.time()

    if SaveMatrix:  # Saving data in matrix
        num_r = 11
        if GenMatrix:  # To save the data in matrix by generations
            num_salto = 1
            num_c = ngen+1
            Matrix = np.empty((num_c, num_r,))
            vector = np.arange(0, num_c, num_salto)
        else:  # To save the data in matrix by function evaluations
            num_c = (cont_evalf/num_salto) + 1
            Matrix = np.empty((num_c, num_r,))
            vector = np.arange(1, cont_evalf+num_salto, num_salto)
        for i in range(len(vector)):
            Matrix[i, 0] = vector[i]
        Matrix[:, 6] = 0.

    begin_sp = time.time()  # Saving data: speciaton init timing
    if neat_alg:  # before version 3
        for ind in population:
            level_info = level_node(ind)
            ind.nodefeat_set(level_info)

        species(population, neat_h, version)

    end_sp = time.time() # Saving data: speciaton end timing
    time_specie.write('\n%s;%s;%s;%s' % (0, begin_sp, end_sp, str(round(end_sp - begin_sp, 2))))  # Saving data in file

    specie_file.write('\n-------')
    for ind in population:
        specie_file.write('\n%s;%s;%s' % (ind.get_specie(),ind,version))

    # Adding parameters for the LS algorithms
    if funcEval.LS_flag:
        for ind in population:
            sizep = len(ind)+2
            param_ones = np.ones(sizep)
            param_ones[0] = 0
            ind.params_set(param_ones)

    # Evaluate the individuals with an invalid fitness
    invalid_ind = [ind for ind in population if not ind.fitness.valid]
    fitnesses = toolbox.map(toolbox.evaluate, invalid_ind)
    for ind, fit in zip(invalid_ind, fitnesses):
        funcEval.cont_evalp += 1
        ind.fitness.values = fit

    best_ind = best_pop(population)  # best individual of the population

    if testing:
        fitnesst_best = toolbox.map(toolbox.evaluate_test, [best_ind])
        best_ind.fitness_test.values = fitnesst_best[0]

        best.write('\n%s;%s;%s;%s;%s;%s' % (0, funcEval.cont_evalp, best_ind.fitness_test.values[0], best_ind.fitness.values[0], len(best_ind), avg_nodes(population)))
    else:
        best.write('\n%s;%s;%s;%s;%s;%s' % (
        0, funcEval.cont_evalp, None, best_ind.fitness.values[0], len(best_ind),
        avg_nodes(population)))

    data_pop = avg_nodes(population)

    if SaveMatrix:
        idx = 0
        Matrix[idx, 1] = best_ind.fitness.values[0]                           # best ind train data
        if testing:
            Matrix[idx, 2] = best_ind.fitness_test.values[0]                  # best ind test data
        else:
            Matrix[idx, 2] = None
        Matrix[idx, 3] = len(best_ind)                                        # best individual size (number of nodes)
        Matrix[idx, 4] = data_pop[0]                                          # avg size of the population
        Matrix[idx, 5] = 0.
        Matrix[idx, 6] = 1                                                    # id to know if the current row is full
        Matrix[idx, 7] = data_pop[1]                                          # max size of the population
        Matrix[idx, 8] = data_pop[2]                                          # min size of the population
        Matrix[idx, 9] = best_ind.get_specie()                                # specie of the best individual
        Matrix[idx, 10] = max(ind_specie(population), key=lambda x: x[0])[0]  # max number of species on current pop

        np.savetxt('./Matrix/%s/idx_%d_%d.txt' % (problem,num_p, n_corr), Matrix, delimiter=",", fmt="%s")

    if neat_alg:
        SpeciesPunishment(population,params,neat_h)

    if funcEval.LS_flag == 1:
        strg = best_ind.__str__()
        l_strg = add_subt_cf(strg, args=[])
        c = tree2f()
        cd = c.convert(l_strg)
        out.write('\n%s;%s;%s;%s;%s;%s' % (0, len(best_ind), best_ind.LS_applied_get(), best_ind.get_params(), cd, best_ind))
    else:
        out.write('\n%s;%s;%s' % (0, len(best_ind), best_ind))

    for ind in population:
        pop_file.write('\n%s;%s'%(ind.fitness.values[0], ind))

    ls_type = ''
    if LS_select == 1:
        ls_type = 'LSHS'
    elif LS_select == 2:
        ls_type = 'Best-Sp'
    elif LS_select == 3:
        ls_type = 'LSHS-Sp'
    elif LS_select == 4:
        ls_type = 'Best-Pop'
    elif LS_select == 5:
        ls_type = 'All-Pop'
    elif LS_select == 6:
        ls_type = 'LSHS-test'
    elif LS_select == 7:
        ls_type = 'Best set'
    elif LS_select == 8:
        ls_type = 'Random set'
    elif LS_select == 9:
        ls_type = 'Best-Random set'

    ################################
    # Printing data
    ################################
    print '---- Generation %d -----' % (0)
    print 'Problem: ', problem
    print 'Problem No.: ', num_p
    print 'Run No.: ', n_corr
    print 'neat-GP:', neat_alg
    print 'neat-cx:', neat_cx
    print 'Local Search:', funcEval.LS_flag
    if funcEval.LS_flag:
        print 'Local Search Heuristic: %s (%s)' % (LS_select,ls_type)
    print 'Best Ind.:', best_ind
    print 'Best Fitness:', best_ind.fitness.values[0]
    if testing:
        print 'Test fitness:', best_ind.fitness_test.values[0]
    print 'Avg Nodes:', avg_nodes(population)
    print 'Evaluations: ', funcEval.cont_evalp
    ##################################

    end_t = time.time()

    specie_statis.write('\n%s;%s' % (0, count_species(population)))

    if testing:
        time_file.write('\n%s;%s;%s;%s;%s;%s' % (0, begin,end_t, str(round(end_t - begin, 2)), best_ind.fitness.values[0], best_ind.fitness_test.values[0]))
    else:
        time_file.write('\n%s;%s;%s;%s;%s' % (
        0, begin, end_t, str(round(end_t - begin, 2)), best_ind.fitness.values[0]))

    # Begin the generational process
    for gen in range(1, ngen+1):

        if not GenMatrix:
            if funcEval.cont_evalp > cont_evalf:
                break

        begin = time.time()
        print '---- Generation %d -----' % (gen)
        print 'Problem: ', problem
        print 'Problem No.: ', num_p
        print 'Run No.: ', n_corr
        print 'neat-GP:', neat_alg
        print 'neat-cx:', neat_cx
        print 'Local Search:', funcEval.LS_flag
        if funcEval.LS_flag:
            print 'Local Search Heuristic: %s (%s)' % (LS_select, ls_type)

        best_ind = copy.deepcopy(best_pop(population))

        if neat_alg:
            parents = p_selection(population, neat_pelit)
        else:
            parents = toolbox.select(population, len(population))

        begin_cx = time.time()
        if neat_cx and neat_alg:  # neat-Crossover
            n = len(parents)
            mut = 1
            cx = 1
            offspring = neatGP(toolbox, parents, cxpb, mutpb, n, mut, cx, neat_pelit)
        else:
            offspring = varOr(parents, toolbox, cxpb, mutpb)
        end_cx = time.time()
        time_specie.write('\n%s;%s;%s;%s' % (0, begin_cx, end_cx, str(round(end_sp - begin_sp, 2))))

        if neat_alg:  # Speciation of the descendants and fitness evaluation
            begin_sp = time.time()
            specie_parents_child(parents,offspring, neat_h, version)
            offspring[:] = parents+offspring

            invalid_ind = [ind for ind in offspring if not ind.fitness.valid]
            fitnesses = toolbox.map(toolbox.evaluate, invalid_ind)

            for ind, fit in zip(invalid_ind, fitnesses):
                funcEval.cont_evalp += 1
                ind.fitness.values = fit

            end_sp = time.time()
            time_specie.write('\n%s;%s;%s;%s' % (gen, begin_sp, end_sp, str(round(end_sp - begin_sp, 2))))

        else:  # evaluation of descendants
            invalid_ind = [ind for ind in offspring]
            if funcEval.LS_flag:  # evaluating with the parameters set
                new_invalid_ind = []
                for ind in invalid_ind:
                    strg = ind.__str__()
                    l_strg = add_subt(strg, ind)
                    c = tree2f()
                    cd = c.convert(l_strg)
                    new_invalid_ind.append(cd)
                fitness_ls = toolbox.map(toolbox.evaluate, new_invalid_ind)
                for ind, ls_fit in zip(invalid_ind, fitness_ls):
                    funcEval.cont_evalp += 1
                    ind.fitness.values = ls_fit
            else:  # normal evaluation
                fitnesses = toolbox.map(toolbox.evaluate, invalid_ind)
                for ind, fit in zip(invalid_ind, fitnesses):
                    funcEval.cont_evalp += 1
                    ind.fitness.values = fit

            orderbyfit = sorted(offspring, key=lambda ind:ind.fitness.values)
            #print len(orderbyfit),len(best_ind)

            # Applying elitism: Best individual survives
            if best_ind.fitness.values[0] <= orderbyfit[0].fitness.values[0]:
                offspring[:] = [best_ind]+orderbyfit[:len(population)-1]

        if neat_alg:
            SpeciesPunishment(offspring, params, neat_h)

        specie_statis.write('\n%s;%s' % (gen, ind_specie(offspring)))
        population[:] = offspring  # population update

        specie_file.write('\n%s--------------------------------' % gen)
        if neat_alg:
            for ind in population:
                specie_file.write('\n%s;%s;%s' % (ind.get_specie(), ind, version))

        cond_ind = 0
        cont_better = 0

        if funcEval.LS_flag:
            for ind in population:
                ind.LS_applied_set(0)

            if   LS_select == 1:
                trees_h(population, num_p, n_corr,  pset, direccion, problem)
            elif LS_select == 2:
                best_specie(population, num_p, n_corr, pset, direccion, problem)
            elif LS_select == 3:
                specie_h(population, num_p, n_corr, pset, direccion, problem)
            elif LS_select == 4:
                best_pop_ls(population, num_p, n_corr, pset, direccion, problem)
            elif LS_select == 5:
                all_pop(population, num_p, n_corr, pset, direccion, problem)
            elif LS_select == 6:
                trees_h_wo(population, num_p, n_corr, pset, direccion, problem)
            elif LS_select == 7:
                ls_bestset(population, num_p, n_corr, pset, direccion, problem)
            elif LS_select == 8:
                ls_random(population, num_p, n_corr, pset, direccion, problem)
            elif LS_select == 9:
                ls_randbestset(population, num_p, n_corr, pset, direccion, problem)


            invalid_ind = [ind for ind in population]
            new_invalid_ind = []
            for ind in population:
                strg = ind.__str__()
                l_strg = add_subt(strg, ind)
                c = tree2f()
                cd = c.convert(l_strg)
                new_invalid_ind.append(cd)
            fitness_ls = toolbox.map(toolbox.evaluate, new_invalid_ind)
            print 'Fitness comp.:',

            # Information about how goes the optimization
            for ind, ls_fit in zip(invalid_ind, fitness_ls):
                if ind.LS_applied_get() == 1:
                    cond_ind += 1
                    if ind.fitness.values[0] < ls_fit:
                        print '-',
                    elif ind.fitness.values[0] > ls_fit:
                        cont_better += 1
                        print '+',
                    elif ind.fitness.values[0] == ls_fit:
                        print '=',
                funcEval.cont_evalp += 1
                ind.fitness.values = ls_fit
            print ''

            pop_file.write('\n----------------------------------------%s'%(gen))
            for ind in population:
                pop_file.write('\n%s;%s;%s;%s'%(ind.LS_applied_get(),ind.fitness.values[0], ind, [x for x in ind.get_params()]))
        else:
            pop_file.write('\n----------------------------------------')
            for ind in population:
                pop_file.write('\n%s;%s;%s;%s;%s;%s'%(ind.LS_applied_get(),ind.LS_story_get(),ind.off_cx_get(),ind.off_mut_get(),ind.fitness.values[0], ind))

        best_ind = best_pop(population)
        if funcEval.LS_flag:
            strg = best_ind.__str__()
            l_strg = add_subt(strg, best_ind)
            c = tree2f()
            cd = c.convert(l_strg)
            new_invalid_ind.append(cd)
            bestind.write('\n%s;%s;%s' % (gen, best_ind.fitness.values[0],cd))
            if testing:
                fit_best = toolbox.map(toolbox.evaluate_test, [cd])
                best_ind.fitness_test.values = fit_best[0]

                best.write('\n%s;%s;%s;%s;%s;%s;%s' % (gen, funcEval.cont_evalp,  best_ind.fitness.values[0], best_ind.LS_fitness_get(), best_ind.fitness_test.values[0], len(best_ind), avg_nodes(population)))
            else:
                best.write('\n%s;%s;%s;%s;%s;%s;%s' % (
                gen, funcEval.cont_evalp, best_ind.fitness.values[0], best_ind.LS_fitness_get(),
                None, len(best_ind), avg_nodes(population)))
            out.write('\n%s;%s;%s;%s;%s;%s' % (gen, len(best_ind), best_ind.LS_applied_get(), best_ind.get_params(), cd, best_ind))

        else:
            if testing:
                fitnesses_test = toolbox.map(toolbox.evaluate_test, [best_ind])
                best_ind.fitness_test.values = fitnesses_test[0]
                best.write('\n%s;%s;%s;%s;%s;%s' % (gen, funcEval.cont_evalp, best_ind.fitness_test.values[0], best_ind.fitness.values[0], len(best_ind), avg_nodes(population)))
            else:
                best.write('\n%s;%s;%s;%s;%s;%s' % (
                gen, funcEval.cont_evalp, None, best_ind.fitness.values[0], len(best_ind),
                avg_nodes(population)))

            bestind.write('\n%s;%s;%s' % (gen, best_ind.fitness.values[0], best_ind))
            out.write('\n%s;%s;%s' % (gen, len(best_ind), best_ind))

        if funcEval.LS_flag:
            print 'Num. LS:', cond_ind
            print 'Ind. Improvement:', cont_better
            print 'Best Ind. LS:', best_ind.LS_applied_get()

        print 'Best Ind.:', best_ind
        print 'Best Fitness:', best_ind.fitness.values[0]
        if testing:
            print 'Test fitness:',best_ind.fitness_test.values[0]
        print 'Avg Nodes:', avg_nodes(population)
        print 'Evaluations: ', funcEval.cont_evalp

        if SaveMatrix:
            data_pop=avg_nodes(population)
            if GenMatrix:
                idx_aux = np.searchsorted(Matrix[:, 0], gen)
                Matrix[idx_aux, 1] = best_ind.fitness.values[0]
                if testing:
                    Matrix[idx_aux, 2] = best_ind.fitness_test.values[0]
                else:
                    Matrix[idx_aux, 2] = None
                Matrix[idx_aux, 3] = len(best_ind)
                Matrix[idx_aux, 4] = data_pop[0]
                Matrix[idx_aux, 5] = gen
                Matrix[idx_aux, 6] = 1
                Matrix[idx_aux, 7] = data_pop[1]  # max nodes
                Matrix[idx_aux, 8] = data_pop[2]  # min nodes
                if neat_alg:
                    Matrix[idx_aux, 9] = best_ind.get_specie()
                    Matrix[idx_aux, 10] = max(ind_specie(population), key=lambda x: x[0])[0]
            else:
                if funcEval.cont_evalp >= cont_evalf:
                    num_c -= 1
                    idx_aux = num_c
                    Matrix[num_c, 1] = best_ind.fitness.values[0]
                    if testing:
                        Matrix[num_c, 2] = best_ind.fitness_test.values[0]
                    else:
                        Matrix[num_c, 2] = None
                    Matrix[num_c, 3] = len(best_ind)
                    Matrix[num_c, 4] = data_pop[0]
                    Matrix[num_c, 5] = gen
                    Matrix[num_c, 6] = 1
                    Matrix[num_c, 7] = data_pop[1]  # max_nodes
                    Matrix[num_c, 8] = data_pop[2]  # min nodes
                    if neat_alg:
                        Matrix[num_c, 9] = max(ind_specie(population), key=lambda x: x[0])[0]

                else:
                    idx_aux = np.searchsorted(Matrix[:, 0], funcEval.cont_evalp)
                    Matrix[idx_aux, 1] = best_ind.fitness.values[0]
                    if testing:
                        Matrix[idx_aux, 2] = best_ind.fitness_test.values[0]
                    else:
                        Matrix[idx_aux, 2] = None
                    Matrix[idx_aux, 3] = len(best_ind)
                    Matrix[idx_aux, 4] = data_pop[0]
                    Matrix[idx_aux, 5] = gen
                    Matrix[idx_aux, 6] = 1
                    Matrix[idx_aux, 7] = data_pop[1]  # max nodes
                    Matrix[idx_aux, 8] = data_pop[2]  # min nodes
                    if neat_alg:
                        Matrix[idx_aux, 9] = max(ind_specie(population), key=lambda x: x[0])[0]

                id_it = idx_aux-1
                id_beg = 0
                flag = True
                flag2 = False
                while flag:
                    if Matrix[id_it, 6] == 0:
                        id_it -= 1
                        flag2 = True
                    else:
                        id_beg = id_it
                        flag = False
                if flag2:
                    x = Matrix[id_beg, 1:8]
                    Matrix[id_beg:idx_aux, 1:] = Matrix[id_beg, 1:]

            np.savetxt('./Matrix/%s/idx_%d_%d.txt' % (problem, num_p, n_corr), Matrix, delimiter=",", fmt="%s")

        end_t = time.time()
        if testing:
            time_file.write('\n%s;%s;%s;%s;%s;%s' % (gen, begin, end_t, str(round(end_t - begin, 2)), best_ind.fitness.values[0], best_ind.fitness_test.values[0]))
        else:
            time_file.write('\n%s;%s;%s;%s;%s' % (
            gen, begin, end_t, str(round(end_t - begin, 2)), best_ind.fitness.values[0]))
    return population, logbook
Ejemplo n.º 11
0
def evo_neat_GP_LS(population, toolbox, cxpb, mutpb, ngen, neat_alg, neat_cx, neat_h,neat_pelit, LS_flag, LS_select, cont_evalf, num_salto, SaveMatrix, GenMatrix, pset,n_corr, num_p, params, direccion, problem,stats=None,
             halloffame=None, verbose=__debug__):
    """This algorithm reproduce the simplest evolutionary algorithm as
    presented in chapter 7 of [Back2000]_.

    :param population: A list of individuals.
    :param toolbox: A :class:`~deap.base.Toolbox` that contains the evolution
                    operators.
    :param cxpb: The probability of mating two individuals.
    :param mutpb: The probability of mutating an individual.
    :param ngen: The number of generation.
    :param neat_alg: wheter or not to use species stuff.
    :param neat_cx: wheter or not to use neatGP cx
    :param neat_h: indicate the distance allowed between each specie
    :param neat_pelit: probability of being elitist, it's used in the neat cx and mutation
    :param LS_flag: wheter or not to use LocalSearchGP
    :param LS_select: indicate the kind of selection to use the LSGP on the population.
    :param cont_evalf: contador maximo del numero de evaluaciones
    :param n_corr: run number just to wirte the txt file
    :param p: problem number just to wirte the txt file
    :param params:indicate the params for the fitness sharing, the diffetent
                    options are:
                    -DontPenalize(str): 'best_specie' or 'best_of_each_specie'
                    -Penalization_method(int):
                        1.without penalization
                        2.penalization fitness sharing
                        3.new penalization
                    -ShareFitness(str): 'yes' or 'no'
    :param stats: A :class:`~deap.tools.Statistics` object that is updated
                  inplace, optional.
    :param halloffame: A :class:`~deap.tools.HallOfFame` object that will
                       contain the best individuals, optional.
    :param verbose: Whether or not to log the statistics.
    :returns: The final population.

    It uses :math:`\lambda = \kappa = \mu` and goes as follow.
    It first initializes the population (:math:`P(0)`) by evaluating
    every individual presenting an invalid fitness. Then, it enters the
    evolution loop that begins by the selection of the :math:`P(g+1)`
    population. Then the crossover operator is applied on a proportion of
    :math:`P(g+1)` according to the *cxpb* probability, the resulting and the
    untouched individuals are placed in :math:`P'(g+1)`. Thereafter, a
    proportion of :math:`P'(g+1)`, determined by *mutpb*, is
    mutated and placed in :math:`P''(g+1)`, the untouched individuals are
    transferred :math:`P''(g+1)`. Finally, those new individuals are evaluated
    and the evolution loop continues until *ngen* generations are completed.
    Briefly, the operators are applied in the following order ::

        evaluate(population)
        for i in range(ngen):
            offspring = select(population)
            offspring = mate(offspring)
            offspring = mutate(offspring)
            evaluate(offspring)
            population = offspring

    This function expects :meth:`toolbox.mate`, :meth:`toolbox.mutate`,
    :meth:`toolbox.select` and :meth:`toolbox.evaluate` aliases to be
    registered in the toolbox.

    .. [Back2000] Back, Fogel and Michalewicz, "Evolutionary Computation 1 :
       Basic Algorithms and Operators", 2000.
    """

    logbook = tools.Logbook()
    logbook.header = ['gen', 'nevals'] + (stats.fields if stats else [])
    pop_file = open('./Results/%s/pop_file_%d_%d.txt' % (problem, num_p, n_corr), 'a')

    if SaveMatrix:  # Saving data in matrix
        num_r = 9
        if GenMatrix:
            num_salto=1
            num_c=ngen+1
            Matrix= np.empty((num_c, num_r,))
            vector = np.arange(0, num_c, num_salto)
        else:
            num_c = (cont_evalf/num_salto) + 1
            Matrix = np.empty((num_c, num_r,))
            vector = np.arange(1, cont_evalf+num_salto, num_salto)
        for i in range(len(vector)):
            Matrix[i, 0] = vector[i]
            #num_r-1
        Matrix[:, 6] = 0.

    #Creation of the species
    # if neat_alg:
    #     species(population,neat_h)
    #     ind_specie(population)

    if funcEval.LS_flag:
        for ind in population:
            sizep = len(ind)+2
            param_ones = np.ones(sizep)
            param_ones[0] = 0
            ind.params_set(param_ones)

    # Evaluate the individuals with an invalid fitness
    invalid_ind = [ind for ind in population if not ind.fitness.valid]
    fitnesses = toolbox.map(toolbox.evaluate, invalid_ind)
    for ind, fit in zip(invalid_ind, fitnesses):
        funcEval.cont_evalp += 1
        ind.fitness.values = fit

    best = open('./Results/%s/bestind_%d_%d.txt' % (problem, num_p, n_corr), 'a')  # save data

    best_ind = best_pop(population)  # best individual of the population
    fitnesst_best = toolbox.map(toolbox.evaluate_test, [best_ind])
    best_ind.fitness_test.values = fitnesst_best[0]
    best.write('\n%s;%s;%s;%s;%s;%s' % (0, funcEval.cont_evalp, best_ind.fitness_test.values[0], best_ind.fitness.values[0], len(best_ind), avg_nodes(population)))
    data_pop=avg_nodes(population)
    if SaveMatrix:
        idx = 0
        Matrix[idx, 1] = best_ind.fitness.values[0]
        Matrix[idx, 2] = best_ind.fitness_test.values[0]
        Matrix[idx, 3] = len(best_ind)
        Matrix[idx, 4] = data_pop[0]
        Matrix[idx, 5] = 0.
        Matrix[idx, 6] = 1  # just an id to know if the current row is full
        Matrix[idx, 7] = data_pop[1]  # max size
        Matrix[idx, 8] = data_pop[2]  # min size

        np.savetxt('./Matrix/%s/idx_%d_%d.txt' % (problem,num_p, n_corr), Matrix, delimiter=",", fmt="%s")

    if neat_alg:
        SpeciesPunishment(population,params,neat_h)

    out = open('./Results/%s/bestind_str_%d_%d.txt' % (problem, num_p, n_corr), 'a')

    if funcEval.LS_flag == 1:
        strg = best_ind.__str__()
        l_strg = add_subt_cf(strg, args=[])
        c = tree2f()
        cd = c.convert(l_strg)
        out.write('\n%s;%s;%s;%s;%s;%s' % (0, len(best_ind), best_ind.LS_applied_get(), best_ind.get_params(), cd, best_ind))
    else:
        out.write('\n%s;%s;%s' % (0, len(best_ind), best_ind))

    for ind in population:
        pop_file.write('\n%s;%s'%(ind.fitness.values[0], ind))

    ls_type = ''
    if LS_select == 1:
        ls_type = 'LSHS'
    elif LS_select == 2:
        ls_type = 'Best-Sp'
    elif LS_select == 3:
        ls_type = 'LSHS-Sp'
    elif LS_select == 4:
        ls_type = 'Best-Pop'
    elif LS_select == 5:
        ls_type = 'All-Pop'
    elif LS_select == 6:
        ls_type = 'LSHS-test'
    elif LS_select == 7:
        ls_type = 'Best set'
    elif LS_select == 8:
        ls_type = 'Random set'
    elif LS_select == 9:
        ls_type = "Best-Random set"

    print '---- Generation %d -----' % (0)
    print 'Problem: ', problem
    print 'Problem No.: ', num_p
    print 'Run No.: ', n_corr
    print 'neat-GP:', neat_alg
    print 'neat-cx:', neat_cx
    print 'Local Search:', funcEval.LS_flag
    if funcEval.LS_flag:
        print 'Local Search Heuristic: %s (%s)' % (LS_select,ls_type)
    print 'Best Ind.:', best_ind
    print 'Best Fitness:', best_ind.fitness.values[0]
    print 'Test fitness:',best_ind.fitness_test.values[0]
    print 'Avg Nodes:', avg_nodes(population)
    print 'Evaluations: ', funcEval.cont_evalp

    # Begin the generational process
    for gen in range(1, ngen+1):

        if funcEval.cont_evalp > cont_evalf:
            break

        print '---- Generation %d -----' % (gen)
        print 'Problem: ', problem
        print 'Problem No.: ', num_p
        print 'Run No.: ', n_corr
        print 'neat-GP:', neat_alg
        print 'neat-cx:', neat_cx
        print 'Local Search:', funcEval.LS_flag
        if funcEval.LS_flag:
            print 'Local Search Heuristic: %s (%s)' % (LS_select, ls_type)

        best_ind = copy.deepcopy(best_pop(population))
        if neat_alg:
            parents = p_selection(population, gen)
        else:
            parents = toolbox.select(population, len(population))

        if neat_cx:
            n = len(parents)
            mut = 1
            cx = 1
            offspring = neatGP(toolbox, parents, cxpb, mutpb, n, mut, cx, neat_pelit)
        else:
            offspring = varOr(parents, toolbox, cxpb, mutpb)

        if neat_alg:
            specie_parents_child(parents,offspring, neat_h)
            offspring[:] = parents+offspring
            ind_specie(offspring)
            invalid_ind = [ind for ind in offspring if not ind.fitness.valid]
            fitnesses = toolbox.map(toolbox.evaluate, invalid_ind)
            for ind, fit in zip(invalid_ind, fitnesses):
                funcEval.cont_evalp += 1
                ind.fitness.values = fit
        else:
            invalid_ind = [ind for ind in offspring]
            if funcEval.LS_flag:
                new_invalid_ind = []
                for ind in invalid_ind:
                    strg = ind.__str__()
                    l_strg = add_subt(strg, ind)
                    c = tree2f()
                    cd = c.convert(l_strg)
                    new_invalid_ind.append(cd)
                fitness_ls = toolbox.map(toolbox.evaluate, new_invalid_ind)
                for ind, ls_fit in zip(invalid_ind, fitness_ls):
                    funcEval.cont_evalp += 1
                    ind.fitness.values = ls_fit
            else:
                fitnesses = toolbox.map(toolbox.evaluate, invalid_ind)
                for ind, fit in zip(invalid_ind, fitnesses):
                    funcEval.cont_evalp += 1
                    ind.fitness.values = fit

            orderbyfit = sorted(offspring, key=lambda ind:ind.fitness.values)
            print len(orderbyfit),len(best_ind)

            if best_ind.fitness.values[0] <= orderbyfit[0].fitness.values[0]:
                offspring[:] = [best_ind]+orderbyfit[:len(population)-1]

        if neat_alg:
            SpeciesPunishment(offspring, params, neat_h)

        population[:] = offspring  # population update

        cond_ind = 0
        cont_better=0
        if funcEval.LS_flag:
            for ind in population:
                ind.LS_applied_set(0)

            if   LS_select == 1:
                trees_h(population, num_p, n_corr,  pset, direccion)
            elif LS_select == 2:
                best_specie(population, num_p, n_corr, pset, direccion)
            elif LS_select == 3:
                specie_h(population, num_p, n_corr, pset, direccion)
            elif LS_select == 4:
                best_pop_ls(population, num_p, n_corr, pset, direccion)
            elif LS_select == 5:
                all_pop(population, num_p, n_corr, pset, direccion)
            elif LS_select == 6:
                trees_h_wo(population, num_p, n_corr, pset, direccion)
            elif LS_select == 7:
                ls_bestset(population, num_p, n_corr, pset, direccion)
            elif LS_select == 8:
                ls_random(population, num_p, n_corr, pset, direccion)
            elif LS_select == 9:
                ls_randbestset(population, num_p, n_corr, pset, direccion)
            #
            invalid_ind = [ind for ind in population]
            new_invalid_ind = []
            for ind in population:
                strg = ind.__str__()
                l_strg = add_subt(strg, ind)
                c = tree2f()
                cd = c.convert(l_strg)
                new_invalid_ind.append(cd)
            fitness_ls = toolbox.map(toolbox.evaluate, new_invalid_ind)
            print 'Fitness comp.:',
            for ind, ls_fit in zip(invalid_ind, fitness_ls):
                if ind.LS_applied_get() == 1:
                    cond_ind += 1
                    if ind.fitness.values[0] < ls_fit:
                        print '-',
                    elif ind.fitness.values[0] > ls_fit:
                        cont_better += 1
                        print '+',
                    elif ind.fitness.values[0] == ls_fit:
                        print '=',
                funcEval.cont_evalp += 1
                ind.fitness.values = ls_fit
            print ''

            pop_file.write('\n----------------------------------------%s'%(gen))
            for ind in population:
                pop_file.write('\n%s;%s;%s;%s'%(ind.LS_applied_get(),ind.fitness.values[0], ind, [x for x in ind.get_params()]))
        else:
            pop_file.write('\n----------------------------------------')
            for ind in population:
                pop_file.write('\n%s;%s;%s;%s;%s;%s'%(ind.LS_applied_get(),ind.LS_story_get(),ind.off_cx_get(),ind.off_mut_get(),ind.fitness.values[0], ind))

        best_ind = best_pop(population)
        if funcEval.LS_flag:
            strg = best_ind.__str__()
            l_strg = add_subt(strg, best_ind)
            c = tree2f()
            cd = c.convert(l_strg)
            new_invalid_ind.append(cd)
            fit_best = toolbox.map(toolbox.evaluate_test, [cd])
            best_ind.fitness_test.values = fit_best[0]
            best.write('\n%s;%s;%s;%s;%s;%s;%s' % (gen, funcEval.cont_evalp,  best_ind.fitness.values[0], best_ind.LS_fitness_get(), best_ind.fitness_test.values[0], len(best_ind), avg_nodes(population)))
            out.write('\n%s;%s;%s;%s;%s;%s' % (gen, len(best_ind), best_ind.LS_applied_get(), best_ind.get_params(), cd, best_ind))
        else:
            fitnesses_test = toolbox.map(toolbox.evaluate_test, [best_ind])
            best_ind.fitness_test.values = fitnesses_test[0]
            best.write('\n%s;%s;%s;%s;%s;%s' % (gen, funcEval.cont_evalp, best_ind.fitness_test.values[0], best_ind.fitness.values[0], len(best_ind), avg_nodes(population)))
            out.write('\n%s;%s;%s' % (gen, len(best_ind), best_ind))

        if funcEval.LS_flag:
            print 'Num. LS:', cond_ind
            print 'Ind. Improvement:', cont_better
            print 'Best Ind. LS:', best_ind.LS_applied_get()

        print 'Best Ind.:', best_ind
        print 'Best Fitness:', best_ind.fitness.values[0]
        print 'Test fitness:',best_ind.fitness_test.values[0]
        print 'Avg Nodes:', avg_nodes(population)
        print 'Evaluations: ', funcEval.cont_evalp

        if SaveMatrix:
            data_pop=avg_nodes(population)
            if GenMatrix:
                idx_aux = np.searchsorted(Matrix[:, 0], gen)
                Matrix[idx_aux, 1] = best_ind.fitness.values[0]
                Matrix[idx_aux, 2] = best_ind.fitness_test.values[0]
                Matrix[idx_aux, 3] = len(best_ind)
                Matrix[idx_aux, 4] = data_pop[0]
                Matrix[idx_aux, 5] = gen
                Matrix[idx_aux, 6] = 1
                Matrix[idx_aux, 7] = data_pop[1]  # max nodes
                Matrix[idx_aux, 8] = data_pop[2]  # min nodes
            else:
                if funcEval.cont_evalp >= cont_evalf:
                    num_c -= 1
                    idx_aux=num_c
                    Matrix[num_c, 1] = best_ind.fitness.values[0]
                    Matrix[num_c, 2] = best_ind.fitness_test.values[0]
                    Matrix[num_c, 3] = len(best_ind)
                    Matrix[num_c, 4] = data_pop[0]
                    Matrix[num_c, 5] = gen
                    Matrix[num_c, 6] = 1
                    Matrix[num_c, 7] = data_pop[1]  #max_nodes
                    Matrix[num_c, 8] = data_pop[2]  #min nodes
                else:
                    idx_aux = np.searchsorted(Matrix[:, 0], funcEval.cont_evalp)
                    Matrix[idx_aux, 1] = best_ind.fitness.values[0]
                    Matrix[idx_aux, 2] = best_ind.fitness_test.values[0]
                    Matrix[idx_aux, 3] = len(best_ind)
                    Matrix[idx_aux, 4] = data_pop[0]
                    Matrix[idx_aux, 5] = gen
                    Matrix[idx_aux, 6] = 1
                    Matrix[idx_aux, 7] = data_pop[1]  #max nodes
                    Matrix[idx_aux, 8] = data_pop[2]  #min nodes

                id_it = idx_aux-1
                id_beg = 0
                flag = True
                flag2 = False
                while flag:
                    if Matrix[id_it, 6] == 0:
                        id_it -= 1
                        flag2 = True
                    else:
                        id_beg = id_it
                        flag = False
                if flag2:
                    x = Matrix[id_beg, 1:8]
                    Matrix[id_beg:idx_aux, 1:] = Matrix[id_beg, 1:]

            np.savetxt('./Matrix/%s/idx_%d_%d.txt' % (problem, num_p, n_corr), Matrix, delimiter=",", fmt="%s")

    return population, logbook