LAPACK Archives

### [Lapack] bug in dstemr.f when N = 2

 ```I think we have a bug in dstemr.f for the case when N is 2. The following block: 392: IF (CS.NE.ZERO) THEN 393: ISUPPZ(2*M-1) = 1 394: ISUPPZ(2*M-1) = 2 395: ELSE 396: ISUPPZ(2*M-1) = 1 397: ISUPPZ(2*M-1) = 1 398: END IF should really be (note lines 394, 397): 392: IF (CS.NE.ZERO) THEN 393: ISUPPZ(2*M-1) = 1 394: ISUPPZ(2*M) = 2 395: ELSE 396: ISUPPZ(2*M-1) = 1 397: ISUPPZ(2*M) = 1 398: END IF similarly the block: 415: IF (SN.NE.ZERO) THEN 416: IF (CS.NE.ZERO) THEN 417: ISUPPZ(2*M-1) = 1 418: ISUPPZ(2*M-1) = 2 419: ELSE 420: ISUPPZ(2*M-1) = 1 421: ISUPPZ(2*M-1) = 1 422: END IF 423: ELSE should really be (note lines 418, 421): 415: IF (SN.NE.ZERO) THEN 416: IF (CS.NE.ZERO) THEN 417: ISUPPZ(2*M-1) = 1 418: ISUPPZ(2*M) = 2 419: ELSE 420: ISUPPZ(2*M-1) = 1 421: ISUPPZ(2*M) = 1 422: END IF 423: ELSE -- regards Ashutosh Mahajan http://coral.ie.lehigh.edu/~asm4 ```
 Current Thread [Lapack] bug in dstemr.f when N = 2, Ashutosh Mahajan <= [Lapack] bug in dstemr.f when N = 2, Julien Langou

For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or