From c57182219ee505a7456bd91061b202b80924623d Mon Sep 17 00:00:00 2001
From: igor175 <igor175@8a072113-8704-0410-8d35-dd094bca7971>
Date: Mon, 22 Apr 2013 08:57:05 +0000
Subject: [PATCH] changed comments in (c,z)hetf2.f

---
 SRC/chetf2.f |  2 +-
 SRC/zhetf2.f | 26 +++++++++++++++++++++++---
 2 files changed, 24 insertions(+), 4 deletions(-)

diff --git a/SRC/chetf2.f b/SRC/chetf2.f
index 27295bc..17d1321 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 2fb4090..14c5bba 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
 *
-- 
2.7.4