sgbbrd(3) on Linux

The text 29788ea71b8d, as Debian 2.2 (Potato), Debian 2.1 (Slink), Debian 2.0 (Hamm) shipped it.

reduce a real general m-by-n band matrix A to upper bidiagonal form B by an orthogonal transformation