用推力执行傅里叶变换

时间:2012-12-02 11:52:04

标签: cuda fft thrust

Thrust是开始编程CUDA的惊人包装器。 我想知道有什么东西可以用推力封装NVIDIA CUFFT,或者我们需要实现自己吗?

1 个答案:

答案 0 :(得分:6)

这是一个非常晚的答案,只是为了从未答复的清单中删除这个问题。

使用带有推力的cuFFT应该非常简单,唯一要做的就是将thrust::device_vector强制转换为原始指针。下面报告了一个非常简单的例子:

#include <iostream>
#include <cufft.h>
#include <stdlib.h>
#include <thrust/host_vector.h>
#include <thrust/device_vector.h>
#include <thrust/generate.h>
#include <thrust/transform.h>

int main(void){

    int N=4;

    // --- Setting up input device vector    
    thrust::device_vector<cuFloatComplex> d_in(N,make_cuComplex(1.f,2.f)), d_out(N);

    cufftHandle plan;
    cufftPlan1d(&plan, N, CUFFT_C2C, 1);

    cufftExecC2C(plan, thrust::raw_pointer_cast(d_in.data()), thrust::raw_pointer_cast(d_out.data()), CUFFT_FORWARD);

    // --- Setting up output host vector    
    thrust::host_vector<cuFloatComplex> h_out(d_out);

    for (int i=0; i<N; i++) printf("Element #%i; Real part = %f; Imaginary part: %f\n",i,h_out[i].x,h_out[i].y);

    getchar();
}