Codebase list octave-iso2mesh / debian/1.9.6+ds-3 removedupnodes.m
debian/1.9.6+ds-3

Tree @debian/1.9.6+ds-3 (Download .tar.gz)

removedupnodes.m @debian/1.9.6+ds-3 — raw · history · blame 

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
function [newnode,newelem]=removedupnodes(node,elem,tol)
%
% [newnode,newelem]=removedupnodes(node,elem)
%
% removing the duplicated nodes from a mesh
%
% author: Qianqian Fang <q.fang at neu.edu>
%
% input:
%   elem: integer array with dimensions of NE x 4, each row contains
%         the indices of all the nodes for each tetrahedron
%   node: node coordinates, 3 columns for x, y and z respectively
%
% output:
%   newnode: nodes without duplicates
%   newelem: elements with only the unique nodes
%
%
% -- this function is part of iso2mesh toolbox (http://iso2mesh.sf.net)
%

if(nargin>=3 && tol~=0)
    node=round(node/tol)*tol;
end
[newnode,I,J]=unique(node,'rows');
if(iscell(elem))
    newelem=cellfun(@(x) J(x)', elem,'UniformOutput',false);
else
    newelem=J(elem);
end