jaccard_indicator function

Compute indicator matrix of pairwise distances smaller than threshold