rpl/lapack/lapack/dlagv2.f - diff

Return to dlagv2.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.4, 2010/08/06 15:32:26	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
*	*