From eb41c88b5f1eb7621e26d92a071e05f983a1701c Mon Sep 17 00:00:00 2001
From: sh-zheng <2294474733@qq.com>
Date: Wed, 13 Mar 2024 23:41:28 +0800
Subject: [PATCH] Fix bug for dlakyf and dkytf2

---
 SRC/dkytf2.f | 8 ++++----
 SRC/dlakyf.f | 8 ++++----
 2 files changed, 8 insertions(+), 8 deletions(-)

diff --git a/SRC/dkytf2.f b/SRC/dkytf2.f
index 69670b2ea..d24e21516 100644
--- 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
index e715a63c9..1f69c5ed6 100644
--- 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