|Re: [eigen] Iterators with dense matrices 2|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] Iterators with dense matrices 2
- From: Benoit Jacob <jacob.benoit.1@xxxxxxxxx>
- Date: Sun, 21 Feb 2010 19:27:29 -0500
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:in-reply-to:references :date:message-id:subject:from:to:content-type; bh=0GH1VUvXB58weWxbD2F4YInNyNblj3kEdAH3grE06RE=; b=tnq7F3Uhx/s7WRzf56nx655qLxK9f/Rfg1ZYKTRkxRui21icszuW87NCZprmTm1wom B90jc+NNauu7C28UI7n4pPT4CrbwD/lfkDJCDW1USPAMFWqozp2aOnzDAJEHSbdWtLU7 O6H2aVUmYv3z+HhhA2quZNMtgt8evykWVEEXU=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type; b=ZFa9PESQEH3BVD3ZsbHS1U1KSWzOHfRtuq5K1Y9O886Wpqe5ZHx53QKnDsKVybh+cz FpiN9tCNV9Kh/AZHHx+fVxSJyL6eSZ/sdum01vDwTM50Q8yL+tdImFCPZl4dZ8OpIeZW PY2XBiI5mttuTqwZINXW8bzlGxS+oXWVkLYkU=
To summarize the current situation:
- we agree on the potential usefulness of iterators and would be OK
to have them in Eigen.
- we already renamed the old end() function to tail() so that the
name "end()" is now available to be used for iterator stuff (and
begin() is available too).
- but we are currently focusing on finishing Eigen 3.0 and this is
not a priority for us as this is a feature addition that can be made
at a later date.
- patches welcome :) but anyone trying to write such a patch would
have to pay attention to consistency with the Sparse module which
already has iterators.
> Still, iterator support might be helpful to re-use algorithms that take
> iterators as input. But there is a much easier solution to achieve that
> (at least that's how I solved this problem): just use both interfaces in
> parallel. For example, initialise your matrix as uBlas matrix and then
> initialise an Eigen2 matrix using exactly the same block of memory
> occupied by the original uBlas matrix.
Indeed, you can do that by mapping the memory buffer using Map<MatrixType>.