a very minor correction, consider
?gejsv(joba, jobu, jobv, jobr, jobt, jobp, m, n, a, lda, sva, u, ldu, v, ldv, work, lwork,iwork, info)
in description of work the manual indicates that the size of vector "work" is "max(4, m+n)"
but for description of "lwork" which is the length of work it gives lwork as "max(6,m+n)"
All the best,