Bug report sent by Alexander V. Kobotov (from Intel) on Mon, 6 Apr 2009 to "lapack...

"(d/s)tgsen.f: iwork(1) always referenced: line 455: iwork( 1 ) = LIWMIN,
documentation says that if IJOB=0 it shouldn't, so NULL pointer causes a
sigfault."

There is indeed the same problem for array work.

I have changed the header of the routines ctgsen.f, dtgsen.f, stgsen.f, and
ztgsen.f. No matter what an array needs to be of size at least 1. So
IWORK is of size at least 1. It was written in the header of the routine:
"IWORK   (workspace/output) INTEGER array, dimension (MAX(1,LIWORK))"
So since IWORK is of size at least 1, there is no reason not to reference it.
When there is a workspace query with IJOB.EQ.0, IWORK(1) is set to 1.

diff --git a/SRC/ctgsen.f b/SRC/ctgsen.f
index fb3930ec8cf223477c3d030e84d521502a2e50db..68fd961467dc9dd5637ada384d31559a3f1ab584 100644 (file)
--- a/SRC/ctgsen.f
+++ b/SRC/ctgsen.f
@@ -169,8 +169,7 @@
 *          message related to LWORK is issued by XERBLA.
 *
 *  IWORK   (workspace/output) INTEGER array, dimension (MAX(1,LIWORK))
-*          IF IJOB = 0, IWORK is not referenced.  Otherwise,
-*          on exit, if INFO = 0, IWORK(1) returns the optimal LIWORK.
+*          On exit, if INFO = 0, IWORK(1) returns the optimal LIWORK.
 *
 *  LIWORK  (input) INTEGER
 *          The dimension of the array IWORK. LIWORK >= 1.

diff --git a/SRC/dtgsen.f b/SRC/dtgsen.f
index 03a9e35b2dc2ac4f3c683ca4f69d02d9b950c9f1..06a16a302b5ba9449ee583de27e9108dfbf013da 100644 (file)
--- a/SRC/dtgsen.f
+++ b/SRC/dtgsen.f
@@ -179,8 +179,7 @@
 *          message related to LWORK is issued by XERBLA.
 *
 *  IWORK   (workspace/output) INTEGER array, dimension (MAX(1,LIWORK))
-*          IF IJOB = 0, IWORK is not referenced.  Otherwise,
-*          on exit, if INFO = 0, IWORK(1) returns the optimal LIWORK.
+*          On exit, if INFO = 0, IWORK(1) returns the optimal LIWORK.
 *
 *  LIWORK  (input) INTEGER
 *          The dimension of the array IWORK. LIWORK >= 1.

diff --git a/SRC/stgsen.f b/SRC/stgsen.f
index 58a4ba9acd6f4799d3a6397a91f6bb7b495b1f6c..fec207b1a02dad2940939f0d709a604c83c57dc3 100644 (file)
--- a/SRC/stgsen.f
+++ b/SRC/stgsen.f
@@ -178,8 +178,7 @@
 *          message related to LWORK is issued by XERBLA.
 *
 *  IWORK   (workspace/output) INTEGER array, dimension (MAX(1,LIWORK))
-*          IF IJOB = 0, IWORK is not referenced.  Otherwise,
-*          on exit, if INFO = 0, IWORK(1) returns the optimal LIWORK.
+*          On exit, if INFO = 0, IWORK(1) returns the optimal LIWORK.
 *
 *  LIWORK  (input) INTEGER
 *          The dimension of the array IWORK. LIWORK >= 1.

diff --git a/SRC/ztgsen.f b/SRC/ztgsen.f
index 399a1d9075cc872b27fb787870114f4a82b69607..ff7dcae7827802dddb7996639f65acf3b8581e95 100644 (file)
--- a/SRC/ztgsen.f
+++ b/SRC/ztgsen.f
@@ -169,8 +169,7 @@
 *          message related to LWORK is issued by XERBLA.
 *
 *  IWORK   (workspace/output) INTEGER array, dimension (MAX(1,LIWORK))
-*          IF IJOB = 0, IWORK is not referenced.  Otherwise,
-*          on exit, if INFO = 0, IWORK(1) returns the optimal LIWORK.
+*          On exit, if INFO = 0, IWORK(1) returns the optimal LIWORK.
 *
 *  LIWORK  (input) INTEGER
 *          The dimension of the array IWORK. LIWORK >= 1.