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
31
32
33
34
|
function eid = getintersecttri(tmppath)
%
% eid=getintersecttri(tmppath)
%
% get the IDs of self-intersecting elements from tetgen
% call this when tetgen complains about self-intersection
%
% author: Qianqian Fang, <q.fang at neu.edu>
%
% input:
% tmppath: working dir, use mwpath('') in most cases
%
% output:
% eid: an array of all intersecting surface elements,
% one can read the corresponding node/elem by
% [no,el]=readoff(mwpath('post_vmesh.off'));
%
% -- this function is part of iso2mesh toolbox (http://iso2mesh.sf.net)
%
exesuff = getexeext;
exesuff = fallbackexeext(exesuff, 'tetgen');
[status, str] = system(['"' mcpath('tetgen') exesuff '" -d "' ...
tmppath 'post_vmesh.poly"']);
eid = [];
if (status == 0)
id = regexp(str, ' #([0-9]+) ', 'tokens');
for j = 1:length(id)
eid(end + 1) = str2num(id{j}{1});
end
end
eid = unique(eid);
|
