Avoid out-of-bounds accesses in LAPACK EIG tests
see https://github.com/Reference-LAPACK/lapack/issues/333
This commit is contained in:
parent
32c7063cb0
commit
2aad88d5b9
|
@ -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
|
||||||
|
|
Loading…
Reference in New Issue