From 5db45a2c4ea2c317b8f3e1b575f14d1ca66dba19 Mon Sep 17 00:00:00 2001 From: Weslley S Pereira Date: Tue, 23 May 2023 14:51:32 -0600 Subject: [PATCH] Fix 103 thanks to @elivanova --- SRC/dtgsna.f | 2 +- SRC/stgsna.f | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) 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 *