Avoid out-of-bounds accesses in LAPACK EIG teststags/v0.3.6^2
@@ -304,7 +304,8 @@ | |||||
10 CONTINUE | 10 CONTINUE | ||||
* | * | ||||
IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN | IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN | ||||
DO 20 J = 1, N - 1 | |||||
CMK DO 20 J = 1, N - 1 | |||||
DO 20 J = 2, N - 1 | |||||
CALL CHER2( CUPLO, N, -CMPLX( E( J ) ), U( 1, J ), 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 | 20 CONTINUE | ||||
@@ -323,7 +323,7 @@ | |||||
10 CONTINUE | 10 CONTINUE | ||||
* | * | ||||
IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN | IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN | ||||
DO 20 J = 1, N - 1 | |||||
DO 20 J = 2, N - 1 | |||||
CALL CHPR2( CUPLO, N, -CMPLX( E( J ) ), U( 1, J ), 1, | CALL CHPR2( CUPLO, N, -CMPLX( E( J ) ), U( 1, J ), 1, | ||||
$ U( 1, J-1 ), 1, WORK ) | $ U( 1, J-1 ), 1, WORK ) | ||||
20 CONTINUE | 20 CONTINUE | ||||
@@ -304,7 +304,8 @@ | |||||
10 CONTINUE | 10 CONTINUE | ||||
* | * | ||||
IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN | IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN | ||||
DO 20 J = 1, N - 1 | |||||
CMK DO 20 J = 1, N - 1 | |||||
DO 20 J = 2, N - 1 | |||||
CALL ZHER2( CUPLO, N, -DCMPLX( E( J ) ), U( 1, J ), 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 | 20 CONTINUE | ||||
@@ -323,7 +323,8 @@ | |||||
10 CONTINUE | 10 CONTINUE | ||||
* | * | ||||
IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN | IF( N.GT.1 .AND. KBAND.EQ.1 ) THEN | ||||
DO 20 J = 1, N - 1 | |||||
CMK DO 20 J = 1, N - 1 | |||||
DO 20 J = 2, N - 1 | |||||
CALL ZHPR2( CUPLO, N, -DCMPLX( E( J ) ), U( 1, J ), 1, | CALL ZHPR2( CUPLO, N, -DCMPLX( E( J ) ), U( 1, J ), 1, | ||||
$ U( 1, J-1 ), 1, WORK ) | $ U( 1, J-1 ), 1, WORK ) | ||||
20 CONTINUE | 20 CONTINUE | ||||