Re: [eigen] FFT for Eigen |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] FFT for Eigen
- From: Márton Danóczy <marton78@xxxxxxxxx>
- Date: Tue, 19 May 2009 14:28:29 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:in-reply-to:references :from:date:message-id:subject:to:content-type :content-transfer-encoding; bh=fliLcYSARXjRPn69ND0uZJ3RFlEl5j19OndJZvOxiNA=; b=jMGZSt0SXzxLt8FF8cy315mbOnixbjkFvQS7NiElDwBP0cNlGymeirfnS3IW9JjDeT lm6NINpYTu8zoETM1ApIgUFnTb1A4NLxO5f40/Ul6wsvg5NahGmE+6K5VQiVQTUtWdf4 Yn5K/YVDLsVsrSo9y5HBuVCm7nfZFUVpzWZng=
- 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=C7tSTn5HUo144Mwcr6ECAW6CUyvXM1lHDfA1Clg8lys04oqAAPvCiv4djfYpAC2gBC zR4pmGftQohIYLWLG/k5ZjQkkQLFLEHSNI8sGHQJmMqE+32GBYlB9GbqI5rssIqa4tvC 2ZA5VYiTpLop88nhpAA1NqDimcEOZX3mHA2UU=
> 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