version 1.13, 2012/12/14 14:22:57
|
version 1.19, 2018/05/29 07:18:42
|
Line 2
|
Line 2
|
* |
* |
* =========== DOCUMENTATION =========== |
* =========== DOCUMENTATION =========== |
* |
* |
* Online html documentation available at |
* Online html documentation available at |
* http://www.netlib.org/lapack/explore-html/ |
* http://www.netlib.org/lapack/explore-html/ |
* |
* |
*> \htmlonly |
*> \htmlonly |
*> Download ZUNMLQ + dependencies |
*> Download ZUNMLQ + dependencies |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zunmlq.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zunmlq.f"> |
*> [TGZ]</a> |
*> [TGZ]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zunmlq.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zunmlq.f"> |
*> [ZIP]</a> |
*> [ZIP]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zunmlq.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zunmlq.f"> |
*> [TXT]</a> |
*> [TXT]</a> |
*> \endhtmlonly |
*> \endhtmlonly |
* |
* |
* Definition: |
* Definition: |
* =========== |
* =========== |
* |
* |
* SUBROUTINE ZUNMLQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, |
* SUBROUTINE ZUNMLQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, |
* WORK, LWORK, INFO ) |
* WORK, LWORK, INFO ) |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
* CHARACTER SIDE, TRANS |
* CHARACTER SIDE, TRANS |
* INTEGER INFO, K, LDA, LDC, LWORK, M, N |
* INTEGER INFO, K, LDA, LDC, LWORK, M, N |
Line 28
|
Line 28
|
* .. Array Arguments .. |
* .. Array Arguments .. |
* COMPLEX*16 A( LDA, * ), C( LDC, * ), TAU( * ), WORK( * ) |
* COMPLEX*16 A( LDA, * ), C( LDC, * ), TAU( * ), WORK( * ) |
* .. |
* .. |
* |
* |
* |
* |
*> \par Purpose: |
*> \par Purpose: |
* ============= |
* ============= |
Line 136
|
Line 136
|
*> The dimension of the array WORK. |
*> The dimension of the array WORK. |
*> If SIDE = 'L', LWORK >= max(1,N); |
*> If SIDE = 'L', LWORK >= max(1,N); |
*> if SIDE = 'R', LWORK >= max(1,M). |
*> if SIDE = 'R', LWORK >= max(1,M). |
*> For optimum performance LWORK >= N*NB if SIDE 'L', and |
*> For good performance, LWORK should generally be larger. |
*> LWORK >= M*NB if SIDE = 'R', where NB is the optimal |
|
*> blocksize. |
|
*> |
*> |
*> If LWORK = -1, then a workspace query is assumed; the routine |
*> If LWORK = -1, then a workspace query is assumed; the routine |
*> only calculates the optimal size of the WORK array, returns |
*> only calculates the optimal size of the WORK array, returns |
Line 156
|
Line 154
|
* Authors: |
* Authors: |
* ======== |
* ======== |
* |
* |
*> \author Univ. of Tennessee |
*> \author Univ. of Tennessee |
*> \author Univ. of California Berkeley |
*> \author Univ. of California Berkeley |
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date November 2011 |
*> \date December 2016 |
* |
* |
*> \ingroup complex16OTHERcomputational |
*> \ingroup complex16OTHERcomputational |
* |
* |
Line 169
|
Line 167
|
SUBROUTINE ZUNMLQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, |
SUBROUTINE ZUNMLQ( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, |
$ WORK, LWORK, INFO ) |
$ WORK, LWORK, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.4.0) -- |
* -- LAPACK computational routine (version 3.7.0) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* November 2011 |
* December 2016 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER SIDE, TRANS |
CHARACTER SIDE, TRANS |
Line 185
|
Line 183
|
* ===================================================================== |
* ===================================================================== |
* |
* |
* .. Parameters .. |
* .. Parameters .. |
INTEGER NBMAX, LDT |
INTEGER NBMAX, LDT, TSIZE |
PARAMETER ( NBMAX = 64, LDT = NBMAX+1 ) |
PARAMETER ( NBMAX = 64, LDT = NBMAX+1, |
|
$ TSIZE = LDT*NBMAX ) |
* .. |
* .. |
* .. Local Scalars .. |
* .. Local Scalars .. |
LOGICAL LEFT, LQUERY, NOTRAN |
LOGICAL LEFT, LQUERY, NOTRAN |
CHARACTER TRANST |
CHARACTER TRANST |
INTEGER I, I1, I2, I3, IB, IC, IINFO, IWS, JC, LDWORK, |
INTEGER I, I1, I2, I3, IB, IC, IINFO, IWT, JC, LDWORK, |
$ LWKOPT, MI, NB, NBMIN, NI, NQ, NW |
$ LWKOPT, MI, NB, NBMIN, NI, NQ, NW |
* .. |
* .. |
* .. Local Arrays .. |
|
COMPLEX*16 T( LDT, NBMAX ) |
|
* .. |
|
* .. External Functions .. |
* .. External Functions .. |
LOGICAL LSAME |
LOGICAL LSAME |
INTEGER ILAENV |
INTEGER ILAENV |
Line 246
|
Line 242
|
* |
* |
IF( INFO.EQ.0 ) THEN |
IF( INFO.EQ.0 ) THEN |
* |
* |
* Determine the block size. NB may be at most NBMAX, where NBMAX |
* Compute the workspace requirements |
* is used to define the local array T. |
|
* |
* |
NB = MIN( NBMAX, ILAENV( 1, 'ZUNMLQ', SIDE // TRANS, M, N, K, |
NB = MIN( NBMAX, ILAENV( 1, 'ZUNMLQ', SIDE // TRANS, M, N, K, |
$ -1 ) ) |
$ -1 ) ) |
LWKOPT = MAX( 1, NW )*NB |
LWKOPT = MAX( 1, NW )*NB + TSIZE |
WORK( 1 ) = LWKOPT |
WORK( 1 ) = LWKOPT |
END IF |
END IF |
* |
* |
Line 272
|
Line 267
|
NBMIN = 2 |
NBMIN = 2 |
LDWORK = NW |
LDWORK = NW |
IF( NB.GT.1 .AND. NB.LT.K ) THEN |
IF( NB.GT.1 .AND. NB.LT.K ) THEN |
IWS = NW*NB |
IF( LWORK.LT.NW*NB+TSIZE ) THEN |
IF( LWORK.LT.IWS ) THEN |
NB = (LWORK-TSIZE) / LDWORK |
NB = LWORK / LDWORK |
|
NBMIN = MAX( 2, ILAENV( 2, 'ZUNMLQ', SIDE // TRANS, M, N, K, |
NBMIN = MAX( 2, ILAENV( 2, 'ZUNMLQ', SIDE // TRANS, M, N, K, |
$ -1 ) ) |
$ -1 ) ) |
END IF |
END IF |
ELSE |
|
IWS = NW |
|
END IF |
END IF |
* |
* |
IF( NB.LT.NBMIN .OR. NB.GE.K ) THEN |
IF( NB.LT.NBMIN .OR. NB.GE.K ) THEN |
Line 292
|
Line 284
|
* |
* |
* Use blocked code |
* Use blocked code |
* |
* |
|
IWT = 1 + NW*NB |
IF( ( LEFT .AND. NOTRAN ) .OR. |
IF( ( LEFT .AND. NOTRAN ) .OR. |
$ ( .NOT.LEFT .AND. .NOT.NOTRAN ) ) THEN |
$ ( .NOT.LEFT .AND. .NOT.NOTRAN ) ) THEN |
I1 = 1 |
I1 = 1 |
Line 324
|
Line 317
|
* H = H(i) H(i+1) . . . H(i+ib-1) |
* H = H(i) H(i+1) . . . H(i+ib-1) |
* |
* |
CALL ZLARFT( 'Forward', 'Rowwise', NQ-I+1, IB, A( I, I ), |
CALL ZLARFT( 'Forward', 'Rowwise', NQ-I+1, IB, A( I, I ), |
$ LDA, TAU( I ), T, LDT ) |
$ LDA, TAU( I ), WORK( IWT ), LDT ) |
IF( LEFT ) THEN |
IF( LEFT ) THEN |
* |
* |
* H or H**H is applied to C(i:m,1:n) |
* H or H**H is applied to C(i:m,1:n) |
Line 342
|
Line 335
|
* Apply H or H**H |
* Apply H or H**H |
* |
* |
CALL ZLARFB( SIDE, TRANST, 'Forward', 'Rowwise', MI, NI, IB, |
CALL ZLARFB( SIDE, TRANST, 'Forward', 'Rowwise', MI, NI, IB, |
$ A( I, I ), LDA, T, LDT, C( IC, JC ), LDC, WORK, |
$ A( I, I ), LDA, WORK( IWT ), LDT, |
$ LDWORK ) |
$ C( IC, JC ), LDC, WORK, LDWORK ) |
10 CONTINUE |
10 CONTINUE |
END IF |
END IF |
WORK( 1 ) = LWKOPT |
WORK( 1 ) = LWKOPT |