powered by
Fix sign ambiguity of eigen-vectors by making U positive diagonal
make.pos.diag(U)
target matrix