diff --git a/SRC/dtgsna.f b/SRC/dtgsna.f index b9b3ad8afb..013dc91bd6 100644 --- a/SRC/dtgsna.f +++ b/SRC/dtgsna.f @@ -632,8 +632,8 @@ SUBROUTINE DTGSNA( JOB, HOWMNY, SELECT, N, A, LDA, B, LDB, VL, C1 = TWO*( ALPHAR*ALPHAR+ALPHAI*ALPHAI+BETA*BETA ) C2 = FOUR*BETA*BETA*ALPHAI*ALPHAI ROOT1 = C1 + SQRT( C1*C1-4.0D0*C2 ) - ROOT2 = C2 / ROOT1 ROOT1 = ROOT1 / TWO + ROOT2 = C2 / ROOT1 COND = MIN( SQRT( ROOT1 ), SQRT( ROOT2 ) ) END IF * diff --git a/SRC/stgsna.f b/SRC/stgsna.f index 40f822915f..430f3c4b7c 100644 --- a/SRC/stgsna.f +++ b/SRC/stgsna.f @@ -632,8 +632,8 @@ SUBROUTINE STGSNA( JOB, HOWMNY, SELECT, N, A, LDA, B, LDB, VL, C1 = TWO*( ALPHAR*ALPHAR+ALPHAI*ALPHAI+BETA*BETA ) C2 = FOUR*BETA*BETA*ALPHAI*ALPHAI ROOT1 = C1 + SQRT( C1*C1-4.0*C2 ) - ROOT2 = C2 / ROOT1 ROOT1 = ROOT1 / TWO + ROOT2 = C2 / ROOT1 COND = MIN( SQRT( ROOT1 ), SQRT( ROOT2 ) ) END IF *