getintersecttri.m

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);