updatebd.m

来自「C++编译指导」· M 代码 · 共 26 行

function bdEdge = updatebd(bdEdge,marker2,d2p)
% UPDATEDBD refine the boundary edges
%
% USAGE
%    bdEdge = updatebd(bdEdge,marker,d2p)
%
% INPUT
%   bdEdge:  set of boundary edges
%   marker:  new node index for marked edge
%      d2p:  index mapping from dual edge to primary edge
% 
% OUTPUT
%   bdEdge:  set of refined boundary edges
%

NB = size(bdEdge,1);
for k = 1:NB
     i = bdEdge(k,1); j = bdEdge(k,2);
     if (marker2(d2p(i,j))>0)
         bdEdge(k,:) = [i,marker2(d2p(i,j))];
         bdEdge(size(bdEdge,1)+1,:) = [marker2(d2p(i,j)),j];
     end
end
%--------------------------------------------------------------------------
% End of function UPDATEDBD
%--------------------------------------------------------------------------