small tweak for the last column of the panel

iyamazaki · iyamazaki · commit 0aa9a6491719 · 2017-11-03T12:07:18.000-04:00
diff --git a/SRC/clahef_aa.f b/SRC/clahef_aa.f
@@ -166,7 +166,7 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
       PARAMETER    ( ZERO = (0.0E+0, 0.0E+0), ONE = (1.0E+0, 0.0E+0) )
 *
 *     .. Local Scalars ..
-      INTEGER      J, K, K1, I1, I2
+      INTEGER      J, K, K1, I1, I2, MJ
       COMPLEX      PIV, ALPHA
 *     ..
 *     .. External Functions ..
@@ -205,6 +205,14 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,
 *
          K = J1+J-1
+         IF( J.EQ.M ) THEN
+*
+*            Only need to compute T(J, J)
+*
+             MJ = 1
+         ELSE
+             MJ = M-J+1
+         END IF
 *
 *        H(J:N, J) := A(J, J:N) - H(J:N, 1:(J-1)) * L(J1:(J-1), J),
 *         where H(J:N, J) has been initialized to be A(J, J:N)
@@ -218,7 +226,7 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *           first column
 *
             CALL CLACGV( J-K1, A( 1, J ), 1 )
-            CALL CGEMV( 'No transpose', M-J+1, J-K1,
+            CALL CGEMV( 'No transpose', MJ, J-K1,
      $                 -ONE, H( J, K1 ), LDH,
      $                       A( 1, J ), 1,
      $                  ONE, H( J, J ), 1 )
@@ -227,15 +235,15 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *
 *        Copy H(i:n, i) into WORK
 *
-         CALL CCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )
+         CALL CCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )
 *
          IF( J.GT.K1 ) THEN
 *
 *           Compute WORK := WORK - L(J-1, J:N) * T(J-1,J),
 *            where A(J-1, J) stores T(J-1, J) and A(J-2, J:N) stores U(J-1, J:N)
 *
             ALPHA = -CONJG( A( K-1, J ) )
-            CALL CAXPY( M-J+1, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )
+            CALL CAXPY( MJ, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )
          END IF
 *
 *        Set A(J, J) = T(J, J)
@@ -349,6 +357,14 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,
 *
          K = J1+J-1
+         IF( J.EQ.M ) THEN
+*
+*            Only need to compute T(J, J)
+*
+             MJ = 1
+         ELSE
+             MJ = M-J+1
+         END IF
 *
 *        H(J:N, J) := A(J:N, J) - H(J:N, 1:(J-1)) * L(J, J1:(J-1))^T,
 *         where H(J:N, J) has been initialized to be A(J:N, J)
@@ -362,7 +378,7 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *           first column
 *
             CALL CLACGV( J-K1, A( J, 1 ), LDA )
-            CALL CGEMV( 'No transpose', M-J+1, J-K1,
+            CALL CGEMV( 'No transpose', MJ, J-K1,
      $                 -ONE, H( J, K1 ), LDH,
      $                       A( J, 1 ), LDA,
      $                  ONE, H( J, J ), 1 )
@@ -371,15 +387,15 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *
 *        Copy H(J:N, J) into WORK
 *
-         CALL CCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )
+         CALL CCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )
 *
          IF( J.GT.K1 ) THEN
 *
 *           Compute WORK := WORK - L(J:N, J-1) * T(J-1,J),
 *            where A(J-1, J) = T(J-1, J) and A(J, J-2) = L(J, J-1)
 *
             ALPHA = -CONJG( A( J, K-1 ) )
-            CALL CAXPY( M-J+1, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )
+            CALL CAXPY( MJ, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )
          END IF
 *
 *        Set A(J, J) = T(J, J)
diff --git a/SRC/clasyf_aa.f b/SRC/clasyf_aa.f
@@ -166,7 +166,7 @@ SUBROUTINE CLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
       PARAMETER          ( ZERO = 0.0E+0, ONE = 1.0E+0 )
 *
 *     .. Local Scalars ..
-      INTEGER            J, K, K1, I1, I2
+      INTEGER            J, K, K1, I1, I2, MJ
       COMPLEX            PIV, ALPHA
 *     ..
 *     .. External Functions ..
@@ -205,6 +205,14 @@ SUBROUTINE CLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,
 *
          K = J1+J-1
+         IF( J.EQ.M ) THEN
+*
+*            Only need to compute T(J, J)
+*
+             MJ = 1
+         ELSE
+             MJ = M-J+1
+         END IF
 *
 *        H(J:M, J) := A(J, J:M) - H(J:M, 1:(J-1)) * L(J1:(J-1), J),
 *         where H(J:M, J) has been initialized to be A(J, J:M)
@@ -217,23 +225,23 @@ SUBROUTINE CLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, K is J+1, skipping only the
 *           first column
 *
-            CALL CGEMV( 'No transpose', M-J+1, J-K1,
+            CALL CGEMV( 'No transpose', MJ, J-K1,
      $                 -ONE, H( J, K1 ), LDH,
      $                       A( 1, J ), 1,
      $                  ONE, H( J, J ), 1 )
          END IF
 *
 *        Copy H(i:M, i) into WORK
 *
-         CALL CCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )
+         CALL CCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )
 *
          IF( J.GT.K1 ) THEN
 *
 *           Compute WORK := WORK - L(J-1, J:M) * T(J-1,J),
 *            where A(J-1, J) stores T(J-1, J) and A(J-2, J:M) stores U(J-1, J:M)
 *
             ALPHA = -A( K-1, J )
-            CALL CAXPY( M-J+1, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )
+            CALL CAXPY( MJ, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )
          END IF
 *
 *        Set A(J, J) = T(J, J)
@@ -345,6 +353,14 @@ SUBROUTINE CLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,
 *
          K = J1+J-1
+         IF( J.EQ.M ) THEN
+*
+*            Only need to compute T(J, J)
+*
+             MJ = 1
+         ELSE
+             MJ = M-J+1
+         END IF
 *
 *        H(J:M, J) := A(J:M, J) - H(J:M, 1:(J-1)) * L(J, J1:(J-1))^T,
 *         where H(J:M, J) has been initialized to be A(J:M, J)
@@ -357,23 +373,23 @@ SUBROUTINE CLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, K is J+1, skipping only the
 *           first column
 *
-            CALL CGEMV( 'No transpose', M-J+1, J-K1,
+            CALL CGEMV( 'No transpose', MJ, J-K1,
      $                 -ONE, H( J, K1 ), LDH,
      $                       A( J, 1 ), LDA,
      $                  ONE, H( J, J ), 1 )
          END IF
 *
 *        Copy H(J:M, J) into WORK
 *
-         CALL CCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )
+         CALL CCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )
 *
          IF( J.GT.K1 ) THEN
 *
 *           Compute WORK := WORK - L(J:M, J-1) * T(J-1,J),
 *            where A(J-1, J) = T(J-1, J) and A(J, J-2) = L(J, J-1)
 *
             ALPHA = -A( J, K-1 )
-            CALL CAXPY( M-J+1, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )
+            CALL CAXPY( MJ, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )
          END IF
 *
 *        Set A(J, J) = T(J, J)
diff --git a/SRC/dlasyf_aa.f b/SRC/dlasyf_aa.f
@@ -166,7 +166,7 @@ SUBROUTINE DLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
       PARAMETER          ( ZERO = 0.0D+0, ONE = 1.0D+0 )
 *
 *     .. Local Scalars ..
-      INTEGER            J, K, K1, I1, I2
+      INTEGER            J, K, K1, I1, I2, MJ
       DOUBLE PRECISION   PIV, ALPHA
 *     ..
 *     .. External Functions ..
@@ -205,6 +205,14 @@ SUBROUTINE DLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,
 *
          K = J1+J-1
+         IF( J.EQ.M ) THEN
+*
+*            Only need to compute T(J, J)
+*
+             MJ = 1
+         ELSE
+             MJ = M-J+1
+         END IF
 *
 *        H(J:M, J) := A(J, J:M) - H(J:M, 1:(J-1)) * L(J1:(J-1), J),
 *         where H(J:M, J) has been initialized to be A(J, J:M)
@@ -217,23 +225,23 @@ SUBROUTINE DLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, K is J+1, skipping only the
 *           first column
 *
-            CALL DGEMV( 'No transpose', M-J+1, J-K1,
+            CALL DGEMV( 'No transpose', MJ, J-K1,
      $                 -ONE, H( J, K1 ), LDH,
      $                       A( 1, J ), 1,
      $                  ONE, H( J, J ), 1 )
          END IF
 *
 *        Copy H(i:M, i) into WORK
 *
-         CALL DCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )
+         CALL DCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )
 *
          IF( J.GT.K1 ) THEN
 *
 *           Compute WORK := WORK - L(J-1, J:M) * T(J-1,J),
 *            where A(J-1, J) stores T(J-1, J) and A(J-2, J:M) stores U(J-1, J:M)
 *
             ALPHA = -A( K-1, J )
-            CALL DAXPY( M-J+1, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )
+            CALL DAXPY( MJ, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )
          END IF
 *
 *        Set A(J, J) = T(J, J)
@@ -345,6 +353,14 @@ SUBROUTINE DLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,
 *
          K = J1+J-1
+         IF( J.EQ.M ) THEN
+*
+*            Only need to compute T(J, J)
+*
+             MJ = 1
+         ELSE
+             MJ = M-J+1
+         END IF
 *
 *        H(J:M, J) := A(J:M, J) - H(J:M, 1:(J-1)) * L(J, J1:(J-1))^T,
 *         where H(J:M, J) has been initialized to be A(J:M, J)
@@ -357,23 +373,23 @@ SUBROUTINE DLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, K is J+1, skipping only the
 *           first column
 *
-            CALL DGEMV( 'No transpose', M-J+1, J-K1,
+            CALL DGEMV( 'No transpose', MJ, J-K1,
      $                 -ONE, H( J, K1 ), LDH,
      $                       A( J, 1 ), LDA,
      $                  ONE, H( J, J ), 1 )
          END IF
 *
 *        Copy H(J:M, J) into WORK
 *
-         CALL DCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )
+         CALL DCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )
 *
          IF( J.GT.K1 ) THEN
 *
 *           Compute WORK := WORK - L(J:M, J-1) * T(J-1,J),
 *            where A(J-1, J) = T(J-1, J) and A(J, J-2) = L(J, J-1)
 *
             ALPHA = -A( J, K-1 )
-            CALL DAXPY( M-J+1, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )
+            CALL DAXPY( MJ, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )
          END IF
 *
 *        Set A(J, J) = T(J, J)
diff --git a/SRC/slasyf_aa.f b/SRC/slasyf_aa.f
@@ -166,7 +166,7 @@ SUBROUTINE SLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
       PARAMETER          ( ZERO = 0.0E+0, ONE = 1.0E+0 )
 *
 *     .. Local Scalars ..
-      INTEGER            J, K, K1, I1, I2
+      INTEGER            J, K, K1, I1, I2, MJ
       REAL               PIV, ALPHA
 *     ..
 *     .. External Functions ..
@@ -205,6 +205,14 @@ SUBROUTINE SLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,
 *
          K = J1+J-1
+         IF( J.EQ.M ) THEN
+*
+*            Only need to compute T(J, J)
+*
+             MJ = 1
+         ELSE
+             MJ = M-J+1
+         END IF
 *
 *        H(J:M, J) := A(J, J:M) - H(J:M, 1:(J-1)) * L(J1:(J-1), J),
 *         where H(J:M, J) has been initialized to be A(J, J:M)
@@ -217,23 +225,23 @@ SUBROUTINE SLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, K is J+1, skipping only the
 *           first column
 *
-            CALL SGEMV( 'No transpose', M-J+1, J-K1,
+            CALL SGEMV( 'No transpose', MJ, J-K1,
      $                 -ONE, H( J, K1 ), LDH,
      $                       A( 1, J ), 1,
      $                  ONE, H( J, J ), 1 )
          END IF
 *
 *        Copy H(i:M, i) into WORK
 *
-         CALL SCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )
+         CALL SCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )
 *
          IF( J.GT.K1 ) THEN
 *
 *           Compute WORK := WORK - L(J-1, J:M) * T(J-1,J),
 *            where A(J-1, J) stores T(J-1, J) and A(J-2, J:M) stores U(J-1, J:M)
 *
             ALPHA = -A( K-1, J )
-            CALL SAXPY( M-J+1, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )
+            CALL SAXPY( MJ, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )
          END IF
 *
 *        Set A(J, J) = T(J, J)
@@ -345,6 +353,14 @@ SUBROUTINE SLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,
 *
          K = J1+J-1
+         IF( J.EQ.M ) THEN
+*
+*            Only need to compute T(J, J)
+*
+             MJ = 1
+         ELSE
+             MJ = M-J+1
+         END IF
 *
 *        H(J:M, J) := A(J:M, J) - H(J:M, 1:(J-1)) * L(J, J1:(J-1))^T,
 *         where H(J:M, J) has been initialized to be A(J:M, J)
@@ -357,23 +373,23 @@ SUBROUTINE SLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,
 *         > for the rest of the columns, K is J+1, skipping only the
 *           first column
 *
-            CALL SGEMV( 'No transpose', M-J+1, J-K1,
+            CALL SGEMV( 'No transpose', MJ, J-K1,
      $                 -ONE, H( J, K1 ), LDH,
      $                       A( J, 1 ), LDA,
      $                  ONE, H( J, J ), 1 )
          END IF
 *
 *        Copy H(J:M, J) into WORK
 *
-         CALL SCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )
+         CALL SCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )
 *
          IF( J.GT.K1 ) THEN
 *
 *           Compute WORK := WORK - L(J:M, J-1) * T(J-1,J),
 *            where A(J-1, J) = T(J-1, J) and A(J, J-2) = L(J, J-1)
 *
             ALPHA = -A( J, K-1 )
-            CALL SAXPY( M-J+1, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )
+            CALL SAXPY( MJ, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )
          END IF
 *
 *        Set A(J, J) = T(J, J)
diff --git a/SRC/zlahef_aa.f b/SRC/zlahef_aa.f
diff --git a/SRC/zlasyf_aa.f b/SRC/zlasyf_aa.f

Original file line number	Diff line number	Diff line change
`@@ -166,7 +166,7 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`166`	`166`	`PARAMETER ( ZERO = (0.0E+0, 0.0E+0), ONE = (1.0E+0, 0.0E+0) )`
`167`	`167`	`*`
`168`	`168`	`* .. Local Scalars ..`
`169`		`- INTEGER J, K, K1, I1, I2`
	`169`	`+ INTEGER J, K, K1, I1, I2, MJ`
`170`	`170`	`COMPLEX PIV, ALPHA`
`171`	`171`	`* ..`
`172`	`172`	`* .. External Functions ..`
`@@ -205,6 +205,14 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`205`	`205`	`* > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,`
`206`	`206`	`*`
`207`	`207`	`K = J1+J-1`
	`208`	`+ IF( J.EQ.M ) THEN`
	`209`	`+*`
	`210`	`+* Only need to compute T(J, J)`
	`211`	`+*`
	`212`	`+ MJ = 1`
	`213`	`+ ELSE`
	`214`	`+ MJ = M-J+1`
	`215`	`+ END IF`
`208`	`216`	`*`
`209`	`217`	`* H(J:N, J) := A(J, J:N) - H(J:N, 1:(J-1)) * L(J1:(J-1), J),`
`210`	`218`	`* where H(J:N, J) has been initialized to be A(J, J:N)`
`@@ -218,7 +226,7 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`218`	`226`	`* first column`
`219`	`227`	`*`
`220`	`228`	`CALL CLACGV( J-K1, A( 1, J ), 1 )`
`221`		`- CALL CGEMV( 'No transpose', M-J+1, J-K1,`
	`229`	`+ CALL CGEMV( 'No transpose', MJ, J-K1,`
`222`	`230`	`$ -ONE, H( J, K1 ), LDH,`
`223`	`231`	`$ A( 1, J ), 1,`
`224`	`232`	`$ ONE, H( J, J ), 1 )`
`@@ -227,15 +235,15 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`227`	`235`	`*`
`228`	`236`	`* Copy H(i:n, i) into WORK`
`229`	`237`	`*`
`230`		`- CALL CCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )`
	`238`	`+ CALL CCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )`
`231`	`239`	`*`
`232`	`240`	`IF( J.GT.K1 ) THEN`
`233`	`241`	`*`
`234`	`242`	`* Compute WORK := WORK - L(J-1, J:N) * T(J-1,J),`
`235`	`243`	`* where A(J-1, J) stores T(J-1, J) and A(J-2, J:N) stores U(J-1, J:N)`
`236`	`244`	`*`
`237`	`245`	`ALPHA = -CONJG( A( K-1, J ) )`
`238`		`- CALL CAXPY( M-J+1, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )`
	`246`	`+ CALL CAXPY( MJ, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )`
`239`	`247`	`END IF`
`240`	`248`	`*`
`241`	`249`	`* Set A(J, J) = T(J, J)`
`@@ -349,6 +357,14 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`349`	`357`	`* > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,`
`350`	`358`	`*`
`351`	`359`	`K = J1+J-1`
	`360`	`+ IF( J.EQ.M ) THEN`
	`361`	`+*`
	`362`	`+* Only need to compute T(J, J)`
	`363`	`+*`
	`364`	`+ MJ = 1`
	`365`	`+ ELSE`
	`366`	`+ MJ = M-J+1`
	`367`	`+ END IF`
`352`	`368`	`*`
`353`	`369`	`* H(J:N, J) := A(J:N, J) - H(J:N, 1:(J-1)) * L(J, J1:(J-1))^T,`
`354`	`370`	`* where H(J:N, J) has been initialized to be A(J:N, J)`
`@@ -362,7 +378,7 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`362`	`378`	`* first column`
`363`	`379`	`*`
`364`	`380`	`CALL CLACGV( J-K1, A( J, 1 ), LDA )`
`365`		`- CALL CGEMV( 'No transpose', M-J+1, J-K1,`
	`381`	`+ CALL CGEMV( 'No transpose', MJ, J-K1,`
`366`	`382`	`$ -ONE, H( J, K1 ), LDH,`
`367`	`383`	`$ A( J, 1 ), LDA,`
`368`	`384`	`$ ONE, H( J, J ), 1 )`
`@@ -371,15 +387,15 @@ SUBROUTINE CLAHEF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`371`	`387`	`*`
`372`	`388`	`* Copy H(J:N, J) into WORK`
`373`	`389`	`*`
`374`		`- CALL CCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )`
	`390`	`+ CALL CCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )`
`375`	`391`	`*`
`376`	`392`	`IF( J.GT.K1 ) THEN`
`377`	`393`	`*`
`378`	`394`	`* Compute WORK := WORK - L(J:N, J-1) * T(J-1,J),`
`379`	`395`	`* where A(J-1, J) = T(J-1, J) and A(J, J-2) = L(J, J-1)`
`380`	`396`	`*`
`381`	`397`	`ALPHA = -CONJG( A( J, K-1 ) )`
`382`		`- CALL CAXPY( M-J+1, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )`
	`398`	`+ CALL CAXPY( MJ, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )`
`383`	`399`	`END IF`
`384`	`400`	`*`
`385`	`401`	`* Set A(J, J) = T(J, J)`
Original file line number	Diff line number	Diff line change
`@@ -166,7 +166,7 @@ SUBROUTINE CLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`166`	`166`	`PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 )`
`167`	`167`	`*`
`168`	`168`	`* .. Local Scalars ..`
`169`		`- INTEGER J, K, K1, I1, I2`
	`169`	`+ INTEGER J, K, K1, I1, I2, MJ`
`170`	`170`	`COMPLEX PIV, ALPHA`
`171`	`171`	`* ..`
`172`	`172`	`* .. External Functions ..`
`@@ -205,6 +205,14 @@ SUBROUTINE CLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`205`	`205`	`* > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,`
`206`	`206`	`*`
`207`	`207`	`K = J1+J-1`
	`208`	`+ IF( J.EQ.M ) THEN`
	`209`	`+*`
	`210`	`+* Only need to compute T(J, J)`
	`211`	`+*`
	`212`	`+ MJ = 1`
	`213`	`+ ELSE`
	`214`	`+ MJ = M-J+1`
	`215`	`+ END IF`
`208`	`216`	`*`
`209`	`217`	`* H(J:M, J) := A(J, J:M) - H(J:M, 1:(J-1)) * L(J1:(J-1), J),`
`210`	`218`	`* where H(J:M, J) has been initialized to be A(J, J:M)`
`@@ -217,23 +225,23 @@ SUBROUTINE CLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`217`	`225`	`* > for the rest of the columns, K is J+1, skipping only the`
`218`	`226`	`* first column`
`219`	`227`	`*`
`220`		`- CALL CGEMV( 'No transpose', M-J+1, J-K1,`
	`228`	`+ CALL CGEMV( 'No transpose', MJ, J-K1,`
`221`	`229`	`$ -ONE, H( J, K1 ), LDH,`
`222`	`230`	`$ A( 1, J ), 1,`
`223`	`231`	`$ ONE, H( J, J ), 1 )`
`224`	`232`	`END IF`
`225`	`233`	`*`
`226`	`234`	`* Copy H(i:M, i) into WORK`
`227`	`235`	`*`
`228`		`- CALL CCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )`
	`236`	`+ CALL CCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )`
`229`	`237`	`*`
`230`	`238`	`IF( J.GT.K1 ) THEN`
`231`	`239`	`*`
`232`	`240`	`* Compute WORK := WORK - L(J-1, J:M) * T(J-1,J),`
`233`	`241`	`* where A(J-1, J) stores T(J-1, J) and A(J-2, J:M) stores U(J-1, J:M)`
`234`	`242`	`*`
`235`	`243`	`ALPHA = -A( K-1, J )`
`236`		`- CALL CAXPY( M-J+1, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )`
	`244`	`+ CALL CAXPY( MJ, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )`
`237`	`245`	`END IF`
`238`	`246`	`*`
`239`	`247`	`* Set A(J, J) = T(J, J)`
`@@ -345,6 +353,14 @@ SUBROUTINE CLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`345`	`353`	`* > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,`
`346`	`354`	`*`
`347`	`355`	`K = J1+J-1`
	`356`	`+ IF( J.EQ.M ) THEN`
	`357`	`+*`
	`358`	`+* Only need to compute T(J, J)`
	`359`	`+*`
	`360`	`+ MJ = 1`
	`361`	`+ ELSE`
	`362`	`+ MJ = M-J+1`
	`363`	`+ END IF`
`348`	`364`	`*`
`349`	`365`	`* H(J:M, J) := A(J:M, J) - H(J:M, 1:(J-1)) * L(J, J1:(J-1))^T,`
`350`	`366`	`* where H(J:M, J) has been initialized to be A(J:M, J)`
`@@ -357,23 +373,23 @@ SUBROUTINE CLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`357`	`373`	`* > for the rest of the columns, K is J+1, skipping only the`
`358`	`374`	`* first column`
`359`	`375`	`*`
`360`		`- CALL CGEMV( 'No transpose', M-J+1, J-K1,`
	`376`	`+ CALL CGEMV( 'No transpose', MJ, J-K1,`
`361`	`377`	`$ -ONE, H( J, K1 ), LDH,`
`362`	`378`	`$ A( J, 1 ), LDA,`
`363`	`379`	`$ ONE, H( J, J ), 1 )`
`364`	`380`	`END IF`
`365`	`381`	`*`
`366`	`382`	`* Copy H(J:M, J) into WORK`
`367`	`383`	`*`
`368`		`- CALL CCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )`
	`384`	`+ CALL CCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )`
`369`	`385`	`*`
`370`	`386`	`IF( J.GT.K1 ) THEN`
`371`	`387`	`*`
`372`	`388`	`* Compute WORK := WORK - L(J:M, J-1) * T(J-1,J),`
`373`	`389`	`* where A(J-1, J) = T(J-1, J) and A(J, J-2) = L(J, J-1)`
`374`	`390`	`*`
`375`	`391`	`ALPHA = -A( J, K-1 )`
`376`		`- CALL CAXPY( M-J+1, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )`
	`392`	`+ CALL CAXPY( MJ, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )`
`377`	`393`	`END IF`
`378`	`394`	`*`
`379`	`395`	`* Set A(J, J) = T(J, J)`
Original file line number	Diff line number	Diff line change
`@@ -166,7 +166,7 @@ SUBROUTINE DLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`166`	`166`	`PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )`
`167`	`167`	`*`
`168`	`168`	`* .. Local Scalars ..`
`169`		`- INTEGER J, K, K1, I1, I2`
	`169`	`+ INTEGER J, K, K1, I1, I2, MJ`
`170`	`170`	`DOUBLE PRECISION PIV, ALPHA`
`171`	`171`	`* ..`
`172`	`172`	`* .. External Functions ..`
`@@ -205,6 +205,14 @@ SUBROUTINE DLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`205`	`205`	`* > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,`
`206`	`206`	`*`
`207`	`207`	`K = J1+J-1`
	`208`	`+ IF( J.EQ.M ) THEN`
	`209`	`+*`
	`210`	`+* Only need to compute T(J, J)`
	`211`	`+*`
	`212`	`+ MJ = 1`
	`213`	`+ ELSE`
	`214`	`+ MJ = M-J+1`
	`215`	`+ END IF`
`208`	`216`	`*`
`209`	`217`	`* H(J:M, J) := A(J, J:M) - H(J:M, 1:(J-1)) * L(J1:(J-1), J),`
`210`	`218`	`* where H(J:M, J) has been initialized to be A(J, J:M)`
`@@ -217,23 +225,23 @@ SUBROUTINE DLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`217`	`225`	`* > for the rest of the columns, K is J+1, skipping only the`
`218`	`226`	`* first column`
`219`	`227`	`*`
`220`		`- CALL DGEMV( 'No transpose', M-J+1, J-K1,`
	`228`	`+ CALL DGEMV( 'No transpose', MJ, J-K1,`
`221`	`229`	`$ -ONE, H( J, K1 ), LDH,`
`222`	`230`	`$ A( 1, J ), 1,`
`223`	`231`	`$ ONE, H( J, J ), 1 )`
`224`	`232`	`END IF`
`225`	`233`	`*`
`226`	`234`	`* Copy H(i:M, i) into WORK`
`227`	`235`	`*`
`228`		`- CALL DCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )`
	`236`	`+ CALL DCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )`
`229`	`237`	`*`
`230`	`238`	`IF( J.GT.K1 ) THEN`
`231`	`239`	`*`
`232`	`240`	`* Compute WORK := WORK - L(J-1, J:M) * T(J-1,J),`
`233`	`241`	`* where A(J-1, J) stores T(J-1, J) and A(J-2, J:M) stores U(J-1, J:M)`
`234`	`242`	`*`
`235`	`243`	`ALPHA = -A( K-1, J )`
`236`		`- CALL DAXPY( M-J+1, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )`
	`244`	`+ CALL DAXPY( MJ, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )`
`237`	`245`	`END IF`
`238`	`246`	`*`
`239`	`247`	`* Set A(J, J) = T(J, J)`
`@@ -345,6 +353,14 @@ SUBROUTINE DLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`345`	`353`	`* > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,`
`346`	`354`	`*`
`347`	`355`	`K = J1+J-1`
	`356`	`+ IF( J.EQ.M ) THEN`
	`357`	`+*`
	`358`	`+* Only need to compute T(J, J)`
	`359`	`+*`
	`360`	`+ MJ = 1`
	`361`	`+ ELSE`
	`362`	`+ MJ = M-J+1`
	`363`	`+ END IF`
`348`	`364`	`*`
`349`	`365`	`* H(J:M, J) := A(J:M, J) - H(J:M, 1:(J-1)) * L(J, J1:(J-1))^T,`
`350`	`366`	`* where H(J:M, J) has been initialized to be A(J:M, J)`
`@@ -357,23 +373,23 @@ SUBROUTINE DLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`357`	`373`	`* > for the rest of the columns, K is J+1, skipping only the`
`358`	`374`	`* first column`
`359`	`375`	`*`
`360`		`- CALL DGEMV( 'No transpose', M-J+1, J-K1,`
	`376`	`+ CALL DGEMV( 'No transpose', MJ, J-K1,`
`361`	`377`	`$ -ONE, H( J, K1 ), LDH,`
`362`	`378`	`$ A( J, 1 ), LDA,`
`363`	`379`	`$ ONE, H( J, J ), 1 )`
`364`	`380`	`END IF`
`365`	`381`	`*`
`366`	`382`	`* Copy H(J:M, J) into WORK`
`367`	`383`	`*`
`368`		`- CALL DCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )`
	`384`	`+ CALL DCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )`
`369`	`385`	`*`
`370`	`386`	`IF( J.GT.K1 ) THEN`
`371`	`387`	`*`
`372`	`388`	`* Compute WORK := WORK - L(J:M, J-1) * T(J-1,J),`
`373`	`389`	`* where A(J-1, J) = T(J-1, J) and A(J, J-2) = L(J, J-1)`
`374`	`390`	`*`
`375`	`391`	`ALPHA = -A( J, K-1 )`
`376`		`- CALL DAXPY( M-J+1, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )`
	`392`	`+ CALL DAXPY( MJ, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )`
`377`	`393`	`END IF`
`378`	`394`	`*`
`379`	`395`	`* Set A(J, J) = T(J, J)`
Original file line number	Diff line number	Diff line change
`@@ -166,7 +166,7 @@ SUBROUTINE SLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`166`	`166`	`PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 )`
`167`	`167`	`*`
`168`	`168`	`* .. Local Scalars ..`
`169`		`- INTEGER J, K, K1, I1, I2`
	`169`	`+ INTEGER J, K, K1, I1, I2, MJ`
`170`	`170`	`REAL PIV, ALPHA`
`171`	`171`	`* ..`
`172`	`172`	`* .. External Functions ..`
`@@ -205,6 +205,14 @@ SUBROUTINE SLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`205`	`205`	`* > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,`
`206`	`206`	`*`
`207`	`207`	`K = J1+J-1`
	`208`	`+ IF( J.EQ.M ) THEN`
	`209`	`+*`
	`210`	`+* Only need to compute T(J, J)`
	`211`	`+*`
	`212`	`+ MJ = 1`
	`213`	`+ ELSE`
	`214`	`+ MJ = M-J+1`
	`215`	`+ END IF`
`208`	`216`	`*`
`209`	`217`	`* H(J:M, J) := A(J, J:M) - H(J:M, 1:(J-1)) * L(J1:(J-1), J),`
`210`	`218`	`* where H(J:M, J) has been initialized to be A(J, J:M)`
`@@ -217,23 +225,23 @@ SUBROUTINE SLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`217`	`225`	`* > for the rest of the columns, K is J+1, skipping only the`
`218`	`226`	`* first column`
`219`	`227`	`*`
`220`		`- CALL SGEMV( 'No transpose', M-J+1, J-K1,`
	`228`	`+ CALL SGEMV( 'No transpose', MJ, J-K1,`
`221`	`229`	`$ -ONE, H( J, K1 ), LDH,`
`222`	`230`	`$ A( 1, J ), 1,`
`223`	`231`	`$ ONE, H( J, J ), 1 )`
`224`	`232`	`END IF`
`225`	`233`	`*`
`226`	`234`	`* Copy H(i:M, i) into WORK`
`227`	`235`	`*`
`228`		`- CALL SCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )`
	`236`	`+ CALL SCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )`
`229`	`237`	`*`
`230`	`238`	`IF( J.GT.K1 ) THEN`
`231`	`239`	`*`
`232`	`240`	`* Compute WORK := WORK - L(J-1, J:M) * T(J-1,J),`
`233`	`241`	`* where A(J-1, J) stores T(J-1, J) and A(J-2, J:M) stores U(J-1, J:M)`
`234`	`242`	`*`
`235`	`243`	`ALPHA = -A( K-1, J )`
`236`		`- CALL SAXPY( M-J+1, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )`
	`244`	`+ CALL SAXPY( MJ, ALPHA, A( K-2, J ), LDA, WORK( 1 ), 1 )`
`237`	`245`	`END IF`
`238`	`246`	`*`
`239`	`247`	`* Set A(J, J) = T(J, J)`
`@@ -345,6 +353,14 @@ SUBROUTINE SLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`345`	`353`	`* > for the rest of the columns, J1 is 2, and J1+J-1 is J+1,`
`346`	`354`	`*`
`347`	`355`	`K = J1+J-1`
	`356`	`+ IF( J.EQ.M ) THEN`
	`357`	`+*`
	`358`	`+* Only need to compute T(J, J)`
	`359`	`+*`
	`360`	`+ MJ = 1`
	`361`	`+ ELSE`
	`362`	`+ MJ = M-J+1`
	`363`	`+ END IF`
`348`	`364`	`*`
`349`	`365`	`* H(J:M, J) := A(J:M, J) - H(J:M, 1:(J-1)) * L(J, J1:(J-1))^T,`
`350`	`366`	`* where H(J:M, J) has been initialized to be A(J:M, J)`
`@@ -357,23 +373,23 @@ SUBROUTINE SLASYF_AA( UPLO, J1, M, NB, A, LDA, IPIV,`
`357`	`373`	`* > for the rest of the columns, K is J+1, skipping only the`
`358`	`374`	`* first column`
`359`	`375`	`*`
`360`		`- CALL SGEMV( 'No transpose', M-J+1, J-K1,`
	`376`	`+ CALL SGEMV( 'No transpose', MJ, J-K1,`
`361`	`377`	`$ -ONE, H( J, K1 ), LDH,`
`362`	`378`	`$ A( J, 1 ), LDA,`
`363`	`379`	`$ ONE, H( J, J ), 1 )`
`364`	`380`	`END IF`
`365`	`381`	`*`
`366`	`382`	`* Copy H(J:M, J) into WORK`
`367`	`383`	`*`
`368`		`- CALL SCOPY( M-J+1, H( J, J ), 1, WORK( 1 ), 1 )`
	`384`	`+ CALL SCOPY( MJ, H( J, J ), 1, WORK( 1 ), 1 )`
`369`	`385`	`*`
`370`	`386`	`IF( J.GT.K1 ) THEN`
`371`	`387`	`*`
`372`	`388`	`* Compute WORK := WORK - L(J:M, J-1) * T(J-1,J),`
`373`	`389`	`* where A(J-1, J) = T(J-1, J) and A(J, J-2) = L(J, J-1)`
`374`	`390`	`*`
`375`	`391`	`ALPHA = -A( J, K-1 )`
`376`		`- CALL SAXPY( M-J+1, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )`
	`392`	`+ CALL SAXPY( MJ, ALPHA, A( J, K-2 ), 1, WORK( 1 ), 1 )`
`377`	`393`	`END IF`
`378`	`394`	`*`
`379`	`395`	`* Set A(J, J) = T(J, J)`