I am using the default FFT in eigen:

Code: Select all

#include <unsupported/Eigen/FFT>
using namespace Eigen;
FFT<double> fft;
fft.SetFlag(fft.HalfSpectrum);
fft.fwd(*A,*B);


This works when A is defined as std::vector<complex<double> >.
but doesn't work if I define A as vector of eigen: VectorXcd
What do I have to change?
thanks

From a quick look at the source code, my guess would be:

Code: Select all

fft.fwd(A,B)

jitseniesen wrote:From a quick look at the source code, my guess would be:

Code: Select all

fft.fwd(A,B)

if the type is VectorXcd * then the right code is

Code: Select all

fft.fwd((*A),(*B))

and not

Code: Select all

fft.fwd(*A,*B)