?(sy/he)tri2 function could pass NB>N to ?(sy/he)tri2x

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

?(sy/he)tri2 function could pass NB>N to ?(sy/he)tri2x

Postby nmozarto » Fri Mar 04, 2011 1:51 am

Problem in ?(sy/he)tri2 was found.

For matrices with small N subroutine pass parameter NBMAX to ?(sy/he)tri2x which could be large than N.
In test for this functionality this is hide. Local test subroutine ilaenv.f pass proper values of NB. But if use ordinary ilaenv we will get mistake.

If we implemented

IF ( NBMAX .GE. N ) THEN
MINSIZE = N
ELSE
MINSIZE = (N+NBMAX+1)*(NBMAX+3)
END IF

and will call

IF( NBMAX .GE. N ) THEN
CALL SSYTRI( UPLO, N, A, LDA, IPIV, WORK, INFO )
ELSE
CALL SSYTRI2X( UPLO, N, A, LDA, IPIV, WORK, NBMAX, INFO )
END IF

the problem will be solved.
nmozarto
 
Posts: 5
Joined: Thu Jan 20, 2011 1:05 am
Location: Intel

Re: ?(sy/he)tri2 function could pass NB>N to ?(sy/he)tri2x

Postby admin » Fri Mar 04, 2011 4:49 pm

Indeed!
Thanks for reporting the problem and submitting a fix.
This has been corrected in the LAPACK repository and added in the 3.3.0 Errata (http://www.netlib.org/lapack/Errata/index2.html)
Julie
admin
Site Admin
 
Posts: 501
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