Re: [eigen] [PATCH] Toeplitz matrix specialization

[Benoit Jacob <jacob.benoit.1@xxxxxxxxx>]

Hi,

Thanks for your patch.
We are in the process of rethinking completely how we deal with
special matrices.
The first kind of such matrices in Eigen was diagonal matrices.
Then, Mauro came by with a patch for tridiagonal matrices, showing
that we need to rethink that completely.
We're going to decide for a general approach to this subject very
soon, and then it'll be a perfect time to add Toeplitz matrices.
But it's very good that you sent this now as it gives us one more
example to keep in mind when deciding for an approach for special
matrices.

Cheers,
Benoit

2009/1/27 Ricard Marxer Piñón <email@xxxxxxxxxxxxxxxx>:
> Here is a first draft of a Toeplitz Matrix specialization.
> ( http://en.wikipedia.org/wiki/Toeplitz_matrix )
>
> A couple of things:
>   - I added a few methods to MatrixBase to easily construct Toeplitz
> matrices the same way as Diagonal matrices (asToeplitz(), isToeplitz())
>   - I added a flag ConstantDiagsBit (maybe it would be better to call it
> ToeplitzBit, I don't know)
>
> Let me know if I'm going on the right direction and what I did wrong.
>
> Toeplitz matrices allow for specialized solver Levinson-Durbin recursion
> which exploits its structure.
>
> --
> ricard
> http://www.ricardmarxer.com
> http://www.caligraft.com
>