if L is None:
        L = np.eye(n)