The formal numerical scaling for solving Eq. 98c is O(N 4 ), due to the N 2 dimensionality of both H j and Ω k . It is thus ideal to find alternative but equivalent EOMs that reduce this scaling. Below, we derive three sets of linearized EOMs that are equivalent to the EOMs in Eq. 98c.