version 1.7, 2012/08/22 09:48:42
|
version 1.15, 2018/05/29 07:18:41
|
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 ZUNCSD + dependencies |
*> Download ZUNCSD + dependencies |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zuncsd.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zuncsd.f"> |
*> [TGZ]</a> |
*> [TGZ]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zuncsd.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zuncsd.f"> |
*> [ZIP]</a> |
*> [ZIP]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zuncsd.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zuncsd.f"> |
*> [TXT]</a> |
*> [TXT]</a> |
*> \endhtmlonly |
*> \endhtmlonly |
* |
* |
* Definition: |
* Definition: |
* =========== |
* =========== |
Line 24
|
Line 24
|
* U1, LDU1, U2, LDU2, V1T, LDV1T, V2T, |
* U1, LDU1, U2, LDU2, V1T, LDV1T, V2T, |
* LDV2T, WORK, LWORK, RWORK, LRWORK, |
* LDV2T, WORK, LWORK, RWORK, LRWORK, |
* IWORK, INFO ) |
* IWORK, INFO ) |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
* CHARACTER JOBU1, JOBU2, JOBV1T, JOBV2T, SIGNS, TRANS |
* CHARACTER JOBU1, JOBU2, JOBV1T, JOBV2T, SIGNS, TRANS |
* INTEGER INFO, LDU1, LDU2, LDV1T, LDV2T, LDX11, LDX12, |
* INTEGER INFO, LDU1, LDU2, LDV1T, LDV2T, LDX11, LDX12, |
Line 39
|
Line 39
|
* $ X12( LDX12, * ), X21( LDX21, * ), X22( LDX22, |
* $ X12( LDX12, * ), X21( LDX21, * ), X22( LDX22, |
* $ * ) |
* $ * ) |
* .. |
* .. |
* |
* |
* |
* |
*> \par Purpose: |
*> \par Purpose: |
* ============= |
* ============= |
Line 188
|
Line 188
|
*> |
*> |
*> \param[out] U1 |
*> \param[out] U1 |
*> \verbatim |
*> \verbatim |
*> U1 is COMPLEX*16 array, dimension (P) |
*> U1 is COMPLEX*16 array, dimension (LDU1,P) |
*> If JOBU1 = 'Y', U1 contains the P-by-P unitary matrix U1. |
*> If JOBU1 = 'Y', U1 contains the P-by-P unitary matrix U1. |
*> \endverbatim |
*> \endverbatim |
*> |
*> |
Line 201
|
Line 201
|
*> |
*> |
*> \param[out] U2 |
*> \param[out] U2 |
*> \verbatim |
*> \verbatim |
*> U2 is COMPLEX*16 array, dimension (M-P) |
*> U2 is COMPLEX*16 array, dimension (LDU2,M-P) |
*> If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) unitary |
*> If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) unitary |
*> matrix U2. |
*> matrix U2. |
*> \endverbatim |
*> \endverbatim |
Line 215
|
Line 215
|
*> |
*> |
*> \param[out] V1T |
*> \param[out] V1T |
*> \verbatim |
*> \verbatim |
*> V1T is COMPLEX*16 array, dimension (Q) |
*> V1T is COMPLEX*16 array, dimension (LDV1T,Q) |
*> If JOBV1T = 'Y', V1T contains the Q-by-Q matrix unitary |
*> If JOBV1T = 'Y', V1T contains the Q-by-Q matrix unitary |
*> matrix V1**H. |
*> matrix V1**H. |
*> \endverbatim |
*> \endverbatim |
Line 229
|
Line 229
|
*> |
*> |
*> \param[out] V2T |
*> \param[out] V2T |
*> \verbatim |
*> \verbatim |
*> V2T is COMPLEX*16 array, dimension (M-Q) |
*> V2T is COMPLEX*16 array, dimension (LDV2T,M-Q) |
*> If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) unitary |
*> If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) unitary |
*> matrix V2**H. |
*> matrix V2**H. |
*> \endverbatim |
*> \endverbatim |
Line 303
|
Line 303
|
* 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 June 2017 |
* |
* |
*> \ingroup complex16OTHERcomputational |
*> \ingroup complex16OTHERcomputational |
* |
* |
Line 320
|
Line 320
|
$ LDV2T, WORK, LWORK, RWORK, LRWORK, |
$ LDV2T, WORK, LWORK, RWORK, LRWORK, |
$ IWORK, INFO ) |
$ IWORK, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.4.0) -- |
* -- LAPACK computational routine (version 3.7.1) -- |
* -- 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 |
* June 2017 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER JOBU1, JOBU2, JOBV1T, JOBV2T, SIGNS, TRANS |
CHARACTER JOBU1, JOBU2, JOBV1T, JOBV2T, SIGNS, TRANS |
Line 356
|
Line 356
|
$ LBBCSDWORKOPT, LORBDBWORK, LORBDBWORKMIN, |
$ LBBCSDWORKOPT, LORBDBWORK, LORBDBWORKMIN, |
$ LORBDBWORKOPT, LORGLQWORK, LORGLQWORKMIN, |
$ LORBDBWORKOPT, LORGLQWORK, LORGLQWORKMIN, |
$ LORGLQWORKOPT, LORGQRWORK, LORGQRWORKMIN, |
$ LORGLQWORKOPT, LORGQRWORK, LORGQRWORKMIN, |
$ LORGQRWORKOPT, LWORKMIN, LWORKOPT |
$ LORGQRWORKOPT, LWORKMIN, LWORKOPT, P1, Q1 |
LOGICAL COLMAJOR, DEFAULTSIGNS, LQUERY, WANTU1, WANTU2, |
LOGICAL COLMAJOR, DEFAULTSIGNS, LQUERY, WANTU1, WANTU2, |
$ WANTV1T, WANTV2T |
$ WANTV1T, WANTV2T |
INTEGER LRWORKMIN, LRWORKOPT |
INTEGER LRWORKMIN, LRWORKOPT |
LOGICAL LRQUERY |
LOGICAL LRQUERY |
* .. |
* .. |
* .. External Subroutines .. |
* .. External Subroutines .. |
EXTERNAL XERBLA, ZBBCSD, ZLACPY, ZLAPMR, ZLAPMT, ZLASCL, |
EXTERNAL XERBLA, ZBBCSD, ZLACPY, ZLAPMR, ZLAPMT, |
$ ZLASET, ZUNBDB, ZUNGLQ, ZUNGQR |
$ ZUNBDB, ZUNGLQ, ZUNGQR |
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
LOGICAL LSAME |
LOGICAL LSAME |
EXTERNAL LSAME |
EXTERNAL LSAME |
* .. |
* .. |
* .. Intrinsic Functions |
* .. Intrinsic Functions |
INTRINSIC COS, INT, MAX, MIN, SIN |
INTRINSIC INT, MAX, MIN |
* .. |
* .. |
* .. Executable Statements .. |
* .. Executable Statements .. |
* |
* |
Line 392
|
Line 392
|
INFO = -8 |
INFO = -8 |
ELSE IF( Q .LT. 0 .OR. Q .GT. M ) THEN |
ELSE IF( Q .LT. 0 .OR. Q .GT. M ) THEN |
INFO = -9 |
INFO = -9 |
ELSE IF( ( COLMAJOR .AND. LDX11 .LT. MAX(1,P) ) .OR. |
ELSE IF ( COLMAJOR .AND. LDX11 .LT. MAX( 1, P ) ) THEN |
$ ( .NOT.COLMAJOR .AND. LDX11 .LT. MAX(1,Q) ) ) THEN |
INFO = -11 |
INFO = -11 |
ELSE IF (.NOT. COLMAJOR .AND. LDX11 .LT. MAX( 1, Q ) ) THEN |
|
INFO = -11 |
|
ELSE IF (COLMAJOR .AND. LDX12 .LT. MAX( 1, P ) ) THEN |
|
INFO = -13 |
|
ELSE IF (.NOT. COLMAJOR .AND. LDX12 .LT. MAX( 1, M-Q ) ) THEN |
|
INFO = -13 |
|
ELSE IF (COLMAJOR .AND. LDX21 .LT. MAX( 1, M-P ) ) THEN |
|
INFO = -15 |
|
ELSE IF (.NOT. COLMAJOR .AND. LDX21 .LT. MAX( 1, Q ) ) THEN |
|
INFO = -15 |
|
ELSE IF (COLMAJOR .AND. LDX22 .LT. MAX( 1, M-P ) ) THEN |
|
INFO = -17 |
|
ELSE IF (.NOT. COLMAJOR .AND. LDX22 .LT. MAX( 1, M-Q ) ) THEN |
|
INFO = -17 |
ELSE IF( WANTU1 .AND. LDU1 .LT. P ) THEN |
ELSE IF( WANTU1 .AND. LDU1 .LT. P ) THEN |
INFO = -20 |
INFO = -20 |
ELSE IF( WANTU2 .AND. LDU2 .LT. M-P ) THEN |
ELSE IF( WANTU2 .AND. LDU2 .LT. M-P ) THEN |
Line 458
|
Line 471
|
IB22D = IB21E + MAX( 1, Q - 1 ) |
IB22D = IB21E + MAX( 1, Q - 1 ) |
IB22E = IB22D + MAX( 1, Q ) |
IB22E = IB22D + MAX( 1, Q ) |
IBBCSD = IB22E + MAX( 1, Q - 1 ) |
IBBCSD = IB22E + MAX( 1, Q - 1 ) |
CALL ZBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q, 0, |
CALL ZBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q, |
$ 0, U1, LDU1, U2, LDU2, V1T, LDV1T, V2T, LDV2T, 0, |
$ THETA, THETA, U1, LDU1, U2, LDU2, V1T, LDV1T, |
$ 0, 0, 0, 0, 0, 0, 0, RWORK, -1, CHILDINFO ) |
$ V2T, LDV2T, THETA, THETA, THETA, THETA, THETA, |
|
$ THETA, THETA, THETA, RWORK, -1, CHILDINFO ) |
LBBCSDWORKOPT = INT( RWORK(1) ) |
LBBCSDWORKOPT = INT( RWORK(1) ) |
LBBCSDWORKMIN = LBBCSDWORKOPT |
LBBCSDWORKMIN = LBBCSDWORKOPT |
LRWORKOPT = IBBCSD + LBBCSDWORKOPT - 1 |
LRWORKOPT = IBBCSD + LBBCSDWORKOPT - 1 |
Line 474
|
Line 488
|
ITAUQ1 = ITAUP2 + MAX( 1, M - P ) |
ITAUQ1 = ITAUP2 + MAX( 1, M - P ) |
ITAUQ2 = ITAUQ1 + MAX( 1, Q ) |
ITAUQ2 = ITAUQ1 + MAX( 1, Q ) |
IORGQR = ITAUQ2 + MAX( 1, M - Q ) |
IORGQR = ITAUQ2 + MAX( 1, M - Q ) |
CALL ZUNGQR( M-Q, M-Q, M-Q, 0, MAX(1,M-Q), 0, WORK, -1, |
CALL ZUNGQR( M-Q, M-Q, M-Q, U1, MAX(1,M-Q), U1, WORK, -1, |
$ CHILDINFO ) |
$ CHILDINFO ) |
LORGQRWORKOPT = INT( WORK(1) ) |
LORGQRWORKOPT = INT( WORK(1) ) |
LORGQRWORKMIN = MAX( 1, M - Q ) |
LORGQRWORKMIN = MAX( 1, M - Q ) |
IORGLQ = ITAUQ2 + MAX( 1, M - Q ) |
IORGLQ = ITAUQ2 + MAX( 1, M - Q ) |
CALL ZUNGLQ( M-Q, M-Q, M-Q, 0, MAX(1,M-Q), 0, WORK, -1, |
CALL ZUNGLQ( M-Q, M-Q, M-Q, U1, MAX(1,M-Q), U1, WORK, -1, |
$ CHILDINFO ) |
$ CHILDINFO ) |
LORGLQWORKOPT = INT( WORK(1) ) |
LORGLQWORKOPT = INT( WORK(1) ) |
LORGLQWORKMIN = MAX( 1, M - Q ) |
LORGLQWORKMIN = MAX( 1, M - Q ) |
IORBDB = ITAUQ2 + MAX( 1, M - Q ) |
IORBDB = ITAUQ2 + MAX( 1, M - Q ) |
CALL ZUNBDB( TRANS, SIGNS, M, P, Q, X11, LDX11, X12, LDX12, |
CALL ZUNBDB( TRANS, SIGNS, M, P, Q, X11, LDX11, X12, LDX12, |
$ X21, LDX21, X22, LDX22, 0, 0, 0, 0, 0, 0, WORK, |
$ X21, LDX21, X22, LDX22, THETA, THETA, U1, U2, |
$ -1, CHILDINFO ) |
$ V1T, V2T, WORK, -1, CHILDINFO ) |
LORBDBWORKOPT = INT( WORK(1) ) |
LORBDBWORKOPT = INT( WORK(1) ) |
LORBDBWORKMIN = LORBDBWORKOPT |
LORBDBWORKMIN = LORBDBWORKOPT |
LWORKOPT = MAX( IORGQR + LORGQRWORKOPT, IORGLQ + LORGLQWORKOPT, |
LWORKOPT = MAX( IORGQR + LORGQRWORKOPT, IORGLQ + LORGLQWORKOPT, |
Line 551
|
Line 565
|
END IF |
END IF |
IF( WANTV2T .AND. M-Q .GT. 0 ) THEN |
IF( WANTV2T .AND. M-Q .GT. 0 ) THEN |
CALL ZLACPY( 'U', P, M-Q, X12, LDX12, V2T, LDV2T ) |
CALL ZLACPY( 'U', P, M-Q, X12, LDX12, V2T, LDV2T ) |
CALL ZLACPY( 'U', M-P-Q, M-P-Q, X22(Q+1,P+1), LDX22, |
IF( M-P .GT. Q) THEN |
$ V2T(P+1,P+1), LDV2T ) |
CALL ZLACPY( 'U', M-P-Q, M-P-Q, X22(Q+1,P+1), LDX22, |
CALL ZUNGLQ( M-Q, M-Q, M-Q, V2T, LDV2T, WORK(ITAUQ2), |
$ V2T(P+1,P+1), LDV2T ) |
$ WORK(IORGLQ), LORGLQWORK, INFO ) |
END IF |
|
IF( M .GT. Q ) THEN |
|
CALL ZUNGLQ( M-Q, M-Q, M-Q, V2T, LDV2T, WORK(ITAUQ2), |
|
$ WORK(IORGLQ), LORGLQWORK, INFO ) |
|
END IF |
END IF |
END IF |
ELSE |
ELSE |
IF( WANTU1 .AND. P .GT. 0 ) THEN |
IF( WANTU1 .AND. P .GT. 0 ) THEN |
Line 579
|
Line 597
|
$ WORK(IORGQR), LORGQRWORK, INFO ) |
$ WORK(IORGQR), LORGQRWORK, INFO ) |
END IF |
END IF |
IF( WANTV2T .AND. M-Q .GT. 0 ) THEN |
IF( WANTV2T .AND. M-Q .GT. 0 ) THEN |
|
P1 = MIN( P+1, M ) |
|
Q1 = MIN( Q+1, M ) |
CALL ZLACPY( 'L', M-Q, P, X12, LDX12, V2T, LDV2T ) |
CALL ZLACPY( 'L', M-Q, P, X12, LDX12, V2T, LDV2T ) |
CALL ZLACPY( 'L', M-P-Q, M-P-Q, X22(P+1,Q+1), LDX22, |
IF( M .GT. P+Q ) THEN |
$ V2T(P+1,P+1), LDV2T ) |
CALL ZLACPY( 'L', M-P-Q, M-P-Q, X22(P1,Q1), LDX22, |
|
$ V2T(P+1,P+1), LDV2T ) |
|
END IF |
CALL ZUNGQR( M-Q, M-Q, M-Q, V2T, LDV2T, WORK(ITAUQ2), |
CALL ZUNGQR( M-Q, M-Q, M-Q, V2T, LDV2T, WORK(ITAUQ2), |
$ WORK(IORGQR), LORGQRWORK, INFO ) |
$ WORK(IORGQR), LORGQRWORK, INFO ) |
END IF |
END IF |
Line 599
|
Line 621
|
* Permute rows and columns to place identity submatrices in top- |
* Permute rows and columns to place identity submatrices in top- |
* left corner of (1,1)-block and/or bottom-right corner of (1,2)- |
* left corner of (1,1)-block and/or bottom-right corner of (1,2)- |
* block and/or bottom-right corner of (2,1)-block and/or top-left |
* block and/or bottom-right corner of (2,1)-block and/or top-left |
* corner of (2,2)-block |
* corner of (2,2)-block |
* |
* |
IF( Q .GT. 0 .AND. WANTU2 ) THEN |
IF( Q .GT. 0 .AND. WANTU2 ) THEN |
DO I = 1, Q |
DO I = 1, Q |