diff options
author | igor175 <igor175@8a072113-8704-0410-8d35-dd094bca7971> | 2011-12-01 07:22:37 +0000 |
---|---|---|
committer | igor175 <igor175@8a072113-8704-0410-8d35-dd094bca7971> | 2011-12-01 07:22:37 +0000 |
commit | 4bb091197931c4b836fd6eec84b958c2f24a4337 (patch) | |
tree | c60da79f8629d1178f7f404fef6dc4e1a4aa31c0 /SRC | |
parent | bcc9b0d09e5a5471072be7ca6f6697e744ca89d5 (diff) | |
download | lapack-4bb091197931c4b836fd6eec84b958c2f24a4337.tar.gz lapack-4bb091197931c4b836fd6eec84b958c2f24a4337.tar.bz2 lapack-4bb091197931c4b836fd6eec84b958c2f24a4337.zip |
fixes to dlasyf_rook.f
Diffstat (limited to 'SRC')
-rw-r--r-- | SRC/dlasyf_rook.f | 27 |
1 files changed, 16 insertions, 11 deletions
diff --git a/SRC/dlasyf_rook.f b/SRC/dlasyf_rook.f index 986839e2..6c8532df 100644 --- a/SRC/dlasyf_rook.f +++ b/SRC/dlasyf_rook.f @@ -170,9 +170,13 @@ *> *> \verbatim *> -*> November 2011, Igor Kozachenko, University of California Berkeley +*> November 2011, Igor Kozachenko, +*> Computer Science Division, +*> University of California, Berkeley *> -*> September 2007, Craig Lucas, University of Manchester +*> September 2007, Sven Hammarling, Nicholas J. Higham, Craig Lucas, +*> School of Mathematics, +*> University of Manchester *> *> \endverbatim * @@ -320,7 +324,7 @@ * IF( K.LT.N ) $ CALL DGEMV( 'No transpose', K, N-K, -ONE, - $ A( 1, K+1 ), LDA, W( IMAX, KW+1 ), LDW + $ A( 1, K+1 ), LDA, W( IMAX, KW+1 ), LDW, $ ONE,W( 1, KW-1 ), 1 ) * * JMAX is the column-index of the largest off-diagonal @@ -399,13 +403,13 @@ * IF( ( KSTEP.EQ.2 ) .AND. ( P.NE.K ) ) THEN * -* Copy non-updated column KK+1 to column P +* Copy non-updated column K to column P * - CALL DCOPY( P-K, A( K, K ), 1, A( P, K ), LDA ) - CALL DCOPY( N-P+1, A( P, K ), 1, A( P, P ), 1 ) + CALL DCOPY( K-P, A( P+1, K ), 1, A( P, P+1 ), LDA ) + CALL DCOPY( P, A( 1, K ), 1, A( 1, P ), 1 ) * -* Interchange rows K and P in last K columns of A -* and last KK columns of W +* Interchange rows K and P in last N-K+1 columns of A +* and last N-K+2 columns of W * CALL DSWAP( N-K+1, A( K, K ), LDA, A( P, K ), LDA ) CALL DSWAP( N-KK+1, W( K, KKW ), LDW, W( P, KKW ), LDW ) @@ -422,7 +426,8 @@ $ LDA ) CALL DCOPY( KP, A( 1, KK ), 1, A( 1, KP ), 1 ) * -* Interchange rows KK and KP in last KK columns of A and W +* Interchange rows KK and KP in last N-KK+1 columns +* of A and W * CALL DSWAP( N-KK+1, A( KK, KK ), LDA, A( KP, KK ), LDA ) CALL DSWAP( N-KK+1, W( KK, KKW ), LDW, W( KP, KKW ), @@ -712,13 +717,13 @@ * IF( ( KSTEP.EQ.2 ) .AND. ( P.NE.K ) ) THEN * -* Copy non-updated column KK-1 to column P +* Copy non-updated column K to column P * CALL DCOPY( P-K, A( K, K ), 1, A( P, K ), LDA ) CALL DCOPY( N-P+1, A( P, K ), 1, A( P, P ), 1 ) * * Interchange rows K and P in first K columns of A -* and first KK columns of W +* and first K+1 columns of W * CALL DSWAP( K, A( K, 1 ), LDA, A( P, 1 ), LDA ) CALL DSWAP( KK, W( K, 1 ), LDW, W( P, 1 ), LDW ) |