Thrust: Operator «*» не поддерживается

Question

я следующую функцию для заполнения вектора т с шагом от -time/2 времени/2 и размер шага DT:

#define THRUST_PREC thrust::complex<double> 
__host__ void generate_time(thrust::device_vector<THRUST_PREC> *t, const double dt, const double time) 
{ 
    THRUST_PREC start = -time/2.0; 
    THRUST_PREC step = dt; 
    thrust::sequence((*t).begin(), (*t).end(), start, step); 
} 

При компиляции, я получаю error : no operator "*" matches these operands , Зачем? Есть ли способ заполнить вектор так, как я, или я должен заполнить его по-старому (aka loop)?

Edit: Полная ошибка: Error 1 error : no operator "*" matches these operands C:\Program Files\NVIDIA GPU Computing Toolkit\CUDA\v7.5\include\thrust\system\detail\generic\sequence.inl 48 1 FFT_test

Answer 1

Это похоже на ошибку в thrust::complex. Операция умножения между const thrust:complex<double> и signed long не определена.

/usr/local/cuda/bin/../targets/x86_64-linux/include/thrust/system/detail/generic/sequence.inl(48): 
error: no operator "*" matches these operands 
      operand types are: const thrust::complex<double> * signed long 
      detected during: 
      .... 

Но как ни странно, вы могли бы использовать вместо thrust::transform. Следующий код работает.

#define THRUST_PREC thrust::complex<double> 
__host__ void generate_time(thrust::device_vector<THRUST_PREC> *t, const double dt, const double time) 
{ 
    THRUST_PREC start = -time/2.0; 
    THRUST_PREC step = dt; 
    thrust::transform(thrust::make_counting_iterator(0), 
        thrust::make_counting_iterator(0) + t->size(), 
        t->begin(), 
        start + step * _1); 
} 

В любом случае, внутренняя реализация использует индекс (типа signed long в thrust::sequence), чтобы вычислить требуемую последовательность с выражением

start + step * index; 

Дело предотвращает thrust::sequence от работы является то, что operator *(...) является не очень перегружен.

thrust::complex<double> a(1,1); 
double double_b = 4; 
float float_b = 4; 
int int_b = 4; 
long long_b = 4; 

a *= double_b; // ok 
a *= float_b; // ok 
a *= int_b; // ok 
a *= long_b; // ok 

std::cout << a * double_b << std::endl; // ok 
std::cout << a * float_b << std::endl; // error 
std::cout << a * int_b << std::endl; // error 
std::cout << a * long_b << std::endl; // error