diff options
| author | langou <julien.langou@ucdenver.edu> | 2017-06-12 05:16:16 +0000 |
|---|---|---|
| committer | GitHub <noreply@github.com> | 2017-06-12 05:16:16 +0000 |
| commit | a26a82ff4a30a9240767621e8f0fe30390c55446 (patch) | |
| tree | 9cb06867074c60a318b468c6c13d7c7044eaef21 /SRC/zlasyf_aa.f | |
| parent | 5c392680ef6362c234c2f2ab3dadd9a9148ffb8f (diff) | |
| parent | fd5cae5d069bf8544067a105e6c46d7c02598865 (diff) | |
| download | lapack-a26a82ff4a30a9240767621e8f0fe30390c55446.tar.gz lapack-a26a82ff4a30a9240767621e8f0fe30390c55446.tar.bz2 lapack-a26a82ff4a30a9240767621e8f0fe30390c55446.zip | |
Merge pull request #166 from iyamazaki/master
fixing the definition of "info" in Aasen's
Diffstat (limited to 'SRC/zlasyf_aa.f')
| -rw-r--r-- | SRC/zlasyf_aa.f | 44 |
1 files changed, 4 insertions, 40 deletions
diff --git a/SRC/zlasyf_aa.f b/SRC/zlasyf_aa.f index 3e1b7639..b6c229c7 100644 --- a/SRC/zlasyf_aa.f +++ b/SRC/zlasyf_aa.f @@ -19,11 +19,11 @@ * =========== * * SUBROUTINE ZLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV, -* H, LDH, WORK, INFO ) +* H, LDH, WORK ) * * .. Scalar Arguments .. * CHARACTER UPLO -* INTEGER J1, M, NB, LDA, LDH, INFO +* INTEGER J1, M, NB, LDA, LDH * .. * .. Array Arguments .. * INTEGER IPIV( * ) @@ -127,16 +127,6 @@ *> WORK is COMPLEX*16 workspace, dimension (M). *> \endverbatim *> -*> \param[out] INFO -*> \verbatim -*> INFO is INTEGER -*> = 0: successful exit -*> < 0: if INFO = -i, the i-th argument had an illegal value -*> > 0: if INFO = i, D(i,i) is exactly zero. The factorization -*> has been completed, but the block diagonal matrix D is -*> exactly singular, and division by zero will occur if it -*> is used to solve a system of equations. -*> \endverbatim * * Authors: * ======== @@ -152,7 +142,7 @@ * * ===================================================================== SUBROUTINE ZLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV, - $ H, LDH, WORK, INFO ) + $ H, LDH, WORK ) * * -- LAPACK computational routine (version 3.7.1) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- @@ -163,7 +153,7 @@ * * .. Scalar Arguments .. CHARACTER UPLO - INTEGER M, NB, J1, LDA, LDH, INFO + INTEGER M, NB, J1, LDA, LDH * .. * .. Array Arguments .. INTEGER IPIV( * ) @@ -192,7 +182,6 @@ * .. * .. Executable Statements .. * - INFO = 0 J = 1 * * K1 is the first column of the panel to be factorized @@ -315,12 +304,6 @@ * Set A(J, J+1) = T(J, J+1) * A( K, J+1 ) = WORK( 2 ) - IF( (A( K, J ).EQ.ZERO ) .AND. (A( K, J+1 ).EQ.ZERO) .AND. - $ ((K.EQ.1) .OR. (A( K-1, J ).EQ.ZERO)) ) THEN - IF(INFO .EQ. 0) THEN - INFO = J - ENDIF - END IF * IF( J.LT.NB ) THEN * @@ -341,13 +324,6 @@ CALL ZLASET( 'Full', 1, M-J-1, ZERO, ZERO, $ A( K, J+2 ), LDA) END IF - ELSE - IF( (A( K, J ).EQ.ZERO) .AND. - $ ((K.EQ.1) .OR. (A( J-1, J ).EQ.ZERO)) ) THEN - IF (INFO.EQ.0) THEN - INFO = J - END IF - END IF END IF J = J + 1 GO TO 10 @@ -468,11 +444,6 @@ * Set A(J+1, J) = T(J+1, J) * A( J+1, K ) = WORK( 2 ) - IF( (A( J, K ).EQ.ZERO) .AND. (A( J+1, K ).EQ.ZERO) .AND. - $ ((K.EQ.1) .OR. (A( J, K-1 ).EQ.ZERO)) ) THEN - IF (INFO .EQ. 0) - $ INFO = J - END IF * IF( J.LT.NB ) THEN * @@ -493,13 +464,6 @@ CALL ZLASET( 'Full', M-J-1, 1, ZERO, ZERO, $ A( J+2, K ), LDA ) END IF - ELSE - IF( (A( J, K ).EQ.ZERO) .AND. - $ ((K.EQ.1) .OR. (A( J, K-1 ).EQ.ZERO)) ) THEN - IF (INFO.EQ.0) THEN - INFO = J - END IF - END IF END IF J = J + 1 GO TO 30 |