CUDA не работает, как ожидалось?

Question

Я запрограммировал код CUDA.

unsigned long mask_buffer;
int s;
off_t p,

for(p=0;p!=5000;p++)
{
    for(s=start;s!=end;s++)
    {
        ref_off = *(((unsigned int*)(idx_base)) + p);

        if((int)(first_indexes[s-start_sequence] % 8 - ref_off % 8) < 0)
        {
            int shamt2 = (ref_off % 8 - first_indexes[s-start_sequence] % 8);
            mask_buffer = *((unsigned long *)(msk_base + (ref_off - first_indexes[s-start_sequence])/8)) >> shamt2;

            if( ( (*(unsigned long *)(seqmaskc + 16 * (s-start_sequence))) ^ mask_buffer ) << shamt2) 
                continue;
        }

        else if((int)(first_indexes[s-start_sequence] % 8 - ref_off % 8) == 0)
        {
            mask_buffer = *((unsigned long *)(msk_base + (ref_off)/8));

            if( (*(unsigned long *)(seqmaskc + 16 * (s-start_sequence)) ^ mask_buffer))
                continue;
        }

        else
        {
            int shamt2 = 8 - (first_indexes[s-start_sequence] % 8 - ref_off % 8);
            mask_buffer = *((unsigned long *)(msk_base + (ref_off/8- first_indexes[s-start_sequence]/8) - 1)) >> shamt2;

            if( ( (*(unsigned long *)(seqmaskc + 16 * (s-start_sequence))) ^ mask_buffer ) << shamt2) 
                continue;
        }

        int shamt = (ref_off % 4 - first_indexes[s-start_sequence] % 4) * 2;

        memcpy(reference_blk, ref_base + ref_off / 4 - first_indexes[s-start_sequence] / 4, sequence_bytes);

        for (rp = last_rp ; rp != (unsigned long *) reference_blk ; rp--) 
        {
            unsigned long tmp = ((*rp) & ((1 << shamt) - 1)) << (8 * sizeof(unsigned long) - shamt);
            *rp = (*rp >> shamt) | shifted_in;
            shifted_in = tmp;
        }

        *rp = (*rp >> shamt) | shifted_in;

        if (sequence_length & 0x3)
            reference_blk[sequence_length >> 2] &= (1 << ((sequence_length & 0x3) << 1)) - 1;

        for ( i = sequence_length >> 2 ; i & (SEQUENCE_ALIGN - 1) ; i++ )
            reference_blk[i] = 0;

        //-- instead of memcmp --//
        int v = 0;
        char *p1 = (char *)sequence;
        char *p2 = (char *)reference_blk;
        int tmp_asd = sequence_bytes;

        while(tmp_asd!=0)
        {
            v = *(p1++) - *(p2++);

            if(v!=0)
                break;

            tmp_asd--;
        }

        if(v == 0)
        {
            mat_count[s - (int)start_sequence]++;      /* Maintain count */
            mat_position[s - (int)start_sequence] = ref_off-first_indexes[s-start_sequence]; /* Record latest position */
        }

    }
}

Этот цикл for является основной функцией моего кода.Но проблема в том, что переменная «p» никогда не увеличивается больше 5 или 6. У меня GT530 на моем компьютере, и моя версия драйвера CUDA и версия времени выполнения также 4.0.В чем проблема в этом коде ???

Answer 1

Предполагая, что это действительно ядро, о котором говорилось в вашем исходном посте до того, как вы отредактировали новый код, вы не можете вызывать стандартные функции C-библиотеки, такие как memcpy внутри ядра CUDA (вы делаете это примерно на полпутичерез код).Вы можете вызывать только __device__ функции в ядре.Так что, если это не повторная реализация memcpy в CUDA, вызов этой функции внутри вашего ядра не будет работать ...

Кроме того, если вы собираетесь пытаться писать ядра CUDA, которые должнычтобы проверить результаты того, что написали другие потоки в группе, например, что вы делаете в разделе сравнения памяти вашего ядра, вам нужно будет поместить некоторые точки синхронизации в ваш код, чтобы убедиться, что все потоки достиглиопределенный момент, прежде чем начать проверять результаты из этих потоков.Если вы пытаетесь проверить, что другой блок мог записать в основную память, для этого также есть примитивы синхронизации, которые гарантируют, что результаты любого предыдущего блока потока будут записаны в основную память.