Codebase list octave-iso2mesh / lintian-fixes/main getintersecttri.m
lintian-fixes/main

Tree @lintian-fixes/main (Download .tar.gz)

getintersecttri.m @lintian-fixes/main — 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
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);