Page 1 of 1
?gejsv accesses V when JOBV='n'
Posted:
Thu Feb 04, 2016 1:33 amby dbakshee
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.
Re: ?gejsv accesses V when JOBV='n'
Posted:
Thu Feb 04, 2016 2:46 amby dbakshee
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'.
Re: ?gejsv accesses V when JOBV='n'
Posted:
Thu Feb 04, 2016 11:02 pmby admin
thank you Dimitri.
I added your suggestion to ?gejsv in our SVN repository - rev 1686
Julie