A block Lanczos method for computing the singular values and corresponding singular vectors of a matrix (1981)

by G H Golub, F T Luk, M L Overton
Venue:ACM Transactions on Mathematical Software