def extract(self, N, F): L = np.zeros((self.tv_dim, self.tv_dim)) L[self.itril] = N.dot(self.T_iS_Tt) L += np.tril(L, -1).T + self.Im Cxx = inv(L) B = self.T_iS.dot(F) Ex = Cxx.dot(B) return Ex