--- rpl/lapack/lapack/dlanv2.f 2011/11/21 20:42:56 1.9
+++ rpl/lapack/lapack/dlanv2.f 2020/05/21 21:45:59 1.19
@@ -1,29 +1,29 @@
-*> \brief \b DLANV2
+*> \brief \b DLANV2 computes the Schur factorization of a real 2-by-2 nonsymmetric matrix in standard form.
*
* =========== DOCUMENTATION ===========
*
-* Online html documentation available at
-* http://www.netlib.org/lapack/explore-html/
+* Online html documentation available at
+* http://www.netlib.org/lapack/explore-html/
*
*> \htmlonly
-*> Download DLANV2 + dependencies
-*>
-*> [TGZ]
-*>
-*> [ZIP]
-*>
+*> Download DLANV2 + dependencies
+*>
+*> [TGZ]
+*>
+*> [ZIP]
+*>
*> [TXT]
-*> \endhtmlonly
+*> \endhtmlonly
*
* Definition:
* ===========
*
* SUBROUTINE DLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )
-*
+*
* .. Scalar Arguments ..
* DOUBLE PRECISION A, B, C, CS, D, RT1I, RT1R, RT2I, RT2R, SN
* ..
-*
+*
*
*> \par Purpose:
* =============
@@ -104,12 +104,12 @@
* Authors:
* ========
*
-*> \author Univ. of Tennessee
-*> \author Univ. of California Berkeley
-*> \author Univ. of Colorado Denver
-*> \author NAG Ltd.
+*> \author Univ. of Tennessee
+*> \author Univ. of California Berkeley
+*> \author Univ. of Colorado Denver
+*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date December 2016
*
*> \ingroup doubleOTHERauxiliary
*
@@ -127,10 +127,10 @@
* =====================================================================
SUBROUTINE DLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )
*
-* -- LAPACK auxiliary routine (version 3.4.0) --
+* -- LAPACK auxiliary routine (version 3.7.0) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* December 2016
*
* .. Scalar Arguments ..
DOUBLE PRECISION A, B, C, CS, D, RT1I, RT1R, RT2I, RT2R, SN
@@ -161,7 +161,6 @@
IF( C.EQ.ZERO ) THEN
CS = ONE
SN = ZERO
- GO TO 10
*
ELSE IF( B.EQ.ZERO ) THEN
*
@@ -174,12 +173,12 @@
A = TEMP
B = -C
C = ZERO
- GO TO 10
+*
ELSE IF( ( A-D ).EQ.ZERO .AND. SIGN( ONE, B ).NE.SIGN( ONE, C ) )
$ THEN
CS = ONE
SN = ZERO
- GO TO 10
+*
ELSE
*
TEMP = A - D
@@ -207,6 +206,7 @@
SN = C / TAU
B = B - C
C = ZERO
+*
ELSE
*
* Complex eigenvalues, or real (almost) equal eigenvalues.
@@ -269,8 +269,6 @@
*
END IF
*
- 10 CONTINUE
-*
* Store eigenvalues in (RT1R,RT1I) and (RT2R,RT2I).
*
RT1R = A