lp = np.flatnonzero(np.tril(ones(BIG.shape)))