%% repair the solution that violates the boundary constraints

function x = Repair(x, LU)

xL = LU(1,:);

xU = LU(2,:);

for i = 1 : size(x, 1)

indexLower = find(x(i, :) < xL);

indexUper = find(x(i, :) > xU);

x(i, indexLower) = min(xU(indexLower), 2 * xL(indexLower) - x(i, indexLower));

x(i, indexUper) = max(xL(indexUper), 2 * xU(indexUper) - x(i, indexUper));

% x(i, indexLower) = xL(indexLower);

% x(i, indexUper) = xU(indexUper);

end

end