LAPACK Archives

[Lapack] possible problem in pxstein (scalapack)

Hi!

Not sure if this is a real problem or not...
pxstein are doing the following

BN = ISPLIT( NBLK )
ONENRM = ABS( D( B1 ) ) + ABS( E( B1 ) )
ONENRM = MAX( ONENRM, ABS( D( BN ) )+ABS( E( BN-1 ) ) )

That last E(bn-1) will sometimes ref e(0)
and the following do loop will ref e(-1) when bn=1
DO 120 J = B1 + 1, BN - 1
ONENRM = MAX( ONENRM, ABS( D( J ) )+ABS( E( J-1 ) )+ ...

the xcsep test triggers this.

-- 
Ake Sandgren, HPC2N, Umea University, S-90187 Umea, Sweden
Internet: ake@Domain.Removed   Phone: +46 90 7866134 Fax: +46 90 7866126
Mobile: +46 70 7716134 WWW: http://www.hpc2n.umu.se


<Prev in Thread] Current Thread [Next in Thread>


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