From ef12723763577e826f6b0b62767d38ab4ae6f017 Mon Sep 17 00:00:00 2001 From: EduardFedorenkov Date: Fri, 17 Feb 2023 15:32:00 +0700 Subject: [PATCH] fix bug in complex precision tests (c|z)het21 --- TESTING/EIG/chet21.f | 4 ++-- TESTING/EIG/zhet21.f | 4 ++-- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/TESTING/EIG/chet21.f b/TESTING/EIG/chet21.f index a274681c98..1bd35cb41c 100644 --- a/TESTING/EIG/chet21.f +++ b/TESTING/EIG/chet21.f @@ -304,9 +304,9 @@ SUBROUTINE CHET21( ITYPE, UPLO, N, KBAND, A, LDA, D, E, U, LDU, V, 10 CONTINUE * IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN - DO 20 J = 2, N - 1 + DO 20 J = 1, N - 1 CALL CHER2( CUPLO, N, -CMPLX( E( J ) ), U( 1, J ), 1, - $ U( 1, J-1 ), 1, WORK, N ) + $ U( 1, J+1 ), 1, WORK, N ) 20 CONTINUE END IF WNORM = CLANHE( '1', CUPLO, N, WORK, N, RWORK ) diff --git a/TESTING/EIG/zhet21.f b/TESTING/EIG/zhet21.f index d254f85e9c..b927a502fa 100644 --- a/TESTING/EIG/zhet21.f +++ b/TESTING/EIG/zhet21.f @@ -304,9 +304,9 @@ SUBROUTINE ZHET21( ITYPE, UPLO, N, KBAND, A, LDA, D, E, U, LDU, V, 10 CONTINUE * IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN - DO 20 J = 2, N - 1 + DO 20 J = 1, N - 1 CALL ZHER2( CUPLO, N, -DCMPLX( E( J ) ), U( 1, J ), 1, - $ U( 1, J-1 ), 1, WORK, N ) + $ U( 1, J+1 ), 1, WORK, N ) 20 CONTINUE END IF WNORM = ZLANHE( '1', CUPLO, N, WORK, N, RWORK )