I'm trying to do an SVD of a tridiagonal matrix by using first dgbbrd
(bidiagonalization) and then dbdsdc (divide-and-conquer SVD). However,
I'm using clapack and can't get it to work. Probably I haven't done the
indexing correctly. I have tried to look at the comments and code in the
source files, but this is just confusing to me. (Looks like fortran indexing...)
Does anyone have any experience with this, and know how to properly index banded matrices?