diff options
| author | igor175 <igor175@8a072113-8704-0410-8d35-dd094bca7971> | 2013-07-09 04:02:51 +0000 |
|---|---|---|
| committer | igor175 <igor175@8a072113-8704-0410-8d35-dd094bca7971> | 2013-07-09 04:02:51 +0000 |
| commit | 127f04d2f09abfb7e2fe87ee1452d7299d13afe3 (patch) | |
| tree | 81d641ed684b282b48dfaef6cb53793c04827280 /SRC/zhetf2_rook.f | |
| parent | 2742b5208903229322750eb81c52435dadeb7282 (diff) | |
Modified new 'rook' pivoting routines and drivers for Hermitian indefinite matrices, Complex precision: zhetf2_rook.f zlahef_rook.f zhetri_rook.f
Diffstat (limited to 'SRC/zhetf2_rook.f')
| -rw-r--r-- | SRC/zhetf2_rook.f | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/SRC/zhetf2_rook.f b/SRC/zhetf2_rook.f index ec8da924..bddb8194 100644 --- a/SRC/zhetf2_rook.f +++ b/SRC/zhetf2_rook.f @@ -26,7 +26,7 @@ * .. * .. Array Arguments .. * INTEGER IPIV( * ) -* COMPLEX*16 A( LDA, * ) +* COMPLEX*16 A( LDA, * ) * .. * * @@ -360,7 +360,7 @@ * * Case(2) * Equivalent to testing for -* ABS( DBLE( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX +* ABS( REAL( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX * (used to handle NaN and Inf) * IF( .NOT.( ABS( DBLE( A( IMAX, IMAX ) ) ) @@ -408,7 +408,7 @@ * * KK is the column of A where pivoting step stopped * - KK = K + KSTEP - 1 + KK = K - KSTEP + 1 * * For only a 2x2 pivot, interchange rows and columns K and P * in the leading submatrix A(1:k,1:k) @@ -492,7 +492,7 @@ * * Store U(k) in column k * - CALL ZSCAL( K-1, D11, A( 1, K ), 1 ) + CALL ZDSCAL( K-1, D11, A( 1, K ), 1 ) ELSE * * Store L(k) in column K @@ -672,7 +672,7 @@ * * Case(2) * Equivalent to testing for -* ABS( DBLE( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX +* ABS( REAL( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX * (used to handle NaN and Inf) * IF( .NOT.( ABS( DBLE( A( IMAX, IMAX ) ) ) |