Browse Source
Merge pull request #2442 from martin-frbg/lapackpr390
Apply fix from Reference-LAPACK PR 390
tags/v0.3.9
Martin Kroeker
GitHub
6 years ago
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
2 changed files with
4 additions and
0 deletions
-
lapack-netlib/SRC/dcombssq.f
-
lapack-netlib/SRC/scombssq.f
|
|
|
@@ -80,6 +80,8 @@ |
|
|
|
IF( V1( 1 ).GE.V2( 1 ) ) THEN |
|
|
|
IF( V1( 1 ).NE.ZERO ) THEN |
|
|
|
V1( 2 ) = V1( 2 ) + ( V2( 1 ) / V1( 1 ) )**2 * V2( 2 ) |
|
|
|
ELSE |
|
|
|
V1( 2 ) = V1( 2 ) + V2( 2 ) |
|
|
|
END IF |
|
|
|
ELSE |
|
|
|
V1( 2 ) = V2( 2 ) + ( V1( 1 ) / V2( 1 ) )**2 * V1( 2 ) |
|
|
|
|
|
|
|
@@ -80,6 +80,8 @@ |
|
|
|
IF( V1( 1 ).GE.V2( 1 ) ) THEN |
|
|
|
IF( V1( 1 ).NE.ZERO ) THEN |
|
|
|
V1( 2 ) = V1( 2 ) + ( V2( 1 ) / V1( 1 ) )**2 * V2( 2 ) |
|
|
|
ELSE |
|
|
|
V1( 2 ) = V1( 2 ) + V2( 2 ) |
|
|
|
END IF |
|
|
|
ELSE |
|
|
|
V1( 2 ) = V2( 2 ) + ( V1( 1 ) / V2( 1 ) )**2 * V1( 2 ) |
|
|
|
|
