version 1.18, 2018/05/29 07:18:17
|
version 1.19, 2023/08/07 08:39:21
|
Line 177
|
Line 177
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date December 2016 |
|
* |
|
*> \ingroup complex16OTHEReigen |
*> \ingroup complex16OTHEReigen |
* |
* |
* ===================================================================== |
* ===================================================================== |
SUBROUTINE ZGGGLM( N, M, P, A, LDA, B, LDB, D, X, Y, WORK, LWORK, |
SUBROUTINE ZGGGLM( N, M, P, A, LDA, B, LDB, D, X, Y, WORK, LWORK, |
$ INFO ) |
$ INFO ) |
* |
* |
* -- LAPACK driver routine (version 3.7.0) -- |
* -- LAPACK driver routine -- |
* -- 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..-- |
* December 2016 |
|
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER INFO, LDA, LDB, LWORK, M, N, P |
INTEGER INFO, LDA, LDB, LWORK, M, N, P |
Line 271
|
Line 268
|
* |
* |
* Quick return if possible |
* Quick return if possible |
* |
* |
IF( N.EQ.0 ) |
IF( N.EQ.0 ) THEN |
$ RETURN |
DO I = 1, M |
|
X(I) = CZERO |
|
END DO |
|
DO I = 1, P |
|
Y(I) = CZERO |
|
END DO |
|
RETURN |
|
END IF |
* |
* |
* Compute the GQR factorization of matrices A and B: |
* Compute the GQR factorization of matrices A and B: |
* |
* |
Line 285
|
Line 289
|
* |
* |
CALL ZGGQRF( N, M, P, A, LDA, WORK, B, LDB, WORK( M+1 ), |
CALL ZGGQRF( N, M, P, A, LDA, WORK, B, LDB, WORK( M+1 ), |
$ WORK( M+NP+1 ), LWORK-M-NP, INFO ) |
$ WORK( M+NP+1 ), LWORK-M-NP, INFO ) |
LOPT = WORK( M+NP+1 ) |
LOPT = INT( WORK( M+NP+1 ) ) |
* |
* |
* Update left-hand-side vector d = Q**H*d = ( d1 ) M |
* Update left-hand-side vector d = Q**H*d = ( d1 ) M |
* ( d2 ) N-M |
* ( d2 ) N-M |