C-сегментная ошибка сегментации - PullRequest
Новая
       13

C-сегментная ошибка сегментации

2 голосов
/ 13 августа 2011

, поэтому я пытался создать эмулятор GPGPU с c & pthreads, но столкнулся с довольно странной проблемой, которую я не знаю, почему она возникла. Код как ниже: 

#include <stdlib.h>
#include <stdio.h>
#include <pthread.h>
#include <assert.h>

// simplifies malloc
#define MALLOC(a) (a *)malloc(sizeof(a))

// Index of x/y coordinate
#define x (0)
#define y (1)

// Defines size of a block
#define BLOCK_DIM_X (3)
#define BLOCK_DIM_Y (2)

// Defines size of the grid, i.e., how many blocks
#define GRID_DIM_X (5)
#define GRID_DIM_Y (7)

// Defines the number of threads in the grid
#define GRID_SIZE (BLOCK_DIM_X * BLOCK_DIM_Y * GRID_DIM_X * GRID_DIM_Y)

// execution environment for the kernel
typedef struct exec_env {
   int threadIdx[2];  // thread location 
   int blockIdx[2];
   int blockDim[2];
   int gridDim[2];

   float *A,*B;       // parameters for the thread
   float *C;
} exec_env;

// kernel 
void *kernel(void *arg)
{
    exec_env *env = (exec_env *) arg;

    // compute number of threads in a block
    int sz = env->blockDim[x] * env->blockDim[y]; 

    // compute the index of the first thread in the block
    int k = sz * (env->blockIdx[y]*env->gridDim[x] + env->blockIdx[x]);

    // compute the index of a thread inside a block
    k = k + env->threadIdx[y]*env->blockDim[x] + env->threadIdx[x]; 

    // check whether it is in range
    assert(k >= 0 && k < GRID_SIZE && "Wrong index computation"); 

    // print coordinates in block and grid and computed index
    /*printf("tx:%d ty:%d bx:%d by:%d idx:%d\n",env->threadIdx[x],
                                              env->threadIdx[y],
                                              env->blockIdx[x],
                                              env->blockIdx[y], k);
    */

    // retrieve two operands 
    float *A = &env->A[k];
    float *B = &env->B[k]; 
    printf("%f %f \n",*A, *B);
    // retrieve pointer to result 
    float *C = &env->C[k]; 

    // do actual computation here !!!

    // For assignment replace the following line with 
    // the code to do matrix addition and multiplication.
    *C = *A + *B; 

    // free execution environment (not needed anymore)
    free(env);

    return NULL;
}


// main function
int main(int argc, char **argv)
{
   float A[GRID_SIZE] = {-1}; 
   float B[GRID_SIZE] = {-1};
   float C[GRID_SIZE] = {-1}; 

   pthread_t threads[GRID_SIZE]; 

   int i=0, bx, by, tx, ty; 
   //Error location
   /*for (i = 0; i < GRID_SIZE;i++){
        A[i] = i;
        B[i] = i+1;
        printf("%f %f\n ", A[i], B[i]);
   }*/

   // Step 1: create execution environment for threads and create thread
   for (bx=0;bx<GRID_DIM_X;bx++) {
      for (by=0;by<GRID_DIM_Y;by++) {
         for (tx=0;tx<BLOCK_DIM_X;tx++) {
            for (ty=0;ty<BLOCK_DIM_Y;ty++) { 

               exec_env *e = MALLOC(exec_env);
               assert(e != NULL && "memory exhausted"); 

               e->threadIdx[x]=tx;
               e->threadIdx[y]=ty;

               e->blockIdx[x]=bx;
               e->blockIdx[y]=by; 

               e->blockDim[x]=BLOCK_DIM_X;
               e->blockDim[y]=BLOCK_DIM_Y;

               e->gridDim[x]=GRID_DIM_X;
               e->gridDim[y]=GRID_DIM_Y;

               // set parameters
               e->A = A; 
               e->B = B; 
               e->C = C;

               // create thread
               pthread_create(&threads[i++],NULL,kernel,(void *)e); 
            }
         }
      }
   }

   // Step 2: wait for completion of all threads
   for (i=0;i<GRID_SIZE;i++) { 
      pthread_join(threads[i], NULL); 
   }

   // Step 3: print result       
   for (i=0;i<GRID_SIZE;i++) { 
      printf("%f ",C[i]); 
   }
   printf("\n");

   return 0;
}

Хорошо, здесь этот код работает нормально, но как только я раскомментирую «Местоположение ошибки» (для цикла, который назначает A [i] = i и B [i] = i + 1, я получаю ошибку из-за ошибки сегментации в Unix, и по этим случайным 0 в C в Cygwin. Я должен признать, что мои основы в C довольно плохие, поэтому весьма вероятно, что я что-то упустил. Если кто-то может дать представление о том, что идет не так, это было бы очень признательно Спасибо.

1 Ответ

5 голосов
/ 13 августа 2011

Это работает, когда вы комментируете это, потому что i по-прежнему равно 0, когда начинаются 4 вложенных цикла.

У вас есть это: 

for (i = 0; i < GRID_SIZE;i++){
    A[i] = i;
    B[i] = i+1;
    printf("%f %f\n ", A[i], B[i]);
}

/* What value is `i` now ? */

А потом 

pthread_create(&threads[i++],NULL,kernel,(void *)e);
                        ^

Так что pthread_create попытается получить доступ к некоторым интересным индексам.

...