SortBoundary.m

Go to the documentation of this file.

function [bet,belt]=SortBoundary(be,bel)
% function [bet,belt]=SortBoundary(be,bel)
%   Sort boundaries mesh

% start search
Label=unique(bel);
  
  %Connect=zeros(2,Th.nq,'int32');
Connect([1 2],be(1,:))=[be(2,:);bel];
bet=0*be;
belt=0*bel;
i=1;
for l=Label % On trie le bord numero l
  I=find(Connect(2,:)==l);
  while (~isempty(I))
    % [I;Connect(1,I)]
    ideb=setdiff(I,Connect(1,I));
    if isempty(ideb)
      ideb=I(1);
    end
    for ii=ideb
      j=ii;
      while (Connect(2,j) == l)
        bet(:,i)=[j;Connect(1,j)];
        belt(i)=l;
        Connect(2,j)=-l;
        j=Connect(1,j);
        i=i+1;
      end
    end
    I=find(Connect(2,:)==l);
  end
end