--- rpl/lapack/lapack/dlagv2.f 2010/04/21 13:45:16 1.2 +++ rpl/lapack/lapack/dlagv2.f 2011/07/22 07:38:06 1.9 @@ -1,10 +1,10 @@ SUBROUTINE DLAGV2( A, LDA, B, LDB, ALPHAR, ALPHAI, BETA, CSL, SNL, $ 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, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* November 2006 +* June 2010 * * .. Scalar Arguments .. INTEGER LDA, LDB @@ -142,6 +142,7 @@ SNR = ZERO A( 2, 1 ) = ZERO B( 2, 1 ) = ZERO + WI = ZERO * * Check if B is singular * @@ -154,6 +155,7 @@ A( 2, 1 ) = ZERO B( 1, 1 ) = ZERO B( 2, 1 ) = ZERO + WI = ZERO * ELSE IF( ABS( B( 2, 2 ) ).LE.ULP ) THEN CALL DLARTG( A( 2, 2 ), A( 2, 1 ), CSR, SNR, T ) @@ -165,6 +167,7 @@ A( 2, 1 ) = ZERO B( 2, 1 ) = ZERO B( 2, 2 ) = ZERO + WI = ZERO * ELSE * @@ -239,7 +242,7 @@ CALL DLASV2( B( 1, 1 ), B( 1, 2 ), B( 2, 2 ), R, T, SNR, $ 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 * CALL DROT( 2, A( 1, 1 ), LDA, A( 2, 1 ), LDA, CSL, SNL )