changed comments in (c,z)hetf2.f

diff --git a/SRC/chetf2.f b/SRC/chetf2.f
index 27295bc77c3076cf6c946243cd746c2ae74fdb12..17d13216ca5591a26b95526ba9c189c017227930 100644 (file)
--- 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 2fb40902514d3570fb0dea0c2e4aafcf30270f18..14c5bba54cb9c399529fd2c0bf22805119bad46c 100644 (file)
--- 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,7 +346,10 @@
                   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,7 +541,10 @@
                   KP = IMAX
                   KSTEP = 2
                END IF
+*
             END IF
+*
+*           ============================================================
 *
             KK = K + KSTEP - 1
             IF( KP.NE.KK ) THEN