| 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: Rodrigo Benenson <rodrigo.benenson@xxxxxxxxx>
- Date: Tue, 25 Jan 2011 14:48:24 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:in-reply-to:references:from:date :message-id:subject:to:content-type:content-transfer-encoding; bh=F6Tg40WMc02NuAomGiXKaGAbZHOTIV274M0GVGErEsg=; b=kw8LkUh1N3WFRZHg8TOXOVslIWSGrLkTY9UIRyOKfrTlnx3CgqJSqb9E0O24HVsEbY VM24MVIu8HcZnNs3Req8INN801BnTlKUUB8nXAE0gfuNZLDGVbD5TqEfw1Czk4JLevjh huGQ/57I9nugUMPnHK64r5W/Ws5s0712CcQZQ=
- 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=CspyfYOWcjZYJuOBK8n8TgWvKvF3MuWMANtQwT0UTu3sYJ5B5h+XRNAdskhDDPvzRw OR9fSdS3zRhRwPmYVvCxSzTKU0yI0rElCRLibMF9edwV7JHcXRh69lglroF3ATIYMeAv 5F6ycJS9olATLG/w0Fo/dOEX+xET4eZAX9YQk=
It is there a recommended "Eigen like" library for 3d arrays manipulation ?
I am looking for something to could be compared performance-wise with Eigen....
Regards,
rodrigob.
On Tue, Jun 29, 2010 at 2:50 PM, Mani chandra <mchandra@xxxxxxxxxx> wrote:
> 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
>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>
>>>
>>>
>>>
>>
>>
>
>
>
>