version 1.11, 2012/07/31 11:06:36
|
version 1.18, 2017/06/17 11:06:29
|
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 DORMQL + dependencies |
*> Download DORMQL + dependencies |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dormql.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dormql.f"> |
*> [TGZ]</a> |
*> [TGZ]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dormql.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dormql.f"> |
*> [ZIP]</a> |
*> [ZIP]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dormql.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dormql.f"> |
*> [TXT]</a> |
*> [TXT]</a> |
*> \endhtmlonly |
*> \endhtmlonly |
* |
* |
* Definition: |
* Definition: |
* =========== |
* =========== |
* |
* |
* SUBROUTINE DORMQL( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, |
* SUBROUTINE DORMQL( 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 .. |
* DOUBLE PRECISION A( LDA, * ), C( LDC, * ), TAU( * ), WORK( * ) |
* DOUBLE PRECISION 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 doubleOTHERcomputational |
*> \ingroup doubleOTHERcomputational |
* |
* |
Line 169
|
Line 167
|
SUBROUTINE DORMQL( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, |
SUBROUTINE DORMQL( 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 |
INTEGER I, I1, I2, I3, IB, IINFO, IWS, LDWORK, LWKOPT, |
INTEGER I, I1, I2, I3, IB, IINFO, IWT, LDWORK, LWKOPT, |
$ MI, NB, NBMIN, NI, NQ, NW |
$ MI, NB, NBMIN, NI, NQ, NW |
* .. |
* .. |
* .. Local Arrays .. |
|
DOUBLE PRECISION T( LDT, NBMAX ) |
|
* .. |
|
* .. External Functions .. |
* .. External Functions .. |
LOGICAL LSAME |
LOGICAL LSAME |
INTEGER ILAENV |
INTEGER ILAENV |
Line 239
|
Line 235
|
INFO = -7 |
INFO = -7 |
ELSE IF( LDC.LT.MAX( 1, M ) ) THEN |
ELSE IF( LDC.LT.MAX( 1, M ) ) THEN |
INFO = -10 |
INFO = -10 |
|
ELSE IF( LWORK.LT.NW .AND. .NOT.LQUERY ) THEN |
|
INFO = -12 |
END IF |
END IF |
* |
* |
IF( INFO.EQ.0 ) THEN |
IF( INFO.EQ.0 ) THEN |
|
* |
|
* Compute the workspace requirements |
|
* |
IF( M.EQ.0 .OR. N.EQ.0 ) THEN |
IF( M.EQ.0 .OR. N.EQ.0 ) THEN |
LWKOPT = 1 |
LWKOPT = 1 |
ELSE |
ELSE |
* |
|
* Determine the block size. NB may be at most NBMAX, where |
|
* NBMAX is used to define the local array T. |
|
* |
|
NB = MIN( NBMAX, ILAENV( 1, 'DORMQL', SIDE // TRANS, M, N, |
NB = MIN( NBMAX, ILAENV( 1, 'DORMQL', SIDE // TRANS, M, N, |
$ K, -1 ) ) |
$ K, -1 ) ) |
LWKOPT = NW*NB |
LWKOPT = NW*NB + TSIZE |
END IF |
END IF |
WORK( 1 ) = LWKOPT |
WORK( 1 ) = LWKOPT |
* |
|
IF( LWORK.LT.NW .AND. .NOT.LQUERY ) THEN |
|
INFO = -12 |
|
END IF |
|
END IF |
END IF |
* |
* |
IF( INFO.NE.0 ) THEN |
IF( INFO.NE.0 ) THEN |
Line 276
|
Line 269
|
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, 'DORMQL', SIDE // TRANS, M, N, K, |
NBMIN = MAX( 2, ILAENV( 2, 'DORMQL', 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 296
|
Line 286
|
* |
* |
* 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 320
|
Line 311
|
* H = H(i+ib-1) . . . H(i+1) H(i) |
* H = H(i+ib-1) . . . H(i+1) H(i) |
* |
* |
CALL DLARFT( 'Backward', 'Columnwise', NQ-K+I+IB-1, IB, |
CALL DLARFT( 'Backward', 'Columnwise', NQ-K+I+IB-1, IB, |
$ A( 1, I ), LDA, TAU( I ), T, LDT ) |
$ A( 1, I ), LDA, TAU( I ), WORK( IWT ), LDT ) |
IF( LEFT ) THEN |
IF( LEFT ) THEN |
* |
* |
* H or H**T is applied to C(1:m-k+i+ib-1,1:n) |
* H or H**T is applied to C(1:m-k+i+ib-1,1:n) |
Line 336
|
Line 327
|
* Apply H or H**T |
* Apply H or H**T |
* |
* |
CALL DLARFB( SIDE, TRANS, 'Backward', 'Columnwise', MI, NI, |
CALL DLARFB( SIDE, TRANS, 'Backward', 'Columnwise', MI, NI, |
$ IB, A( 1, I ), LDA, T, LDT, C, LDC, WORK, |
$ IB, A( 1, I ), LDA, WORK( IWT ), LDT, C, LDC, |
$ LDWORK ) |
$ WORK, LDWORK ) |
10 CONTINUE |
10 CONTINUE |
END IF |
END IF |
WORK( 1 ) = LWKOPT |
WORK( 1 ) = LWKOPT |