vb_neighbor_index

 nearest neighbor index of each vertex
  xxF = vb_neighbor_index(V,F)
 --- Input
 V : vertex of surface
 F : patch index
 --- Output
 xxF{n} : neighboring vertex index of n-th vertex
        = [ 2nd-vertex of j-th patch, 3rd-vertex of j-th patch, j]
          [(# of patch connected to n-th vertex) x 3]
 Ver 1.0  by M. Sato  2004-2-10

 xxF{n} : 頂点 n に隣接する面の頂点番号と面番号

 Copyright (C) 2011, ATR All Rights Reserved.
 License : New BSD License(see VBMEG_LICENSE.txt)

0001 function    xxF = vb_neighbor_index(V,F)
0002 % nearest neighbor index of each vertex
0003 %  xxF = vb_neighbor_index(V,F)
0004 % --- Input
0005 % V : vertex of surface
0006 % F : patch index
0007 % --- Output
0008 % xxF{n} : neighboring vertex index of n-th vertex
0009 %        = [ 2nd-vertex of j-th patch, 3rd-vertex of j-th patch, j]
0010 %          [(# of patch connected to n-th vertex) x 3]
0011 % Ver 1.0  by M. Sato  2004-2-10
0012 %
0013 % xxF{n} : 頂点 n に隣接する面の頂点番号と面番号
0014 %
0015 % Copyright (C) 2011, ATR All Rights Reserved.
0016 % License : New BSD License(see VBMEG_LICENSE.txt)
0017 
0018 Npoint = size(V,1);          % 頂点数
0019 Npatch = size(F,1);          % 三角面数
0020 xxF    = cell(Npoint,1);     % 隣接頂点インデックスリスト
0021 
0022 % ３角面に関するループ
0023 for j=1:Npatch,
0024     % ３角面の頂点インデックス
0025     k1 = F(j,1); 
0026     k2 = F(j,2); 
0027     k3 = F(j,3); 
0028     
0029     % 各頂点の隣接点リストに他の頂点と面番号を加える
0030     xxF{k1} = [ xxF{k1} ; k2 , k3 ,j];
0031     xxF{k2} = [ xxF{k2} ; k3 , k1 ,j];
0032     xxF{k3} = [ xxF{k3} ; k1 , k2 ,j];
0033 end;