tors used to
          compute the reduction.
          On exit with LTRA = .TRUE., the leading N-by-(K+N) part of
          this array contains the matrix Aout and in the zero parts
          information about the elementary reflectors.

  LDA     INTEGER
          The leading dimension of the array A.
          LDA &gt;= MAX(1,K+N),  if LTRA = .FALSE.;
          LDA &gt;= MAX(1,N),    if LTRA = .TRUE..

  B       (input/output) DOUBLE PRECISION array, dimension
                  (LDB,K+N)   if LTRB = .FAL