Changed comment in (sdcz)lasyf.f and (cz)lahef.f about interchanges for 2x2 pivot
authorigor175 <igor175@8a072113-8704-0410-8d35-dd094bca7971>
Fri, 2 Nov 2012 03:02:56 +0000 (03:02 +0000)
committerigor175 <igor175@8a072113-8704-0410-8d35-dd094bca7971>
Fri, 2 Nov 2012 03:02:56 +0000 (03:02 +0000)
SRC/clahef.f
SRC/clasyf.f
SRC/dlasyf.f
SRC/slasyf.f
SRC/zlahef.f
SRC/zlasyf.f

index 004ce603e95c4923cf5bbec8636c090bb6374e9b..f6c84a85638457a2478e266440793486748dc1e5 100644 (file)
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K-1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K-1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = REAL( A( KK, KK ) )
                CALL CCOPY( KK-1-KP, A( KP+1, KK ), 1, A( KP, KP+1 ),
      $            CALL CCOPY( KP-1, A( 1, KK ), 1, A( 1, KP ), 1 )
 *
 *              Interchange rows KK and KP in last K+1 to N columns of A
-*              (columns K (or K and K-1) of A will be later overwritten)
-*              and interchange rows KK and KP in last KKW to NB columns
-*              of W
+*              (columns K (or K and K-1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in last KKW to NB columns of W.
 *
                IF( K.LT.N )
      $            CALL CSWAP( N-K, A( KK, K+1 ), LDA, A( KP, K+1 ),
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K+1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K+1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = REAL( A( KK, KK ) )
                CALL CCOPY( KP-KK-1, A( KK+1, KK ), 1, A( KP, KK+1 ),
      $            CALL CCOPY( N-KP, A( KP+1, KK ), 1, A( KP+1, KP ), 1 )
 *
 *              Interchange rows KK and KP in first K-1 columns of A
-*              (columns K (or K and K+1) of A will be later overwritten)
-*              and interchange rows KK and KP in first KK columns of W
+*              (columns K (or K and K+1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in first KK columns of W.
 *
                IF( K.GT.1 )
      $            CALL CSWAP( K-1, A( KK, 1 ), LDA, A( KP, 1 ), LDA )
index 5c4e2bc4969d6d1d7e22b87f6dac31e035fa6614..e5312f190210b65cdb015b896a9c94b2249fcf0e 100644 (file)
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K-1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K-1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = A( KK, KK )
                CALL CCOPY( KK-1-KP, A( KP+1, KK ), 1, A( KP, KP+1 ),
      $            CALL CCOPY( KP-1, A( 1, KK ), 1, A( 1, KP ), 1 )
 *
 *              Interchange rows KK and KP in last K+1 to N columns of A
-*              (columns K (or K and K-1) of A will be later overwritten)
-*              and interchange rows KK and KP in last KKW to NB columns
-*              of W
+*              (columns K (or K and K-1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in last KKW to NB columns of W.
 *
                IF( K.LT.N )
      $            CALL CSWAP( N-K, A( KK, K+1 ), LDA, A( KP, K+1 ),
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K+1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K+1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = A( KK, KK )
                CALL CCOPY( KP-KK-1, A( KK+1, KK ), 1, A( KP, KK+1 ),
      $            CALL CCOPY( N-KP, A( KP+1, KK ), 1, A( KP+1, KP ), 1 )
 *
 *              Interchange rows KK and KP in first K-1 columns of A
-*              (columns K (or K and K+1) of A will be later overwritten)
-*              and interchange rows KK and KP in first KK columns of W
+*              (columns K (or K and K+1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in first KK columns of W.
 *
                IF( K.GT.1 )
      $            CALL CSWAP( K-1, A( KK, 1 ), LDA, A( KP, 1 ), LDA )
index c02aa1d243708820cfffe1139f3f8ce59ba61b55..57cc0fe60be849dd82e7b93b2316e06fd87643bd 100644 (file)
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K-1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K-1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = A( KK, KK )
                CALL DCOPY( KK-1-KP, A( KP+1, KK ), 1, A( KP, KP+1 ),
      $            CALL DCOPY( KP-1, A( 1, KK ), 1, A( 1, KP ), 1 )
 *
 *              Interchange rows KK and KP in last K+1 to N columns of A
-*              (columns K (or K and K-1) of A will be later overwritten)
-*              and interchange rows KK and KP in last KKW to NB columns
-*              of W
+*              (columns K (or K and K-1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in last KKW to NB columns of W.
 *
                IF( K.LT.N )
      $            CALL DSWAP( N-K, A( KK, K+1 ), LDA, A( KP, K+1 ),
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K+1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K+1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = A( KK, KK )
                CALL DCOPY( KP-KK-1, A( KK+1, KK ), 1, A( KP, KK+1 ),
      $            CALL DCOPY( N-KP, A( KP+1, KK ), 1, A( KP+1, KP ), 1 )
 *
 *              Interchange rows KK and KP in first K-1 columns of A
-*              (columns K (or K and K+1) of A will be later overwritten)
-*              and interchange rows KK and KP in first KK columns of W
+*              (columns K (or K and K+1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in first KK columns of W.
 *
                IF( K.GT.1 )
      $            CALL DSWAP( K-1, A( KK, 1 ), LDA, A( KP, 1 ), LDA )
index 7c0840f7a849bf4944e496cc31ac9286964cac80..9a96fe430af99b9719ecfb580b0697c0558e5d8e 100644 (file)
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K-1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K-1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = A( KK, KK )
                CALL SCOPY( KK-1-KP, A( KP+1, KK ), 1, A( KP, KP+1 ),
      $            CALL SCOPY( KP-1, A( 1, KK ), 1, A( 1, KP ), 1 )
 *
 *              Interchange rows KK and KP in last K+1 to N columns of A
-*              (columns K (or K and K-1) of A will be later overwritten)
-*              and interchange rows KK and KP in last KKW to NB columns
-*              of W
+*              (columns K (or K and K-1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in last KKW to NB columns of W.
 *
                IF( K.LT.N )
      $            CALL SSWAP( N-K, A( KK, K+1 ), LDA, A( KP, K+1 ),
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K+1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K+1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = A( KK, KK )
                CALL SCOPY( KP-KK-1, A( KK+1, KK ), 1, A( KP, KK+1 ),
      $            CALL SCOPY( N-KP, A( KP+1, KK ), 1, A( KP+1, KP ), 1 )
 *
 *              Interchange rows KK and KP in first K-1 columns of A
-*              (columns K (or K and K+1) of A will be later overwritten)
-*              and interchange rows KK and KP in first KK columns of W
+*              (columns K (or K and K+1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in first KK columns of W.
 *
                IF( K.GT.1 )
      $            CALL SSWAP( K-1, A( KK, 1 ), LDA, A( KP, 1 ), LDA )
index 2566c7d380676eb995a02d2b23c0e05c9c9de4ff..c9623bd75d0c1626a14447d8a6f143506af6a327 100644 (file)
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K-1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K-1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = DBLE( A( KK, KK ) )
                CALL ZCOPY( KK-1-KP, A( KP+1, KK ), 1, A( KP, KP+1 ),
      $            CALL ZCOPY( KP-1, A( 1, KK ), 1, A( 1, KP ), 1 )
 *
 *              Interchange rows KK and KP in last K+1 to N columns of A
-*              (columns K (or K and K-1) of A will be later overwritten)
-*              and interchange rows KK and KP in last KKW to NB columns
-*              of W
+*              (columns K (or K and K-1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in last KKW to NB columns of W.
 *
                IF( K.LT.N )
      $            CALL ZSWAP( N-K, A( KK, K+1 ), LDA, A( KP, K+1 ),
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K+1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K+1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = DBLE( A( KK, KK ) )
                CALL ZCOPY( KP-KK-1, A( KK+1, KK ), 1, A( KP, KK+1 ),
      $            CALL ZCOPY( N-KP, A( KP+1, KK ), 1, A( KP+1, KP ), 1 )
 *
 *              Interchange rows KK and KP in first K-1 columns of A
-*              (columns K (or K and K+1) of A will be later overwritten)
-*              and interchange rows KK and KP in first KK columns of W
+*              (columns K (or K and K+1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in first KK columns of W.
 *
                IF( K.GT.1 )
      $            CALL ZSWAP( K-1, A( KK, 1 ), LDA, A( KP, 1 ), LDA )
index 3249e649c771440715e7f40c5363c2dd2845f318..f31d152fae1da2cf55e807517d7e490422c34414 100644 (file)
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K-1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K-1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = A( KK, KK )
                CALL ZCOPY( KK-1-KP, A( KP+1, KK ), 1, A( KP, KP+1 ),
      $            CALL ZCOPY( KP-1, A( 1, KK ), 1, A( 1, KP ), 1 )
 *
 *              Interchange rows KK and KP in last K+1 to N columns of A
-*              (columns K (or K and K-1) of A will be later overwritten)
-*              and interchange rows KK and KP in last KKW to NB columns
-*              of W
+*              (columns K (or K and K-1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in last KKW to NB columns of W.
 *
                IF( K.LT.N )
      $            CALL ZSWAP( N-K, A( KK, K+1 ), LDA, A( KP, K+1 ),
 *
 *              Copy non-updated column KK to column KP of submatrix A
 *              at step K. No need to copy element into column K
-*              (or K and K+1) of A, since these columns will be later
-*              overwritten.
+*              (or K and K+1 for 2x2 pivot) of A, since these columns
+*              will be later overwritten.
 *
                A( KP, KP ) = A( KK, KK )
                CALL ZCOPY( KP-KK-1, A( KK+1, KK ), 1, A( KP, KK+1 ),
      $            CALL ZCOPY( N-KP, A( KP+1, KK ), 1, A( KP+1, KP ), 1 )
 *
 *              Interchange rows KK and KP in first K-1 columns of A
-*              (columns K (or K and K+1) of A will be later overwritten)
-*              and interchange rows KK and KP in first KK columns of W
+*              (columns K (or K and K+1 for 2x2 pivot) of A will be
+*              later overwritten). Interchange rows KK and KP
+*              in first KK columns of W.
 *
                IF( K.GT.1 )
      $            CALL ZSWAP( K-1, A( KK, 1 ), LDA, A( KP, 1 ), LDA )