FFT example

fft-example.cpp

#include <vexcl/vexcl.hpp>

int main() {
    const size_t n = 1024;

    vex::Context ctx(
            vex::Filter::Env             &&
            vex::Filter::DoublePrecision &&
            vex::Filter::Count(1)
            );
    std::cout << ctx << std::endl;

    vex::FFT<double, cl_double2> fft(ctx, n);                     // Takes double, returns double2 (complex)
    vex::FFT<cl_double2, double> ifft(ctx, n, vex::fft::inverse); // Takes double2, returns double

    vex::vector<double> rhs(ctx, n), u(ctx, n), K(ctx, n);
    rhs = 0;

    u = ifft( K * fft(rhs) );
}

Makefile

fft-example: fft-example.cpp
	g++ -o $@ $^ -std=c++11 -I${HOME}/work/vexcl -lOpenCL -lboost_system