Re: [eigen] A complex FFT for Eigen |

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

*To*: eigen@xxxxxxxxxxxxxxxxxxx*Subject*: Re: [eigen] A complex FFT for Eigen*From*: Tim Molteno <tim@xxxxxxxxxxxxxxxxxxx>*Date*: Fri, 28 Nov 2008 09:22:32 +1300

Benoit Jacob wrote:

- your code doesn't leverage Eigen as much as it could to enable high performance. For example, in apply(), replacing the for loop by single operations on vectors would let Eigen handle loop unrolling and vectorization. Well, except that Eigen doesn't vectorize std::complex currently, so I wonder if your code should be adapted to that or if this means that we should consider vectorizing std::complex. (So far we haven't seen many uses for std::complex anyway so we didn't optimize it as much as we maybe could).

data[1:N] += data[N+1:2N]

- the factorize() function here is not vectorizable at all because it doesn't do any contiguous access. Aren't there other, more hardware-friendly algorithms?

Cheers Tim ---

**Follow-Ups**:**Re: [eigen] A complex FFT for Eigen***From:*Benoit Jacob

**References**:**[eigen] A complex FFT for Eigen***From:*Tim Molteno

**Re: [eigen] A complex FFT for Eigen***From:*Benoit Jacob

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] projects using Eigen2 ?** - Next by Date:
**Re: [eigen] A complex FFT for Eigen** - Previous by thread:
**Re: [eigen] A complex FFT for Eigen** - Next by thread:
**Re: [eigen] A complex FFT for Eigen**

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