Find Holes in Triangular Mesh

function [holeCellArray,bounding_triangles,holeLengths] = findTriMeshHoles(faces,vertices)
% Finds holes in a triangular mesh
% Note: Does not work if a hole shares more than one vertex with other holes
% Input:
% faces = M x 3
% vertices = N x 3 (optional if you want the hole lengths)
% Output:
% holeCellArray = P x 1 cell array containing a list of holes, which are
% traced in consecutive order (list of scalar indices)
% bounding_triangles = Q x 3 list of faces that contain a bounding edge (does
% not contain triangles that only has a single bounding vertex)
% holeLengths = P x 1 vector containing the perimeter of each hole
Please let me know if you find any errors. Thanks!