diff options
author | Julien Langou <julien.langou@ucdenver.edu> | 2016-11-03 08:48:54 +0100 |
---|---|---|
committer | Julien Langou <julien.langou@ucdenver.edu> | 2016-11-03 08:48:54 +0100 |
commit | bbff7393714b29a6ff70e8c1565784cb16a0e746 (patch) | |
tree | 146d4bc89db9148abc4b19e2e453623f1a6c7bbe /SRC/zgelq.f | |
parent | bd47060bcb3a470520622de69ac1426ca4186f5e (diff) |
Lots of trailing whitespaces in the files of Syd. Cleaning this. No big deal.
Diffstat (limited to 'SRC/zgelq.f')
-rw-r--r-- | SRC/zgelq.f | 78 |
1 files changed, 39 insertions, 39 deletions
diff --git a/SRC/zgelq.f b/SRC/zgelq.f index 2e188df9..33125b3d 100644 --- a/SRC/zgelq.f +++ b/SRC/zgelq.f @@ -1,26 +1,26 @@ -* +* * Definition: * =========== * -* SUBROUTINE CGELQ( M, N, A, LDA, WORK1, LWORK1, WORK2, LWORK2, +* SUBROUTINE CGELQ( M, N, A, LDA, WORK1, LWORK1, WORK2, LWORK2, * INFO) -* +* * .. Scalar Arguments .. * INTEGER INFO, LDA, M, N, LWORK1, LWORK2 * .. * .. Array Arguments .. * COMPLEX*16 A( LDA, * ), WORK1( * ), WORK2( * ) * .. -* +* * *> \par Purpose: * ============= *> *> \verbatim -*> -*> ZGELQ computes an LQ factorization of an M-by-N matrix A, -*> using ZLASWLQ when A is short and wide -*> (N sufficiently greater than M), and otherwise ZGELQT: +*> +*> ZGELQ computes an LQ factorization of an M-by-N matrix A, +*> using ZLASWLQ when A is short and wide +*> (N sufficiently greater than M), and otherwise ZGELQT: *> A = L * Q . *> \endverbatim * @@ -43,10 +43,10 @@ *> \verbatim *> A is COMPLEX*16 array, dimension (LDA,N) *> On entry, the M-by-N matrix A. -*> On exit, the elements on and below the diagonal of the array -*> contain the M-by-min(M,N) lower trapezoidal matrix L +*> On exit, the elements on and below the diagonal of the array +*> contain the M-by-min(M,N) lower trapezoidal matrix L *> (L is lower triangular if M <= N); -*> the elements above the diagonal are the rows of +*> the elements above the diagonal are the rows of *> blocked V representing Q (see Further Details). *> \endverbatim *> @@ -60,13 +60,13 @@ *> \verbatim *> WORK1 is COMPLEX*16 array, dimension (MAX(1,LWORK1)) *> WORK1 contains part of the data structure used to store Q. -*> WORK1(1): algorithm type = 1, to indicate output from +*> WORK1(1): algorithm type = 1, to indicate output from *> ZLASWLQ or ZGELQT *> WORK1(2): optimum size of WORK1 *> WORK1(3): minimum size of WORK1 *> WORK1(4): horizontal block size *> WORK1(5): vertical block size -*> WORK1(6:LWORK1): data structure needed for Q, computed by +*> WORK1(6:LWORK1): data structure needed for Q, computed by *> ZLASWLQ or ZGELQT *> \endverbatim *> @@ -74,25 +74,25 @@ *> \verbatim *> LWORK1 is INTEGER *> The dimension of the array WORK1. -*> If LWORK1 = -1, then a query is assumed. In this case the +*> If LWORK1 = -1, then a query is assumed. In this case the *> routine calculates the optimal size of WORK1 and -*> returns this value in WORK1(2), and calculates the minimum -*> size of WORK1 and returns this value in WORK1(3). -*> No error message related to LWORK1 is issued by XERBLA when +*> returns this value in WORK1(2), and calculates the minimum +*> size of WORK1 and returns this value in WORK1(3). +*> No error message related to LWORK1 is issued by XERBLA when *> LWORK1 = -1. *> \endverbatim *> *> \param[out] WORK2 *> \verbatim *> (workspace) COMPLEX*16 array, dimension (MAX(1,LWORK2)) -*> +*> *> \endverbatim *> \param[in] LWORK2 *> \verbatim *> LWORK2 is INTEGER *> The dimension of the array WORK2. *> If LWORK2 = -1, then a query is assumed. In this case the -*> routine calculates the optimal size of WORK2 and +*> routine calculates the optimal size of WORK2 and *> returns this value in WORK2(1), and calculates the minimum *> size of WORK2 and returns this value in WORK2(2). *> No error message related to LWORK2 is issued by XERBLA when @@ -121,19 +121,19 @@ *> Depending on the matrix dimensions M and N, and row and column *> block sizes MB and NB returned by ILAENV, GELQ will use either *> LASWLQ(if the matrix is short-and-wide) or GELQT to compute -*> the LQ decomposition. +*> the LQ decomposition. *> The output of LASWLQ or GELQT representing Q is stored in A and in -*> array WORK1(6:LWORK1) for later use. -*> WORK1(2:5) contains the matrix dimensions M,N and block sizes MB, NB -*> which are needed to interpret A and WORK1(6:LWORK1) for later use. -*> WORK1(1)=1 indicates that the code needed to take WORK1(2:5) and +*> array WORK1(6:LWORK1) for later use. +*> WORK1(2:5) contains the matrix dimensions M,N and block sizes MB, NB +*> which are needed to interpret A and WORK1(6:LWORK1) for later use. +*> WORK1(1)=1 indicates that the code needed to take WORK1(2:5) and *> decide whether LASWLQ or GELQT was used is the same as used below in -*> GELQ. For a detailed description of A and WORK1(6:LWORK1), see +*> GELQ. For a detailed description of A and WORK1(6:LWORK1), see *> Further Details in LASWLQ or GELQT. *> \endverbatim *> * ===================================================================== - SUBROUTINE ZGELQ( M, N, A, LDA, WORK1, LWORK1, WORK2, LWORK2, + SUBROUTINE ZGELQ( M, N, A, LDA, WORK1, LWORK1, WORK2, LWORK2, $ INFO) * * -- LAPACK computational routine (version 3.5.0) -- @@ -175,8 +175,8 @@ * LQUERY = ( LWORK1.EQ.-1 .OR. LWORK2.EQ.-1 ) * -* Determine the block size -* +* Determine the block size +* IF ( MIN(M,N).GT.0 ) THEN MB = ILAENV( 1, 'ZGELQ ', ' ', M, N, 1, -1) NB = ILAENV( 1, 'ZGELQ ', ' ', M, N, 2, -1) @@ -198,18 +198,18 @@ END IF * * Determine if the workspace size satisfies minimum size -* - LMINWS = .FALSE. +* + LMINWS = .FALSE. IF((LWORK1.LT.MAX(1,MB*M*NBLCKS+5) $ .OR.(LWORK2.LT.MB*M)).AND.(LWORK2.GE.M).AND.(LWORK1.GE.M+5) $ .AND.(.NOT.LQUERY)) THEN IF (LWORK1.LT.MAX(1,MB*M*NBLCKS+5)) THEN LMINWS = .TRUE. MB = 1 - END IF + END IF IF (LWORK1.LT.MAX(1,M*NBLCKS+5)) THEN LMINWS = .TRUE. - NB = N + NB = N END IF IF (LWORK2.LT.MB*M) THEN LMINWS = .TRUE. @@ -223,13 +223,13 @@ INFO = -2 ELSE IF( LDA.LT.MAX( 1, M ) ) THEN INFO = -4 - ELSE IF( LWORK1.LT.MAX( 1, MB*M*NBLCKS+5 ) + ELSE IF( LWORK1.LT.MAX( 1, MB*M*NBLCKS+5 ) $ .AND.(.NOT.LQUERY).AND. (.NOT.LMINWS)) THEN INFO = -6 ELSE IF( (LWORK2.LT.MAX(1,M*MB)).AND.(.NOT.LQUERY) $ .AND.(.NOT.LMINWS) ) THEN - INFO = -8 - END IF + INFO = -8 + END IF * IF( INFO.EQ.0) THEN WORK1(1) = 1 @@ -257,12 +257,12 @@ * IF((N.LE.M).OR.(NB.LE.M).OR.(NB.GE.N)) THEN CALL ZGELQT( M, N, MB, A, LDA, WORK1(6), MB, WORK2, INFO) - ELSE - CALL ZLASWLQ( M, N, MB, NB, A, LDA, WORK1(6), MB, WORK2, + ELSE + CALL ZLASWLQ( M, N, MB, NB, A, LDA, WORK1(6), MB, WORK2, $ LWORK2, INFO) END IF RETURN -* +* * End of ZGELQ * - END
\ No newline at end of file + END |