From 2aad88d5b9ded514d65c257cea818165447e5b78 Mon Sep 17 00:00:00 2001
From: Martin Kroeker <martin@ruby.chemie.uni-freiburg.de>
Date: Sat, 27 Apr 2019 23:01:49 +0200
Subject: [PATCH] Avoid out-of-bounds accesses in LAPACK EIG tests

see https://github.com/Reference-LAPACK/lapack/issues/333
---
 lapack-netlib/TESTING/EIG/chet21.f | 3 ++-
 lapack-netlib/TESTING/EIG/chpt21.f | 2 +-
 lapack-netlib/TESTING/EIG/zhet21.f | 3 ++-
 lapack-netlib/TESTING/EIG/zhpt21.f | 3 ++-
 4 files changed, 7 insertions(+), 4 deletions(-)

diff --git a/lapack-netlib/TESTING/EIG/chet21.f b/lapack-netlib/TESTING/EIG/chet21.f
index 8dbdb521e..5aff64904 100644
--- a/lapack-netlib/TESTING/EIG/chet21.f
+++ b/lapack-netlib/TESTING/EIG/chet21.f
@@ -304,7 +304,8 @@
    10    CONTINUE
 *
          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,
      $                     U( 1, J-1 ), 1, WORK, N )
    20       CONTINUE
diff --git a/lapack-netlib/TESTING/EIG/chpt21.f b/lapack-netlib/TESTING/EIG/chpt21.f
index 4b9279470..e151a8bd8 100644
--- a/lapack-netlib/TESTING/EIG/chpt21.f
+++ b/lapack-netlib/TESTING/EIG/chpt21.f
@@ -323,7 +323,7 @@
    10    CONTINUE
 *
          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,
      $                     U( 1, J-1 ), 1, WORK )
    20       CONTINUE
diff --git a/lapack-netlib/TESTING/EIG/zhet21.f b/lapack-netlib/TESTING/EIG/zhet21.f
index 32a09741e..f6cb2d70a 100644
--- a/lapack-netlib/TESTING/EIG/zhet21.f
+++ b/lapack-netlib/TESTING/EIG/zhet21.f
@@ -304,7 +304,8 @@
    10    CONTINUE
 *
          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,
      $                     U( 1, J-1 ), 1, WORK, N )
    20       CONTINUE
diff --git a/lapack-netlib/TESTING/EIG/zhpt21.f b/lapack-netlib/TESTING/EIG/zhpt21.f
index f9268661a..ef9e4418d 100644
--- a/lapack-netlib/TESTING/EIG/zhpt21.f
+++ b/lapack-netlib/TESTING/EIG/zhpt21.f
@@ -323,7 +323,8 @@
    10    CONTINUE
 *
          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,
      $                     U( 1, J-1 ), 1, WORK )
    20       CONTINUE