TSVD.jl icon indicating copy to clipboard operation
TSVD.jl copied to clipboard

a problem with TSVD.teig - anything I did wrong?

Open jamesjcai opened this issue 4 years ago • 2 comments

using LinearAlgebra, TSVD X=rand(4,5) X=X'X F=TSVD.teig(X) X*F[2] ≈ F[2]*diagm(F[1]) false

jamesjcai avatar Mar 17 '20 16:03 jamesjcai

Just repeating my reply from Slack here for the record. So teig isn't documented and exported because it's not complete although it's pretty close. That means that F[2] doesn't contain the eigenvectors of X but the eigenvectors of the tridiagonal reduction of X. F[4] contains the Lanczos vectors so it should be possible to construct the eigenvectors of X. Let's keep this issue open until I've fixed this.

andreasnoack avatar Mar 17 '20 20:03 andreasnoack

If I want to use TSVD.jl to compute truncated svd, what function should I call?

Vilin97 avatar Jun 27 '20 16:06 Vilin97