Skip to content
Prev 9459 / 63421 Next

prcomp need Conj(t(s$vt)) (PR#2924)

In prcomp

        s <- La.svd(x, nu = 0)
        s$v <- t(s$vt)

the second above line should be 

        s$v <- Conj(t(s$vt))

(to cover complex cases).

Paul Gilbert