innersurf

 outface=innersurf(node,face,outface)

 extract the interior triangles (shared by two enclosed compartments) of a complex surface

 author: Qianqian Fang, <q.fang at neu.edu>

 input:
    node:  node coordinates
    face:  surface triangle list
    outface: (optional) the exterior triangle list, if not given, will
           be computed using outersurf().

 output:
    inface: the collection of interior triangles of the surface mesh

 -- this function is part of iso2mesh toolbox (http://iso2mesh.sf.net)

0001 function [inface, outface]=innersurf(node,face,outface)
0002 %
0003 % outface=innersurf(node,face,outface)
0004 %
0005 % extract the interior triangles (shared by two enclosed compartments) of a complex surface
0006 %
0007 % author: Qianqian Fang, <q.fang at neu.edu>
0008 %
0009 % input:
0010 %    node:  node coordinates
0011 %    face:  surface triangle list
0012 %    outface: (optional) the exterior triangle list, if not given, will
0013 %           be computed using outersurf().
0014 %
0015 % output:
0016 %    inface: the collection of interior triangles of the surface mesh
0017 %
0018 % -- this function is part of iso2mesh toolbox (http://iso2mesh.sf.net)
0019 %
0020 
0021 if(nargin<3)
0022   outface=outersurf(node,face);
0023 end
0024 
0025 [I,J]=ismember(sort(face,2),sort(outface,2),'rows');
0026 
0027 inface=face(I==0,:);