updatebd.m

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

function bdedge = updatebd(bdedge,marker,d2p)
% UPDATEDBD refine the boundary edges
%
% bdedge = updatebd(marker,bdedge,d2p)
%
% Input:
%   marker: See function mark;
%   bdedge: boundary edge set
%   d2p: index mapping from dual edge to primary edge
% 
% Output:
%   bdedge: refined boundary edges set
%
% L. Chen 03-26-2006

ND = size(bdedge,1);
for k = 1:ND
    i = bdedge(k,1); j = bdedge(k,2);
    if marker(d2p(i,j)) >0
       bdedge(k,:) = [i,marker(d2p(i,j))];
       bdedge(size(bdedge,1)+1,:) = [j,marker(d2p(i,j))];
    end
end