115 SUBROUTINE cung2l( M, N, K, A, LDA, TAU, WORK, INFO )
123 INTEGER INFO, K, LDA, M, N
126 COMPLEX A( lda, * ), TAU( * ), WORK( * )
133 parameter ( one = ( 1.0e+0, 0.0e+0 ),
134 $ zero = ( 0.0e+0, 0.0e+0 ) )
152 ELSE IF( n.LT.0 .OR. n.GT.m )
THEN
154 ELSE IF( k.LT.0 .OR. k.GT.n )
THEN
156 ELSE IF( lda.LT.max( 1, m ) )
THEN
160 CALL xerbla(
'CUNG2L', -info )
183 a( m-n+ii, ii ) = one
184 CALL clarf(
'Left', m-n+ii, ii-1, a( 1, ii ), 1, tau( i ), a,
186 CALL cscal( m-n+ii-1, -tau( i ), a( 1, ii ), 1 )
187 a( m-n+ii, ii ) = one - tau( i )
191 DO 30 l = m - n + ii + 1, m
subroutine cscal(N, CA, CX, INCX)
CSCAL
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine cung2l(M, N, K, A, LDA, TAU, WORK, INFO)
CUNG2L generates all or part of the unitary matrix Q from a QL factorization determined by cgeqlf (un...
subroutine clarf(SIDE, M, N, V, INCV, TAU, C, LDC, WORK)
CLARF applies an elementary reflector to a general rectangular matrix.