version 1.2, 2010/04/21 13:45:16
|
version 1.5, 2010/08/07 13:18:06
|
Line 1
|
Line 1
|
SUBROUTINE DLAGV2( A, LDA, B, LDB, ALPHAR, ALPHAI, BETA, CSL, SNL, |
SUBROUTINE DLAGV2( A, LDA, B, LDB, ALPHAR, ALPHAI, BETA, CSL, SNL, |
$ CSR, SNR ) |
$ CSR, SNR ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.2) -- |
* -- LAPACK auxiliary routine (version 3.2.2) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* November 2006 |
* June 2010 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER LDA, LDB |
INTEGER LDA, LDB |
Line 142
|
Line 142
|
SNR = ZERO |
SNR = ZERO |
A( 2, 1 ) = ZERO |
A( 2, 1 ) = ZERO |
B( 2, 1 ) = ZERO |
B( 2, 1 ) = ZERO |
|
WI = ZERO |
* |
* |
* Check if B is singular |
* Check if B is singular |
* |
* |
Line 154
|
Line 155
|
A( 2, 1 ) = ZERO |
A( 2, 1 ) = ZERO |
B( 1, 1 ) = ZERO |
B( 1, 1 ) = ZERO |
B( 2, 1 ) = ZERO |
B( 2, 1 ) = ZERO |
|
WI = ZERO |
* |
* |
ELSE IF( ABS( B( 2, 2 ) ).LE.ULP ) THEN |
ELSE IF( ABS( B( 2, 2 ) ).LE.ULP ) THEN |
CALL DLARTG( A( 2, 2 ), A( 2, 1 ), CSR, SNR, T ) |
CALL DLARTG( A( 2, 2 ), A( 2, 1 ), CSR, SNR, T ) |
Line 165
|
Line 167
|
A( 2, 1 ) = ZERO |
A( 2, 1 ) = ZERO |
B( 2, 1 ) = ZERO |
B( 2, 1 ) = ZERO |
B( 2, 2 ) = ZERO |
B( 2, 2 ) = ZERO |
|
WI = ZERO |
* |
* |
ELSE |
ELSE |
* |
* |