def l1_norm_off_diag(A): "convenience method for l1 norm, excluding diagonal" # let's speed this up later # assume A symmetric, sparse too return np.linalg.norm(A - np.diag(A.diagonal()), ord=1)