Use short-name (K instead of J1+J-1) for clarity.

diff --git a/SRC/clahef_aa.f b/SRC/clahef_aa.f
index cd07f32..01e8f98 100644 (file)
--- a/SRC/clahef_aa.f
+++ b/SRC/clahef_aa.f
@@ -258,7 +258,7 @@
 *           Compute WORK(2:N) = T(J, J) L(J, (J+1):N)
 *            where A(J, J) stores T(J, J) and A(J-1, (J+1):N) stores U(J, (J+1):N)
 *
-            IF( (J1+J-1).GT.1 ) THEN
+            IF( K.GT.1 ) THEN
                ALPHA = -A( K, J )
                CALL CAXPY( M-J, ALPHA, A( K-1, J+1 ), LDA,
      $                                 WORK( 2 ), 1 )
@@ -412,7 +412,7 @@
 *           Compute WORK(2:N) = T(J, J) L((J+1):N, J)
 *            where A(J, J) = T(J, J) and A((J+1):N, J-1) = L((J+1):N, J)
 *
-            IF( (J1+J-1).GT.1 ) THEN
+            IF( K.GT.1 ) THEN
                ALPHA = -A( J, K )
                CALL CAXPY( M-J, ALPHA, A( J+1, K-1 ), 1,
      $                                 WORK( 2 ), 1 )

diff --git a/SRC/dlasyf_aa.f b/SRC/dlasyf_aa.f
index 5139dd2..cc0b80f 100644 (file)
--- a/SRC/dlasyf_aa.f
+++ b/SRC/dlasyf_aa.f
@@ -256,7 +256,7 @@
 *           Compute WORK(2:N) = T(J, J) L(J, (J+1):N)
 *            where A(J, J) stores T(J, J) and A(J-1, (J+1):N) stores U(J, (J+1):N)
 *
-            IF( (J1+J-1).GT.1 ) THEN
+            IF( K.GT.1 ) THEN
                ALPHA = -A( K, J )
                CALL DAXPY( M-J, ALPHA, A( K-1, J+1 ), LDA,
      $                                 WORK( 2 ), 1 )
@@ -406,7 +406,7 @@
 *           Compute WORK(2:N) = T(J, J) L((J+1):N, J)
 *            where A(J, J) = T(J, J) and A((J+1):N, J-1) = L((J+1):N, J)
 *
-            IF( (J1+J-1).GT.1 ) THEN
+            IF( K.GT.1 ) THEN
                ALPHA = -A( J, K )
                CALL DAXPY( M-J, ALPHA, A( J+1, K-1 ), 1,
      $                                 WORK( 2 ), 1 )

diff --git a/SRC/slasyf_aa.f b/SRC/slasyf_aa.f
index 5a20ef8..953d574 100644 (file)
--- a/SRC/slasyf_aa.f
+++ b/SRC/slasyf_aa.f
@@ -256,7 +256,7 @@
 *           Compute WORK(2:N) = T(J, J) L(J, (J+1):N)
 *            where A(J, J) stores T(J, J) and A(J-1, (J+1):N) stores U(J, (J+1):N)
 *
-            IF( (J1+J-1).GT.1 ) THEN
+            IF( K.GT.1 ) THEN
                ALPHA = -A( K, J )
                CALL SAXPY( M-J, ALPHA, A( K-1, J+1 ), LDA,
      $                                 WORK( 2 ), 1 )
@@ -406,7 +406,7 @@
 *           Compute WORK(2:N) = T(J, J) L((J+1):N, J)
 *            where A(J, J) = T(J, J) and A((J+1):N, J-1) = L((J+1):N, J)
 *
-            IF( (J1+J-1).GT.1 ) THEN
+            IF( K.GT.1 ) THEN
                ALPHA = -A( J, K )
                CALL SAXPY( M-J, ALPHA, A( J+1, K-1 ), 1,
      $                                 WORK( 2 ), 1 )

diff --git a/SRC/zlahef_aa.f b/SRC/zlahef_aa.f
index 7a74f66..45d1b67 100644 (file)
--- a/SRC/zlahef_aa.f
+++ b/SRC/zlahef_aa.f
@@ -258,7 +258,7 @@
 *           Compute WORK(2:N) = T(J, J) L(J, (J+1):N)
 *            where A(J, J) stores T(J, J) and A(J-1, (J+1):N) stores U(J, (J+1):N)
 *
-            IF( (J1+J-1).GT.1 ) THEN
+            IF( K.GT.1 ) THEN
                ALPHA = -A( K, J )
                CALL ZAXPY( M-J, ALPHA, A( K-1, J+1 ), LDA,
      $                                 WORK( 2 ), 1 )
@@ -412,7 +412,7 @@
 *           Compute WORK(2:N) = T(J, J) L((J+1):N, J)
 *            where A(J, J) = T(J, J) and A((J+1):N, J-1) = L((J+1):N, J)
 *
-            IF( (J1+J-1).GT.1 ) THEN
+            IF( K.GT.1 ) THEN
                ALPHA = -A( J, K )
                CALL ZAXPY( M-J, ALPHA, A( J+1, K-1 ), 1,
      $                                 WORK( 2 ), 1 )