ndsort.m

改进的AD算法用于重排节点编号以优化带宽！

% 根据各节点跨距对节点重新排序，若跨距相等则考虑相邻单元内节点编号总和平均值,
% 若再相同考虑含最小节点编号的编号将小，含最大节点编号的编号将大;
function refsh=ndsort(ndm,nnum)
refsh=zeros(nnum,1);
k=1;
while k<=nnum
    ndmin=min( ndm(:,1) );                        % 查找当前节点跨距最小值;
    j=0;
    for i=1:nnum
        if ndm(i,1)==ndmin
            j=j+1;
            min_ndist(j)=i;                       % 记录当前最小节点跨距值的出现位置min_ndist;
        end
    end
    for i=1:j
        ndm( min_ndist(i),1 )=inf;                % 将min_ndist处的节点跨距值置大数，下次作比较时将直接跳过;
    end
    if j>=2                                       % 若跨距相等，则比较其相邻单元节点编号总和平均值;
        jj=1;
        while jj<=j
            nsmin=min( ndm(min_ndist,2) );        % 寻找拥有当前最小节点跨距值的节点中其相邻单元节点编号总和平均值的最小值;
            l=0;
            for i=1:j
                if ndm( min_ndist(i),2 )==nsmin
                    l=l+1;
                    min_nsum(l)=min_ndist(i);     % 记录当前最小跨距相等的节点相邻单元节点编号总和平均值的最小值的出现位置min_nsum;
                end
            end
            for i=1:l
                ndm( min_nsum(i),2 )=inf;         % 将min_nsum处的节点相邻单元节点编号总和平均值置大数，下次作比较时将直接跳过;
            end
            if l>=2                               % 若相邻单元节点编号总和平均值相等，则比较相邻单元节点编号最小值;
                jjj=1;
                while jjj<=l
                    nmmin=min( ndm(min_nsum,3) );     % 寻找拥有当前相邻单元最小节点编号值;
                    m=0;
                    for i=1:l
                        if ndm( min_nsum(i),3 )==nmmin;
                            m=m+1;
                            min_nmin(m)=min_nsum(i);  % 记录当前相邻单元节点编号最小值的最小值的出现位置min_nmin;
                        end
                    end
                    for i=1:m
                        ndm( min_nmin(i),3 )=inf;     % 将min_nmin处的节点相邻单元节点编号最小值置大数，下次作比较时将直接跳过;
                    end
                    if m>=2
                        for i=1:m
                            refsh(k)=min_nmin(i);
                            k=k+1;
                            jj=jj+1;
                            jjj=jjj+1;
                        end
                    else
                        refsh(k)=min_nmin(m);
                        k=k+1;
                        jj=jj+1;
                        jjj=jjj+1;
                    end
                    clear min_nmin;
                end
            else                                  % 否则将此节点位置进行排序到refsh中;
                refsh(k)=min_nsum(l);
                k=k+1;
                jj=jj+1;
            end
            clear min_nsum;
        end
    else
        refsh(k)=min_ndist(j);
        k=k+1;
    end
    clear min_ndist;
end