Re: [eigen] Skyline matrix |

[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]

*To*: eigen@xxxxxxxxxxxxxxxxxxx*Subject*: Re: [eigen] Skyline matrix*From*: guillaume saupin <guillaume.saupin@xxxxxx>*Date*: Tue, 27 Oct 2009 10:40:19 +0100

Gael Guennebaud a écrit :

this is not an easy question because we are still unsure how to manageefficiently and with minimal code all these kind of "special" matrices.Of course an easy solution would be to make it inheritSparseMatrixBase and let the generic algorithms based on sparseiterators do the job. However, this approach won't be optimal becausethe iterators won't take advantage of the specificity of aSkylineMatrix which, e.g., should allow vectorization.Actually skyline matrices are very similar to banded matrices, and inparticular it seems to me that they are similar enough to share thesame algorithms. So for instance these two kind of matrices couldinherit the same matrix base and the algorithms could be implementedin a generic way via the concept of "range inner vector" which is acolumn (or a row) vector with an associated start index... We alsoneed efficient basic operators dealing with such vector. E.g. theaddition of two such vector can be efficiently done by first copyingthe non overlapping part (start and end), and then sum the overlappingpart using Eigen's core module feature.

guillaume

These are just some initial thoughts and the discussion is very open! gaelOn Mon, Oct 19, 2009 at 12:55 PM, guillaume saupin<guillaume.saupin@xxxxxx <mailto:guillaume.saupin@xxxxxx>> wrote:Hello, We are planning to use your library in our projects, but we need a skyline matrix. Therefore I'd like to implement one for eigen, but I don't now where to start. Is there a specific class that can be a good starting point / skeleton to start with ? The SparseMatrixBase might be a good choice. Should this SkylineMatrix inherit from SparseMatrixBase, or be a separate class ? Thanks, guillaume

**Follow-Ups**:**Re: [eigen] Skyline matrix***From:*Gael Guennebaud

**References**:**[eigen] Skyline matrix***From:*guillaume saupin

**Re: [eigen] Skyline matrix***From:*Gael Guennebaud

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] Mapping array of scalars into quaternions** - Next by Date:
**Re: [eigen] Skyline matrix** - Previous by thread:
**Re: [eigen] Skyline matrix** - Next by thread:
**Re: [eigen] Skyline matrix**

Mail converted by MHonArc 2.6.19+ | http://listengine.tuxfamily.org/ |