updated comments in (s,d,c,z)lasyf.f and (c,z)lahef.f

6 files changed, 34 insertions, 24 deletions

diff --git a/SRC/clahef.f b/SRC/clahef.f
index 5007928e..9ba5694b 100644
--- a/SRC/clahef.f
+++ b/SRC/clahef.f
@@ -269,7 +269,8 @@
          ABSAKK = ABS( REAL( W( K, KW ) ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+*        column K, and COLMAX is its absolute value.
+*        Determine both COLMAX and IMAX.
 *
          IF( K.GT.1 ) THEN
             IMAX = ICAMAX( K-1, W( 1, KW ), 1 )
@@ -280,7 +281,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K
@@ -584,7 +585,8 @@
          ABSAKK = ABS( REAL( W( K, K ) ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+*        column K, and COLMAX is its absolute value.
+*        Determine both COLMAX and IMAX.
 *
          IF( K.LT.N ) THEN
             IMAX = K + ICAMAX( N-K, W( K+1, K ), 1 )
@@ -595,7 +597,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K
diff --git a/SRC/clasyf.f b/SRC/clasyf.f
index c091dfd5..6dd8c07d 100644
--- a/SRC/clasyf.f
+++ b/SRC/clasyf.f
@@ -266,7 +266,8 @@
          ABSAKK = CABS1( W( K, KW ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+*        column K, and COLMAX is its absolute value.
+*        Determine both COLMAX and IMAX.
 *
          IF( K.GT.1 ) THEN
             IMAX = ICAMAX( K-1, W( 1, KW ), 1 )
@@ -277,7 +278,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K
@@ -561,7 +562,8 @@
          ABSAKK = CABS1( W( K, K ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+*        column K, and COLMAX is its absolute value.
+*        Determine both COLMAX and IMAX.
 *
          IF( K.LT.N ) THEN
             IMAX = K + ICAMAX( N-K, W( K+1, K ), 1 )
@@ -572,7 +574,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K
diff --git a/SRC/dlasyf.f b/SRC/dlasyf.f
index 84a4a300..3f1d65fe 100644
--- a/SRC/dlasyf.f
+++ b/SRC/dlasyf.f
@@ -257,7 +257,8 @@
          ABSAKK = ABS( W( K, KW ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+*        column K, and COLMAX is its absolute value.
+*        Determine both COLMAX and IMAX.
 *
          IF( K.GT.1 ) THEN
             IMAX = IDAMAX( K-1, W( 1, KW ), 1 )
@@ -268,7 +269,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K
@@ -552,7 +553,8 @@
          ABSAKK = ABS( W( K, K ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+*        column K, and COLMAX is its absolute value.
+*        Determine both COLMAX and IMAX.
 *
          IF( K.LT.N ) THEN
             IMAX = K + IDAMAX( N-K, W( K+1, K ), 1 )
@@ -563,7 +565,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K
diff --git a/SRC/slasyf.f b/SRC/slasyf.f
index d6fdfd3f..d5fd6b5b 100644
--- a/SRC/slasyf.f
+++ b/SRC/slasyf.f
@@ -257,7 +257,8 @@
          ABSAKK = ABS( W( K, KW ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+*        column K, and COLMAX is its absolute value.
+*        Determine both COLMAX and IMAX.
 *
          IF( K.GT.1 ) THEN
             IMAX = ISAMAX( K-1, W( 1, KW ), 1 )
@@ -268,7 +269,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K
@@ -552,7 +553,8 @@
          ABSAKK = ABS( W( K, K ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+*        column K, and COLMAX is its absolute value.
+*        Determine both COLMAX and IMAX.
 *
          IF( K.LT.N ) THEN
             IMAX = K + ISAMAX( N-K, W( K+1, K ), 1 )
@@ -563,7 +565,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K
diff --git a/SRC/zlahef.f b/SRC/zlahef.f
index c552401c..b7d5e47a 100644
--- a/SRC/zlahef.f
+++ b/SRC/zlahef.f
@@ -269,7 +269,8 @@
          ABSAKK = ABS( DBLE( W( K, KW ) ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+*        column K, and COLMAX is its absolute value.
+*        Determine both COLMAX and IMAX.
 *
          IF( K.GT.1 ) THEN
             IMAX = IZAMAX( K-1, W( 1, KW ), 1 )
@@ -280,7 +281,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K
@@ -592,7 +593,8 @@
          ABSAKK = ABS( DBLE( W( K, K ) ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+*        column K, and COLMAX is its absolute value.
+*        Determine both COLMAX and IMAX.
 *
          IF( K.LT.N ) THEN
             IMAX = K + IZAMAX( N-K, W( K+1, K ), 1 )
@@ -603,7 +605,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K
diff --git a/SRC/zlasyf.f b/SRC/zlasyf.f
index a50b8cf1..4cd32298 100644
--- a/SRC/zlasyf.f
+++ b/SRC/zlasyf.f
@@ -266,7 +266,7 @@
          ABSAKK = CABS1( W( K, KW ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+
 *
          IF( K.GT.1 ) THEN
             IMAX = IZAMAX( K-1, W( 1, KW ), 1 )
@@ -277,7 +277,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K
@@ -561,7 +561,7 @@
          ABSAKK = CABS1( W( K, K ) )
 *
 *        IMAX is the row-index of the largest off-diagonal element in
-*        column K, and COLMAX is its absolute value
+
 *
          IF( K.LT.N ) THEN
             IMAX = K + IZAMAX( N-K, W( K+1, K ), 1 )
@@ -572,7 +572,7 @@
 *
          IF( MAX( ABSAKK, COLMAX ).EQ.ZERO ) THEN
 *
-*           Column K is zero: set INFO and continue
+*           Column K is zero or underflow: set INFO and continue
 *
             IF( INFO.EQ.0 )
      $         INFO = K