Testing My Adaptive Restart Genetic Algorithm (Matlab Code)

function y = ObjectiveFunction(x)
% To customize this GA code for your problems, you only need 2 things.
% First, add objective function on the blank space on m.file named "ObjectiveFunction.m". 
% Second, define your problem on the blank space on m.file named "GA.m"
% The rest will be automatically handled by the programme.
y=(sum(sin(x) .^2, 2) - exp(-sum(x .^ 2, 2))) .* exp(-sum(sin(sqrt(abs(x))) .^2, 2));

function Y=population(p,nv,lb,ub)
% nv = number of variables
% p = population size
% lb = Lower bound
% ub= Upper bound
for i = 1:p
    for j = 1:nv
        Y(i,j)=(ub(j)-lb(j))*rand+lb(j);
    end
end
Y = Y;

function Y=crossover(P,n)
% P = population
% n = pair of chromosomes to be crossovered
[x1 y1]=size(P);
Z=zeros(2*n,y1);
for i = 1:n
    r1=randi(x1,1,2);
    while r1(1)==r1(2)
        r1=randi(x1,1,2);
    end
    A1=P(r1(1),:);
    A2=P(r1(2),:);
.
.
.

function Y=mutation(P,n)
% P = population
% n = pair of chromosomes to be mutated
[x1 y1]=size(P);
Z=zeros(2*n,y1);
for i = 1:n
    r1=randi(x1,1,2);
    while r1(1)==r1(2)
        r1=randi(x1,1,2);
    end
    A1=P(r1(1),:);
    A2=P(r1(2),:);
.
.
.

function Y=local_search(X,s,lb,ub)
% X = current best solution
% step size
[x y]=size(X);
A=ones(2*y,1)*X;
j=1;
for i=1:y
        L1=X(1,j)+s*rand;
        if L1 > ub(i)
            L1 = ub(i);
        end
.
.
.

function Y=evaluation(P,ot)
[x1 y1]=size(P);
H=zeros(x1,1);
for i = 1:x1
   H(i,1)= ObjectiveFunction(P(i,:)); 
end
if ot == 1
    Y = H;
else
    Y=10^6-H;
end

function [YY1 YY2] = selection(P1,B,p,s)
% P1 - population
% B - fitness value 
% p - population size
% s = select top s chromsomes
%-------------------------------------------------------------------------
% Top selection operation
B=B';
for i =1:s
    [r1 c1]=find(B==max(B));
    Y1(i,:)=P1(max(c1),:);
    Fn(i)=B(max(c1));
    P1(max(c1),:)=[];
    B(max(c1))=[];
    clear r1 c1
.
.
.

clear all
clc
close all
tic
%--------------------------------------------------------------------------
% To customize this GA code for your problems, you only need 2 things.
% First, add objective function on the blank space on m.file named "ObjectiveFunction.m". 
% Second, define your problem on the blank space on m.file named "GA.m"
% The rest will be automatically handled by the programme.
nv = 2;                    % number of variables
lb = [-10 -10];            % lower bound
ub = [10 10];              % upper bound
ot = -1;                   % minimization ot = -1; maximization ot = 1
t =  3;                    % computing time (s)
%--------------------------------------------------------------------------
% Maximize performance of the GA (optional)
p=50; % population size
c=10;  % crossover rate
m=10; % mutation rate
s=20;  % adaptive restart search process
g=3;  % keep top chromosomes
r=3; % number of chromosomes in initial guess
ms = 0.001; % max step size for local search
%-------------------------------------------------------------------
% Stoping criteria
tg=10000000; % number of generattion - set be be large to use computing time
%--------------------------------------------------------------------------
P1=population(r, nv, lb, ub); % Initial guess
w=1;
 
for j = 1:tg
    P=population(p-r, nv, lb, ub);
    P(p-r+1:p,:)=P1;
    for i=1:tg   
        % Extended population
        P(p+1:p+2*c,:)=crossover(P,c);
.
.
.
Sorry! This is only a half of the code.