diff options
| author | igor175 <igor175@8a072113-8704-0410-8d35-dd094bca7971> | 2013-04-22 08:57:05 +0000 |
|---|---|---|
| committer | igor175 <igor175@8a072113-8704-0410-8d35-dd094bca7971> | 2013-04-22 08:57:05 +0000 |
| commit | c57182219ee505a7456bd91061b202b80924623d (patch) | |
| tree | c91e60c477bcff8492a8706859560fb085dead0e | |
| parent | 9523660b3a2574b95a022006eb62a51a558443e2 (diff) | |
| download | lapack-c57182219ee505a7456bd91061b202b80924623d.tar.gz lapack-c57182219ee505a7456bd91061b202b80924623d.tar.bz2 lapack-c57182219ee505a7456bd91061b202b80924623d.zip | |
changed comments in (c,z)hetf2.f
| -rw-r--r-- | SRC/chetf2.f | 2 | ||||
| -rw-r--r-- | SRC/zhetf2.f | 26 |
2 files changed, 24 insertions, 4 deletions
diff --git a/SRC/chetf2.f b/SRC/chetf2.f index 27295bc7..17d13216 100644 --- a/SRC/chetf2.f +++ b/SRC/chetf2.f @@ -1,4 +1,4 @@ -*> \brief \b CHETF2 computes the factorization of a complex Hermitian matrix, using the diagonal pivoting method (unblocked algorithm). +*> \brief \b CHETF2 computes the factorization of a complex Hermitian matrix, using the diagonal pivoting method (unblocked algorithm calling Level 2 BLAS). * * =========== DOCUMENTATION =========== * diff --git a/SRC/zhetf2.f b/SRC/zhetf2.f index 2fb40902..14c5bba5 100644 --- a/SRC/zhetf2.f +++ b/SRC/zhetf2.f @@ -1,4 +1,4 @@ -*> \brief \b ZHETF2 computes the factorization of a complex Hermitian matrix, using the diagonal pivoting method (unblocked algorithm). +*> \brief \b ZHETF2 computes the factorization of a complex Hermitian matrix, using the diagonal pivoting method (unblocked algorithm, calling Level 2 BLAS). * * =========== DOCUMENTATION =========== * @@ -302,6 +302,11 @@ KP = K A( K, K ) = DBLE( A( K, K ) ) ELSE +* +* ============================================================ +* +* Test for interchange +* IF( ABSAKK.GE.ALPHA*COLMAX ) THEN * * no interchange, use 1-by-1 pivot block @@ -310,7 +315,8 @@ ELSE * * JMAX is the column-index of the largest off-diagonal -* element in row IMAX, and ROWMAX is its absolute value +* element in row IMAX, and ROWMAX is its absolute value. +* Determine only ROWMAX. * JMAX = IMAX + IZAMAX( K-IMAX, A( IMAX, IMAX+1 ), LDA ) ROWMAX = CABS1( A( IMAX, JMAX ) ) @@ -324,6 +330,7 @@ * no interchange, use 1-by-1 pivot block * KP = K +* ELSE IF( ABS( DBLE( A( IMAX, IMAX ) ) ).GE.ALPHA*ROWMAX ) $ THEN * @@ -339,8 +346,11 @@ KP = IMAX KSTEP = 2 END IF +* END IF * +* ============================================================ +* KK = K - KSTEP + 1 IF( KP.NE.KK ) THEN * @@ -487,6 +497,11 @@ KP = K A( K, K ) = DBLE( A( K, K ) ) ELSE +* +* ============================================================ +* +* Test for interchange +* IF( ABSAKK.GE.ALPHA*COLMAX ) THEN * * no interchange, use 1-by-1 pivot block @@ -495,7 +510,8 @@ ELSE * * JMAX is the column-index of the largest off-diagonal -* element in row IMAX, and ROWMAX is its absolute value +* element in row IMAX, and ROWMAX is its absolute value. +* Determine only ROWMAX. * JMAX = K - 1 + IZAMAX( IMAX-K, A( IMAX, K ), LDA ) ROWMAX = CABS1( A( IMAX, JMAX ) ) @@ -509,6 +525,7 @@ * no interchange, use 1-by-1 pivot block * KP = K +* ELSE IF( ABS( DBLE( A( IMAX, IMAX ) ) ).GE.ALPHA*ROWMAX ) $ THEN * @@ -524,8 +541,11 @@ KP = IMAX KSTEP = 2 END IF +* END IF * +* ============================================================ +* KK = K + KSTEP - 1 IF( KP.NE.KK ) THEN * |