Mauch MIREX 2010

function song = integratesmallparts(song)

parts = song.parts;

nParts = numel(parts);

%% break up parts (every part will have one instance)

iNewpart = 0;

partindices = [];

for iPart = 1:nParts

    parts(iPart).nInd = length(parts(iPart).indices);

    for iInstance = 1:parts(iPart).nInd

        iNewpart = iNewpart + 1;

        newparts(iNewpart) = parts(iPart);

        newparts(iNewpart).indices = parts(iPart).indices(iInstance);

        newparts(iNewpart).nInd = 1;

        partindices = [partindices; parts(iPart).indices(iInstance)];

    end

end

nNewparts = numel(newparts);

%% sort the parts in order of occurrence

[sortind, sortind] = sort(partindices);

for iNewpart = 1:nNewparts 

    sortednewparts(iNewpart) = newparts(sortind(iNewpart));

end

%%

eraselist = [];

for iNewpart = 2:nNewparts

    if sortednewparts(iNewpart).n < 12

        eraselist = [eraselist; iNewpart];

        sortednewparts(iNewpart-1).n = ...

            sortednewparts(iNewpart-1).n + sortednewparts(iNewpart).n;

    end

end

sortednewparts(eraselist) = [];

song.integratedparts = sortednewparts;