?gejsv accesses V when JOBV='n'

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

?gejsv accesses V when JOBV='n'

Postby dbakshee » Thu Feb 04, 2016 1:33 am

call cgejsv (joba='C', jobu='U', jobv='N', jobr='N', jobt='T', jobp='P', m=5, n=5, ...)

The function attempts to access array V and produces segmentation fault.
dbakshee
 
Posts: 9
Joined: Fri Sep 18, 2015 5:44 am

Re: ?gejsv accesses V when JOBV='n'

Postby dbakshee » Thu Feb 04, 2016 2:46 am

This is probably not a bug, but documentation is confusing. The descriptions of U and V say:
If JOBU = 'N' U is not referenced.
If JOBV = 'N' V is not referenced.

However, the description of JOBT says for JOBT='T":
the caller should provide both U and V, as one of the matrices is used as workspace if the matrix A is transposed

Should the documentation be improved? Perhaps thus:
If JOBU = 'N' U is not referenced, unless JOBT='T'.
If JOBV = 'N' V is not referenced, unless JOBT='T'.
dbakshee
 
Posts: 9
Joined: Fri Sep 18, 2015 5:44 am

Re: ?gejsv accesses V when JOBV='n'

Postby admin » Thu Feb 04, 2016 11:02 pm

thank you Dimitri.
I added your suggestion to ?gejsv in our SVN repository - rev 1686
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