Bug in the new ?syswapr function

Post here if you want to report a bug to the LAPACK team

Bug in the new ?syswapr function

Postby nmozarto » Mon Jan 24, 2011 3:28 am

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
nmozarto
 
Posts: 5
Joined: Thu Jan 20, 2011 1:05 am
Location: Intel

Re: Bug in the new ?syswapr function

Postby admin » Mon Jan 24, 2011 4:24 pm

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
admin
Site Admin
 
Posts: 486
Joined: Wed Dec 08, 2004 7:07 pm

Re: Bug in the new ?syswapr function

Postby admin » Wed Jan 26, 2011 11:01 am

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
admin
Site Admin
 
Posts: 486
Joined: Wed Dec 08, 2004 7:07 pm


Return to Bug report

Who is online

Users browsing this forum: No registered users and 1 guest