Fix bug for dlakyf and dkytf2

sh-zheng · Mar 13, 2024 · eb41c88 · eb41c88
1 parent 3e3f8f9
commit eb41c88
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 8 deletions.
diff --git a/SRC/dkytf2.f b/SRC/dkytf2.f
@@ -295,9 +295,9 @@ SUBROUTINE DKYTF2( UPLO, N, A, LDA, IPIV, INFO )
 *        COLMAX1 and COLMAX2 are their absolute values.
 *
          IF(K.GT.2) THEN
-            IMAX1 = ISAMAX( K-2, A( 1, K ), 1 )
+            IMAX1 = IDAMAX( K-2, A( 1, K ), 1 )
             COLMAX1 = ABS( A( IMAX1, K ) )
-            IMAX2 = ISAMAX( K-2, A( 1, K-1 ), 1 )
+            IMAX2 = IDAMAX( K-2, A( 1, K-1 ), 1 )
             COLMAX2 = ABS( A( IMAX2, K-1 ) )
          ELSE
             IMAX1 = 0
@@ -452,9 +452,9 @@ SUBROUTINE DKYTF2( UPLO, N, A, LDA, IPIV, INFO )
 *        COLMAX1 and COLMAX2 are their absolute values.
 *
          IF(K.LT.N-1) THEN
-            IMAX1 = K+1 + ISAMAX( N-K-1, A( K+2, K ), 1 )
+            IMAX1 = K+1 + IDAMAX( N-K-1, A( K+2, K ), 1 )
             COLMAX1 = ABS( A( IMAX1, K ) )
-            IMAX2 = K+1 + ISAMAX( N-K-1, A( K+2, K+1 ), 1 )
+            IMAX2 = K+1 + IDAMAX( N-K-1, A( K+2, K+1 ), 1 )
             COLMAX2 = ABS( A( IMAX2, K+1 ) )
          ELSE
             IMAX1 = 0

diff --git a/SRC/dlakyf.f b/SRC/dlakyf.f
@@ -293,9 +293,9 @@ SUBROUTINE DLAKYF( UPLO, N, NB, KB, A, LDA, IPIV, W, LDW, INFO )
 *        COLMAX1 and COLMAX2 are their absolute values.
 *
          IF(K.GT.2) THEN
-            IMAX1 = ISAMAX( K-2, W( 1, KW ), 1 )
+            IMAX1 = IDAMAX( K-2, W( 1, KW ), 1 )
             COLMAX1 = ABS( W( IMAX1, KW ) )
-            IMAX2 = ISAMAX( K-2, W( 1, KW-1 ), 1 )
+            IMAX2 = IDAMAX( K-2, W( 1, KW-1 ), 1 )
             COLMAX2 = ABS( W( IMAX2, KW-1 ) )
          ELSE
             IMAX1 = 0
@@ -604,9 +604,9 @@ SUBROUTINE DLAKYF( UPLO, N, NB, KB, A, LDA, IPIV, W, LDW, INFO )
 *        COLMAX1 and COLMAX2 are their absolute values.
 *
          IF(K.LT.N-1) THEN
-            IMAX1 = K+1 + ISAMAX( N-K-1, W( K+2, K ), 1 )
+            IMAX1 = K+1 + IDAMAX( N-K-1, W( K+2, K ), 1 )
             COLMAX1 = ABS( W( IMAX1, K ) )
-            IMAX2 = K+1 + ISAMAX( N-K-1, W( K+2, K+1 ), 1 )
+            IMAX2 = K+1 + IDAMAX( N-K-1, W( K+2, K+1 ), 1 )
             COLMAX2 = ABS( W( IMAX2, K+1 ) )
          ELSE
             IMAX1 = 0