CUDA優化問題

Question

這裏有一個簡單的程序：

void multiply(const int* v_in, const int* w_in, int n_v, int n_w, int* w_out) 
    { 
     for(int i=0; i<n_w; i++) 
     { 
     int sum=0; 
     for(int j=0; j<n_v; j++) 
      sum += (w_in[i]*v_in[j])>>1; 
     w_out[i]=sum; 
     } 
    } 

設定n_v，n_w〜10^6。很顯然，在CUDA中至少有十幾種等同的方法可以實現，用不同的方式將（n_v * n_w）操作細分爲線程，有或沒有共享內存......從理論上講，哪種方法最快？

Answer 1

簡單：

void multiply(const int* v_in, const int* w_in, int n_v, int n_w, int* w_out) 
    { 
     int *v = shared; // dynamic 
     for(int i = block.rank; i < n_w; i += block.size) 
     { 
     int w = w_in[i]; // coalesced 
     int sum=0; 
     for(int j=0; j<n_v; j += block.size) { // assumption 
      v[block.rank] = v_in[j+block.rank]; 
      __synch(); 
      for(int k = 0; k < block.size; ++k) 
       sum += (w*v[k])>>1; // 
      __synch(); // ouch 
     } 
     w_out[i] = sum; // ditto 
     } 
    }