rpl/lapack/lapack/dlanv2.f - diff

Return to dlanv2.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.18, 2018/05/29 07:17:58	version 1.19, 2020/05/21 21:45:59
Line 161	Line 161
IF( C.EQ.ZERO ) THEN	IF( C.EQ.ZERO ) THEN
CS = ONE	CS = ONE
SN = ZERO	SN = ZERO
GO TO 10
*	*
ELSE IF( B.EQ.ZERO ) THEN	ELSE IF( B.EQ.ZERO ) THEN
*	*
Line 174	Line 173
A = TEMP	A = TEMP
B = -C	B = -C
C = ZERO	C = ZERO
GO TO 10	*
ELSE IF( ( A-D ).EQ.ZERO .AND. SIGN( ONE, B ).NE.SIGN( ONE, C ) )	ELSE IF( ( A-D ).EQ.ZERO .AND. SIGN( ONE, B ).NE.SIGN( ONE, C ) )
$ THEN	$ THEN
CS = ONE	CS = ONE
SN = ZERO	SN = ZERO
GO TO 10	*
ELSE	ELSE
*	*
TEMP = A - D	TEMP = A - D
Line 207	Line 206
SN = C / TAU	SN = C / TAU
B = B - C	B = B - C
C = ZERO	C = ZERO
	*
ELSE	ELSE
*	*
* Complex eigenvalues, or real (almost) equal eigenvalues.	* Complex eigenvalues, or real (almost) equal eigenvalues.
Line 269	Line 269
*	*
END IF	END IF
*	*
10 CONTINUE
*
* Store eigenvalues in (RT1R,RT1I) and (RT2R,RT2I).	* Store eigenvalues in (RT1R,RT1I) and (RT2R,RT2I).
*	*
RT1R = A	RT1R = A