I need to use global variable in nested functions.
On running the below code, I'm getting an error
The GLOBAL or PERSISTENT declaration of "initial_flag" appears in a nested function, but should be in the outermost function where it is used.
Here is the logic of the below code:
There is a global variable named initial_flag which is defined in "demoEA.m" and "benchmark_func.m" and also in the nested function within "benchmark_func.m". The purpose of this variable is to avoid running code
if (func_num == 1) fhd = @f1;
elseif (func_num == 2) fhd = @f2;
end
multiple times in "benchmark_func.m" and avoid loading "datafiles/f01.mat" multiple times which contains the definition of (xopt, lb, ub) used in the nested functions.
Here is the Code.
"demoEA.m"
clear; clc; close all;
lb = -100;
ub = 100;
D = 1000;
NP = 100;
runs = 2;
Max_FEs = 1000;
global initial_flag;
format short e;
for func_num = 1:2
for run = 1:runs
initial_flag = 0;
pop = lb + (ub-lb)*rand(D, NP);
val = benchmark_func(pop, func_num);
FEs = NP;
while (FEs <= Max_FEs-NP)
newpop = lb + (ub-lb)*rand(D, NP);
newval = benchmark_func(newpop, func_num);
FEs = FEs + NP;
index = (newval < val);
val(index) = newval(index);
pop(:, index) = newpop(:, index);
fprintf(1, 'func_num = %d, run = %d, FEs = %d\n', func_num, run, FEs);
fprintf(1, 'min(val) = %e\n\n', min(val));
end
end
end
"benchmark_func.m"
function fit = benchmark_func(x, func_num)
global initial_flag
persistent fhd
if (initial_flag == 0)
if (func_num == 1) fhd = @f1;
elseif (func_num == 2) fhd = @f2;
end
end
fit = fhd(x);
function fit = elliptic(x)
[D ps] = size(x);
condition = 1e+6;
coefficients = condition .^ linspace(0, 1, D);
fit = coefficients * T_irreg(x).^2;
end
function fit = rosenbrock(x)
[D ps] = size(x);
fit = sum(100.*(x(1:D-1,:).^2-x(2:D, :)).^2+(x(1:D-1, :)-1).^2);
end
function fit = f1(x)
global initial_flag
persistent xopt lb ub
[D ps] = size(x);
if (initial_flag == 0)
load 'datafiles/f01.mat';
initial_flag = 1;
end
idx = checkBounds(x, lb, ub);
x = x - repmat(xopt, 1, ps);
fit = elliptic(x);
fit(idx) = NaN;
end
function fit = f2(x)
global initial_flag
persistent xopt lb ub
[D ps] = size(x);
if (initial_flag == 0)
load 'datafiles/f02.mat';
initial_flag = 1;
end
idx = checkBounds(x, lb, ub);
x = x - repmat(xopt, 1, ps);
fit = rosenbrock(x);
fit(idx) = NaN;
end
function g = T_irreg(f)
a = 0.1;
g = f;
idx = (f > 0);
g(idx) = log(f(idx))/a;
g(idx) = exp(g(idx) + 0.49*(sin(g(idx)) + sin(0.79*g(idx)))).^a;
idx = (f < 0);
g(idx) = log(-f(idx))/a;
g(idx) = -exp(g(idx) + 0.49*(sin(0.55*g(idx)) + sin(0.31*g(idx)))).^a;
end
function indices = checkBounds(x, lb, ub)
indices = find(sum(x > ub | x < lb) > 0);
end
end
Please help me to solve this issue. For complete code please refer attachments.
