214 SUBROUTINE zggrqf( M, P, N, A, LDA, TAUA, B, LDB, TAUB, WORK,
223 INTEGER INFO, LDA, LDB, LWORK, M, N, P
226 COMPLEX*16 A( lda, * ), B( ldb, * ), TAUA( * ), TAUB( * ),
234 INTEGER LOPT, LWKOPT, NB, NB1, NB2, NB3
244 INTRINSIC int, max, min
251 nb1 = ilaenv( 1,
'ZGERQF',
' ', m, n, -1, -1 )
252 nb2 = ilaenv( 1,
'ZGEQRF',
' ', p, n, -1, -1 )
253 nb3 = ilaenv( 1,
'ZUNMRQ',
' ', m, n, p, -1 )
254 nb = max( nb1, nb2, nb3 )
255 lwkopt = max( n, m, p )*nb
257 lquery = ( lwork.EQ.-1 )
260 ELSE IF( p.LT.0 )
THEN
262 ELSE IF( n.LT.0 )
THEN
264 ELSE IF( lda.LT.max( 1, m ) )
THEN
266 ELSE IF( ldb.LT.max( 1, p ) )
THEN
268 ELSE IF( lwork.LT.max( 1, m, p, n ) .AND. .NOT.lquery )
THEN
272 CALL xerbla(
'ZGGRQF', -info )
274 ELSE IF( lquery )
THEN
280 CALL zgerqf( m, n, a, lda, taua, work, lwork, info )
285 CALL zunmrq(
'Right',
'Conjugate Transpose', p, n, min( m, n ),
286 $ a( max( 1, m-n+1 ), 1 ), lda, taua, b, ldb, work,
288 lopt = max( lopt, int( work( 1 ) ) )
292 CALL zgeqrf( p, n, b, ldb, taub, work, lwork, info )
293 work( 1 ) = max( lopt, int( work( 1 ) ) )
subroutine zggrqf(M, P, N, A, LDA, TAUA, B, LDB, TAUB, WORK, LWORK, INFO)
ZGGRQF
subroutine zunmrq(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
ZUNMRQ
subroutine zgerqf(M, N, A, LDA, TAU, WORK, LWORK, INFO)
ZGERQF
subroutine zgeqrf(M, N, A, LDA, TAU, WORK, LWORK, INFO)
ZGEQRF VARIANT: left-looking Level 3 BLAS of the algorithm.
subroutine xerbla(SRNAME, INFO)
XERBLA