fix #302

langou · langou · commit eecc00bb1cc8 · 2018-11-26T08:43:43.000-07:00
diff --git a/SRC/dlanv2.f b/SRC/dlanv2.f
@@ -161,7 +161,6 @@ SUBROUTINE DLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )
       IF( C.EQ.ZERO ) THEN
          CS = ONE
          SN = ZERO
-         GO TO 10
 *
       ELSE IF( B.EQ.ZERO ) THEN
 *
@@ -174,12 +173,12 @@ SUBROUTINE DLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )
          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 @@ SUBROUTINE DLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )
             SN = C / TAU
             B = B - C
             C = ZERO
+*
          ELSE
 *
 *           Complex eigenvalues, or real (almost) equal eigenvalues.
@@ -268,8 +268,6 @@ SUBROUTINE DLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )
          END IF
 *
       END IF
-*
-   10 CONTINUE
 *
 *     Store eigenvalues in (RT1R,RT1I) and (RT2R,RT2I).
 *
diff --git a/SRC/slanv2.f b/SRC/slanv2.f
@@ -161,7 +161,6 @@ SUBROUTINE SLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )
       IF( C.EQ.ZERO ) THEN
          CS = ONE
          SN = ZERO
-         GO TO 10
 *
       ELSE IF( B.EQ.ZERO ) THEN
 *
@@ -174,12 +173,12 @@ SUBROUTINE SLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )
          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 @@ SUBROUTINE SLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )
             SN = C / TAU
             B = B - C
             C = ZERO
+*
          ELSE
 *
 *           Complex eigenvalues, or real (almost) equal eigenvalues.
@@ -268,8 +268,6 @@ SUBROUTINE SLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )
          END IF
 *
       END IF
-*
-   10 CONTINUE
 *
 *     Store eigenvalues in (RT1R,RT1I) and (RT2R,RT2I).
 *

Original file line number	Diff line number	Diff line change
`@@ -161,7 +161,6 @@ SUBROUTINE DLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )`
`161`	`161`	`IF( C.EQ.ZERO ) THEN`
`162`	`162`	`CS = ONE`
`163`	`163`	`SN = ZERO`
`164`		`- GO TO 10`
`165`	`164`	`*`
`166`	`165`	`ELSE IF( B.EQ.ZERO ) THEN`
`167`	`166`	`*`
`@@ -174,12 +173,12 @@ SUBROUTINE DLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )`
`174`	`173`	`A = TEMP`
`175`	`174`	`B = -C`
`176`	`175`	`C = ZERO`
`177`		`- GO TO 10`
	`176`	`+*`
`178`	`177`	`ELSE IF( ( A-D ).EQ.ZERO .AND. SIGN( ONE, B ).NE.SIGN( ONE, C ) )`
`179`	`178`	`$ THEN`
`180`	`179`	`CS = ONE`
`181`	`180`	`SN = ZERO`
`182`		`- GO TO 10`
	`181`	`+*`
`183`	`182`	`ELSE`
`184`	`183`	`*`
`185`	`184`	`TEMP = A - D`
`@@ -207,6 +206,7 @@ SUBROUTINE DLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )`
`207`	`206`	`SN = C / TAU`
`208`	`207`	`B = B - C`
`209`	`208`	`C = ZERO`
	`209`	`+*`
`210`	`210`	`ELSE`
`211`	`211`	`*`
`212`	`212`	`* Complex eigenvalues, or real (almost) equal eigenvalues.`
`@@ -268,8 +268,6 @@ SUBROUTINE DLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )`
`268`	`268`	`END IF`
`269`	`269`	`*`
`270`	`270`	`END IF`
`271`		`-*`
`272`		`- 10 CONTINUE`
`273`	`271`	`*`
`274`	`272`	`* Store eigenvalues in (RT1R,RT1I) and (RT2R,RT2I).`
`275`	`273`	`*`
Original file line number	Diff line number	Diff line change
`@@ -161,7 +161,6 @@ SUBROUTINE SLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )`
`161`	`161`	`IF( C.EQ.ZERO ) THEN`
`162`	`162`	`CS = ONE`
`163`	`163`	`SN = ZERO`
`164`		`- GO TO 10`
`165`	`164`	`*`
`166`	`165`	`ELSE IF( B.EQ.ZERO ) THEN`
`167`	`166`	`*`
`@@ -174,12 +173,12 @@ SUBROUTINE SLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )`
`174`	`173`	`A = TEMP`
`175`	`174`	`B = -C`
`176`	`175`	`C = ZERO`
`177`		`- GO TO 10`
	`176`	`+*`
`178`	`177`	`ELSE IF( (A-D).EQ.ZERO .AND. SIGN( ONE, B ).NE.`
`179`	`178`	`$ SIGN( ONE, C ) ) THEN`
`180`	`179`	`CS = ONE`
`181`	`180`	`SN = ZERO`
`182`		`- GO TO 10`
	`181`	`+*`
`183`	`182`	`ELSE`
`184`	`183`	`*`
`185`	`184`	`TEMP = A - D`
`@@ -207,6 +206,7 @@ SUBROUTINE SLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )`
`207`	`206`	`SN = C / TAU`
`208`	`207`	`B = B - C`
`209`	`208`	`C = ZERO`
	`209`	`+*`
`210`	`210`	`ELSE`
`211`	`211`	`*`
`212`	`212`	`* Complex eigenvalues, or real (almost) equal eigenvalues.`
`@@ -268,8 +268,6 @@ SUBROUTINE SLANV2( A, B, C, D, RT1R, RT1I, RT2R, RT2I, CS, SN )`
`268`	`268`	`END IF`
`269`	`269`	`*`
`270`	`270`	`END IF`
`271`		`-*`
`272`		`- 10 CONTINUE`
`273`	`271`	`*`
`274`	`272`	`* Store eigenvalues in (RT1R,RT1I) and (RT2R,RT2I).`
`275`	`273`	`*`