, Gregorio Quintana-Orti.
    </p>

    <p>
      The distinguishing features of these routines include the fact
      that
      <ul>
        <li>
          they accept the banded matrix in the compact "GB"
          ("General Banded") format
        </li>
        <li>
          the QR factors are also stored compactly, by overwriting
          the input matrix storage.
        </li>
        <li>
          large problems are handled efficiently by taking advantage
          of the blocked version of the al