function [e,loop] = separateLoops(e)
% remove loops from e; make a list of the number of loops leaving
% each node.
loop = zeros(size(e,1),1);
for i=1:size(e,1),
	loops = e(i,:)==i;
	loop(i) = sum(loops)>0;
	e(i,loops) = -1;
end