0
__global__ void HYPER (int tFast, int tFastLenth, int kilo, int lenPrzvFast, double eps, int AF,double *arrINTLighFast, int *arrPrzvFncFst, int dv_ptr)
{
for(int j = 0;j<(tFast*tFastLenth);j++)
{ arrINTLighFast[j]=0;
}
for(int j = 0;j<(kilo);j++) arrPrzvFncFst[j]=0;
for(int j = 1;j<(tFast*tFastLenth);j++)
{ arrINTLighFast[j]= arrINTLighFast[j-1] + AF*exp(-j/(eps+tFast)); }
for(int j = 0;j<(tFast*tFastLenth-1);j++)
{
for(int i=(arrINTLighFast[j]);i< (arrINTLighFast[j+1]);i++)
{arrPrzvFncFst[i]=j;}
}
for(int j = 0;j<lenPrzvFast;j++)
{ devvecPrzvFncFst61440Slow149998[j]= arrPrzvFncFst[j] ;}
}
int main (void)
{
const int tFast = 9;
const int tFastLenth = 6;
double arrINTLighFast[tFast*tFastLenth];
int arrPrzvFncFst[61500];
int AF = 1000;
int kilo = 1024;
int kilo150 = 149998;
const double eps=0.0000001;
const int lenPrzvFast=61500;
thrust::host_vector<int> vecPrzvFncFst61440Slow149998;
int Len_vecPrzv=(lenPrzvFast+kilo150);
for (int j=0;j<Len_vecPrzv;j++) vecPrzvFncFst61440Slow149998.push_back(0);
for (int j=0;j<Len_vecPrzv;j++) vecPrzvFncFst61440Slow149998 [j] = 0;
thrust::device_vector<int> devvecPrzvFncFst61440Slow149998 = vecPrzvFncFst61440Slow149998;
int *dv_ptr = thrust::raw_pointer_cast(devvecPrzvFncFst61440Slow149998.data());
HYPER <<<blocks, threads>>>(tFast, tFastLenth, kilo, lenPrzvFast, eps, AF, arrINTLighFast, arrPrzvFncFst, dv_ptr);
thrust::host_vector<int> HostvecPrzvFncFst61440Slow149998 = devvecPrzvFncFst61440Slow149998;
std::cout << "Device vector is: " << std::endl;
for(int j = 0; j<vecPrzvFncFst61440Slow149998.size(); j++)
std::cout << "vecPrzvFncFst61440Slow149998[" << j << "] = " << HostvecPrzvFncFst61440Slow149998[j] << std::endl;
return 0;
}
還有就是我不能在函數中使用向量問題發揮作用,所以我決定用推力:: raw_pointer_cast。然而,我有問題:在編譯期間,我有一個錯誤:標識符「devvecPrzvFncFst61440Slow149998」未定義。第二個是我不能明確地找出如何將int * dv_ptr傳遞給函數和原型,這裏有一個錯誤:類型爲「int *」的參數與類型爲「int」的參數不兼容。我看着在互聯網中,但沒有解決如何順利解決我上面怎麼投推力:: device_vector <int>與原始指針
謝謝您的時間