Page 1 of 1

[DZSTEMR] NSPLIT unitialized; N=2 code wrong

PostPosted: Fri Jul 26, 2013 5:38 am
by lawrence mulholland
1. NSPLIT unitialized

line 740 in zstemr, line 726 in dstemr
IF( NSPLIT.GT.1 .OR. N.EQ.2 ) THEN

If n==2 then NSPLIT is unitialized

2. ZSTEMR
In n==2 code there are four instances of the same bug
eg line 527
ISUPPZ(2*M-1) = 1
ISUPPZ(2*M-1) = 2

should be
ISUPPZ(2*M-1) = 1
ISUPPZ(2*M) = 2