function ALL = allEquid2d(ARRAY);

% allequid(<array>) -- Find nearest neighbors for all points in an array.
%
% The allnn2d() function takes the array <array>, and finds all the individuals
% in the array's nearest neighbors.  If more than one point lies exactly the same
% distance from the argument, however, allnn2d() will still only return one point.  
% In other words, with allnn2d it is possible to have at MOST one nearest neighbor.
%
% With allEquid, ALL points that are equidistant from a given point are returned.
% This is an N^2 algorithm.    
 
DIST = [];  						% initialize the DISTance matrix

N = size(ARRAY,1);			  

for n = 1:N							% loop over every point and find distance
   for m = 1:N						% for each other point in the group
      if (m == n)
         DIST(n,m) = inf;
      else
         x1 = ARRAY(n,2);
         y1 = ARRAY(n,3);
         
         x2 = ARRAY(m,2);
         y2 = ARRAY(m,3);
         
         DIST(n,m) = dist2d(x1,y1,x2,y2);
      end		
   end
end

currentlow = inf;
longlist = [];
zerowidth = 0;
q = 0;

% this is all just to find out how wide to make our zeros array

for i = 1:size(DIST,1)
   for j = 1:size(DIST,2)  
      if (DIST(i,j) < currentlow)   % if the jth element of the ith row is less than
         currentlow = DIST(i,j);		% our current lowest element, change currentlow
      end
   end
   for j = 1:size(DIST,2)
      if (DIST(i,j) == currentlow)  % if the cell holds the currentlow value
         q = q+1;							% increment our counter		
      end
      if (q > zerowidth)
         zerowidth = q;					% update the width of our final array
      end
      q = 0;
   end
   currentlow = inf;
end

final = zeros(size(ARRAY,1), zerowidth); % build our zeros matrix

for i = 1:size(DIST,1)
   for j = 1:size(DIST,2)
      
      if (DIST(i,j) < currentlow)
         currentlow = DIST(i,j);
      end
   end
   available = 1;
   for j = 1:size(DIST,2)
      if (DIST(i,j) == currentlow)	% put the neighbors          
         final (i,available) = j;	% into their proper place in the final array.
         available = available + 1;
      end
   end
   
   currentlow = inf;
end


ALL = final;