Re: [eigen] banded matrices in Eigen |
[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]
LauraThank you,Hi Rhys,Thinking of another approach, could somebody tell if it is possible to use the ThreadedConjugateGradient solver with a custom matrix (that is to say, a matrix class defined by myself)? Or maybe with a matrix which is derived from the Diagonal class?
Yes, I meant a sparse matrix with that kind of structure (symmetric and multi-diagonal). I was using a SparseMatrix, but maybe it is better to use something more similar to the Diagonal class.
2014-02-07 4:06 GMT+01:00 Rhys Ulerich <rhys.ulerich@xxxxxxxxx>:
Hi Laura,
Just to be clear, do you mean dense, symmetric, banded a la BLAS SBMV?
> I would also be interested in sparse multi-diagonal symmetric matrices.. I am
> dealing with sparse matrices that are composed only of a certain number of
> diagonals, and, since they are symmetric, some of this diagonals are
> identical.
Or do you mean a genuinely sparse matrix (so compressed storage of
some sort) that happens to exhibit the same structure?
I'm not aware of any but would eagerly hear news. Banded matrices are
> Has there been by any chance some improvement since the last question
> about this topic?
the ugly stepchild of linear algebra libraries.
- Rhys
Mail converted by MHonArc 2.6.19+ | http://listengine.tuxfamily.org/ |