function newindices = singleindices(oldindices,m)

% This function solves the problem of performing calculations such 
% as computing the spectrum of the Laplacian on SG when our chosen
% data storage scheme is redundant. We take a set of indices and 
% level number and create a set of indices which does not include
% both indices for any point which has more than one.

newindices = [] ;

for k = oldindices
  l = alternateindex(k,m) ;
  
  if(sum(l == newindices) == 0)
    newindices = [newindices k] ;
  end
end
