0001 function [vertices, faces] = mergeMeshVertices(vertices, faces, vertexInds, varargin)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018 newPos = vertices(vertexInds(1), :);
0019 if nargin > 3
0020 newPos = varargin{1};
0021 end
0022
0023 vertices(vertexInds(1), :) = newPos;
0024 vertices(vertexInds(2:end), :) = NaN;
0025
0026
0027 faces(ismember(faces, vertexInds)) = vertexInds(1);
0028
0029
0030 nFaces = size(faces, 1);
0031 dgnFaces = false(nFaces, 1);
0032 dims = [1 2;1 3;2 3];
0033 for i = 1:3
0034 dgnFaces = dgnFaces | faces(:,dims(i,1)) == faces(:,dims(i,2));
0035 end
0036
0037
0038 faces(dgnFaces, :) = [];