Question

Matlab-How can i exit the for loop when counts exceeds the max value, then continue to the last loop and plot accepted values.

counts = 0;
for k = 2:length(partDiameter)
overlaps = 1;
while overlaps == 1
counts = counts + 1;
if counts >= 10000
disp('Max particles packed')
inCoord = [inCoord;[x,y]];
return;
end
x = (partDiameter(k)-L)*rand(1,1) + (L-partDiameter(k)/2);
y = (partDiameter(k)-L)*rand(1,1) + (L-partDiameter(k)/2);

% check if this x,y postion overlaps with previously packed
% partices
for m = 1:length(inCoord(:,1))
% distance between particles
realDist = sqrt((x-inCoord(m,1))^2 + (y-inCoord(m,2))^2);
overDist = partDiameter(m)/2 + partDiameter(k)/2;
if counts >= 10000
disp('Max particles packed')
inCoord = [inCoord;[x,y]];
continue;
end
if (realDist < overDist)
overlaps = 1;
break;
elseif (realDist >= overDist && m == length(inCoord(:,1)))
overlaps = 0;
inCoord = [inCoord;[x,y]];
solidFrac = sum(pi*(partDiameter(m)/2).^2)/L;
  
end
end
end
  
end
% Plot particles
figure(2)
hold on
for k = 1:length(partDiameter)
rectangle('Position',...
[inCoord(k,1)-partDiameter(k)/2, inCoord(k,2)-partDiameter(k)/2,...
partDiameter(k),partDiameter(k)],...
'Curvature', [1,1],...
'FaceColor','b',...
'EdgeColor','k')
end
xlim([0,L])
ylim([0,L])
axis square

Answer 1

`Hey,

Note: Brother in case of any queries, just comment in box I would be very happy to assist all your queries

Brother use break statement to get out of the loop immediately. It will exit the first for loop and will execute other remaining statements.

counts = 0;
for k = 2:length(partDiameter)
overlaps = 1;
while overlaps == 1
counts = counts + 1;
if counts >= 10000
disp('Max particles packed')
inCoord = [inCoord;[x,y]];
break;
end
x = (partDiameter(k)-L)*rand(1,1) + (L-partDiameter(k)/2);
y = (partDiameter(k)-L)*rand(1,1) + (L-partDiameter(k)/2);

% check if this x,y postion overlaps with previously packed
% partices
for m = 1:length(inCoord(:,1))
% distance between particles
realDist = sqrt((x-inCoord(m,1))^2 + (y-inCoord(m,2))^2);
overDist = partDiameter(m)/2 + partDiameter(k)/2;
if counts >= 10000
disp('Max particles packed')
inCoord = [inCoord;[x,y]];
break;
end
if (realDist < overDist)
overlaps = 1;
break;
elseif (realDist >= overDist && m == length(inCoord(:,1)))
overlaps = 0;
inCoord = [inCoord;[x,y]];
solidFrac = sum(pi*(partDiameter(m)/2).^2)/L;
  
end
end
end
  
end
% Plot particles
figure(2)
hold on
for k = 1:length(partDiameter)
rectangle('Position',...
[inCoord(k,1)-partDiameter(k)/2, inCoord(k,2)-partDiameter(k)/2,...
partDiameter(k),partDiameter(k)],...
'Curvature', [1,1],...
'FaceColor','b',...
'EdgeColor','k')
end
xlim([0,L])
ylim([0,L])
axis square

Kindly revert for any queries

Thanks.