fix bug in complex precision tests (c|z)het21

EduardFedorenkov · EduardFedorenkov · commit ef1272376357 · 2023-02-17T15:32:00.000+07:00
diff --git 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
@@ -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 )
