Boost logo

Ublas :

Subject: [ublas] Help with new matrix adaptor
From: Arnaldo Leon (arnaldo.leon_at_[hidden])
Date: 2009-02-07 14:54:53


I've been wanting to implement a read-only matrix adaptor for lagged
matrices, which are similar to Hankel matrices, but without the
square-matrix constraint, that is:

a(i,j) = a(i+1,j-1)

This would allow me to store the data in memory in a very compact
form. In my current implementation (where I copy the data) a matrix
with 100 lags consumes 100x more memory than necessary.

I've looked at symmetric_adaptor and I think it could be used as an
example of what I want to do, but the documentation is extremely

What do the find1() and find2() methods do? What is the "rank" variable?

Is there an easier way to set up const_iterator1 and const_iterator2?
Is there any sample code or tutorials somewhere? What is index1() and

Has anybody else implemented such an adaptor?

Is my approach correct? or should I be using some other template?

The only operations I will need to be using on this matrix is prod()
and column(), so I only need read-only access.

Thanks in advance...