function d = ABMS_close_pts(pts)
%

[num_pts,~] = size(pts);
d = zeros(num_pts,num_pts);

for p1 = 1:num_pts
    for p2 = 1:num_pts
        d(p1,p2) = norm(pts(p1,:)-pts(p2,:));
    end
end
