Re: [eigen] 3 dimensional arrays

[ Thread Index | Date Index | More Archives ]

We do have a GPU port of our code, where about 70% of the computation is done on the GPU. We're getting nice speed ups but we can do much better. But at the moment it runs only on one tesla card. We want to change to Eigen since all development on Blitz++ has stopped.


On 06/29/2010 06:03 PM, Rohit Garg wrote:
If it has a lot of array based computation aka element wise stuff, may
be should look into gpu's as well.

On Tue, Jun 29, 2010 at 4:58 PM, Mani chandra<mchandra@xxxxxxxxxx>  wrote:
  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?

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).


On 06/29/2010 01:12 PM, Hauke Heibel wrote:
On Tue, Jun 29, 2010 at 9:23 AM, Mani chandra<mchandra@xxxxxxxxxx>
  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+