function indexes = h3d_filter_results(m, minimum_apart_deg)
%function indexes = h3d_filter_results(m, minimum_apart_deg)
	% m is an ordered 3xn matrix containing elements of S^2.
	
	indexes(1) = 1;
	k = 2;
	
	for i=2:size(m,2)
		v_old = m(:,indexes(k-1));
		v_new = m(:,i);
		
		if apart(v_old,v_new,minimum_apart_deg)
			indexes(k) = i;
			k = k + 1;
		end
	end
	
	
function res = apart(v1,v2,minimum_apart_deg)

	angle = acos(v1' * v2);
	
	if abs(angle) > deg2rad(minimum_apart_deg)
		res = 1;
	else
		res = 0;
	end