Skip to content

Commit 87ac1ce

Browse files
authored
Apply fix from Reference-LAPACK PR390, NaN not propagating
1 parent 9e40c08 commit 87ac1ce

1 file changed

Lines changed: 2 additions & 0 deletions

File tree

lapack-netlib/SRC/scombssq.f

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -80,6 +80,8 @@ SUBROUTINE SCOMBSSQ( V1, V2 )
8080
IF( V1( 1 ).GE.V2( 1 ) ) THEN
8181
IF( V1( 1 ).NE.ZERO ) THEN
8282
V1( 2 ) = V1( 2 ) + ( V2( 1 ) / V1( 1 ) )**2 * V2( 2 )
83+
ELSE
84+
V1( 2 ) = V1( 2 ) + V2( 2 )
8385
END IF
8486
ELSE
8587
V1( 2 ) = V2( 2 ) + ( V1( 1 ) / V2( 1 ) )**2 * V1( 2 )

0 commit comments

Comments
 (0)