Several issues in LAPACK 3.4.2 and ScaLAPACK 2.0.2

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

Several issues in LAPACK 3.4.2 and ScaLAPACK 2.0.2

Postby ivan_latkin » Thu Sep 20, 2012 7:49 am


Please find my comments regarding issues in LAPACK 3.4.2 and ScaLAPACK 2.0.2 below:

1. LAPACKE interfaces:
• lapacke_?tprfb.c – ‘myldwork’ parameter is redundant;
• lapacke_?tprfb_work.c – ‘mywork’ and ‘myldwork’ should be renamed to ‘work’ and ‘ldwork’;
• lapacke_ctpqrt.c, lapacke_ctpqrt_work.c – wrong order of parameters (‘t,b,ldb,ldt‘ instead of ‘b,ldb,t,ldt’);
• lapacke_?tpqrt2.c – parameter ‘l’ between ‘n’ and ‘a’ is absent;
2. Tests:
• ?drvhex.f(617), ?drvsyx.f(623), c and z precisions: IWORK(N+1) was changed to RWORK, but RWORK( 2*NRHS+1 ) should be used here, otherwise BERR and FERR ( which are stored from RWORK(1) to RWORK( 2*NRHS )) are destroyed before check, resulting 10^308 residual errors.
• ?drvge.f(d,s –line 555, c,z – line 556), ?drvgb.f (651) – add check ‘.AND. INFO.LE.N’ to compute RPVGRW correctly
3. Comments and source files :
• ?sbevd.f(159): In description comment ‘The dimension of array LIWORK...’ – rename LIWORK with IWORK
• chesvxx.f, zhesvxx.f: Size of WORK array should be 5*N (not 2*N), see line with call of ?hetrf

• p?laqr5
- STOP statement is used in library function - it may crash program without any message. It should be changed to RETURN, i believe.

Best regards,
Ivan Latkin
Posts: 3
Joined: Thu Jan 12, 2012 1:28 am
Location: Intel Co.

Re: Several issues in LAPACK 3.4.2 and ScaLAPACK 2.0.2

Postby admin » Fri Sep 21, 2012 11:29 pm

Dear Ivan,
Thank you very much for reporting those issues.
We corrected them all in our repositories.
Site Admin
Posts: 612
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