See post on the forum from Pauli Virtanen.

	http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4560

Documentation of xSTEVR says that VL and VU are not referenced if RANGE = 'A'
or 'I', however the code does reference VL and VU for any value of RANGE:

      VLL = VL
      VUU = VU

So simply protect these two lines of code:

      IF( VALEIG ) THEN
         VLL = VL
         VUU = VU
      END IF

Julien.

diff --git a/SRC/dstevr.f b/SRC/dstevr.f
index 076ab92e9dc2fe9e7e31e298b95db1e791cb1a92..f393cfd6bdbc3c2b6394b360a8bd3ed7349c9687 100644 (file)
--- a/SRC/dstevr.f
+++ b/SRC/dstevr.f
@@ -435,8 +435,10 @@
 *     Scale matrix to allowable range, if necessary.
 *
       ISCALE = 0
-      VLL = VL
-      VUU = VU
+      IF( VALEIG ) THEN
+         VLL = VL
+         VUU = VU
+      END IF
 *
       TNRM = DLANST( 'M', N, D, E )
       IF( TNRM.GT.ZERO .AND. TNRM.LT.RMIN ) THEN

diff --git a/SRC/sstevr.f b/SRC/sstevr.f
index 556a78b6ed22425e2e1b10d018e989fb618cb1b7..74157cead082d88e8e8a57e03abc6020cc928f50 100644 (file)
--- a/SRC/sstevr.f
+++ b/SRC/sstevr.f
@@ -436,8 +436,10 @@
 *     Scale matrix to allowable range, if necessary.
 *
       ISCALE = 0
-      VLL = VL
-      VUU = VU
+      IF( VALEIG ) THEN
+         VLL = VL
+         VUU = VU
+      END IF
 *
       TNRM = SLANST( 'M', N, D, E )
       IF( TNRM.GT.ZERO .AND. TNRM.LT.RMIN ) THEN