|Re: [eigen] A complex FFT for Eigen|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] A complex FFT for Eigen
- From: "Benoit Jacob" <jacob.benoit.1@xxxxxxxxx>
- Date: Sun, 30 Nov 2008 17:14:18 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:to :subject:in-reply-to:mime-version:content-type :content-transfer-encoding:content-disposition:references; bh=5S9usDD5ki+5KlvFZ5+MGJQRVYyk2xSGZVjVAaQMeHI=; b=ljDefb5Lx6kD/8nuWYPe1bGNzXyEBghTK55X+rAs7SX0o0lSAZuJtVjBGM7/zTrMpf xx8xFu7gzL4EuJLGHP3zkPkA1jS+4FzlKnqr6LZAO2O2ehUzUA7yOS/avOnyDM1F2gYS 29lN+Lckh7Mod+jfpbhq6vPkv+ke3uaC0UTdY=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:in-reply-to:mime-version :content-type:content-transfer-encoding:content-disposition :references; b=oCueQTXm0w1gSoPwDO2k1ddg2wvnmaFXTmklkFIcMtcyjs+4d9BKuFgh1PjaQRbCkd 6FjNDJuQQ+pkFZ1kCE/qn0KIDqHg1B/RWe/6JizBUseEqMy8OysGzUv27L46f8QvJQ19 TeFQDTkQZoIt+XW4l6V5VS/Drzyf+XbAQv3KE=
I think we now know where we are going, so we can start adding an
experimental FFT module to Eigen.
Let's summarize the goals:
- generalization to large/dynamic sizes, with an optimized fixed-side kernel
- allow optional FFTW backend. This must be thought of early in the
process as it impacts the API
- replace for loops by Eigen expressions as much as possible
Tim, are you willing to do this yourself? You could use an existing
Eigen module, like Regression or LU, as a model. We could help you
with any Eigen-specific aspect. One important thing is that the API,
the naming of variables, the coding style are uniform with the rest of
Eigen. If you do it yourself, send your first patch to this list
(generated by svn diff) and then you'll be able to apply for your own
KDE SVN account. I'd like the above 3 goals to be met before I merge
this new module.
Otherwise, i'll do it when I have time, but that may be not soon as
there are other things to do.
2008/11/30 Tim Molteno <tim@xxxxxxxxxxxxxxxxxxx>:
> On Sun, 30 Nov 2008, Benoit Jacob wrote:
>> Hehe yeah, it's a very strange aspect of C++ syntax that we too
>> struggled with when we first encountered it...
>> it the faulty line
>> the problem is that data itself is of a templated type, so when you
>> call the templated method segment<int>(int) on it, you need to add the
>> "template" keyword, like this:
>> decimation.apply(data.template segment<N>(0),-1);
>> Otherwise the C++ compiler doesn't understand that the < sign means
>> template parameter, and tries to interprete it as operator<.
> Excellent! and thank you for your help. I have it working now and have
> attached the most recent version to this email. I have a small standalone
> testharness here as well.