LAPACK/ScaLAPACK Development

Bug in new ?syswapr function was found.
The LDA is not passed to ?SYSWAPR (from ?SYTRI2X suboutine) where A is A(N,N), this causes trashing of the matrix.

For example :
SUBROUTINE SSYSWAPR( UPLO, N, A, I1, I2)
* .. Scalar Arguments ..
CHARACTER UPLO
INTEGER I1, I2, N
* .. Array Arguments ..
REAL A(N,N)

The problem appears when we use matrix A(LDA,N). The case is not checked by testsuite, where in ?DRVSY always LDA = MAX( 1, N ).
Please pay your attention on this problem

Hi,
Indeed, we need to modify and add this testing.
This is listed as bug0067 in LAPACK Errata: http://www.netlib.org/lapack/Errata/index2.html
Thank you
Julie

Hi again,
we just have corrected the problem in our SVN repository (rev 893 and 894).
our SVN repository is open for read-only

Code: Select all

svn co http://icl.cs.utk.edu/svn/lapack-dev/lapack/trunk


This has not yet been extensively tested.
Please let us know if this fixes things for you.
Thank again for reporting the problem
Julie