summaryrefslogtreecommitdiff
path: root/SRC
diff options
context:
space:
mode:
Diffstat (limited to 'SRC')
-rw-r--r--SRC/clahef.f10
-rw-r--r--SRC/clasyf.f10
-rw-r--r--SRC/dlasyf.f10
-rw-r--r--SRC/slasyf.f10
-rw-r--r--SRC/zlahef.f10
-rw-r--r--SRC/zlasyf.f8
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