Я пишу на CUDA c++ 10.0. У меня есть массив на девайсе, который мне нужно регулярно заполнять новыми случайными числами. Для этого я написал kernel, с использованием параллелизации и прочего. Для генерации случайных чисел я решил использовать библиотеку thrust::random
.
Допустим, мне надо заполнить массив из 100 чисел, и я это разбиваю на условные 10 блоков по 10 потоков. Получается вот такой кернель:
__global__ void gen(float* a, thrust::minstd_rand rng)
{
int index = threadIdx.x + blockIdx.x * blockDim.x;
thrust::uniform_real_distribution<float> dist(-10,10);
a[index] = dist(gen);
}
Сам вызов функции выглядит так:
int main()
{
// инициализация devA на девайсе
// ...
thrust::minstd_rand rgen(time(0)); // нужно каждый раз разное, поэтому инициализируем временем системы
gen<<<bloks, threads>>>(devA, rgen); // собственно вызов
// ...
}
Сначала вопрос, где размещаются переменные, объявленные в кернеле, сразу на девайсе? thrust::uniform_real_distribution<float> dist
и int index
уже на девайсе, или как то там внутри бегают между CPU/GPU?
И, собственно, сам вопрос. Данная функция уже написана неправильно, так каждый вызов этой функции копирует себе состояние rgen, соответственно получаем много повторяющихся значений. Так же не уверен, что создание dist
внутри кернеля целесообразно. Как вариант избежать копирований - rng
передавать в виде ссылки на этот самый генератор, но тогда он должен находиться на девайсе, а еще, скорее всего, появится race condition, когда все потоки обратятся к этому генератору чуть ли не одновременно. Помогите разобраться, что с этим можно сделать? Или может быть я вообще не в ту сторону это делаю?
Скажу уточню, в итоге мне нужно генерировать огромные массивы случайных чисел, постоянно в ходе жизненного цикла программы, поэтому избавляться от параллелизма куды и запускать все в единственном блоке\потоке - очень нежелательно.