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

by G Golub, F T Luk, M L Overton
Venue:ACM Trans. Math. Software