Re: [eigen] 3 dimensional arrays

[Mani chandra <mchandra@xxxxxxxxxx>]

 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.

Mani

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?
> >
> > 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