Skip to content

Commit 23587aa

Browse files
authored
Merge pull request #798 from EduardFedorenkov/master
fix bug in complex precision tests (c|z)het21
2 parents eb8f5fa + ef12723 commit 23587aa

File tree

2 files changed

+4
-4
lines changed

2 files changed

+4
-4
lines changed

TESTING/EIG/chet21.f

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -304,9 +304,9 @@ SUBROUTINE CHET21( ITYPE, UPLO, N, KBAND, A, LDA, D, E, U, LDU, V,
304304
10 CONTINUE
305305
*
306306
IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN
307-
DO 20 J = 2, N - 1
307+
DO 20 J = 1, N - 1
308308
CALL CHER2( CUPLO, N, -CMPLX( E( J ) ), U( 1, J ), 1,
309-
$ U( 1, J-1 ), 1, WORK, N )
309+
$ U( 1, J+1 ), 1, WORK, N )
310310
20 CONTINUE
311311
END IF
312312
WNORM = CLANHE( '1', CUPLO, N, WORK, N, RWORK )

TESTING/EIG/zhet21.f

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -304,9 +304,9 @@ SUBROUTINE ZHET21( ITYPE, UPLO, N, KBAND, A, LDA, D, E, U, LDU, V,
304304
10 CONTINUE
305305
*
306306
IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN
307-
DO 20 J = 2, N - 1
307+
DO 20 J = 1, N - 1
308308
CALL ZHER2( CUPLO, N, -DCMPLX( E( J ) ), U( 1, J ), 1,
309-
$ U( 1, J-1 ), 1, WORK, N )
309+
$ U( 1, J+1 ), 1, WORK, N )
310310
20 CONTINUE
311311
END IF
312312
WNORM = ZLANHE( '1', CUPLO, N, WORK, N, RWORK )

0 commit comments

Comments
 (0)