Correct a bug found by Ashutosh Mahajan from Mathematics and Computer Science Divisio...
authorlangou <langou@users.noreply.github.com>
Wed, 16 Dec 2009 01:51:40 +0000 (01:51 +0000)
committerlangou <langou@users.noreply.github.com>
Wed, 16 Dec 2009 01:51:40 +0000 (01:51 +0000)
commitecf41ea0c0ac83c20b0f4268bc4ca1e4b3ff9663
treec567032b88624794141dc9d81e75d1965bbfc0ae
parentbe91663bbc3e38d7debd10ea3dda660f85ff69fe
Correct a bug found by Ashutosh Mahajan from Mathematics and Computer Science Division of Argonne National Labs.
Email sent on Fri, 11 Dec 2009 14:51:07 to lapack@cs.utk.edu.

Subject: [Lapack] bug in dsyevr.f when N = 1

when dsyevr is called with N=1, then it returns without setting up ISUPPZ values.

331:       IF( N.EQ.1 ) THEN
...
342:          IF( WANTZ )
343:      $      Z( 1, 1 ) = ONE
344:          RETURN
345:       END IF

It should rather do

IF( WANTZ )
    Z( 1, 1 ) = ONE
    ISUPPZ(1) = 1
    ISUPPZ(2) = 1
END IF
SRC/cheevr.f
SRC/dsyevr.f
SRC/ssyevr.f
SRC/zheevr.f