LAPACK_?syconv interface problem

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

LAPACK_?syconv interface problem

Postby Victor_K » Thu Jul 09, 2015 12:41 am

FORTRAN interface
SUBROUTINE DSYCONV( UPLO, WAY, N, A, LDA, IPIV, WORK, INFO )
has argument WORK which contains some useful information on exit:

DSYCONV convert A given by TRF into L and D and vice-versa.
Get Non-diag elements of D (returned in workspace) and
apply or reverse permutation done in TRF.

Look at
lapack_int LAPACKE_dsyconv( int matrix_layout, char uplo, char way, lapack_int n,
double* a, lapack_int lda, const lapack_int* ipiv )

The argument work is omitted and there is no way to get non-diag elements of D.


This is a bug because work array is used to pass off-diagonal elements of matrix D:
Victor_K
 
Posts: 5
Joined: Sat Nov 27, 2010 9:48 am

Re: LAPACK_?syconv interface problem

Postby admin » Thu Jul 09, 2015 10:58 pm

Thank you Victor,
This bug has been listed as bug 129 in http://www.netlib.org/lapack/Errata/index2.html
Julie
admin
Site Admin
 
Posts: 609
Joined: Wed Dec 08, 2004 7:07 pm


Return to Bug report

Who is online

Users browsing this forum: No registered users and 3 guests

cron