Browse Source
Merge pull request #3087 from martin-frbg/lapack477
Apply Reference-LAPACK PR 477 for convergence problems in CHGEQZ/ZHGEQZtags/v0.3.14^2
Martin Kroeker
GitHub
4 years ago
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
2 changed files with 16 additions and 4 deletions
Unified View
Diff Options
-
+8 -2lapack-netlib/SRC/chgeqz.f
-
+8 -2lapack-netlib/SRC/zhgeqz.f
+ 8
- 2
lapack-netlib/SRC/chgeqz.f
View File
| @@ -743,8 +743,14 @@ | |||||
| * | * | ||||
| * Exceptional shift. Chosen for no particularly good reason. | * Exceptional shift. Chosen for no particularly good reason. | ||||
| * | * | ||||
| ESHIFT = ESHIFT + (ASCALE*H(ILAST,ILAST-1))/ | |||||
| $ (BSCALE*T(ILAST-1,ILAST-1)) | |||||
| IF( ( IITER / 20 )*20.EQ.IITER .AND. | |||||
| $ BSCALE*ABS1(T( ILAST, ILAST )).GT.SAFMIN ) THEN | |||||
| ESHIFT = ESHIFT + ( ASCALE*H( ILAST, | |||||
| $ ILAST ) )/( BSCALE*T( ILAST, ILAST ) ) | |||||
| ELSE | |||||
| ESHIFT = ESHIFT + ( ASCALE*H( ILAST, | |||||
| $ ILAST-1 ) )/( BSCALE*T( ILAST-1, ILAST-1 ) ) | |||||
| END IF | |||||
| SHIFT = ESHIFT | SHIFT = ESHIFT | ||||
| END IF | END IF | ||||
| * | * | ||||
+ 8
- 2
lapack-netlib/SRC/zhgeqz.f
View File
| @@ -744,8 +744,14 @@ | |||||
| * | * | ||||
| * Exceptional shift. Chosen for no particularly good reason. | * Exceptional shift. Chosen for no particularly good reason. | ||||
| * | * | ||||
| ESHIFT = ESHIFT + (ASCALE*H(ILAST,ILAST-1))/ | |||||
| $ (BSCALE*T(ILAST-1,ILAST-1)) | |||||
| IF( ( IITER / 20 )*20.EQ.IITER .AND. | |||||
| $ BSCALE*ABS1(T( ILAST, ILAST )).GT.SAFMIN ) THEN | |||||
| ESHIFT = ESHIFT + ( ASCALE*H( ILAST, | |||||
| $ ILAST ) )/( BSCALE*T( ILAST, ILAST ) ) | |||||
| ELSE | |||||
| ESHIFT = ESHIFT + ( ASCALE*H( ILAST, | |||||
| $ ILAST-1 ) )/( BSCALE*T( ILAST-1, ILAST-1 ) ) | |||||
| END IF | |||||
| SHIFT = ESHIFT | SHIFT = ESHIFT | ||||
| END IF | END IF | ||||
| * | * | ||||