docs: def kscale(matrix, k=7, dists=None):
      """ Returns the local scale based on the k-th nearest neighbour """
      dists = (kdists(matrix, k=k) if dists is None else dists)
      scale = dists.dot(dists.T)
      return scale