dlabrd(3) on Linux
The text 9831b230f0f0, as Debian 3.0 (Woody) shipped it.
reduce the first NB rows and columns of a real general m by n matrix A to upper or lower bidiagonal form by an orthogonal transformation Q' * A * P, and returns the matrices X and Y which are needed to apply the transformation to the unreduced part of A