Dear Sir/Madam,
I am a Ph.D. candidate in TU, Munich. I have a question about rankk
updating. Making use of the subroutine 'dsyrk', one can update the symmetric
matrix with only half of the matrix is updated. But with this subroutine,
the other half part, which is not referenced, is still in the memory, of the
case which I want to avoid.
My question is, if I want to update a matrix like C=alfa*A'*A + C, and C is
a mbym matrix in packed symmetric storage, A is a n by m matrix, how can I
implement this operation with BLAS or Lapack library, with high performance.
I know that one can use 'dspr' to realize this updating, but since the Blas
level 1 operation is less efficient than level 3, I hope there are some ways
to improve the computation speed.
Best regards
Weiyong

Weiyong YI
D80290 Muenchen
Tel: (+49)8928923185
email: chxyike@Domain.Removed
 next part 
An HTML attachment was scrubbed...
URL:
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100304/124aa839/attachment.html
