PHT = np.dot(P, H.T)