fix scaling (Reference-LAPACK PR 830)

2023-05-24 16:52:56 +02:00
parent e2f4adfc60
commit 912bfd9bfb
1 changed files with 1 additions and 1 deletions
--- a/lapack-netlib/SRC/stgsna.f
+++ b/lapack-netlib/SRC/stgsna.f
@@ -632,8 +632,8 @@
               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
 *