Re: [eigen] 3 dimensional arrays |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] 3 dimensional arrays
- From: Rohit Garg <rpg.314@xxxxxxxxx>
- Date: Tue, 29 Jun 2010 18:03:34 +0530
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:mime-version:received:in-reply-to :references:from:date:message-id:subject:to:content-type :content-transfer-encoding; bh=42pPIws/wJAxHNzHpJh+fhlIHt6/Az7yIwx79B3372I=; b=AqxfGKdLe+1HDkOx78s5s3tkdea/QnCn2Pck2lWh1Lb8jLdySoDIH4PF/FcOwyjmZX LEcBCoDKkWIz2DPrJ2LMytQdjAJjuWYUQ7P7X11aelqExPeeUpSDw4GH/Zkz18mWlZLk N+qqwSJERNiLT2aYhbWgFq85biOmQO1xpnO1s=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :content-type:content-transfer-encoding; b=ajPoLql2lyf2AtkjJ5cetOtw6ft0GZmLN+MkO/vArq/S9f2tddhLNQXC60SVYFdeKK 3JTYSPpC5RgJBv7c9P+hpj8Kf31eYhhA50JvjbTTIdyEmsZqKPQdN/WmEG5atN6Rc4qP 4NvbnRbB2GbKeZCFwbihZhAclYSdnEzdd06+I=
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
>>>>
>>>>
>>>>
>>>
>>>
>>>
>>
>>
>
>
>
>
--
Rohit Garg
http://rpg-314.blogspot.com/
Senior Undergraduate
Department of Physics
Indian Institute of Technology
Bombay