Re: [eigen] 3 dimensional arrays

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


Nice. We really want to port our CFD code written using Blitz++ to Eigen and manipulating 3D arrays is essential. So I guess we wait then? With the 3D data structure that you are suggesting, can the arrays then be used just like 3D arrays?

Mani
On 06/29/2010 04:40 PM, Benoit Jacob wrote:
2010/6/29 Mani chandra<mchandra@xxxxxxxxxx>:
  Can one use Eigen to perform operations similar to Blitz++ on
multidimensional arrays?
Not currently. Eigen is definitely focused on the cases of 1D
(vectors) and 2D (matrices).

However, as Hauke mentions (but his example is inside out :) ), we
support using Array as Scalar type, e.g.

     Matrix<ArrayXd, Dynamic, Dynamic>

will give you a 3D data structure. Currently this isn't well
optimized, and we don't offer a tensor API, but both issues can be
overcome (the speed issue hopefully by the time 3.0 or 3.1 is
released, the tensor API would be a nice extension module).

Benoit

Mani
On 06/29/2010 01:12 PM, Hauke Heibel wrote:
On Tue, Jun 29, 2010 at 9:23 AM, Mani chandra<mchandra@xxxxxxxxxx>    wrote:
  Any plans to support multi-dimensional arrays in Eigen in the future?
Only to some extend. It is unlikely that you'll ever be seeing
constructors like

ArrayXXXd m(rows,cols,slices)

but we plan to support (maybe that is already partially working in the
development branch)  so called "arrays of structures", i.e.

Array<VectorXd, Rows, Cols>

Our core devs, Benoit&    Gael know more.

- Hauke












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