Problem with ?SYTRS2 subroutine

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

Problem with ?SYTRS2 subroutine

Postby nmozarto » Thu Jan 27, 2011 5:47 am

Problem in new function ?SYTRS2 was found: wrong answer is got for a matrix A(LDA,N) and for B(LDB,N). If we take LDA=N, LDB=N then everything is OK.

The bug is in calling ?TRSM and ?SCAL functions:
CALL ?TRSM('L','L','T','U',N,NRHS,ONE,A,N,B,N)
CALL ?SCAL( NRHS, ONE / A( I, I ), B( I, 1 ), N )
The calling should be
CALL ?TRSM('L','L','T','U',N,NRHS,ONE,A,LDA,B,N)
CALL ?SCAL( NRHS, ONE / A( I, I ), B( I, 1 ), LDB )
The problem vanishes after these corrections
nmozarto
 
Posts: 5
Joined: Thu Jan 20, 2011 1:05 am
Location: Intel

Re: Problem with ?SYTRS2 subroutine

Postby admin » Thu Jan 27, 2011 10:22 am

Thanks again, I just made the correction
admin
Site Admin
 
Posts: 490
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

cron