function mesh = mmgReadMesh(filename)

fid = fopen(filename, 'r');

if ~fid

    error('Could not open %s for reading.', filename);

end

str = fgets(fid);

while str ~= -1

    if startsWith(str,'Vertices')

        N = str2double(fgets(fid));

        points = fscanf(fid,'%f %f %f %d',[4 N])';

    elseif startsWith(str,'Tetrahedra')

        N = str2double(fgets(fid));

        cells = fscanf(fid,'%d %d %d %d %d',[5 N])';

    end

    str = fgets(fid);

end

fclose(fid);

mesh.points = points(:,1:3);

mesh.cells = int32(cells(:,1:4));

mesh.cellData.class = uint8(cells(:,5));

mesh.cellTypes = repmat(uint8(10), size(cells,1), 1);

end