aboutsummaryrefslogtreecommitdiff
path: root/SRC/claed8.f
diff options
context:
space:
mode:
authorjulie <julielangou@users.noreply.github.com>2010-03-18 19:41:50 +0000
committerjulie <julielangou@users.noreply.github.com>2010-03-18 19:41:50 +0000
commit486f708677442b1a528addbd6696799c95bf2a0d (patch)
treedde81f9f446071a99c2b9543f3aa4eefd443708d /SRC/claed8.f
parent86e2baba6e03e6a204c86152a3e0d8f1246f2b99 (diff)
bug0049 :: output parameter GIVPTR stays uninitialized in case of quick exit (n.eq.0) or if the rank-1 modifier is small enough in *laed8 sent by Alexander Kobotov (Intel) on Wed, 3 Feb 2010 to lapack@cs.utk.edu
Diffstat (limited to 'SRC/claed8.f')
-rw-r--r--SRC/claed8.f8
1 files changed, 7 insertions, 1 deletions
diff --git a/SRC/claed8.f b/SRC/claed8.f
index 07cf124c..79ac0314 100644
--- a/SRC/claed8.f
+++ b/SRC/claed8.f
@@ -174,6 +174,13 @@
RETURN
END IF
*
+* Need to initialize GIVPTR to O here in case of quick exit
+* to prevent an unspecified code behavior (usually sigfault)
+* when IWORK array on entry to *stedc is not zeroed
+* (or at least some IWORK entries which used in *laed7 for GIVPTR).
+*
+ GIVPTR = 0
+*
* Quick return if possible
*
IF( N.EQ.0 )
@@ -241,7 +248,6 @@
* components of Z are zero in this new basis.
*
K = 0
- GIVPTR = 0
K2 = N + 1
DO 60 J = 1, N
IF( RHO*ABS( Z( J ) ).LE.TOL ) THEN