LAPACK Archives

[Lapack] Issue in *la_{po,sy,he}rfsx_extended.f

Julie,

there's one more issue in the testing suite at TESTING/LIN/cebchvxx.f:

line 106:
     $                   CGBSVXX, SLACPY, LSAMEN

should be
     $                   CGBSVXX, CLACPY, LSAMEN

and line 177:
         CALL DLACPY('ALL', KL+KU+1, N, AB, LDAB, ABCOPY, LDAB)

should be
         CALL CLACPY('ALL', KL+KU+1, N, AB, LDAB, ABCOPY, LDAB)

Best regards,

Michael.

From: julie langou [mailto:julie@Domain.Removed]
Sent: Friday, April 02, 2010 11:09 PM
To: Chuvelev, Michael
Cc: lapack@Domain.Removed; Henry, Greg; Kobotov, Alexander V
Subject: Re: [Lapack] Issue in *la_{po,sy,he}rfsx_extended.f

Thank you Michael for the bug report,
I just applied your fix in LAPACK repository  (revision 736).
Regards,
Julie
On Apr 2, 2010, at 6:32 AM, Chuvelev, Michael wrote:


Dear lapackers,

let me notify you of the issue in LAPACK 3.2.1 *la_{po,sy,he}rfsx_extended.f 
sources. For instanse, in sla_syrfsx_extended.f:

DY is described as
*     DY             (input) REAL array, dimension (N)
*     Workspace to hold the intermediate solution.

whereas N-by-NRHS array is referenced at DY due to the following statement 
(line 384):

            CALL SSYTRS( UPLO, N, NRHS, AF, LDAF, IPIV, DY, N, INFO )

This may lead to an outside memory corruption. I assume one RHS should be 
referenced here, because all the RHSs are considered individually in the loop 
by J. That is, substituting NRHS by 1:

            CALL SSYTRS( UPLO, N, 1, AF, LDAF, IPIV, DY, N, INFO )

should fix the issue. This is actual for any subroutine mentioned.

BTW, *la_{ge,gb}rfsx_extended.f source code is free from this issue - only one 
RHS is referenced in the corresponding solves.

Best regards,

Michael.


--------------------------------------------------------------------

Closed Joint Stock Company Intel A/O

Registered legal address: Krylatsky Hills Business Park,

17 Krylatskaya Str., Bldg 4, Moscow 121614,

Russian Federation



This e-mail and any attachments may contain confidential material for

the sole use of the intended recipient(s). Any review or distribution

by others is strictly prohibited. If you are not the intended

recipient, please contact the sender and delete all copies.
_______________________________________________
Lapack mailing list
Lapack@Domain.Removed<mailto:Lapack@Domain.Removed>
http://lists.eecs.utk.edu/mailman/listinfo/lapack

**********************************************
Julie Langou; Research Associate in Computer Science
Innovative Computing Laboratory;
University of Tennessee from Denver, Colorado ;-)
julie@Domain.Removed<mailto:julie@Domain.Removed>; http://www.cs.utk.edu/~julie/








--------------------------------------------------------------------
Closed Joint Stock Company Intel A/O
Registered legal address: Krylatsky Hills Business Park, 
17 Krylatskaya Str., Bldg 4, Moscow 121614, 
Russian Federation

This e-mail and any attachments may contain confidential material for
the sole use of the intended recipient(s). Any review or distribution
by others is strictly prohibited. If you are not the intended
recipient, please contact the sender and delete all copies.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100405/e2d7b89a/attachment.html
 

<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