2 Method Notations We use X>, X1, Tr(X) and vec(X) to denote the transpose, inverse, trace, and column-wise vectorization of amatrixX. We use X Y to represent the Kronecker product