LAPACK Archives

[Lapack] Moment matrix eigenvector update

Hi,

if X is an N x M (N>>M) matrix let C = X' * X and C = VDV'

if w is a Mx1 vector what is the diagonalization of ( C + ww' )?

I understand that the trick is to diagonalize (D+pp') where p=Vw
and that this is the heart of the divide and conquer algorithm
but are there any LAPACK function to do it?

I am also interested to do an SVD on
[ X ]
[ -- ]
[ w']

which is an analogous problem and interesting if C is rank defficient.

Cheers,

Laszlo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: 
http://lists.cs.utk.edu/private/lapack/attachments/20061115/55249aa3/attachment.htm
 

<Prev in Thread] Current Thread [Next in Thread>


For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or