Skip to main content
7 events
when toggle format what by license comment
May 23, 2017 at 12:35 history edited CommunityBot
replaced http://stackoverflow.com/ with https://stackoverflow.com/
Jun 12, 2013 at 14:13 comment added Max Hutchinson I don't know of any. You should probably just explicitly populate it. If the matrix is banded and symmetric (or Hermitian), then you can use BLAS (netlib.org/blas/dsbmv.f of zhbmv.f) and LAPACK (netlib.org/lapack/complex16/dsbev.f or zhbev.f).
Jun 12, 2013 at 11:44 comment added lorniper however, in order to save memory, I want to use banded storage, do you know any routine that is able to format a sparse matrix A in banded form?
Jun 11, 2013 at 21:43 comment added Max Hutchinson With ARPACK, you don't actually pass a matrix to the library. Instead, it uses a "reverse communication interface" to give you a vector, x, and ask you to compute y = Ax for it. In a sense, your 'matrix' is really a linear operator that you can implement however you like. Ref: caam.rice.edu/software/ARPACK/UG/node9.html
Jun 11, 2013 at 20:21 comment added lorniper what did you mean by "store the data however you like". the problem is essentially two steps, 1) how to declare A as a sparse matrix, 2) how to solve Ax=lBy?
Jun 11, 2013 at 18:57 vote accept lorniper
Jun 11, 2013 at 16:16 history answered Max Hutchinson CC BY-SA 3.0