122 SUBROUTINE sgelq2( M, N, A, LDA, TAU, WORK, INFO )
130 INTEGER INFO, LDA, M, N
133 REAL A( lda, * ), TAU( * ), WORK( * )
140 parameter ( one = 1.0e+0 )
159 ELSE IF( n.LT.0 )
THEN
161 ELSE IF( lda.LT.max( 1, m ) )
THEN
165 CALL xerbla(
'SGELQ2', -info )
175 CALL slarfg( n-i+1, a( i, i ), a( i, min( i+1, n ) ), lda,
183 CALL slarf(
'Right', m-i, n-i+1, a( i, i ), lda, tau( i ),
184 $ a( i+1, i ), lda, work )
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine sgelq2(M, N, A, LDA, TAU, WORK, INFO)
SGELQ2 computes the LQ factorization of a general rectangular matrix using an unblocked algorithm...
subroutine slarfg(N, ALPHA, X, INCX, TAU)
SLARFG generates an elementary reflector (Householder matrix).
subroutine slarf(SIDE, M, N, V, INCV, TAU, C, LDC, WORK)
SLARF applies an elementary reflector to a general rectangular matrix.