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