Re: [eigen] FFT for Eigen

[ Thread Index | Date Index | More Archives ]

I can believe that this is probably a very efficient storage scheme.
We could offer this as an option if really it's not too hard to
implement (i didn't start thinking about this).

The default should remain the current for many reasons, but as an
option why not.


2009/5/19 Márton Danóczy <marton78@xxxxxxxxx>:
>> I concur: I don't think that it would be very useful to have complex
>> matrices with the real and imaginary parts stored separately. Most
>> operations -- and the more costly ones -- would run slower in such a
>> scheme. The basic issue here is memory locality.
> What about storing them packet by packet? That is, in case of floats,
> four real parts followed by four imaginary parts. That would not be
> too hard to implement and vectorization of component-wise operations
> would be trivial. And I think even FFTW can handle that using the guru
> interface, by setting up a split fft plan with a stride of
> 2*packetsize.
> Marton

Mail converted by MHonArc 2.6.19+