version 1.3, 2010/08/06 15:28:39
|
version 1.9, 2011/07/22 07:38: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 |
* |
* |
Line 239
|
Line 242
|
CALL DLASV2( B( 1, 1 ), B( 1, 2 ), B( 2, 2 ), R, T, SNR, |
CALL DLASV2( B( 1, 1 ), B( 1, 2 ), B( 2, 2 ), R, T, SNR, |
$ CSR, SNL, CSL ) |
$ CSR, SNL, CSL ) |
* |
* |
* Form (A,B) := Q(A,B)Z' where Q is left rotation matrix and |
* Form (A,B) := Q(A,B)Z**T where Q is left rotation matrix and |
* Z is right rotation matrix computed from DLASV2 |
* Z is right rotation matrix computed from DLASV2 |
* |
* |
CALL DROT( 2, A( 1, 1 ), LDA, A( 2, 1 ), LDA, CSL, SNL ) |
CALL DROT( 2, A( 1, 1 ), LDA, A( 2, 1 ), LDA, CSL, SNL ) |