docs: def noise(matrix):
      """Just a quick function to make a matrix noisy using a standard Poisson
      distribution (contacts are treated as rare events).
      """
  
      D = shortest_path_interpolation(matrix, strict=True)
      return np.random.poisson(lam=D)